A velocity-dependent anomalous radial transport model for (2-D, 2-V) kinetic transport codes
NASA Astrophysics Data System (ADS)
Bodi, Kowsik; Krasheninnikov, Sergei; Cohen, Ron; Rognlien, Tom
2008-11-01
Plasma turbulence constitutes a significant part of radial plasma transport in magnetically confined plasmas. This turbulent transport is modeled in the form of anomalous convection and diffusion coefficients in fluid transport codes. There is a need to model the same in continuum kinetic edge codes [such as the (2-D, 2-V) transport version of TEMPEST, NEO, and the code being developed by the Edge Simulation Laboratory] with non-Maxwellian distributions. We present an anomalous transport model with velocity-dependent convection and diffusion coefficients leading to a diagonal transport matrix similar to that used in contemporary fluid transport models (e.g., UEDGE). Also presented are results of simulations corresponding to radial transport due to long-wavelength ExB turbulence using a velocity-independent diffusion coefficient. A BGK collision model is used to enable comparison with fluid transport codes.
Wall surface temperature calculation in the SolEdge2D-EIRENE transport code
NASA Astrophysics Data System (ADS)
Denis, J.; Pégourié, B.; Bucalossi, J.; Bufferand, H.; Ciraolo, G.; Gardarein, J.-L.; Gaspar, J.; Grisolia, C.; Hodille, E.; Missirlian, M.; Serre, E.; Tamain, P.
2016-02-01
A thermal wall model is developed for the SolEdge2D-EIRENE edge transport code for calculating the surface temperature of the actively-cooled vessel components in interaction with the plasma. This is a first step towards a self-consistent evaluation of the recycling of particles, which depends on the wall surface temperature. The proposed thermal model is built to match both steady-state temperature and time constant of actively-cooled plasma facing components. A benchmark between this model and the Finite Element Modelling code CAST3M is performed in the case of an ITER-like monoblock. An example of application is presented for a SolEdge2D-EIRENE simulation of a medium-power discharge in the WEST tokamak, showing the steady-state wall temperature distribution and the temperature cycling due to an imposed Edge Localised Mode-like event.
Assessment of the 2D MOC solver in MPACT: Michigan parallel characteristics transport code
Collins, B.; Kochunas, B.; Downar, T.
2013-07-01
MPACT (Michigan Parallel Characteristics Transport Code) is a new reactor analysis tool being developed by researchers at the University of Michigan as an advanced pin-resolved transport capability within VERA (Virtual Environment for Reactor Analysis). VERA is the end-user reactor simulation tool being developed by the Consortium for the Advanced Simulation of Light Water Reactors (CASL). The MPACT development project is itself unique for the way it is changing how students perform research to achieve the instructional and research goals of an academic institution, while providing immediate value to the industry. One of the major computational pieces in MPACT is the 2D MOC solver. It is critical that the 2D MOC solver provide an efficient, accurate, and robust solution over a broad range of reactor operating conditions. The C5G7 benchmark is first used to test the accuracy of the method with a fixed set of cross-sections. The VERA Core Physics Progression Problems are then used to compare the accuracy of both the 2D transport solver and also the cross-section treatments. (authors)
NASA Astrophysics Data System (ADS)
Leboeuf, Jean-Noel; Decyk, Viktor; Newman, David; Sanchez, Raul
2012-03-01
The massively parallel, nonlinear, 3D, toroidal, electrostatic, gyrokinetic, PIC, Cartesian geometry UCAN code, with particle ions and adiabatic electrons, has been successfully exercised to identify non-diffusive transport characteristics in DIII-D-like tokamak discharges. The limitation in applying UCAN to larger scale discharges is the 1D domain decomposition in the toroidal (or z-) direction for massively parallel implementation using MPI which has restricted the calculations to a few hundred ion Larmor radii per minor radius. To exceed these sizes, we have implemented 2D domain decomposition in UCAN with the addition of the y-direction to the processor mix. This has been facilitated by use of relevant components in the 2D domain decomposed PLIB2 library of field and particle management routines developed for UCLA's UPIC framework of conventional PIC codes. The gyro-averaging in gyrokinetic codes has necessitated the use of replicated arrays for efficient charge accumulation and particle push. The 2D domain-decomposed UCAN2 code reproduces the original 1D domain results within roundoff. Production calculations at large system sizes have been performed with UCAN2 on 131072 processors of the Cray XE6 at NERSC.
NASA Astrophysics Data System (ADS)
Yamada, Susumu; Kitamura, Akihiro; Kurikami, Hiroshi; Machida, Masahiko
2015-04-01
Fukushima Daiichi Nuclear Power Plant (FDNPP) accident on March 2011 released significant quantities of radionuclides to atmosphere. The most significant nuclide is radioactive cesium isotopes. Therefore, the movement of the cesium is one of the critical issues for the environmental assessment. Since the cesium is strongly sorbed by soil particles, the cesium transport can be regarded as the sediment transport which is mainly brought about by the aquatic system such as a river and a lake. In this research, our target is the sediment transport on Ogaki dam reservoir which is located in about 16 km northwest from FDNPP. The reservoir is one of the principal irrigation dam reservoirs in Fukushima Prefecture and its upstream river basin was heavily contaminated by radioactivity. We simulate the sediment transport on the reservoir using 2-D river simulation code named Nays2D originally developed by Shimizu et al. (The latest version of Nays2D is available as a code included in iRIC (http://i-ric.org/en/), which is a river flow and riverbed variation analysis software package). In general, a 2-D simulation code requires a huge amount of calculation time. Therefore, we parallelize the code and execute it on a parallel computer. We examine the relationship between the behavior of the sediment transport and the height of the reservoir exit. The simulation result shows that almost all the sand that enter into the reservoir deposit close to the entrance of the reservoir for any height of the exit. The amounts of silt depositing within the reservoir slightly increase by raising the height of the exit. However, that of the clay dramatically increases. Especially, more than half of the clay deposits, if the exit is sufficiently high. These results demonstrate that the water level of the reservoir has a strong influence on the amount of the clay discharged from the reservoir. As a result, we conclude that the tuning of the water level has a possibility for controlling the
NASA Astrophysics Data System (ADS)
Kirnev, G.; Fundamenski, W.; Corrigan, G.
2007-06-01
The scrape-off layer (SOL) of the JET tokamak has been modelled using a two-dimensional plasma/neutral code, EDGE2D/NIMBUS, with variable transport coefficients, chosen according to nine candidate theories for radial heat transport in the SOL. Comparison of the radial power width on the outer divertor plates, λq, predicted by modelling and measured experimentally in L-mode and ELM-averaged H-mode at JET is presented. Transport coefficients based on classical and neo-classical ion conduction are found to offer the best agreement with experimentally measured λq magnitude and scaling with target power, upstream density and toroidal field. These results reinforce the findings of an earlier study, based on a simplified model of the SOL (Chankin 1997 Plasma Phys. Control. Fusion 39 1059), and support the earlier estimate of the power width at the entrance of the outer divertor volume in ITER, λq ap 4 mm mapped to the outer mid-plane (Fundamenski et al 2004 Nucl. Fusion 44 20).
Mosleh-Shirazi, Mohammad Amin; Zarrini-Monfared, Zinat; Karbasi, Sareh; Zamani, Ali
2014-01-01
Two-dimensional (2D) arrays of thick segmented scintillators are of interest as X-ray detectors for both 2D and 3D image-guided radiotherapy (IGRT). Their detection process involves ionizing radiation energy deposition followed by production and transport of optical photons. Only a very limited number of optical Monte Carlo simulation models exist, which has limited the number of modeling studies that have considered both stages of the detection process. We present ScintSim1, an in-house optical Monte Carlo simulation code for 2D arrays of scintillation crystals, developed in the MATLAB programming environment. The code was rewritten and revised based on an existing program for single-element detectors, with the additional capability to model 2D arrays of elements with configurable dimensions, material, etc., The code generates and follows each optical photon history through the detector element (and, in case of cross-talk, the surrounding ones) until it reaches a configurable receptor, or is attenuated. The new model was verified by testing against relevant theoretically known behaviors or quantities and the results of a validated single-element model. For both sets of comparisons, the discrepancies in the calculated quantities were all <1%. The results validate the accuracy of the new code, which is a useful tool in scintillation detector optimization. PMID:24600168
Mosleh-Shirazi, Mohammad Amin; Zarrini-Monfared, Zinat; Karbasi, Sareh; Zamani, Ali
2014-01-01
Two-dimensional (2D) arrays of thick segmented scintillators are of interest as X-ray detectors for both 2D and 3D image-guided radiotherapy (IGRT). Their detection process involves ionizing radiation energy deposition followed by production and transport of optical photons. Only a very limited number of optical Monte Carlo simulation models exist, which has limited the number of modeling studies that have considered both stages of the detection process. We present ScintSim1, an in-house optical Monte Carlo simulation code for 2D arrays of scintillation crystals, developed in the MATLAB programming environment. The code was rewritten and revised based on an existing program for single-element detectors, with the additional capability to model 2D arrays of elements with configurable dimensions, material, etc., The code generates and follows each optical photon history through the detector element (and, in case of cross-talk, the surrounding ones) until it reaches a configurable receptor, or is attenuated. The new model was verified by testing against relevant theoretically known behaviors or quantities and the results of a validated single-element model. For both sets of comparisons, the discrepancies in the calculated quantities were all <1%. The results validate the accuracy of the new code, which is a useful tool in scintillation detector optimization. PMID:24600168
2-d Finite Element Code Postprocessor
Energy Science and Technology Software Center (ESTSC)
1996-07-15
ORION is an interactive program that serves as a postprocessor for the analysis programs NIKE2D, DYNA2D, TOPAZ2D, and CHEMICAL TOPAZ2D. ORION reads binary plot files generated by the two-dimensional finite element codes currently used by the Methods Development Group at LLNL. Contour and color fringe plots of a large number of quantities may be displayed on meshes consisting of triangular and quadrilateral elements. ORION can compute strain measures, interface pressures along slide lines, reaction forcesmore » along constrained boundaries, and momentum. ORION has been applied to study the response of two-dimensional solids and structures undergoing finite deformations under a wide variety of large deformation transient dynamic and static problems and heat transfer analyses.« less
A Computer Code for 2-D Transport Calculations in x-y Geometry Using the Interface Current Method.
Energy Science and Technology Software Center (ESTSC)
1990-12-01
Version 00 RICANT performs 2-dimensional neutron transport calculations in x-y geometry using the interface current method. In the interface current method, the angular neutron currents crossing region surfaces are expanded in terms of the Legendre polynomials in the two half-spaces made by the region surfaces.
MAGNUM-2D computer code: user's guide
England, R.L.; Kline, N.W.; Ekblad, K.J.; Baca, R.G.
1985-01-01
Information relevant to the general use of the MAGNUM-2D computer code is presented. This computer code was developed for the purpose of modeling (i.e., simulating) the thermal and hydraulic conditions in the vicinity of a waste package emplaced in a deep geologic repository. The MAGNUM-2D computer computes (1) the temperature field surrounding the waste package as a function of the heat generation rate of the nuclear waste and thermal properties of the basalt and (2) the hydraulic head distribution and associated groundwater flow fields as a function of the temperature gradients and hydraulic properties of the basalt. MAGNUM-2D is a two-dimensional numerical model for transient or steady-state analysis of coupled heat transfer and groundwater flow in a fractured porous medium. The governing equations consist of a set of coupled, quasi-linear partial differential equations that are solved using a Galerkin finite-element technique. A Newton-Raphson algorithm is embedded in the Galerkin functional to formulate the problem in terms of the incremental changes in the dependent variables. Both triangular and quadrilateral finite elements are used to represent the continuum portions of the spatial domain. Line elements may be used to represent discrete conduits. 18 refs., 4 figs., 1 tab.
Validation and testing of the VAM2D computer code
Kool, J.B.; Wu, Y.S. )
1991-10-01
This document describes two modeling studies conducted by HydroGeoLogic, Inc. for the US NRC under contract no. NRC-04089-090, entitled, Validation and Testing of the VAM2D Computer Code.'' VAM2D is a two-dimensional, variably saturated flow and transport code, with applications for performance assessment of nuclear waste disposal. The computer code itself is documented in a separate NUREG document (NUREG/CR-5352, 1989). The studies presented in this report involve application of the VAM2D code to two diverse subsurface modeling problems. The first one involves modeling of infiltration and redistribution of water and solutes in an initially dry, heterogeneous field soil. This application involves detailed modeling over a relatively short, 9-month time period. The second problem pertains to the application of VAM2D to the modeling of a waste disposal facility in a fractured clay, over much larger space and time scales and with particular emphasis on the applicability and reliability of using equivalent porous medium approach for simulating flow and transport in fractured geologic media. Reflecting the separate and distinct nature of the two problems studied, this report is organized in two separate parts. 61 refs., 31 figs., 9 tabs.
Ion Transport in 2-D Graphene Nanochannels
NASA Astrophysics Data System (ADS)
Xie, Quan; Foo, Elbert; Duan, Chuanhua
2015-11-01
Graphene membranes have recently attracted wide attention due to its great potential in water desalination and selective molecular sieving. Further developments of these membranes, including enhancing their mass transport rate and/or molecular selectivity, rely on the understanding of fundamental transport mechanisms through graphene membranes, which has not been studied experimentally before due to fabrication and measurement difficulties. Herein we report the fabrication of the basic constituent of graphene membranes, i.e. 2-D single graphene nanochannels (GNCs) and the study of ion transport in these channels. A modified bonding technique was developed to form GNCs with well-defined geometry and uniform channel height. Ion transport in such GNCs was studied using DC conductance measurement. Our preliminary results showed that the ion transport in GNCs is still governed by surface charge at low concentrations (10-6M to 10-4M). However, GNCs exhibits much higher ionic conductances than silica nanochannels with the same geometries in the surface-charge-governed regime. This conductance enhancement can be attributed to the pre-accumulation of charges on graphene surfaces. The work is supported by the Faculty Startup Fund (Boston University, USA).
In search of a 2-dB coding gain
NASA Technical Reports Server (NTRS)
Yuen, J. H.; Vo, Q. D.
1985-01-01
A recent code search found a (15,1/5), a (14,1/6), and a (15,1/6) convolutional code which, when concatenated with a 10-bit (1023,959) Reed-Solomon (RS) code, achieves a bit-error rate (BER) of 0.000001 at a bit signal-to-noise ratio (SNR) of 0.50 dB, 0.47 dB and 0.42 B, respectively. All of these three codes outperform the Voyager communication system, our baseline, which achieves a BER of 10.000001 at bit SNR of 2.53 db, by more than 2 dB. The 2 dB coding improvement goal was exceeded.
Energy Science and Technology Software Center (ESTSC)
2010-02-01
Neutron transport, calculation of multiplication factor and neutron fluxes in 2-D configurations: cell calculations, 2-D diffusion and transport, and burnup. Preparation of a cross section library for the code BOXER from a basic library in ENDF/B format (ETOBOX).
Recent update of the RPLUS2D/3D codes
NASA Technical Reports Server (NTRS)
Tsai, Y.-L. Peter
1991-01-01
The development of the RPLUS2D/3D codes is summarized. These codes utilize LU algorithms to solve chemical non-equilibrium flows in a body-fitted coordinate system. The motivation behind the development of these codes is the need to numerically predict chemical non-equilibrium flows for the National AeroSpace Plane Program. Recent improvements include vectorization method, blocking algorithms for geometric flexibility, out-of-core storage for large-size problems, and an LU-SW/UP combination for CPU-time efficiency and solution quality.
AnisWave2D: User's Guide to the 2d Anisotropic Finite-DifferenceCode
Toomey, Aoife
2005-01-06
This document describes a parallel finite-difference code for modeling wave propagation in 2D, fully anisotropic materials. The code utilizes a mesh refinement scheme to improve computational efficiency. Mesh refinement allows the grid spacing to be tailored to the velocity model, so that fine grid spacing can be used in low velocity zones where the seismic wavelength is short, and coarse grid spacing can be used in zones with higher material velocities. Over-sampling of the seismic wavefield in high velocity zones is therefore avoided. The code has been implemented to run in parallel over multiple processors and allows large-scale models and models with large velocity contrasts to be simulated with ease.
ELLIPT2D: A Flexible Finite Element Code Written Python
Pletzer, A.; Mollis, J.C.
2001-03-22
The use of the Python scripting language for scientific applications and in particular to solve partial differential equations is explored. It is shown that Python's rich data structure and object-oriented features can be exploited to write programs that are not only significantly more concise than their counter parts written in Fortran, C or C++, but are also numerically efficient. To illustrate this, a two-dimensional finite element code (ELLIPT2D) has been written. ELLIPT2D provides a flexible and easy-to-use framework for solving a large class of second-order elliptic problems. The program allows for structured or unstructured meshes. All functions defining the elliptic operator are user supplied and so are the boundary conditions, which can be of Dirichlet, Neumann or Robbins type. ELLIPT2D makes extensive use of dictionaries (hash tables) as a way to represent sparse matrices.Other key features of the Python language that have been widely used include: operator over loading, error handling, array slicing, and the Tkinter module for building graphical use interfaces. As an example of the utility of ELLIPT2D, a nonlinear solution of the Grad-Shafranov equation is computed using a Newton iterative scheme. A second application focuses on a solution of the toroidal Laplace equation coupled to a magnetohydrodynamic stability code, a problem arising in the context of magnetic fusion research.
CBEAM. 2-D: a two-dimensional beam field code
Dreyer, K.A.
1985-05-01
CBEAM.2-D is a two-dimensional solution of Maxwell's equations for the case of an electron beam propagating through an air medium. Solutions are performed in the beam-retarded time frame. Conductivity is calculated self-consistently with field equations, allowing sophisticated dependence of plasma parameters to be handled. A unique feature of the code is that it is implemented on an IBM PC microcomputer in the BASIC language. Consequently, it should be available to a wide audience.
2D FEM Heat Transfer & E&M Field Code
Energy Science and Technology Software Center (ESTSC)
1992-04-02
TOPAZ and TOPAZ2D are two-dimensional implicit finite element computer codes for heat transfer analysis. TOPAZ2D can also be used to solve electrostatic and magnetostatic problems. The programs solve for the steady-state or transient temperature or electrostatic and magnetostatic potential field on two-dimensional planar or axisymmetric geometries. Material properties may be temperature or potential-dependent and either isotropic or orthotropic. A variety of time and temperature-dependent boundary conditions can be specified including temperature, flux, convection, and radiation.more » By implementing the user subroutine feature, users can model chemical reaction kinetics and allow for any type of functional representation of boundary conditions and internal heat generation. The programs can solve problems of diffuse and specular band radiation in an enclosure coupled with conduction in the material surrounding the enclosure. Additional features include thermal contact resistance across an interface, bulk fluids, phase change, and energy balances.« less
2D FEM Heat Transfer & E&M Field Code
1992-04-02
TOPAZ and TOPAZ2D are two-dimensional implicit finite element computer codes for heat transfer analysis. TOPAZ2D can also be used to solve electrostatic and magnetostatic problems. The programs solve for the steady-state or transient temperature or electrostatic and magnetostatic potential field on two-dimensional planar or axisymmetric geometries. Material properties may be temperature or potential-dependent and either isotropic or orthotropic. A variety of time and temperature-dependent boundary conditions can be specified including temperature, flux, convection, and radiation. By implementing the user subroutine feature, users can model chemical reaction kinetics and allow for any type of functional representation of boundary conditions and internal heat generation. The programs can solve problems of diffuse and specular band radiation in an enclosure coupled with conduction in the material surrounding the enclosure. Additional features include thermal contact resistance across an interface, bulk fluids, phase change, and energy balances.
Modelling RF sources using 2-D PIC codes
Eppley, K.R.
1993-03-01
In recent years, many types of RF sources have been successfully modelled using 2-D PIC codes. Both cross field devices (magnetrons, cross field amplifiers, etc.) and pencil beam devices (klystrons, gyrotrons, TWT`S, lasertrons, etc.) have been simulated. All these devices involve the interaction of an electron beam with an RF circuit. For many applications, the RF structure may be approximated by an equivalent circuit, which appears in the simulation as a boundary condition on the electric field (``port approximation``). The drive term for the circuit is calculated from the energy transfer between beam and field in the drift space. For some applications it may be necessary to model the actual geometry of the structure, although this is more expensive. One problem not entirely solved is how to accurately model in 2-D the coupling to an external waveguide. Frequently this is approximated by a radial transmission line, but this sometimes yields incorrect results. We also discuss issues in modelling the cathode and injecting the beam into the PIC simulation.
Modelling RF sources using 2-D PIC codes
Eppley, K.R.
1993-03-01
In recent years, many types of RF sources have been successfully modelled using 2-D PIC codes. Both cross field devices (magnetrons, cross field amplifiers, etc.) and pencil beam devices (klystrons, gyrotrons, TWT'S, lasertrons, etc.) have been simulated. All these devices involve the interaction of an electron beam with an RF circuit. For many applications, the RF structure may be approximated by an equivalent circuit, which appears in the simulation as a boundary condition on the electric field ( port approximation''). The drive term for the circuit is calculated from the energy transfer between beam and field in the drift space. For some applications it may be necessary to model the actual geometry of the structure, although this is more expensive. One problem not entirely solved is how to accurately model in 2-D the coupling to an external waveguide. Frequently this is approximated by a radial transmission line, but this sometimes yields incorrect results. We also discuss issues in modelling the cathode and injecting the beam into the PIC simulation.
Domino, Stefan; Luketa-Hanlin, Anay; Gallegos, Carlos
2006-10-27
FAA Smoke Transport Code, a physics-based Computational Fluid Dynamics tool, which couples heat, mass, and momentum transfer, has been developed to provide information on smoke transport in cargo compartments with various geometries and flight conditions. The software package contains a graphical user interface for specification of geometry and boundary conditions, analysis module for solving the governing equations, and a post-processing tool. The current code was produced by making substantial improvements and additions to a code obtained from a university. The original code was able to compute steady, uniform, isothermal turbulent pressurization. In addition, a preprocessor and postprocessor were added to arrive at the current software package.
50 CFR Table 2d to Part 679 - Species Codes-Non-FMP Species
Code of Federal Regulations, 2012 CFR
2012-10-01
... 50 Wildlife and Fisheries 13 2012-10-01 2012-10-01 false Species Codes-Non-FMP Species 2d Table 2d to Part 679 Wildlife and Fisheries FISHERY CONSERVATION AND MANAGEMENT, NATIONAL OCEANIC AND ATMOSPHERIC ADMINISTRATION, DEPARTMENT OF COMMERCE (CONTINUED) FISHERIES OF THE EXCLUSIVE ECONOMIC ZONE OFF ALASKA Pt. 679, Table 2d Table 2d to Part...
50 CFR Table 2d to Part 679 - Species Codes-Non-FMP Species
Code of Federal Regulations, 2013 CFR
2013-10-01
... 50 Wildlife and Fisheries 13 2013-10-01 2013-10-01 false Species Codes-Non-FMP Species 2d Table 2d to Part 679 Wildlife and Fisheries FISHERY CONSERVATION AND MANAGEMENT, NATIONAL OCEANIC AND... ALASKA Pt. 679, Table 2d Table 2d to Part 679—Species Codes—Non-FMP Species General use...
50 CFR Table 2d to Part 679 - Species Codes-Non-FMP Species
Code of Federal Regulations, 2014 CFR
2014-10-01
... 50 Wildlife and Fisheries 13 2014-10-01 2014-10-01 false Species Codes-Non-FMP Species 2d Table 2d to Part 679 Wildlife and Fisheries FISHERY CONSERVATION AND MANAGEMENT, NATIONAL OCEANIC AND... ALASKA Pt. 679, Table 2d Table 2d to Part 679—Species Codes—Non-FMP Species General use...
2-D skin-current toroidal-MHD-equilibrium code
Feinberg, B.; Niland, R.A.; Coonrod, J.; Levine, M.A.
1982-09-01
A two-dimensional, toroidal, ideal MHD skin-current equilibrium computer code is described. The code is suitable for interactive implementation on a minicomptuer. Some examples of the use of the code for design and interpretation of toroidal cusp experiments are presented.
Energy Science and Technology Software Center (ESTSC)
2006-10-27
FAA Smoke Transport Code, a physics-based Computational Fluid Dynamics tool, which couples heat, mass, and momentum transfer, has been developed to provide information on smoke transport in cargo compartments with various geometries and flight conditions. The software package contains a graphical user interface for specification of geometry and boundary conditions, analysis module for solving the governing equations, and a post-processing tool. The current code was produced by making substantial improvements and additions to a codemore » obtained from a university. The original code was able to compute steady, uniform, isothermal turbulent pressurization. In addition, a preprocessor and postprocessor were added to arrive at the current software package.« less
Transport Experiments on 2D Correlated Electron Physics in Semiconductors
Tsui, Daniel
2014-03-24
This research project was designed to investigate experimentally the transport properties of the 2D electrons in Si and GaAs, two prototype semiconductors, in several new physical regimes that were previously inaccessible to experiments. The research focused on the strongly correlated electron physics in the dilute density limit, where the electron potential energy to kinetic energy ratio rs>>1, and on the fractional quantum Hall effect related physics in nuclear demagnetization refrigerator temperature range on samples with new levels of purity and controlled random disorder.
A New 2D-Transport, 1D-Diffusion Approximation of the Boltzmann Transport equation
Larsen, Edward
2013-06-17
The work performed in this project consisted of the derivation, implementation, and testing of a new, computationally advantageous approximation to the 3D Boltz- mann transport equation. The solution of the Boltzmann equation is the neutron flux in nuclear reactor cores and shields, but solving this equation is difficult and costly. The new “2D/1D” approximation takes advantage of a special geometric feature of typical 3D reactors to approximate the neutron transport physics in a specific (ax- ial) direction, but not in the other two (radial) directions. The resulting equation is much less expensive to solve computationally, and its solutions are expected to be sufficiently accurate for many practical problems. In this project we formulated the new equation, discretized it using standard methods, developed a stable itera- tion scheme for solving the equation, implemented the new numerical scheme in the MPACT code, and tested the method on several realistic problems. All the hoped- for features of this new approximation were seen. For large, difficult problems, the resulting 2D/1D solution is highly accurate, and is calculated about 100 times faster than a 3D discrete ordinates simulation.
NASA Technical Reports Server (NTRS)
Kapoor, Kamlesh; Anderson, Bernhard H.; Shaw, Robert J.
1994-01-01
A two-dimensional computational code, PRLUS2D, which was developed for the reactive propulsive flows of ramjets and scramjets, was validated for two-dimensional shock-wave/turbulent-boundary-layer interactions. The problem of compression corners at supersonic speeds was solved using the RPLUS2D code. To validate the RPLUS2D code for hypersonic speeds, it was applied to a realistic hypersonic inlet geometry. Both the Baldwin-Lomax and the Chien two-equation turbulence models were used. Computational results showed that the RPLUS2D code compared very well with experimentally obtained data for supersonic compression corner flows, except in the case of large separated flows resulting from the interactions between the shock wave and turbulent boundary layer. The computational results compared well with the experiment results in a hypersonic NASA P8 inlet case, with the Chien two-equation turbulence model performing better than the Baldwin-Lomax model.
TOPAZ2D heat transfer code users manual and thermal property data base
NASA Astrophysics Data System (ADS)
Shapiro, A. B.; Edwards, A. L.
1990-05-01
TOPAZ2D is a two dimensional implicit finite element computer code for heat transfer analysis. This user's manual provides information on the structure of a TOPAZ2D input file. Also included is a material thermal property data base. This manual is supplemented with The TOPAZ2D Theoretical Manual and the TOPAZ2D Verification Manual. TOPAZ2D has been implemented on the CRAY, SUN, and VAX computers. TOPAZ2D can be used to solve for the steady state or transient temperature field on two dimensional planar or axisymmetric geometries. Material properties may be temperature dependent and either isotropic or orthotropic. A variety of time and temperature dependent boundary conditions can be specified including temperature, flux, convection, and radiation. Time or temperature dependent internal heat generation can be defined locally be element or globally by material. TOPAZ2D can solve problems of diffuse and specular band radiation in an enclosure coupled with conduction in material surrounding the enclosure. Additional features include thermally controlled reactive chemical mixtures, thermal contact resistance across an interface, bulk fluid flow, phase change, and energy balances. Thermal stresses can be calculated using the solid mechanics code NIKE2D which reads the temperature state data calculated by TOPAZ2D. A three dimensional version of the code, TOPAZ3D is available.
2D Quantum Transport Modeling in Nanoscale MOSFETs
NASA Technical Reports Server (NTRS)
Svizhenko, Alexei; Anantram, M. P.; Govindan, T. R.; Biegel, Bryan
2001-01-01
With the onset of quantum confinement in the inversion layer in nanoscale MOSFETs, behavior of the resonant level inevitably determines all device characteristics. While most classical device simulators take quantization into account in some simplified manner, the important details of electrostatics are missing. Our work addresses this shortcoming and provides: (a) a framework to quantitatively explore device physics issues such as the source-drain and gate leakage currents, DIBL, and threshold voltage shift due to quantization, and b) a means of benchmarking quantum corrections to semiclassical models (such as density- gradient and quantum-corrected MEDICI). We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions, oxide tunneling and phase-breaking scattering are treated on equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. Quantum simulations are focused on MIT 25, 50 and 90 nm "well- tempered" MOSFETs and compared to classical and quantum corrected models. The important feature of quantum model is smaller slope of Id-Vg curve and consequently higher threshold voltage. These results are quantitatively consistent with I D Schroedinger-Poisson calculations. The effect of gate length on gate-oxide leakage and sub-threshold current has been studied. The shorter gate length device has an order of magnitude smaller current at zero gate bias than the longer gate length device without a significant trade-off in on-current. This should be a device design consideration.
Automatic differentiation of the TACO2D finite element code using ADIFOR
Carle, A.; Fagan, M.
1996-04-01
The need for sensitivities in particular applications is becoming increasingly important in problems such as optimal design or control. In this study, the authors use ADIFOR to generate derivative code for TACO2D, a finite element heat transfer code. The study of TACO2D indicates that ADIFOR-generated derivatives yield accurate derivatives at a fraction of the time requirements of finite difference approximations, and space requirements proportional to the number of variables. The primary focus on TACO2D was for the design of chemical vapor deposition reactors.
TOPAZ2D heat transfer code users manual and thermal property data base
Shapiro, A.B.; Edwards, A.L.
1990-05-01
TOPAZ2D is a two dimensional implicit finite element computer code for heat transfer analysis. This user's manual provides information on the structure of a TOPAZ2D input file. Also included is a material thermal property data base. This manual is supplemented with The TOPAZ2D Theoretical Manual and the TOPAZ2D Verification Manual. TOPAZ2D has been implemented on the CRAY, SUN, and VAX computers. TOPAZ2D can be used to solve for the steady state or transient temperature field on two dimensional planar or axisymmetric geometries. Material properties may be temperature dependent and either isotropic or orthotropic. A variety of time and temperature dependent boundary conditions can be specified including temperature, flux, convection, and radiation. Time or temperature dependent internal heat generation can be defined locally be element or globally by material. TOPAZ2D can solve problems of diffuse and specular band radiation in an enclosure coupled with conduction in material surrounding the enclosure. Additional features include thermally controlled reactive chemical mixtures, thermal contact resistance across an interface, bulk fluid flow, phase change, and energy balances. Thermal stresses can be calculated using the solid mechanics code NIKE2D which reads the temperature state data calculated by TOPAZ2D. A three dimensional version of the code, TOPAZ3D is available. The material thermal property data base, Chapter 4, included in this manual was originally published in 1969 by Art Edwards for use with his TRUMP finite difference heat transfer code. The format of the data has been altered to be compatible with TOPAZ2D. Bob Bailey is responsible for adding the high explosive thermal property data.
PiCode: A New Picture-Embedding 2D Barcode.
Chen, Changsheng; Huang, Wenjian; Zhou, Baojian; Liu, Chenchen; Mow, Wai Ho
2016-08-01
Nowadays, 2D barcodes have been widely used as an interface to connect potential customers and advertisement contents. However, the appearance of a conventional 2D barcode pattern is often too obtrusive for integrating into an aesthetically designed advertisement. Besides, no human readable information is provided before the barcode is successfully decoded. This paper proposes a new picture-embedding 2D barcode, called PiCode, which mitigates these two limitations by equipping a scannable 2D barcode with a picturesque appearance. PiCode is designed with careful considerations on both the perceptual quality of the embedded image and the decoding robustness of the encoded message. Comparisons with the existing beautified 2D barcodes show that PiCode achieves one of the best perceptual qualities for the embedded image, and maintains a better tradeoff between image quality and decoding robustness in various application conditions. PiCode has been implemented in the MATLAB on a PC and some key building blocks have also been ported to Android and iOS platforms. Its practicality for real-world applications has been successfully demonstrated. PMID:27249833
2D MIMO Network Coding with Inter-Route Interference Cancellation
NASA Astrophysics Data System (ADS)
Tran, Gia Khanh; Sakaguchi, Kei; Ono, Fumie; Araki, Kiyomichi
Infrastructure wireless mesh network has been attracting much attention due to the wide range of its application such as public wireless access, sensor network, etc. In recent years, researchers have shown that significant network throughput gain can be achieved by employing network coding in a wireless environment. For further improvement of network throughput in one dimensional (1D) topology, Ono et al. proposed to use multiple antenna technique combined with network coding. In this paper, being inspired by MIMO network coding in 1D topology, the authors establish a novel MIMO network coding algorithm for a 2D topology consisting of two crossing routes. In this algorithm, multiple network coded flows are spatially multiplexed. Owing to the efficient usage of radio resource of network coding and co-channel interference cancellation ability of MIMO, the proposed algorithm shows an 8-fold gain in network capacity compared to conventional methods in the best-case scenario.
Icarus: A 2-D Direct Simulation Monte Carlo (DSMC) Code for Multi-Processor Computers
BARTEL, TIMOTHY J.; PLIMPTON, STEVEN J.; GALLIS, MICHAIL A.
2001-10-01
Icarus is a 2D Direct Simulation Monte Carlo (DSMC) code which has been optimized for the parallel computing environment. The code is based on the DSMC method of Bird[11.1] and models from free-molecular to continuum flowfields in either cartesian (x, y) or axisymmetric (z, r) coordinates. Computational particles, representing a given number of molecules or atoms, are tracked as they have collisions with other particles or surfaces. Multiple species, internal energy modes (rotation and vibration), chemistry, and ion transport are modeled. A new trace species methodology for collisions and chemistry is used to obtain statistics for small species concentrations. Gas phase chemistry is modeled using steric factors derived from Arrhenius reaction rates or in a manner similar to continuum modeling. Surface chemistry is modeled with surface reaction probabilities; an optional site density, energy dependent, coverage model is included. Electrons are modeled by either a local charge neutrality assumption or as discrete simulational particles. Ion chemistry is modeled with electron impact chemistry rates and charge exchange reactions. Coulomb collision cross-sections are used instead of Variable Hard Sphere values for ion-ion interactions. The electro-static fields can either be: externally input, a Langmuir-Tonks model or from a Green's Function (Boundary Element) based Poison Solver. Icarus has been used for subsonic to hypersonic, chemically reacting, and plasma flows. The Icarus software package includes the grid generation, parallel processor decomposition, post-processing, and restart software. The commercial graphics package, Tecplot, is used for graphics display. All of the software packages are written in standard Fortran.
Snapshot 2D tomography via coded aperture x-ray scatter imaging
MacCabe, Kenneth P.; Holmgren, Andrew D.; Tornai, Martin P.; Brady, David J.
2015-01-01
This paper describes a fan beam coded aperture x-ray scatter imaging system which acquires a tomographic image from each snapshot. This technique exploits cylindrical symmetry of the scattering cross section to avoid the scanning motion typically required by projection tomography. We use a coded aperture with a harmonic dependence to determine range, and a shift code to determine cross-range. Here we use a forward-scatter configuration to image 2D objects and use serial exposures to acquire tomographic video of motion within a plane. Our reconstruction algorithm also estimates the angular dependence of the scattered radiance, a step toward materials imaging and identification. PMID:23842254
Field depth extension of 2D barcode scanner based on wavefront coding and projection algorithm
NASA Astrophysics Data System (ADS)
Zhao, Tingyu; Ye, Zi; Zhang, Wenzi; Huang, Weiwei; Yu, Feihong
2008-03-01
Wavefront coding (WFC) used in 2D barcode scanners can extend the depth of field into a great extent with simpler structure compared to the autofocus microscope system. With a cubic phase mask (CPM) employed in the STOP, blurred images will be obtained in charge coupled device (CCD), which can be restored by digital filters. Direct methods are used widely in real-time restoration with good computational efficiency but with details smoothed. Here, the results of direct method are firstly filtered by hard-threshold function. The positions of the steps can be detected by simple differential operators. With the positions corrected by projection algorithm, the exact barcode information is restored. A wavefront coding system with 7mm effective focal length and 6 F-number is designed as an example. Although with the different magnification, images of different object distances can be restored by one point spread function (PSF) with 200mm object distance. A QR code (Quickly Response Code) of 31mm X 27mm is used as a target object. The simulation results showed that the sharp imaging objective distance is from 80mm to 355mm. The 2D barcode scanner with wavefront coding extends field depth with simple structure, low cost and large manufacture tolerance. This combination of the direct filter and projection algorithm proposed here could get the exact 2D barcode information with good computational efficiency.
MULTI2D - a computer code for two-dimensional radiation hydrodynamics
NASA Astrophysics Data System (ADS)
Ramis, R.; Meyer-ter-Vehn, J.; Ramírez, J.
2009-06-01
Simulation of radiation hydrodynamics in two spatial dimensions is developed, having in mind, in particular, target design for indirectly driven inertial confinement energy (IFE) and the interpretation of related experiments. Intense radiation pulses by laser or particle beams heat high-Z target configurations of different geometries and lead to a regime which is optically thick in some regions and optically thin in others. A diffusion description is inadequate in this situation. A new numerical code has been developed which describes hydrodynamics in two spatial dimensions (cylindrical R-Z geometry) and radiation transport along rays in three dimensions with the 4 π solid angle discretized in direction. Matter moves on a non-structured mesh composed of trilateral and quadrilateral elements. Radiation flux of a given direction enters on two (one) sides of a triangle and leaves on the opposite side(s) in proportion to the viewing angles depending on the geometry. This scheme allows to propagate sharply edged beams without ray tracing, though at the price of some lateral diffusion. The algorithm treats correctly both the optically thin and optically thick regimes. A symmetric semi-implicit (SSI) method is used to guarantee numerical stability. Program summaryProgram title: MULTI2D Catalogue identifier: AECV_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AECV_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 151 098 No. of bytes in distributed program, including test data, etc.: 889 622 Distribution format: tar.gz Programming language: C Computer: PC (32 bits architecture) Operating system: Linux/Unix RAM: 2 Mbytes Word size: 32 bits Classification: 19.7 External routines: X-window standard library (libX11.so) and corresponding heading files (X11/*.h) are
NASA Astrophysics Data System (ADS)
Stone, James M.; Norman, Michael L.
1992-06-01
A detailed description of ZEUS-2D, a numerical code for the simulation of fluid dynamical flows including a self-consistent treatment of the effects of magnetic fields and radiation transfer is presented. Attention is given to the hydrodynamic (HD) algorithms which form the foundation for the more complex MHD and radiation HD algorithms. The effect of self-gravity on the flow dynamics is accounted for by an iterative solution of the sparse-banded matrix resulting from discretizing the Poisson equation in multidimensions. The results of an extensive series of HD test problems are presented. A detailed description of the MHD algorithms in ZEUS-2D is presented. A new method of computing the electromotive force is developed using the method of characteristics (MOC). It is demonstrated through the results of an extensive series of MHD test problems that the resulting hybrid MOC-constrained transport method provides for the accurate evolution of all modes of MHD wave families.
Active transport and cluster formation on 2D networks.
Greulich, P; Santen, L
2010-06-01
We introduce a model for active transport on inhomogeneous networks embedded in a diffusive environment which is motivated by vesicular transport on actin filaments. In the presence of a hard-core interaction, particle clusters are observed that exhibit an algebraically decaying distribution in a large parameter regime, indicating the existence of clusters on all scales. The scale-free behavior can be understood by a mechanism promoting preferential attachment of particles to large clusters. The results are compared with a diffusion-limited aggregation model and active transport on a regular network. For both models we observe aggregation of particles to clusters which are characterized by a finite size scale if the relevant time scales and particle densities are considered. PMID:20556462
CAS2D- NONROTATING BLADE-TO-BLADE, STEADY, POTENTIAL TRANSONIC CASCADE FLOW ANALYSIS CODE
NASA Technical Reports Server (NTRS)
Dulikravich, D. S.
1994-01-01
An exact, full-potential-equation model for the steady, irrotational, homoentropic, and homoenergetic flow of a compressible, inviscid fluid through a two-dimensional planar cascade together with its appropriate boundary conditions has been derived. The CAS2D computer program numerically solves an artificially time-dependent form of the actual full-potential-equation, providing a nonrotating blade-to-blade, steady, potential transonic cascade flow analysis code. Comparisons of results with test data and theoretical solutions indicate very good agreement. In CAS2D, the governing equation is discretized by using type-dependent, rotated finite differencing and the finite area technique. The flow field is discretized by providing a boundary-fitted, nonuniform computational mesh. This mesh is generated by using a sequence of conformal mapping, nonorthogonal coordinate stretching, and local, isoparametric, bilinear mapping functions. The discretized form of the full-potential equation is solved iteratively by using successive line over relaxation. Possible isentropic shocks are captured by the explicit addition of an artificial viscosity in a conservative form. In addition, a four-level, consecutive, mesh refinement feature makes CAS2D a reliable and fast algorithm for the analysis of transonic, two-dimensional cascade flows. The results from CAS2D are not directly applicable to three-dimensional, potential, rotating flows through a cascade of blades because CAS2D does not consider the effects of the Coriolis force that would be present in the three-dimensional case. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 200K of 8 bit bytes. The CAS2D program was developed in 1980.
Laure, Chloé; Karamessini, Denise; Milenkovic, Olgica; Charles, Laurence; Lutz, Jean-François
2016-08-26
A 2D approach was studied for the design of polymer-based molecular barcodes. Uniform oligo(alkoxyamine amide)s, containing a monomer-coded binary message, were synthesized by orthogonal solid-phase chemistry. Sets of oligomers with different chain-lengths were prepared. The physical mixture of these uniform oligomers leads to an intentional dispersity (1st dimension fingerprint), which is measured by electrospray mass spectrometry. Furthermore, the monomer sequence of each component of the mass distribution can be analyzed by tandem mass spectrometry (2nd dimension sequencing). By summing the sequence information of all components, a binary message can be read. A 4-bytes extended ASCII-coded message was written on a set of six uniform oligomers. Alternatively, a 3-bytes sequence was written on a set of five oligomers. In both cases, the coded binary information was recovered. PMID:27484303
Hall-Effect Thruster Simulations with 2-D Electron Transport and Hydrodynamic Ions
NASA Technical Reports Server (NTRS)
Mikellides, Ioannis G.; Katz, Ira; Hofer, Richard H.; Goebel, Dan M.
2009-01-01
A computational approach that has been used extensively in the last two decades for Hall thruster simulations is to solve a diffusion equation and energy conservation law for the electrons in a direction that is perpendicular to the magnetic field, and use discrete-particle methods for the heavy species. This "hybrid" approach has allowed for the capture of bulk plasma phenomena inside these thrusters within reasonable computational times. Regions of the thruster with complex magnetic field arrangements (such as those near eroded walls and magnets) and/or reduced Hall parameter (such as those near the anode and the cathode plume) challenge the validity of the quasi-one-dimensional assumption for the electrons. This paper reports on the development of a computer code that solves numerically the 2-D axisymmetric vector form of Ohm's law, with no assumptions regarding the rate of electron transport in the parallel and perpendicular directions. The numerical challenges related to the large disparity of the transport coefficients in the two directions are met by solving the equations in a computational mesh that is aligned with the magnetic field. The fully-2D approach allows for a large physical domain that extends more than five times the thruster channel length in the axial direction, and encompasses the cathode boundary. Ions are treated as an isothermal, cold (relative to the electrons) fluid, accounting for charge-exchange and multiple-ionization collisions in the momentum equations. A first series of simulations of two Hall thrusters, namely the BPT-4000 and a 6-kW laboratory thruster, quantifies the significance of ion diffusion in the anode region and the importance of the extended physical domain on studies related to the impact of the transport coefficients on the electron flow field.
Screening and transport in 2D semiconductor systems at low temperatures
Das Sarma, S.; Hwang, E. H.
2015-01-01
Low temperature carrier transport properties in 2D semiconductor systems can be theoretically well-understood within RPA-Boltzmann theory as being limited by scattering from screened Coulomb disorder arising from random quenched charged impurities in the environment. In this work, we derive a number of analytical formula, supported by realistic numerical calculations, for the relevant density, mobility, and temperature range where 2D transport should manifest strong intrinsic (i.e., arising purely from electronic effects) metallic temperature dependence in different semiconductor materials arising entirely from the 2D screening properties, thus providing an explanation for why the strong temperature dependence of the 2D resistivity can only be observed in high-quality and low-disorder 2D samples and also why some high-quality 2D materials manifest much weaker metallicity than other materials. We also discuss effects of interaction and disorder on the 2D screening properties in this context as well as compare 2D and 3D screening functions to comment why such a strong intrinsic temperature dependence arising from screening cannot occur in 3D metallic carrier transport. Experimentally verifiable predictions are made about the quantitative magnitude of the maximum possible low-temperature metallicity in 2D systems and the scaling behavior of the temperature scale controlling the quantum to classical crossover. PMID:26572738
Multipacting Simulation Study for 56 MHz Quarter Wave Resonator using 2D Code
Naik,D.; Ben-Zvi, I.
2009-01-02
A beam excited 56 MHz Radio Frequency (RF) Niobium Quarter Wave Resonator (QWR) has been proposed to enhance RHIC beam luminosity and bunching. Being a RF cavity, multipacting is expected; therefore an extensive study was carried out with the Multipac 2.1 2D simulation code. The study revealed that multipacting occurs in various bands up to peak surface electric field 50 kV/m and is concentrated mostly above the beam gap and on the outer conductor. To suppress multipacting, a ripple structure was introduced to the outer conductor and the phenomenon was successfully eliminated from the cavity.
Automated design of coupled RF cavities using 2-D and 3-D codes
Smith, Peter; Christiansen, D. W.; Greninger, P. T.; Spalek, G.
2001-01-01
Coupled RF cavities in the Accelerator Production of Tritium Project have been designed using a procedure in which a 2-D code (CCT) searches for a design that meets frequency and coupling requirements, while a 3-D code (HFSS) is used to obtain empirical factors used by CCT to characterize the coupling slot between cavities. Using assumed values of the empirical factors, CCT runs the Superfish code iteratively to solve for a trial cavity design that has a specified frequency and coupling. The frequency shifts and the coupling constant k of the slot are modeled in CCT using a perturbation theory, the results of which are adjusted using the empirical factors. Given a trial design, HFSS is run using periodic boundary conditions to obtain a mode spectrum. The mode spectrum is processed using the DISPER code to obtain values of the coupling and the frequencies with slots. These results are used to calculate a new set of empirical factors, which are fed back into CCT for another design iteration. Cold models have been fabricated and tested to validate the codes, and results will be presented.
A sequential partly iterative approach for multicomponent reactive transport with CORE2D
Samper, J.; Xu, T.; Yang, C.
2008-11-01
CORE{sup 2D} V4 is a finite element code for modeling partly or fully saturated water flow, heat transport and multicomponent reactive solute transport under both local chemical equilibrium and kinetic conditions. It can handle coupled microbial processes and geochemical reactions such as acid-base, aqueous complexation, redox, mineral dissolution/precipitation, gas dissolution/exsolution, ion exchange, sorption via linear and nonlinear isotherms, sorption via surface complexation. Hydraulic parameters may change due to mineral precipitation/dissolution reactions. Coupled transport and chemical equations are solved by using sequential iterative approaches. A sequential partly-iterative approach (SPIA) is presented which improves the accuracy of the traditional sequential noniterative approach (SNIA) and is more efficient than the general sequential iterative approach (SIA). While SNIA leads to a substantial saving of computing time, it introduces numerical errors which are especially large for cation exchange reactions. SPIA improves the efficiency of SIA because the iteration between transport and chemical equations is only performed in nodes with a large mass transfer between solid and liquid phases. The efficiency and accuracy of SPIA are compared to those of SIA and SNIA using synthetic examples and a case study of reactive transport through the Llobregat Delta aquitard in Spain. SPIA is found to be as accurate as SIA while requiring significantly less CPU time. In addition, SPIA is much more accurate than SNIA with only a minor increase in computing time. A further enhancement of the efficiency of SPIA is achieved by improving the efficiency of the Newton-Raphson method used for solving chemical equations. Such an improvement is obtained by working with increments of log-concentrations and ignoring the terms of the Jacobian matrix containing derivatives of activity coefficients. A proof is given for the symmetry and non-singularity of the Jacobian matrix
NASA Astrophysics Data System (ADS)
Stone, James M.; Norman, Michael L.
1992-06-01
In this, the second of a series of three papers, we continue a detailed description of ZEUS-2D, a numerical code for the simulation of fluid dynamical flows in astrophysics including a self-consistent treatment of the effects of magnetic fields and radiation transfer. In this paper, we give a detailed description of the magnetohydrodynamical (MHD) algorithms in ZEUS-2D. The recently developed constrained transport (CT) algorithm is implemented for the numerical evolution of the components of the magnetic field for MHD simulations. This formalism guarantees the numerically evolved field components will satisfy the divergence-free constraint at all times. We find, however, that the method used to compute the electromotive forces must be chosen carefully to propagate accurately all modes of MHD wave families (in particular shear Alfvén waves). A new method of computing the electromotive force is developed using the method of characteristics (MOC). It is demonstrated through the results of an extensive series of MHD test problems that the resulting hybrid MOC-CT method provides for the accurate evolution of all modes of MHD wave families.
Modelling 2001 lahars at Popocatépetl volcano using FLO2D numerical code
NASA Astrophysics Data System (ADS)
Caballero, L.; Capra, L.
2013-12-01
Popocatépetl volcano is located on the central part of the Transmexican Volcanic Belt. It is one of the most active volcanoes in Mexico and endanger more than 25 million people that lives in its surroundings. In the last months, the renewal of its volcanic activity put into alert scientific community. One of the possible scenarios is the 2001 explosive activity, which was characterized by a 8 km eruptive column and the subsequent formation of pumice flows up to 4 km from the crater. Lahars were generated few hours after, remobilizing the new deposits towards NE flank of the volcano, along Huiloac Gorge, almost reaching Santiago Xalitzintla town (Capra et al., 2004). The occurrence of a similar scenario makes very important to reproduce this event to delimitate accurately lahar hazard zones. In this work, 2001 lahar deposit is modeled using FLO2D numerical code. Geophone data is used to reconstruct initial hydrograph and sediment concentration. Sensitivity study of most important parameters used by this code like Manning, and α and β coefficients was conducted in order to achieve a good simulation. Results obtained were compared with field data and demonstrated a good agreement in thickness and flow distribution. A comparison with previously published data with laharZ program (Muñoz-Salinas, 2009) is also made. Additionally, lahars with fluctuating sediment concentrations but with similar volume are simulated to observe the influence of the rheological behavior on lahar distribution.
Embedded morphological dilation coding for 2D and 3D images
NASA Astrophysics Data System (ADS)
Lazzaroni, Fabio; Signoroni, Alberto; Leonardi, Riccardo
2002-01-01
Current wavelet-based image coders obtain high performance thanks to the identification and the exploitation of the statistical properties of natural images in the transformed domain. Zerotree-based algorithms, as Embedded Zerotree Wavelets (EZW) and Set Partitioning In Hierarchical Trees (SPIHT), offer high Rate-Distortion (RD) coding performance and low computational complexity by exploiting statistical dependencies among insignificant coefficients on hierarchical subband structures. Another possible approach tries to predict the clusters of significant coefficients by means of some form of morphological dilation. An example of a morphology-based coder is the Significance-Linked Connected Component Analysis (SLCCA) that has shown performance which are comparable to the zerotree-based coders but is not embedded. A new embedded bit-plane coder is proposed here based on morphological dilation of significant coefficients and context based arithmetic coding. The algorithm is able to exploit both intra-band and inter-band statistical dependencies among wavelet significant coefficients. Moreover, the same approach is used both for two and three-dimensional wavelet-based image compression. Finally we the algorithms are tested on some 2D images and on a medical volume, by comparing the RD results to those obtained with the state-of-the-art wavelet-based coders.
Parallel algorithms for 2-D cylindrical transport equations of Eigenvalue problem
Wei, J.; Yang, S.
2013-07-01
In this paper, aimed at the neutron transport equations of eigenvalue problem under 2-D cylindrical geometry on unstructured grid, the discrete scheme of Sn discrete ordinate and discontinuous finite is built, and the parallel computation for the scheme is realized on MPI systems. Numerical experiments indicate that the designed parallel algorithm can reach perfect speedup, it has good practicality and scalability. (authors)
Acceleration of a Monte Carlo radiation transport code
Hochstedler, R.D.; Smith, L.M.
1996-03-01
Execution time for the Integrated TIGER Series (ITS) Monte Carlo radiation transport code has been reduced by careful re-coding of computationally intensive subroutines. Three test cases for the TIGER (1-D slab geometry), CYLTRAN (2-D cylindrical geometry), and ACCEPT (3-D arbitrary geometry) codes were identified and used to benchmark and profile program execution. Based upon these results, sixteen top time-consuming subroutines were examined and nine of them modified to accelerate computations with equivalent numerical output to the original. The results obtained via this study indicate that speedup factors of 1.90 for the TIGER code, 1.67 for the CYLTRAN code, and 1.11 for the ACCEPT code are achievable. {copyright} {ital 1996 American Institute of Physics.}
Differential Sensitivity Theory applied to the MESA2D code for multi-material problems
Henninger, R.J.; Maudlin, P.J.; Harstad, E.N.
1996-05-01
The technique called Differential Sensitivity Theory (DST) is extended to the multi-component system of equations solved by the MESA2D hydrocode. DST uses adjoint techniques to determine exact sensitivity derivatives, i.e., if R is a calculation result of interest (response R) and {alpha}{sub i} is a calculation input (parameter {alpha}{sub i}), then {partial_derivative}R/{partial_derivative}{alpha}{sub i} is defined as the sensitivity. The advantage of using DST is that for an n-parameter problem {ital all} n sensitivities can be obtained by integrating the solutions from only {ital two} calculations, a MESA calculation and its corresponding adjoint calculation using an Adjoint Continuum Mechanics (ACM) code. Previous papers have described application of the technique to one-dimensional, single-material problems. This work presents the derivation and solution of the additional adjoint equations for the purpose of computing sensitivities for two-dimensional, multi-component problems. As an example, results for a multi-material flyer plate impact problem featuring an oblique impact are given. {copyright} {ital 1996 American Institute of Physics.}
Differential Sensitivity Theory applied to the MESA2D code for multi-material problems
NASA Astrophysics Data System (ADS)
Henninger, R. J.; Maudlin, P. J.; Harstad, E. N.
1996-05-01
The technique called Differential Sensitivity Theory (DST) is extended to the multi-component system of equations solved by the MESA2D hydrocode. DST uses adjoint techniques to determine exact sensitivity derivatives, i.e., if R is a calculation result of interest (response R) and αi is a calculation input (parameter αi), then ∂R/∂αi is defined as the sensitivity. The advantage of using DST is that for an n-parameter problem all n sensitivities can be obtained by integrating the solutions from only two calculations, a MESA calculation and its corresponding adjoint calculation using an Adjoint Continuum Mechanics (ACM) code. Previous papers have described application of the technique to one-dimensional, single-material problems. This work presents the derivation and solution of the additional adjoint equations for the purpose of computing sensitivities for two-dimensional, multi-component problems. As an example, results for a multi-material flyer plate impact problem featuring an oblique impact are given.
Differential sensitivity theory applied to the MESA2D code for multi-material problems
Henninger, R.J.; Maudlin, P.J.; Harstad, E.N.
1995-09-01
The technique called Differential Sensitivity Theory (DST) is extended to the multi-component system of equations solved by the MESA2D hydrocode. DST uses adjoint techniques to determine exact sensitivity derivatives, i.e., if R is a calculation result of interest (response R) and {alpha}{sub i} is a calculation input (parameter {alpha}{sub i}), then {partial_derivative}R/{partial_derivative}{alpha}{sub i} is defined as the sensitivity. The advantage of using DST is that for an n-parameter problem all n sensitivities can be obtained by integrating the solutions from only two calculations, a MESA calculation and its corresponding adjoint calculation using an Adjoint Continuum Mechanics (ACM) code. Previous papers have described application of the technique to one-dimensional, single-material problems. This work presents the derivation and solution of the additional adjoint equations for the purpose of computing sensitivities for two-dimensional, multi-component problems. As an example, results for a multi-material flyer plate impact problem featuring an oblique impact are given.
A 2D/1D coupling neutron transport method based on the matrix MOC and NEM methods
Zhang, H.; Zheng, Y.; Wu, H.; Cao, L.
2013-07-01
A new 2D/1D coupling method based on the matrix MOC method (MMOC) and nodal expansion method (NEM) is proposed for solving the three-dimensional heterogeneous neutron transport problem. The MMOC method, used for radial two-dimensional calculation, constructs a response matrix between source and flux with only one sweep and then solves the linear system by using the restarted GMRES algorithm instead of the traditional trajectory sweeping process during within-group iteration for angular flux update. Long characteristics are generated by using the customization of commercial software AutoCAD. A one-dimensional diffusion calculation is carried out in the axial direction by employing the NEM method. The 2D and ID solutions are coupled through the transverse leakage items. The 3D CMFD method is used to ensure the global neutron balance and adjust the different convergence properties of the radial and axial solvers. A computational code is developed based on these theories. Two benchmarks are calculated to verify the coupling method and the code. It is observed that the corresponding numerical results agree well with references, which indicates that the new method is capable of solving the 3D heterogeneous neutron transport problem directly. (authors)
A novel 2D wavelength-time chaos code in optical CDMA system
NASA Astrophysics Data System (ADS)
Zhang, Qi; Xin, Xiangjun; Wang, Yongjun; Zhang, Lijia; Yu, Chongxiu; Meng, Nan; Wang, Houtian
2012-11-01
Two-dimensional wavelength-time chaos code is proposed and constructed for a synchronous optical code division multiple access system. The access performance is compared between one-dimensional chaos code, WDM/chaos code and the proposed code. Comparison shows that two-dimensional wavelength-time chaos code possesses larger capacity, better spectral efficiency and bit-error ratio than WDM/chaos combinations and one-dimensional chaos code.
Effects of magnetic impurities on transport in 2D topological insulators
NASA Astrophysics Data System (ADS)
Dang, Xiaoqian; Burton, J. D.; Tsymbal, Evgeny
Understanding the transport properties of topological insulators could bring such materials from fundamental research to potential applications. Here we report on the theoretical investigations of the effects of magnetic impurities on transport properties of model two-dimensional (2D) topological insulators (TIs). We utilize the tight-binding form of the Bernevig-Hughes-Zhang model and investigate the transport properties by employing the Landauer-Büttiker formalism. We explore the current distribution in 2D TIs resulting from scattering by a magnetic impurity which breaks time-reversal symmetry. We find that a magnetic impurity could drive anti-resonant behavior of the conductance, as revealed from full backscattering of the electron current flowing at one of the edges of the TI. This phenomenon occurs due to spin-flip scattering when the Fermi energy matches the impurity state and the magnetic moment of the impurity is aligned along the TI edge. Additionally, we explore the effect of an external magnetic gate attached to the system and show that changing the magnetization orientation within the gate allows the control of conductance. This geometric setup could be realized experimentally providing the opportunity to tune transport properties of 2D TIs by a magnetic gate.
2D/1D approximations to the 3D neutron transport equation. I: Theory
Kelley, B. W.; Larsen, E. W.
2013-07-01
A new class of '2D/1D' approximations is proposed for the 3D linear Boltzmann equation. These approximate equations preserve the exact transport physics in the radial directions x and y and diffusion physics in the axial direction z. Thus, the 2D/1D equations are more accurate approximations of the 3D Boltzmann equation than the conventional 3D diffusion equation. The 2D/1D equations can be systematically discretized, to yield accurate simulation methods for 3D reactor core problems. The resulting solutions will be more accurate than 3D diffusion solutions, and less expensive to generate than standard 3D transport solutions. In this paper, we (i) show that the simplest 2D/1D equation has certain desirable properties, (ii) systematically discretize this equation, and (iii) derive a stable iteration scheme for solving the discrete system of equations. In a companion paper [1], we give numerical results that confirm the theoretical predictions of accuracy and iterative stability. (authors)
Ren, X; Domier, C W; Kramer, G; Luhmann, N C; Muscatello, C M; Shi, L; Tobias, B J; Valeo, E
2014-11-01
A synthetic microwave imaging reflectometer (MIR) diagnostic employing the full-wave reflectometer code (FWR2D) has been developed and is currently being used to guide the design of real systems, such as the one recently installed on DIII-D. The FWR2D code utilizes real plasma profiles as input, and it is combined with optical simulation tools for synthetic diagnostic signal generation. A detailed discussion of FWR2D and the process to generate the synthetic signal are presented in this paper. The synthetic signal is also compared to a prescribed density fluctuation spectrum to quantify the imaging quality. An example is presented with H-mode-like plasma profiles derived from a DIII-D discharge, where the MIR focal is located in the pedestal region. It is shown that MIR is suitable for diagnosing fluctuations with poloidal wavenumber up to 2.0 cm(-1) and fluctuation amplitudes less than 5%. PMID:25430276
Impact of high speed civil transports on stratospheric ozone: A 2-D model investigation
Kinnison, D.E.; Connell, P.S.
1996-12-01
This study investigates the effect on stratospheric ozone from a fleet of proposed High Speed Civil Transports (HSCTs). The new LLNL 2-D operator-split chemical-radiative-transport model of the troposphere and stratosphere is used for this HSCT investigation. This model is integrated in a diurnal manner, using an implicit numerical solver. Therefore, rate coefficients are not modified by any sort of diurnal average factor. This model also does not make any assumptions on lumping of chemical species into families. Comparisons to previous model-derived HSCT assessment of ozone change are made, both to the previous LLNL 2-D model and to other models from the international assessment modeling community. The sensitivity to the NO{sub x} emission index and sulfate surface area density is also explored.
2D/1D approximations to the 3D neutron transport equation. II: Numerical comparisons
Kelley, B. W.; Collins, B.; Larsen, E. W.
2013-07-01
In a companion paper [1], (i) several new '2D/1D equations' are introduced as accurate approximations to the 3D Boltzmann transport equation, (ii) the simplest of these approximate equations is systematically discretized, and (iii) a theoretically stable iteration scheme is developed to solve the discrete equations. In this paper, numerical results are presented that confirm the theoretical predictions made in [1]. (authors)
Modeling anomalous radial transport in kinetic transport codes
NASA Astrophysics Data System (ADS)
Bodi, K.; Krasheninnikov, S. I.; Cohen, R. H.; Rognlien, T. D.
2009-11-01
Anomalous transport is typically the dominant component of the radial transport in magnetically confined plasmas, where the physical origin of this transport is believed to be plasma turbulence. A model is presented for anomalous transport that can be used in continuum kinetic edge codes like TEMPEST, NEO and the next-generation code being developed by the Edge Simulation Laboratory. The model can also be adapted to particle-based codes. It is demonstrated that the model with a velocity-dependent diffusion and convection terms can match a diagonal gradient-driven transport matrix as found in contemporary fluid codes, but can also include off-diagonal effects. The anomalous transport model is also combined with particle drifts and a particle/energy-conserving Krook collision operator to study possible synergistic effects with neoclassical transport. For the latter study, a velocity-independent anomalous diffusion coefficient is used to mimic the effect of long-wavelength ExB turbulence.
In-facility transport code review
Spore, J.W.; Boyack, B.E.; Bohl, W.R.
1996-07-01
The following computer codes were reviewed by the In-Facility Transport Working Group for application to the in-facility transport of radioactive aerosols, flammable gases, and/or toxic gases: (1) CONTAIN, (2) FIRAC, (3) GASFLOW, (4) KBERT, and (5) MELCOR. Based on the review criteria as described in this report and the versions of each code available at the time of the review, MELCOR is the best code for the analysis of in-facility transport when multidimensional effects are not significant. When multi-dimensional effects are significant, GASFLOW should be used.
Development of models for the two-dimensional, two-fluid code for sodium boiling NATOF-2D. [LMFBR
Zielinski, R.G.; Kazimi, M.S.
1981-09-01
Several features were incorporated into NATOF-2D, a two-dimensional, two fluid code developed at MIT for the purpose of analysis of sodium boiling transients under LMFBR conditions. They include improved interfacial mass, momentum and energy exchange rate models, and a cell-to-cell radial heat conduction mechanism which was calibrated by simulation of Westinghouse Blanket Heat Transfer Test Program Runs 544 and 545. Finally, a direct method of pressure field solution was implemented into a direct method of pressure field solution was implemented into NATOF-2D, replacing the iterative technique previously available, and resulted in substantially reduced computational costs.
Monte Carlo Ion Transport Analysis Code.
Energy Science and Technology Software Center (ESTSC)
2009-04-15
Version: 00 TRIPOS is a versatile Monte Carlo ion transport analysis code. It has been applied to the treatment of both surface and bulk radiation effects. The media considered is composed of multilayer polyatomic materials.
A 2-D orientation-adaptive prediction filter in lifting structures for image coding.
Gerek, Omer N; Cetin, A Enis
2006-01-01
Lifting-style implementations of wavelets are widely used in image coders. A two-dimensional (2-D) edge adaptive lifting structure, which is similar to Daubechies 5/3 wavelet, is presented. The 2-D prediction filter predicts the value of the next polyphase component according to an edge orientation estimator of the image. Consequently, the prediction domain is allowed to rotate +/-45 degrees in regions with diagonal gradient. The gradient estimator is computationally inexpensive with additional costs of only six subtractions per lifting instruction, and no multiplications are required. PMID:16435541
Monte Carlo Nucleon Meson Transport Code System.
Energy Science and Technology Software Center (ESTSC)
2000-11-17
Version 00 NMTC/JAERI97 is an upgraded version of the code system NMTC/JAERI, which was developed in 1982 at JAERI and is based on the CCC-161/NMTC code system. NMTC/JAERI97 simulates high energy nuclear reactions and nucleon-meson transport processes.
NASA Astrophysics Data System (ADS)
Zhao, Yaqin; Zhong, Xin; Wu, Di; Zhang, Ye; Ren, Guanghui; Wu, Zhilu
2013-09-01
Optical code-division multiple access (OCDMA) systems usually allocate orthogonal or quasi-orthogonal codes to the active users. When transmitting through atmospheric scattering channel, the coding pulses are broadened and the orthogonality of the codes is worsened. In truly asynchronous case, namely both the chips and the bits are asynchronous among each active user, the pulse broadening affects the system performance a lot. In this paper, we evaluate the performance of a 2D asynchronous hard-limiting wireless OCDMA system through atmospheric scattering channel. The probability density function of multiple access interference in truly asynchronous case is given. The bit error rate decreases as the ratio of the chip period to the root mean square delay spread increases and the channel limits the bit rate to different levels when the chip period varies.
STATUS OF THE MCNPX TRANSPORT CODE
Hughes, H.G.; Chadwick, M.B.
2000-10-01
The Monte Carlo particle transport code MCNPX and its associated data have been the focus of a major development effort at Los Alamos for several years. The system has reached a mature state, and has become a significant tool for many intermediate and high-energy particle transport applications. A recent version has been released to the Radiation Safety Information Computational Center (RSICC). A recent report provides an overview of the code and an extensive set of references for the component physics modules used in the code. In this paper we review the status of the developmental version of MCNPX, and describe some important new enhancements, including the use of evaluated nuclear data files for proton transport; the use of photonuclear reaction data; improved elastic and inelastic react ion cross sections for nucleons, antinucleons, pions, and kaons; and two new modes of operation of the code. We also illustrate the use of the new proton and photonuclear data in two representative applications.
Transport studies in 2D transition metal dichalcogenides and black phosphorus.
Du, Yuchen; Neal, Adam T; Zhou, Hong; Ye, Peide D
2016-07-01
Two-dimensional (2D) materials are a new family of materials with interesting physical properties, ranging from insulating hexagonal boron nitride, semiconducting or semi-metallic transition metal dichalcogenides, to gapless metallic graphene. In this review, we provide a brief discussion of transport studies in transition metal dichalcogenides, including both semiconducting and semi-metallic phases, as well as a discussion of the newly emerged narrow bandgap layered material, black phosphorus, in terms of its electrical and quantum transport properties at room and cryogenic temperatures. Ultra-thin layered channel materials with atomic layer thickness in the cross-plane direction, together with relatively high carrier mobility with appropriate passivation techniques, provide the promise for new scientific discoveries and broad device applications. PMID:27187790
Transport studies in 2D transition metal dichalcogenides and black phosphorus
NASA Astrophysics Data System (ADS)
Du, Yuchen; Neal, Adam T.; Zhou, Hong; Ye, Peide D.
2016-07-01
Two-dimensional (2D) materials are a new family of materials with interesting physical properties, ranging from insulating hexagonal boron nitride, semiconducting or semi-metallic transition metal dichalcogenides, to gapless metallic graphene. In this review, we provide a brief discussion of transport studies in transition metal dichalcogenides, including both semiconducting and semi-metallic phases, as well as a discussion of the newly emerged narrow bandgap layered material, black phosphorus, in terms of its electrical and quantum transport properties at room and cryogenic temperatures. Ultra-thin layered channel materials with atomic layer thickness in the cross-plane direction, together with relatively high carrier mobility with appropriate passivation techniques, provide the promise for new scientific discoveries and broad device applications.
Implementation of an anomalous radial transport model for continuum kinetic edge codes
NASA Astrophysics Data System (ADS)
Bodi, K.; Krasheninnikov, S. I.; Cohen, R. H.; Rognlien, T. D.
2007-11-01
Radial plasma transport in magnetic fusion devices is often dominated by plasma turbulence compared to neoclassical collisional transport. Continuum kinetic edge codes [such as the (2d,2v) transport version of TEMPEST and also EGK] compute the collisional transport directly, but there is a need to model the anomalous transport from turbulence for long-time transport simulations. Such a model is presented and results are shown for its implementation in the TEMPEST gyrokinetic edge code. The model includes velocity-dependent convection and diffusion coefficients expressed as a Hermite polynominals in velocity. The specification of the Hermite coefficients can be set, e.g., by specifying the ratio of particle and energy transport as in fluid transport codes. The anomalous transport terms preserve the property of no particle flux into unphysical regions of velocity space. TEMPEST simulations are presented showing the separate control of particle and energy anomalous transport, and comparisons are made with neoclassical transport also included.
Momentum Transport: 2D and 3D Cloud Resolving Model Simulations
NASA Technical Reports Server (NTRS)
Tao, Wei-Kuo
2001-01-01
The major objective of this study is to investigate the momentum budgets associated with several convective systems that developed during the TOGA COARE IOP (west Pacific warm pool region) and GATE (east Atlantic region). The tool for this study is the improved Goddard Cumulas Ensemble (GCE) model which includes a 3-class ice-phase microphysical scheme, explicit cloud radiative interactive processes and air-sea interactive surface processes. The model domain contains 256 x 256 grid points (with 2 km resolution) in the horizontal and 38 grid points (to a depth of 22 km) in the vertical. The 2D domain has 1024 grid points. The simulations were performed over a 7-day time period (December 19-26, 1992, for TOGA COARE and September 1-7, 1994 for GATE). Cyclic literal boundary conditions are required for this type of long-term integration. Two well organized squall systems (TOGA, COARE February 22, 1993, and GATE September 12, 1994) were also simulated using the 3D GCE model. Only 9 h simulations were required to cover the life time of the squall systems. the lateral boundary conditions were open for these two squall systems simulations. the following will be examined: (1) the momentum budgets in the convective and stratiform regions, (2) the relationship between momentum transport and cloud organization (i.e., well organized squall lines versus less organized convective), (3) the differences and similarities in momentum transport between 2D and 3D simulated convective systems, and (4) the differences and similarities in momentum budgets between cloud systems simulated with open and cyclic lateral boundary conditions. Preliminary results indicate that there are only small differences between 2D and 3D simulated momentum budgets. Major differences occur, however, between momentum budgets associated with squall systems simulated using different lateral boundary conditions.
DANTSYS: A diffusion accelerated neutral particle transport code system
Alcouffe, R.E.; Baker, R.S.; Brinkley, F.W.; Marr, D.R.; O`Dell, R.D.; Walters, W.F.
1995-06-01
The DANTSYS code package includes the following transport codes: ONEDANT, TWODANT, TWODANT/GQ, TWOHEX, and THREEDANT. The DANTSYS code package is a modular computer program package designed to solve the time-independent, multigroup discrete ordinates form of the boltzmann transport equation in several different geometries. The modular construction of the package separates the input processing, the transport equation solving, and the post processing (or edit) functions into distinct code modules: the Input Module, one or more Solver Modules, and the Edit Module, respectively. The Input and Edit Modules are very general in nature and are common to all the Solver Modules. The ONEDANT Solver Module contains a one-dimensional (slab, cylinder, and sphere), time-independent transport equation solver using the standard diamond-differencing method for space/angle discretization. Also included in the package are solver Modules named TWODANT, TWODANT/GQ, THREEDANT, and TWOHEX. The TWODANT Solver Module solves the time-independent two-dimensional transport equation using the diamond-differencing method for space/angle discretization. The authors have also introduced an adaptive weighted diamond differencing (AWDD) method for the spatial and angular discretization into TWODANT as an option. The TWOHEX Solver Module solves the time-independent two-dimensional transport equation on an equilateral triangle spatial mesh. The THREEDANT Solver Module solves the time independent, three-dimensional transport equation for XYZ and RZ{Theta} symmetries using both diamond differencing with set-to-zero fixup and the AWDD method. The TWODANT/GQ Solver Module solves the 2-D transport equation in XY and RZ symmetries using a spatial mesh of arbitrary quadrilaterals. The spatial differencing method is based upon the diamond differencing method with set-to-zero fixup with changes to accommodate the generalized spatial meshing.
Parallel Finite Element Electron-Photon Transport Analysis on 2-D Unstructured Mesh
Drumm, C.R.
1999-01-01
A computer code has been developed to solve the linear Boltzmann transport equation on an unstructured mesh of triangles, from a Pro/E model. An arbitriwy arrangement of distinct material regions is allowed. Energy dependence is handled by solving over an arbitrary number of discrete energy groups. Angular de- pendence is treated by Legendre-polynomial expansion of the particle cross sections and a discrete ordinates treatment of the particle fluence. The resulting linear system is solved in parallel with a preconditioned conjugate-gradients method. The solution method is unique, in that the space-angle dependence is solved si- multaneously, eliminating the need for the usual inner iterations. Electron cross sections are obtained from a Goudsrnit-Saunderson modifed version of the CEPXS code. A one-dimensional version of the code has also been develop@ for testing and development purposes.
Incorporating a Turbulence Transport Model into 2-D Hybrid Hall Thruster Simulations
NASA Astrophysics Data System (ADS)
Cha, Eunsun; Cappelli, Mark A.; Fernandez, Eduardo
2014-10-01
2-D hybrid simulations of Hall plasma thrusters that do not resolve cross-field transport-generating fluctuations require a model to capture how electrons migrate across the magnetic field. We describe the results of integrating a turbulent electron transport model into simulations of plasma behavior in a plane spanned by the E and B field vectors. The simulations treat the electrons as a fluid and the heavy species (ions/neutrals) as discrete particles. The transport model assumes that the turbulent eddy cascade in the electron fluid to smaller scales is the primary means of electron energy dissipation. Using this model, we compare simulations to experimental measurements made on a laboratory Hall discharge over a range of discharge voltage. Both the current-voltage trends as well as the plasma properties such as plasma temperature, electron density, and ion velocities seem agree favorably with experiments, where a simple Bohm transport model tends to perform poorly in capturing much of the discharge behavior.
NASA Astrophysics Data System (ADS)
Krause, M.; Camenzind, M.
2001-12-01
In the present paper, we examine the convergence behavior and inter-code reliability of astrophysical jet simulations in axial symmetry. We consider both pure hydrodynamic jets and jets with a dynamically significant magnetic field. The setups were chosen to match the setups of two other publications, and recomputed with the MHD code NIRVANA. We show that NIRVANA and the two other codes give comparable, but not identical results. We explain the differences by the different application of artificial viscosity in the three codes and numerical details, which can be summarized in a resolution effect, in the case without magnetic field: NIRVANA turns out to be a fair code of medium efficiency. It needs approximately twice the resolution as the code by Lind (Lind et al. 1989) and half the resolution as the code by Kössl (Kössl & Müller 1988). We find that some global properties of a hydrodynamical jet simulation, like e.g. the bow shock velocity, converge at 100 points per beam radius (ppb) with NIRVANA. The situation is quite different after switching on the toroidal magnetic field: in this case, global properties converge even at 10 ppb. In both cases, details of the inner jet structure and especially the terminal shock region are still insufficiently resolved, even at our highest resolution of 70 ppb in the magnetized case and 400 ppb for the pure hydrodynamic jet. The magnetized jet even suffers from a fatal retreat of the Mach disk towards the inflow boundary, which indicates that this simulation does not converge, in the end. This is also in definite disagreement with earlier simulations, and challenges further studies of the problem with other codes. In the case of our highest resolution simulation, we can report two new features: first, small scale Kelvin-Helmholtz instabilities are excited at the contact discontinuity next to the jet head. This slows down the development of the long wavelength Kelvin-Helmholtz instability and its turbulent cascade to smaller
Simulation and analysis of solute transport in 2D fracture/pipe networks: The SOLFRAC program
NASA Astrophysics Data System (ADS)
Bodin, Jacques; Porel, Gilles; Delay, Fred; Ubertosi, Fabrice; Bernard, Stéphane; de Dreuzy, Jean-Raynald
2007-01-01
The Time Domain Random Walk (TDRW) method has been recently developed by Delay and Bodin [Delay, F. and Bodin, J., 2001. Time domain random walk method to simulate transport by advection-dispersion and matrix diffusion in fracture networks. Geophys. Res. Lett., 28(21): 4051-4054.] and Bodin et al. [Bodin, J., Porel, G. and Delay, F., 2003c. Simulation of solute transport in discrete fracture networks using the time domain random walk method. Earth Planet. Sci. Lett., 6566: 1-8.] for simulating solute transport in discrete fracture networks. It is assumed that the fracture network can reasonably be represented by a network of interconnected one-dimensional pipes (i.e. flow channels). Processes accounted for are: (1) advection and hydrodynamic dispersion in the channels, (2) matrix diffusion, (3) diffusion into stagnant zones within the fracture planes, (4) sorption reactions onto the fracture walls and in the matrix, (5) linear decay, and (6) mass sharing at fracture intersections. The TDRW method is handy and very efficient in terms of computation costs since it allows for the one-step calculation of the particle residence time in each bond of the network. This method has been programmed in C++, and efforts have been made to develop an efficient and user-friendly software, called SOLFRAC. This program is freely downloadable at the URL http://labo.univ-poitiers.fr/hydrasa/intranet/telechargement.htm. It calculates solute transport into 2D pipe networks, while considering different types of injections and different concepts of local dispersion within each flow channel. Post-simulation analyses are also available, such as the mean velocity or the macroscopic dispersion at the scale of the entire network. The program may be used to evaluate how a given transport mechanism influences the macroscopic transport behaviour of fracture networks. It may also be used, as is the case, e.g., with analytical solutions, to interpret laboratory or field tracer test experiments
50 CFR Table 2d to Part 679 - Species Codes-Non-FMP Species
Code of Federal Regulations, 2011 CFR
2011-10-01
... (unspecified) 625 Lamprey, pacific 600 Lingcod 130 Lumpsucker 216 Pacific flatnose 260 Pacific hagfish 212 Pacific hake 112 Pacific lamprey 600 Pacific saury 220 Pacific tomcod 250 Poacher (Family Algonidae) 219... description Code Arctic char, anadromous 521 Dolly varden, anadromous 531 Eels or eel-like fish 210 Eel,...
TOPAZ - a finite element heat conduction code for analyzing 2-D solids
Shapiro, A.B.
1984-03-01
TOPAZ is a two-dimensional implicit finite element computer code for heat conduction analysis. This report provides a user's manual for TOPAZ and a description of the numerical algorithms used. Sample problems with analytical solutions are presented. TOPAZ has been implemented on the CRAY and VAX computers.
2D simulation of transport and degradation in the River Rhine.
Teichmann, L; Reuschenbach, P; Müller, B; Horn, H
2002-01-01
A simple 2D model has been developed for the simulation of mass transport and degradation of substances in the river Rhine. The model describes mass transport in the flow direction with a convective and a dispersive term. Transversal transport is described by segmenting the river and formulating a transversal exchange coefficient between the segments. Degradation can be formulated with any kinetics from first order to complex enzyme kinetics. The model was verified with monitoring data from the river Rhine. The hydrodynamic parameters such as dispersion coefficients and exchange coefficients were fitted to the conductivity, which was assumed to be non-degradable. The degradation term was fitted to ammonia values. The model was used to simulate measured concentrations of a readily (Aniline) and a poorly biodegradable substance (1,4-Dioxan) 10 m from the left river bank. It was the objective of this research program to develop a model which allows a realistic estimation of the locally and regionally predicted environmental concentration of chemical substances in the EU risk assessment scheme. PMID:12380980
Coupling 2-D cylindrical and 3-D x-y-z transport computations
Abu-Shumays, I.K.; Yehnert, C.E.; Pitcairn, T.N.
1998-06-30
This paper describes a new two-dimensional (2-D) cylindrical geometry to three-dimensional (3-D) rectangular x-y-z splice option for multi-dimensional discrete ordinates solutions to the neutron (photon) transport equation. Of particular interest are the simple transformations developed and applied in order to carry out the required spatial and angular interpolations. The spatial interpolations are linear and equivalent to those applied elsewhere. The angular interpolations are based on a high order spherical harmonics representation of the angular flux. Advantages of the current angular interpolations over previous work are discussed. An application to an intricate streaming problem is provided to demonstrate the advantages of the new method for efficient and accurate prediction of particle behavior in complex geometries.
2-D Circulation Control Airfoil Benchmark Experiments Intended for CFD Code Validation
NASA Technical Reports Server (NTRS)
Englar, Robert J.; Jones, Gregory S.; Allan, Brian G.; Lin, Johb C.
2009-01-01
A current NASA Research Announcement (NRA) project being conducted by Georgia Tech Research Institute (GTRI) personnel and NASA collaborators includes the development of Circulation Control (CC) blown airfoils to improve subsonic aircraft high-lift and cruise performance. The emphasis of this program is the development of CC active flow control concepts for both high-lift augmentation, drag control, and cruise efficiency. A collaboration in this project includes work by NASA research engineers, whereas CFD validation and flow physics experimental research are part of NASA s systematic approach to developing design and optimization tools for CC applications to fixed-wing aircraft. The design space for CESTOL type aircraft is focusing on geometries that depend on advanced flow control technologies that include Circulation Control aerodynamics. The ability to consistently predict advanced aircraft performance requires improvements in design tools to include these advanced concepts. Validation of these tools will be based on experimental methods applied to complex flows that go beyond conventional aircraft modeling techniques. This paper focuses on recent/ongoing benchmark high-lift experiments and CFD efforts intended to provide 2-D CFD validation data sets related to NASA s Cruise Efficient Short Take Off and Landing (CESTOL) study. Both the experimental data and related CFD predictions are discussed.
Efficient simulation of pitch angle collisions in a 2+2-D Eulerian Vlasov code
NASA Astrophysics Data System (ADS)
Banks, Jeff; Berger, R.; Brunner, S.; Tran, T.
2014-10-01
Here we discuss pitch angle scattering collisions in the context of the Eulerian-based kinetic code LOKI that evolves the Vlasov-Poisson system in 2+2-dimensional phase space. The collision operator is discretized using 4th order accurate conservative finite-differencing. The treatment of the Vlasov operator in phase-space uses an approach based on a minimally diffuse, fourth-order-accurate discretization (Banks and Hittinger, IEEE T. Plasma Sci. 39, 2198). The overall scheme is therefore discretely conservative and controls unphysical oscillations. Some details of the numerical scheme will be presented, and the implementation on modern highly concurrent parallel computers will be discussed. We will present results of collisional effects on linear and non-linear Landau damping of electron plasma waves (EPWs). In addition we will present initial results showing the effect of collisions on the evolution of EPWs in two space dimensions. This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344 and funded by the LDRD program at LLNL under project tracking code 12-ERD-061.
FLUKA: A Multi-Particle Transport Code
Ferrari, A.; Sala, P.R.; Fasso, A.; Ranft, J.; /Siegen U.
2005-12-14
This report describes the 2005 version of the Fluka particle transport code. The first part introduces the basic notions, describes the modular structure of the system, and contains an installation and beginner's guide. The second part complements this initial information with details about the various components of Fluka and how to use them. It concludes with a detailed history and bibliography.
Spin-Orbit Interaction and Related Transport Phenomena in 2d Electron and Hole Systems
NASA Astrophysics Data System (ADS)
Khaetskii, A.
Spin-orbit interaction is responsible for many physical phenomena which are under intensive study currently. Here we discuss several of them. The first phenomenon is the edge spin accumulation, which appears due to spin-orbit interaction in 2D mesoscopic structures in the presence of a charge current. We consider the case of a strong spin-orbit-related splitting of the electron spectrum, i.e. a spin precession length is small compared to the mean free path l. The structure can be either in a ballistic regime (when the mean free path is the largest scale in the problem) or quasi-ballistic regime (when l is much smaller than the sample size). We show how physics of edge spin accumulation in different situations should be understood from the point of view of unitarity of boundary scattering. Using transparent method of scattering states, we are able to explain some previous puzzling theoretical results. We clarify the important role of the form of the spin-orbit Hamiltonian, the role of the boundary conditions, etc., and reveal the wrong results obtained in the field by other researchers. The relation between the edge spin density and the bulk spin current in different regimes is discussed. The detailed comparison with the existing theoretical works is presented. Besides, we consider several new transport phenomena which appear in the presence of spin-orbit interaction, for example, magnetotransport phenomena in an external classical magnetic field. In particular, new mechanism of negative magneto-resistance appears which is due to destruction of spin fluxes by the magnetic field, and which can be really pronounced in 2D systems with strong scatterers.
High Energy Particle Transport Code System.
Energy Science and Technology Software Center (ESTSC)
2003-12-17
Version 00 NMTC/JAM is an upgraded version of the code CCC-694/NMTC-JAERI97, which was developed in 1982 at JAERI and is based on the CCC-161/NMTC code system. NMTC/JAM simulates high energy nuclear reactions and nuclear meson transport processes. The applicable energy range of NMTC/JAM was extended in principle up to 200 GeV for nucleons and mesons by introducing the high energy nuclear reaction code Jet-Aa Microscopic (JAM) for the intra-nuclear cascade part. For the evaporation andmore » fission process, a new model, GEM, can be used to describe the light nucleus production from the excited residual nucleus. According to the extension of the applicable energy, the nucleon-nucleus non-elastic, elastic and differential elastic cross section data were upgraded. In addition, the particle transport in a magnetic field was implemented for beam transport calculations. Some new tally functions were added, and the format of input and output of data is more user friendly. These new calculation functions and utilities provide a tool to carry out reliable neutronics study of a large scale target system with complex geometry more accurately and easily than with the previous model. It implements an intranuclear cascade model taking account of the in-medium nuclear effects and the preequilibrium calculation model based on the exciton one. For treating the nucleon transport process, the nucleon-nucleus cross sections are revised to those derived by the systematics of Pearlstein. Moreover, the level density parameter derived by Ignatyuk is included as a new option for particle evaporation calculation. A geometry package based on the Combinatorial Geometry with multi-array system and the importance sampling technique is implemented in the code. Tally function is also employed for obtaining such physical quantities as neutron energy spectra, heat deposition and nuclide yield without editing a history file. The code can simulate both the primary spallation reaction and the
2D Resistive Magnetohydrodynamics Calculations with an Arbitrary Lagrange Eulerian Code
NASA Astrophysics Data System (ADS)
Rousculp, C. L.; Gianakon, T. A.; Lipnikov, K. N.; Nelson, E. M.
2015-11-01
Single fluid resistive MHD is useful for modeling Z-pinch configurations in cylindrical geometry. One such example is thin walled liners for shock physics or HEDP experiments driven by capacitor banks such as the LANL's PHELIX or Sandia-Z. MHD is also useful for modeling high-explosive-driven flux compression generators (FCGs) and their high-current switches. The resistive MHD in our arbitrary Lagrange Eulerian (ALE) code operates in one and two dimensions in both Cartesian and cylindrical geometry. It is implemented as a time-step split operator, which consists of, ideal MHD connected to the explicit hydro momentum and energy equations and a second order mimetic discretization solver for implicit solution of the magnetic diffusion equation. In a staggered grid scheme, a single-component of cell-centered magnetic flux is conserved in the Lagrangian frame exactly, while magnetic forces are accumulated at the nodes. Total energy is conserved to round off. Total flux is conserved under the ALE relaxation and remap. The diffusion solver consistently computes Ohmic heating. Both Neumann and Dirichlet boundary conditions are available with coupling to external circuit models. Example calculations will be shown.
A 2D forward and inverse code for streaming potential problems
NASA Astrophysics Data System (ADS)
Soueid Ahmed, A.; Jardani, A.; Revil, A.
2013-12-01
The self-potential method corresponds to the passive measurement of the electrical field in response to the occurrence of natural sources of current in the ground. One of these sources corresponds to the streaming current associated with the flow of the groundwater. We can therefore apply the self- potential method to recover non-intrusively some information regarding the groundwater flow. We first solve the forward problem starting with the solution of the groundwater flow problem, then computing the source current density, and finally solving a Poisson equation for the electrical potential. We use the finite-element method to solve the relevant partial differential equations. In order to reduce the number of (petrophysical) model parameters required to solve the forward problem, we introduced an effective charge density tensor of the pore water, which can be determined directly from the permeability tensor for neutral pore waters. The second aspect of our work concerns the inversion of the self-potential data using Tikhonov regularization with smoothness and weighting depth constraints. This approach accounts for the distribution of the electrical resistivity, which can be independently and approximately determined from electrical resistivity tomography. A numerical code, SP2DINV, has been implemented in Matlab to perform both the forward and inverse modeling. Three synthetic case studies are discussed.
THE MCNPX MONTE CARLO RADIATION TRANSPORT CODE
WATERS, LAURIE S.; MCKINNEY, GREGG W.; DURKEE, JOE W.; FENSIN, MICHAEL L.; JAMES, MICHAEL R.; JOHNS, RUSSELL C.; PELOWITZ, DENISE B.
2007-01-10
MCNPX (Monte Carlo N-Particle eXtended) is a general-purpose Monte Carlo radiation transport code with three-dimensional geometry and continuous-energy transport of 34 particles and light ions. It contains flexible source and tally options, interactive graphics, and support for both sequential and multi-processing computer platforms. MCNPX is based on MCNP4B, and has been upgraded to most MCNP5 capabilities. MCNP is a highly stable code tracking neutrons, photons and electrons, and using evaluated nuclear data libraries for low-energy interaction probabilities. MCNPX has extended this base to a comprehensive set of particles and light ions, with heavy ion transport in development. Models have been included to calculate interaction probabilities when libraries are not available. Recent additions focus on the time evolution of residual nuclei decay, allowing calculation of transmutation and delayed particle emission. MCNPX is now a code of great dynamic range, and the excellent neutronics capabilities allow new opportunities to simulate devices of interest to experimental particle physics; particularly calorimetry. This paper describes the capabilities of the current MCNPX version 2.6.C, and also discusses ongoing code development.
NASA Technical Reports Server (NTRS)
Fleming, Eric L.; Jackman, Charles H.; Considine, David B.; Stolarski, Richard S.
1999-01-01
In this study, we examine the sensitivity of long lived tracers to changes in the base transport components in our 2-D model. Changes to the strength of the residual circulation in the upper troposphere and stratosphere and changes to the lower stratospheric K(sub zz) had similar effects in that increasing the transport rates decreased the overall stratospheric mean age, and increased the rate of removal of material from the stratosphere. Increasing the stratospheric K(sub yy) increased the mean age due to the greater recycling of air parcels through the middle atmosphere, via the residual circulation, before returning to the troposphere. However, increasing K(sub yy) along with self-consistent increases in the corresponding planetary wave drive, which leads to a stronger residual circulation, more than compensates for the K(sub yy)-effect, and produces significantly younger ages throughout the stratosphere. Simulations with very small tropical stratospheric K(sub yy) decreased the globally averaged age of air by as much as 25% in the middle and upper stratosphere, and resulted in substantially weaker vertical age gradients above 20 km in the extratropics. We found only very small stratospheric tracer sensitivity to the magnitude of the horizontal mixing across the tropopause, and to the strength of the mesospheric gravity wave drag and diffusion used in the model. We also investigated the transport influence on chemically active tracers and found a strong age-tracer correlation, both in concentration and calculated lifetimes. The base model transport gives the most favorable overall comparison with a variety of inert tracer observations, and provides a significant improvement over our previous 1995 model transport. Moderate changes to the base transport were found to provide modest agreement with some of the measurements. Transport scenarios with residence times ranging from moderately shorter to slightly longer relative to the base case simulated N2O lifetimes
Hyde, R.A.
1983-06-08
This paper describes a computer code, RYLLA, which models the deposition of x-rays into thin metal slabs, and transports the resulting photoelectrons, finding the distribution of electrons leaving the slab from both the front and back surfaces. The slab must be homogeneous, but can contain a mixture of up to 5 different elements. Due to the short electron mean free path at low electron energies, RYLLA should be used only for studying thin slabs, roughly < 100 mg/cm/sup 2/ for low Z metals, and < 10 mg/cm/sup 2/ for high Z metals. X-ray energies should be in the range of 1 to 150 keV, as they are deposited only via photoionization and Compton scattering processes. Following photoionization, a hole exists in the electron cloud of the absorbing atom. This fills either by Auger or fluoresence, resulting in lower energy holes which are also filled. Fluoresence photons are transported and absorbed in the same manner as the primary photons, except that they are isotropically produced. Once all photons have been transported and absorbed, and all holes have been filled, a space- and energy-dependent electron source spectrum has been obtained. This is used in a discrete ordinate expansion solution of the 1-D transport equation, which gives the output electron spectra at the two slab surfaces. This paper discusses both the physics and coding of RYLLA. Examples of user input are given, as are some comparisons with other codes.
NASA Astrophysics Data System (ADS)
Velioǧlu, Deniz; Cevdet Yalçıner, Ahmet; Zaytsev, Andrey
2016-04-01
Tsunamis are huge waves with long wave periods and wave lengths that can cause great devastation and loss of life when they strike a coast. The interest in experimental and numerical modeling of tsunami propagation and inundation increased considerably after the 2011 Great East Japan earthquake. In this study, two numerical codes, FLOW 3D and NAMI DANCE, that analyze tsunami propagation and inundation patterns are considered. Flow 3D simulates linear and nonlinear propagating surface waves as well as long waves by solving three-dimensional Navier-Stokes (3D-NS) equations. NAMI DANCE uses finite difference computational method to solve 2D depth-averaged linear and nonlinear forms of shallow water equations (NSWE) in long wave problems, specifically tsunamis. In order to validate these two codes and analyze the differences between 3D-NS and 2D depth-averaged NSWE equations, two benchmark problems are applied. One benchmark problem investigates the runup of long waves over a complex 3D beach. The experimental setup is a 1:400 scale model of Monai Valley located on the west coast of Okushiri Island, Japan. Other benchmark problem is discussed in 2015 National Tsunami Hazard Mitigation Program (NTHMP) Annual meeting in Portland, USA. It is a field dataset, recording the Japan 2011 tsunami in Hilo Harbor, Hawaii. The computed water surface elevation and velocity data are compared with the measured data. The comparisons showed that both codes are in fairly good agreement with each other and benchmark data. The differences between 3D-NS and 2D depth-averaged NSWE equations are highlighted. All results are presented with discussions and comparisons. Acknowledgements: Partial support by Japan-Turkey Joint Research Project by JICA on earthquakes and tsunamis in Marmara Region (JICA SATREPS - MarDiM Project), 603839 ASTARTE Project of EU, UDAP-C-12-14 project of AFAD Turkey, 108Y227, 113M556 and 213M534 projects of TUBITAK Turkey, RAPSODI (CONCERT_Dis-021) of CONCERT
Ensemble Solute Transport in 2-D Operator-Stable Random Fields
NASA Astrophysics Data System (ADS)
Monnig, N. D.; Benson, D. A.
2006-12-01
The heterogeneous velocity field that exists at many scales in an aquifer will typically cause a dissolved solute plume to grow at a rate faster than Fick's Law predicts. Some statistical model must be adopted to account for the aquifer structure that engenders the velocity heterogeneity. A fractional Brownian motion (fBm) model has been shown to create the long-range correlation that can produce continually faster-than-Fickian plume growth. Previous fBm models have assumed isotropic scaling (defined here by a scalar Hurst coefficient). Motivated by field measurements of aquifer hydraulic conductivity, recent techniques were developed to construct random fields with anisotropic scaling with a self-similarity parameter that is defined by a matrix. The growth of ensemble plumes is analyzed for transport through 2-D "operator- stable" fBm hydraulic conductivity (K) fields. Both the longitudinal and transverse Hurst coefficients are important to both plume growth rates and the timing and duration of breakthrough. Smaller Hurst coefficients in the transverse direction lead to more "continuity" or stratification in the direction of transport. The result is continually faster-than-Fickian growth rates, highly non-Gaussian ensemble plumes, and a longer tail early in the breakthrough curve. Contrary to some analytic stochastic theories for monofractal K fields, the plume growth rate never exceeds Mercado's [1967] purely stratified aquifer growth rate of plume apparent dispersivity proportional to mean distance. Apparent super-Mercado growth must be the result of other factors, such as larger plumes corresponding to either a larger initial plume size or greater variance of the ln(K) field.
Validation of comprehensive space radiation transport code
Shinn, J.L.; Simonsen, L.C.; Cucinotta, F.A.
1998-12-01
The HZETRN code has been developed over the past decade to evaluate the local radiation fields within sensitive materials on spacecraft in the space environment. Most of the more important nuclear and atomic processes are now modeled and evaluation within a complex spacecraft geometry with differing material components, including transition effects across boundaries of dissimilar materials, are included. The atomic/nuclear database and transport procedures have received limited validation in laboratory testing with high energy ion beams. The codes have been applied in design of the SAGE-III instrument resulting in material changes to control injurious neutron production, in the study of the Space Shuttle single event upsets, and in validation with space measurements (particle telescopes, tissue equivalent proportional counters, CR-39) on Shuttle and Mir. The present paper reviews the code development and presents recent results in laboratory and space flight validation.
Hallquist, J.O.
1982-02-01
This revised report provides an updated user's manual for DYNA2D, an explicit two-dimensional axisymmetric and plane strain finite element code for analyzing the large deformation dynamic and hydrodynamic response of inelastic solids. A contact-impact algorithm permits gaps and sliding along material interfaces. By a specialization of this algorithm, such interfaces can be rigidly tied to admit variable zoning without the need of transition regions. Spatial discretization is achieved by the use of 4-node solid elements, and the equations-of motion are integrated by the central difference method. An interactive rezoner eliminates the need to terminate the calculation when the mesh becomes too distorted. Rather, the mesh can be rezoned and the calculation continued. The command structure for the rezoner is described and illustrated by an example.
NASA Astrophysics Data System (ADS)
Goble, Nicholas; Watson, John; Manfra, Michael; Gao, Xuan
2014-03-01
Understanding the non-monotonic behavior in the temperature dependent resistance, R(T) , of strongly correlated two-dimensional (2D) carriers in clean semiconductors has been a central issue in the studies of 2D metallic states and metal-insulator transitions. We have studied the transport of high mobility 2D holes in 20nm wide GaAs quantum wells with varying interface roughness by changing the Al fraction x in the AlxGa1-xAs barrier. Prior to this work, no comprehensive study of the non-monotonic resistance peak against controlled barrier characteristics has been conducted. We show that the shape of the electronic contribution to R(T) is qualitatively unchanged throughout all of our measurements, regardless of the percentage of Al in the barrier. It is observed that increasing x or short range interface roughness suppresses both the strength and characteristic temperature scale of the 2D metallicity, pointing to the distinct role of short range versus long range disorder in the 2D metallic transport in this 2D hole system with interaction parameter rs ~ 20. N.G. acknowledges the US DOE GAANN fellowship (P200A090276 & P200A070434). M.J.M. is supported by the Miller Family Foundation and the US DOE, Office of Basic Energy Sciences, DMS (DE-SC0006671). X.P.A.G thanks the NSF for funding support (DMR-0906415).
In situ 2D-extraction of DNA wheels by 3D through-solution transport.
Yonamine, Yusuke; Cervantes-Salguero, Keitel; Nakanishi, Waka; Kawamata, Ibuki; Minami, Kosuke; Komatsu, Hirokazu; Murata, Satoshi; Hill, Jonathan P; Ariga, Katsuhiko
2015-12-28
Controlled transfer of DNA nanowheels from a hydrophilic to a hydrophobic surface was achieved by complexation of the nanowheels with a cationic lipid (2C12N(+)). 2D surface-assisted extraction, '2D-extraction', enabled structure-persistent transfer of DNA wheels, which could not be achieved by simple drop-casting. PMID:26583486
3D Multigroup Sn Neutron Transport Code
Energy Science and Technology Software Center (ESTSC)
2001-02-14
ATTILA is a 3D multigroup transport code with arbitrary order ansotropic scatter. The transport equation is solved in first order form using a tri-linear discontinuous spatial differencing on an arbitrary tetrahedral mesh. The overall solution technique is source iteration with DSA acceleration of the scattering source. Anisotropic boundary and internal sources may be entered in the form of spherical harmonics moments. Alpha and k eigenvalue problems are allowed, as well as fixed source problems. Forwardmore » and adjoint solutions are available. Reflective, vacumn, and source boundary conditions are available. ATTILA can perform charged particle transport calculations using slowing down (CSD) terms. ATTILA can also be used to peform infra-red steady-state calculations for radiative transfer purposes.« less
3D Multigroup Sn Neutron Transport Code
McGee, John; Wareing, Todd; Pautz, Shawn
2001-02-14
ATTILA is a 3D multigroup transport code with arbitrary order ansotropic scatter. The transport equation is solved in first order form using a tri-linear discontinuous spatial differencing on an arbitrary tetrahedral mesh. The overall solution technique is source iteration with DSA acceleration of the scattering source. Anisotropic boundary and internal sources may be entered in the form of spherical harmonics moments. Alpha and k eigenvalue problems are allowed, as well as fixed source problems. Forward and adjoint solutions are available. Reflective, vacumn, and source boundary conditions are available. ATTILA can perform charged particle transport calculations using slowing down (CSD) terms. ATTILA can also be used to peform infra-red steady-state calculations for radiative transfer purposes.
Eyler, L.L.; Budden, M.J.
1985-03-01
The objective of this work is to assess prediction capabilities and features of the MAGNUM-2D computer code in relation to its intended use in the Basalt Waste Isolation Project (BWIP). This objective is accomplished through a code verification and benchmarking task. Results are documented which support correctness of prediction capabilities in areas of intended model application. 10 references, 43 figures, 11 tables.
Space Radiation Transport Code Development: 3DHZETRN
NASA Technical Reports Server (NTRS)
Wilson, John W.; Slaba, Tony C.; Badavi, Francis F.; Reddell, Brandon D.; Bahadori, Amir A.
2015-01-01
The space radiation transport code, HZETRN, has been used extensively for research, vehicle design optimization, risk analysis, and related applications. One of the simplifying features of the HZETRN transport formalism is the straight-ahead approximation, wherein all particles are assumed to travel along a common axis. This reduces the governing equation to one spatial dimension allowing enormous simplification and highly efficient computational procedures to be implemented. Despite the physical simplifications, the HZETRN code is widely used for space applications and has been found to agree well with fully 3D Monte Carlo simulations in many circumstances. Recent work has focused on the development of 3D transport corrections for neutrons and light ions (Z < 2) for which the straight-ahead approximation is known to be less accurate. Within the development of 3D corrections, well-defined convergence criteria have been considered, allowing approximation errors at each stage in model development to be quantified. The present level of development assumes the neutron cross sections have an isotropic component treated within N explicit angular directions and a forward component represented by the straight-ahead approximation. The N = 1 solution refers to the straight-ahead treatment, while N = 2 represents the bi-directional model in current use for engineering design. The figure below shows neutrons, protons, and alphas for various values of N at locations in an aluminum sphere exposed to a solar particle event (SPE) spectrum. The neutron fluence converges quickly in simple geometry with N > 14 directions. The improved code, 3DHZETRN, transports neutrons, light ions, and heavy ions under space-like boundary conditions through general geometry while maintaining a high degree of computational efficiency. A brief overview of the 3D transport formalism for neutrons and light ions is given, and extensive benchmarking results with the Monte Carlo codes Geant4, FLUKA, and
Impurity transport and radiated power loss estimates in NSTX using 2-d USXR arrays
NASA Astrophysics Data System (ADS)
Stutman, Dan; Finkenthal, Michael; Vero, Robert; Roquemore, Lane; Johnson, David; Kaita, Robert
2000-10-01
Large area, low capacitance, absolute photodiodes have been installed in the three ultrasoft X-ray (USXR) arrays operational on NSTX, for the beginning of the high power operation phase. In addition to bandpass filtered measurements in the USXR range, the absolute diodes enable total radiated power measurements and improve the accuracy of the absolute emissivity estimates. The position of the viewing chords with respect to the NSTX vacuum vessel has been measured using a laser pointing technique and the vignetting of the top array by in-vessel structures calibrated using an in-vessel extended light source. The Granetz-Cormack algorithm with Bessel radial functions is used to derive emissivity maps from the measured brightness profiles. The plasma emission in the spectral ranges defined by the bandpass filters on each array (0.3 μm Ti, 10 μm, 100 μm and 500 μm Be) is modeled using impurity line emission data computed with the HULLAC atomic physics package coupled to a 1-d impurity transport code. The evolution of the USXR emissivity, radiated power, and estimated impurity and plasma profiles during MHD phenomena like the Internal Reconnection Event and sawteeth, as well as emission data for ohmic, auxiliary heated and coaxial helicity injection discharges are presented.
Time-Dependent, Parallel Neutral Particle Transport Code System.
Energy Science and Technology Software Center (ESTSC)
2009-09-10
Version 00 PARTISN (PARallel, TIme-Dependent SN) is the evolutionary successor to CCC-547/DANTSYS. The PARTISN code package is a modular computer program package designed to solve the time-independent or dependent multigroup discrete ordinates form of the Boltzmann transport equation in several different geometries. The modular construction of the package separates the input processing, the transport equation solving, and the post processing (or edit) functions into distinct code modules: the Input Module, the Solver Module, and themore » Edit Module, respectively. PARTISN is the evolutionary successor to the DANTSYSTM code system package. The Input and Edit Modules in PARTISN are very similar to those in DANTSYS. However, unlike DANTSYS, the Solver Module in PARTISN contains one, two, and three-dimensional solvers in a single module. In addition to the diamond-differencing method, the Solver Module also has Adaptive Weighted Diamond-Differencing (AWDD), Linear Discontinuous (LD), and Exponential Discontinuous (ED) spatial differencing methods. The spatial mesh may consist of either a standard orthogonal mesh or a block adaptive orthogonal mesh. The Solver Module may be run in parallel for two and three dimensional problems. One can now run 1-D problems in parallel using Energy Domain Decomposition (triggered by Block 5 input keyword npeg>0). EDD can also be used in 2-D/3-D with or without our standard Spatial Domain Decomposition. Both the static (fixed source or eigenvalue) and time-dependent forms of the transport equation are solved in forward or adjoint mode. In addition, PARTISN now has a probabilistic mode for Probability of Initiation (static) and Probability of Survival (dynamic) calculations. Vacuum, reflective, periodic, white, or inhomogeneous boundary conditions are solved. General anisotropic scattering and inhomogeneous sources are permitted. PARTISN solves the transport equation on orthogonal (single level or block-structured AMR) grids in 1-D
2-D Path Corrections for Local and Regional Coda Waves: A Test of Transportability
Mayeda, K M; Malagnini, L; Phillips, W S; Walter, W R; Dreger, D S; Morasca, P
2005-07-13
Reliable estimates of the seismic source spectrum are necessary for accurate magnitude, yield, and energy estimation. In particular, how seismic radiated energy scales with increasing earthquake size has been the focus of recent debate within the community and has direct implications on earthquake source physics studies as well as hazard mitigation. The 1-D coda methodology of Mayeda et al. [2003] has provided the lowest variance estimate of the source spectrum when compared against traditional approaches that use direct S-waves, thus making it ideal for networks that have sparse station distribution. The 1-D coda methodology has been mostly confined to regions of approximately uniform complexity. For larger, more geophysically complicated regions, 2-D path corrections may be required. We will compare performance of 1-D versus 2-D path corrections in a variety of regions. First, the complicated tectonics of the northern California region coupled with high quality broadband seismic data provides for an ideal ''apples-to-apples'' test of 1-D and 2-D path assumptions on direct waves and their coda. Next, we will compare results for the Italian Alps using high frequency data from the University of Genoa. For Northern California, we used the same station and event distribution and compared 1-D and 2-D path corrections and observed the following results: (1) 1-D coda results reduced the amplitude variance relative to direct S-waves by roughly a factor of 8 (800%); (2) Applying a 2-D correction to the coda resulted in up to 40% variance reduction from the 1-D coda results; (3) 2-D direct S-wave results, though better than 1-D direct waves, were significantly worse than the 1-D coda. We found that coda-based moment-rate source spectra derived from the 2-D approach were essentially identical to those from the 1-D approach for frequencies less than {approx}0.7-Hz, however for the high frequencies (0.7 {le} f {le} 8.0-Hz), the 2-D approach resulted in inter-station scatter
Icarus: A 2D direct simulation Monte Carlo (DSMC) code for parallel computers. User`s manual - V.3.0
Bartel, T.; Plimpton, S.; Johannes, J.; Payne, J.
1996-10-01
Icarus is a 2D Direct Simulation Monte Carlo (DSMC) code which has been optimized for the parallel computing environment. The code is based on the DSMC method of Bird and models from free-molecular to continuum flowfields in either cartesian (x, y) or axisymmetric (z, r) coordinates. Computational particles, representing a given number of molecules or atoms, are tracked as they have collisions with other particles or surfaces. Multiple species, internal energy modes (rotation and vibration), chemistry, and ion transport are modelled. A new trace species methodology for collisions and chemistry is used to obtain statistics for small species concentrations. Gas phase chemistry is modelled using steric factors derived from Arrhenius reaction rates. Surface chemistry is modelled with surface reaction probabilities. The electron number density is either a fixed external generated field or determined using a local charge neutrality assumption. Ion chemistry is modelled with electron impact chemistry rates and charge exchange reactions. Coulomb collision cross-sections are used instead of Variable Hard Sphere values for ion-ion interactions. The electrostatic fields can either be externally input or internally generated using a Langmuir-Tonks model. The Icarus software package includes the grid generation, parallel processor decomposition, postprocessing, and restart software. The commercial graphics package, Tecplot, is used for graphics display. The majority of the software packages are written in standard Fortran.
Comparisons of anomalous and collisional radial transport with a continuum kinetic edge code
NASA Astrophysics Data System (ADS)
Bodi, K.; Krasheninnikov, S.; Cohen, R.; Rognlien, T.
2009-05-01
Modeling of anomalous (turbulence-driven) radial transport in controlled-fusion plasmas is necessary for long-time transport simulations. Here the focus is continuum kinetic edge codes such as the (2-D, 2-V) transport version of TEMPEST, NEO, and the code being developed by the Edge Simulation Laboratory, but the model also has wider application. Our previously developed anomalous diagonal transport matrix model with velocity-dependent convection and diffusion coefficients allows contact with typical fluid transport models (e.g., UEDGE). Results are presented that combine the anomalous transport model and collisional transport owing to ion drift orbits utilizing a Krook collision operator that conserves density and energy. Comparison is made of the relative magnitudes and possible synergistic effects of the two processes for typical tokamak device parameters.
Strongly Metallic Electron and Hole 2D Transport in an Ambipolar Si-Vacuum Field Effect Transistor
NASA Astrophysics Data System (ADS)
Hu, Binhui; Yazdanpanah, M. M.; Kane, B. E.; Hwang, E. H.; Das Sarma, S.
2015-07-01
We report experiment and theory on an ambipolar gate-controlled Si(111)-vacuum field effect transistor where we study electron and hole (low-temperature 2D) transport in the same device simply by changing the external gate voltage to tune the system from being a 2D electron system at positive gate voltage to a 2D hole system at negative gate voltage. The electron (hole) conductivity manifests strong (moderate) metallic temperature dependence with the conductivity decreasing by a factor of 8 (2) between 0.3 K and 4.2 K with the peak electron mobility (˜18 m2/V s ) being roughly 20 times larger than the peak hole mobility (in the same sample). Our theory explains the data well using random phase approximation screening of background Coulomb disorder, establishing that the observed metallicity is a direct consequence of the strong temperature dependence of the effective screened disorder.
Electronic Transport Properties of New 2-D Materials GeH and NaSn2As2
NASA Astrophysics Data System (ADS)
He, Bin; Cultrara, Nicholas; Arguilla, Maxx; Goldberger, Joshua; Heremans, Joseph
2-D materials potentially have superior thermoelectric properties compared to traditional 3-D materials due to their layered structure. Here we present electrical and thermoelectric transport properties of 2 types of 2-D materials, GeH and NaSn2As2. GeH is a graphane analog which is prepared using chemical exfoliation of CaGe2 crystals. Intrinsic GeH is proven to be a highly resistive material at room temperature. Resistance and Seebeck coefficient of Ga doped GeH are measured in a cryostat with a gating voltage varying from -100V to 100V. NaSn2As2 is another 2-D system, with Na atom embedded between nearly-2D Sn-As layers. Unlike GeH, NaSn2As2 is a metal based of Hall measurements, with p-type behavior, and with van der Pauw resistances on the order of 5m Ω/square. Thermoelectric transport properties of NaSn2As2 will be reported. This work is support by the NSF EFRI-2DARE project EFRI-1433467.
2D Electrical Resistivity Tomography surveys optimisation of the solutes transports in porous media.
NASA Astrophysics Data System (ADS)
Lekmine, G.; Pessel, M.; Auradou, H.
2009-04-01
Electrical resistivity tomography applied in borehole or cross-borehole is a method often used to follow the invasion process of pollutant [Daily, 1991]. The aim of this work is to test experimentally the electrode arrays and inversion process used to obtain a spatial representation of tracer propagation in porous media. Experiments were conducted in a plexiglas container with glass beads of 166 microns in diameter. The height of the container is 275 mm, its width 85 mm and its thickness 10 mm. 21 electrodes, equally spaced, are placed along each of the lateral sides of the porous medium : these electrodes are used to perform the electrical measurements. The porous medium is lightened from behind and a video camera records the propagation of the fluids. The fluid containing the tracer (i.e the pollutant) is a water solution containing a small amount of dye together with NaCl (0.5g/l up to 2.0g/l). The medium is first saturated by a water solution containing a slight concentration of NaCl so that its density is smaller than the injected fluid. An upward flow is first established, then the denser fluid is injected at the bottom and over the full width of the medium. In this way, the flow is stabilized by gravity avoiding the development of unstable fingers. Still, the fluids are miscible and a mixing front develops during the flow: in the present study, both the determination by optical and electrical imaging of the mean position of the front and its width are of interest. The comparison of the two techniques allows to study the ability of the inversion process to quantify the solute transport. The electrical measurements are acquired by a standard multi electrode system (IRIS Instruments) and the data are inverted with the Res2Dinv software which models the 2D distribution of conductivity contrasts. The obtained bulk conductivity can be related through Archie's law to fluid conductivity by the porosity and the cementation factor which have been experimentally
Li, Yixiang; Qiu, Chunyin; Xu, Shengjun; Ke, Manzhu; Liu, Zhengyou
2015-01-01
Conventional microparticle transports by light or sound are realized along a straight line. Recently, this limit has been overcome in optics as the growing up of the self-accelerating Airy beams, which are featured by many peculiar properties, e.g., bending propagation, diffraction-free and self-healing. However, the bending angles of Airy beams are rather small since they are only paraxial solutions of the two-dimensional (2D) Helmholtz equation. Here we propose a novel micromanipulation by using acoustic Half-Bessel beams, which are strict solutions of the 2D Helmholtz equation. Compared with that achieved by Airy beams, the bending angle of the particle trajectory attained here is much steeper (exceeding 90o). The large-angle bending transport of microparticles, which is robust to complex scattering environment, enables a wide range of applications from the colloidal to biological sciences. PMID:26279478
NASA Astrophysics Data System (ADS)
Li, Yixiang; Qiu, Chunyin; Xu, Shengjun; Ke, Manzhu; Liu, Zhengyou
2015-08-01
Conventional microparticle transports by light or sound are realized along a straight line. Recently, this limit has been overcome in optics as the growing up of the self-accelerating Airy beams, which are featured by many peculiar properties, e.g., bending propagation, diffraction-free and self-healing. However, the bending angles of Airy beams are rather small since they are only paraxial solutions of the two-dimensional (2D) Helmholtz equation. Here we propose a novel micromanipulation by using acoustic Half-Bessel beams, which are strict solutions of the 2D Helmholtz equation. Compared with that achieved by Airy beams, the bending angle of the particle trajectory attained here is much steeper (exceeding 90o). The large-angle bending transport of microparticles, which is robust to complex scattering environment, enables a wide range of applications from the colloidal to biological sciences.
Li, Yixiang; Qiu, Chunyin; Xu, Shengjun; Ke, Manzhu; Liu, Zhengyou
2015-01-01
Conventional microparticle transports by light or sound are realized along a straight line. Recently, this limit has been overcome in optics as the growing up of the self-accelerating Airy beams, which are featured by many peculiar properties, e.g., bending propagation, diffraction-free and self-healing. However, the bending angles of Airy beams are rather small since they are only paraxial solutions of the two-dimensional (2D) Helmholtz equation. Here we propose a novel micromanipulation by using acoustic Half-Bessel beams, which are strict solutions of the 2D Helmholtz equation. Compared with that achieved by Airy beams, the bending angle of the particle trajectory attained here is much steeper (exceeding 90(o)). The large-angle bending transport of microparticles, which is robust to complex scattering environment, enables a wide range of applications from the colloidal to biological sciences. PMID:26279478
Topologically robust transport of entangled photons in a 2D photonic system.
Mittal, Sunil; Orre, Venkata Vikram; Hafezi, Mohammad
2016-07-11
We theoretically study the transport of time-bin entangled photon pairs in a two-dimensional topological photonic system of coupled ring resonators. This system implements the integer quantum Hall model using a synthetic gauge field and exhibits topologically robust edge states. We show that the transport through edge states preserves temporal correlations of entangled photons whereas bulk transport does not preserve these correlations and can lead to significant unwanted temporal bunching or anti-bunching of photons. We study the effect of disorder on the quantum transport properties; while the edge transport remains robust, bulk transport is very susceptible, and in the limit of strong disorder, bulk states become localized. We show that this localization is manifested as an enhanced bunching/anti-bunching of photons. This topologically robust transport of correlations through edge states could enable robust on-chip quantum communication channels and delay lines for information encoded in temporal correlations of photons. PMID:27410836
Los Alamos radiation transport code system on desktop computing platforms
Briesmeister, J.F.; Brinkley, F.W.; Clark, B.A.; West, J.T. )
1990-01-01
The Los Alamos Radiation Transport Code System (LARTCS) consists of state-of-the-art Monte Carlo and discrete ordinates transport codes and data libraries. These codes were originally developed many years ago and have undergone continual improvement. With a large initial effort and continued vigilance, the codes are easily portable from one type of hardware to another. The performance of scientific work-stations (SWS) has evolved to the point that such platforms can be used routinely to perform sophisticated radiation transport calculations. As the personal computer (PC) performance approaches that of the SWS, the hardware options for desk-top radiation transport calculations expands considerably. The current status of the radiation transport codes within the LARTCS is described: MCNP, SABRINA, LAHET, ONEDANT, TWODANT, TWOHEX, and ONELD. Specifically, the authors discuss hardware systems on which the codes run and present code performance comparisons for various machines.
SYMTRAN - A Time-dependent Symmetric Tandem Mirror Transport Code
Hua, D; Fowler, T
2004-06-15
A time-dependent version of the steady-state radial transport model in symmetric tandem mirrors in Ref. [1] has been coded up and first tests performed. Our code, named SYMTRAN, is an adaptation of the earlier SPHERE code for spheromaks, now modified for tandem mirror physics. Motivated by Post's new concept of kinetic stabilization of symmetric mirrors, it is an extension of the earlier TAMRAC rate-equation code omitting radial transport [2], which successfully accounted for experimental results in TMX. The SYMTRAN code differs from the earlier tandem mirror radial transport code TMT in that our code is focused on axisymmetric tandem mirrors and classical diffusion, whereas TMT emphasized non-ambipolar transport in TMX and MFTF-B due to yin-yang plugs and non-symmetric transitions between the plugs and axisymmetric center cell. Both codes exhibit interesting but different non-linear behavior.
Crossover from 3D to 2D Quantum Transport in Bi2Se3/In2Se3 Superlattices
NASA Astrophysics Data System (ADS)
Yanfei, Zhao; Haiwen, Liu; Xin, Guo; Ying, Jiang; Yi, Sun; Huichao, Wang; Yong, Wang; Handong, Li; Maohai, Xie; Xincheng, Xie; Jian, Wang
2015-03-01
The topological insulator/normal insulator (TI/NI) superlattices (SLs) with multiple Dirac channels are predicted to offer great opportunity to design novel materials and investigate new quantum phenomena. Here, we report first transport studies on the SLs composed of TI Bi2Se3 layers sandwiched by NI In2Se3 layers artificially grown by molecular beam epitaxy (MBE). The transport properties of two kinds of SL samples show convincing evidence that the transport dimensionality changes from three-dimensional (3D) to two-dimensional (2D) when decreasing the thickness of building block Bi2Se3 layers, corresponding to the crossover from coherent TI transport to separated TI channels. Our findings provide the possibility to realizing 3D surface states in TI/NI SLs.
Crossover from 3D to 2D quantum transport in Bi2Se3/In2Se3 superlattices.
Zhao, Yanfei; Liu, Haiwen; Guo, Xin; Jiang, Ying; Sun, Yi; Wang, Huichao; Wang, Yong; Li, Han-Dong; Xie, Mao-Hai; Xie, Xin-Cheng; Wang, Jian
2014-09-10
The topological insulator/normal insulator (TI/NI) superlattices (SLs) with multiple Dirac channels are predicted to offer great opportunity to design novel materials and investigate new quantum phenomena. Here, we report first transport studies on the SLs composed of TI Bi2Se3 layers sandwiched by NI In2Se3 layers artificially grown by molecular beam epitaxy (MBE). The transport properties of two kinds of SL samples show convincing evidence that the transport dimensionality changes from three-dimensional (3D) to two-dimensional (2D) when decreasing the thickness of building block Bi2Se3 layers, corresponding to the crossover from coherent TI transport to separated TI channels. Our findings provide the possibility to realizing "3D surface states" in TI/NI SLs. PMID:25102289
Description of Transport Codes for Space Radiation Shielding
NASA Technical Reports Server (NTRS)
Kim, Myung-Hee Y.; Wilson, John W.; Cucinotta, Francis A.
2011-01-01
This slide presentation describes transport codes and their use for studying and designing space radiation shielding. When combined with risk projection models radiation transport codes serve as the main tool for study radiation and designing shielding. There are three criteria for assessing the accuracy of transport codes: (1) Ground-based studies with defined beams and material layouts, (2) Inter-comparison of transport code results for matched boundary conditions and (3) Comparisons to flight measurements. These three criteria have a very high degree with NASA's HZETRN/QMSFRG.
Spin Hall effect and spin transport in graphene and 2D heterostructures
NASA Astrophysics Data System (ADS)
Oezyilmaz, Barbaros
Semiconducting 2D materials offer new opportunities in both alternative technologies and fundamental discoveries by using the spin degree freedom of electrons. One of the main challenges in this field is to identify new materials which allow the control of spin currents by means of the electric field effect. This requires either a sizeable spin-orbit coupling strength or a sizeable bandgap or both. Unfortunately, pristine graphene has a negligibly small spin-orbit coupling strength. Recently we have addressed this problem in three distinct ways. First we have used chemical functionalization to introduce locally sp3 type bonding. Next we used metal ad-atoms to increase spin-orbit coupling via local enhancement of the spin-orbit coupling strength due to resonant scattering. Finally, I will show that the proximity of graphene on transition metal dichalcogenides can also lead to a significant enhancement of the spin-orbit coupling strength. I will complete my talk with a brief discussion on the possibility of all electrical spin injection into complementary 2D crystals such as WS2, MoS2 or black phosphorus. Membership Pending in the abstract Special Instructions field.
High-field and thermal transport in 2D atomic layer devices
NASA Astrophysics Data System (ADS)
Serov, Andrey; Dorgan, Vincent E.; Behnam, Ashkan; English, Chris D.; Li, Zuanyi; Islam, Sharnali; Pop, Eric
2014-06-01
This paper reviews our recent results of high-field electrical and thermal properties of atomically thin two-dimensional materials. We show how self-heating affects velocity saturation in suspended and supported graphene. We also demonstrate that multi-valley transport must be taken into account to describe high-field transport in MoS2. At the same time we characterized thermal properties of suspended and nanoscale graphene samples over a wide range of temperatures. We uncovered the effects of edge scattering and grain boundaries on thermal transport in graphene, and showed how the thermal conductivity varies between diffusive and ballistic heat flow limits.
NASA Astrophysics Data System (ADS)
Iannaccone, G.; Zhang, Q.; Bruzzone, S.; Fiori, G.
2016-01-01
Different proposals of graphene transistors based on off-plane (i.e., vertical) transport, have recently appeared in the literature, exhibiting experimental current modulation of a factor 104-105 at room temperature. These devices overcome the lack of bandgap that undermines the operation of graphene transistors, and positively exploit graphene's ultimate thinness, high conductivity, and low density of states. However, very little is known about vertical transport through graphene and two-dimensional materials, either in terms of experiments or theory. In this paper we will discuss the physics and the electronics of off-plane transport through hetero-structures of graphene and 2D materials. We investigate transport across vertical heterostructures of 2D materials with multi-scale simulations, including first-principle density functional theory and non-equilibrium Green's functions based on NanoTCAD ViDES. We show that unexpected behaviors emerge, which are not observed in the more familiar semiconductor heterostructures based on III-V and II-VI materials systems, and that are not predicted by simplistic physical models. Such properties have a significant impact on the design and performance of transistors for digital or high frequency operations.
Numerical model of water flow and solute accumulation in vertisols using HYDRUS 2D/3D code
NASA Astrophysics Data System (ADS)
Weiss, Tomáš; Dahan, Ofer; Turkeltub, Tuvia
2015-04-01
boundary to the wall of the crack (so that the solute can accumulate due to evaporation on the crack block wall, and infiltrating fresh water can push the solute further down) - in order to do so, HYDRUS 2D/3D code had to be modified by its developers. Unconventionally, the main fitting parameters were: parameter a and n in the soil water retention curve and saturated hydraulic conductivity. The amount of infiltrated water (within a reasonable range), the infiltration function in the crack and the actual evaporation from the crack were also used as secondary fitting parameters. The model supports the previous findings that significant amount (~90%) of water from rain events must infiltrate through the crack. It was also noted that infiltration from the crack has to be increasing with depth and that the highest infiltration rate should be somewhere between 1-3m. This paper suggests a new way how to model vertisols in semi-arid regions. It also supports the previous findings about vertisols: especially, the utmost importance of soil cracks as preferential pathways for water and contaminants and soil cracks as deep evaporators.
NASA Astrophysics Data System (ADS)
Apruzese, J. P.; Giuliani, J. L.; Hansen, S. B.
2012-09-01
In modeling optically thick, high energy density plasmas (HEDP), radiation transport is of comparable importance to atomic physics, hydrodynamics, and other transport processes. Radiation transport theory is the framework for calculating radiation output based on atomic kinetics, and is thus critical to designing and diagnosing experiments in which radiation production is significant. Starting in the 1960s, the astrophysics community established benchmarks for computational radiation transport, based on a parameterized two-level-atom, and mostly applicable to media with very high optical depths and low collisional quenching of radiative transitions. The purpose of the work reported here is to establish a similar computational radiation transport benchmark that is more relevant to a laboratory HEDP environment. Our model consists of 8 levels of mostly K-shell Al ions, ranging from the ground state of Li-like Al to its bare nucleus. Rates connecting the levels are given by well-known, readily reproducible analytic prescriptions. The results presented consist of level populations as a function of position within the cylindrical medium, and emitted line profiles. For the plasma conditions considered, the magnitudes and spatial variations of the populations are a sensitive indicator of the accuracy of the radiation transport model, and are critical in calculating experimentally relevant quantities such as radiative powers and line ratios.
Code System To Analyze Radiological Impact From Radwaste Transportation.
Energy Science and Technology Software Center (ESTSC)
1988-05-01
Version 00 RADSHIP-2 is a computer code system used to analyze the environmental impact of radwaste transportation in Taiwan. The specific transport scheme including the land transport by truck and sea transport by ship or barge were considered in the analysis for normal transport and transport accident conditions. The code combines meteorological, population, health physics, transportation, packaging and material factors and has the capability to obtain the results of the expected annual population radiation exposure,more » the expected number of annual latent cancer fatalities and the annual probability of a given number of early fatalities.« less
Pangolin v1.0, a conservative 2-D transport model for large scale parallel calculation
NASA Astrophysics Data System (ADS)
Praga, A.; Cariolle, D.; Giraud, L.
2014-07-01
To exploit the possibilities of parallel computers, we designed a large-scale bidimensional atmospheric transport model named Pangolin. As the basis for a future chemistry-transport model, a finite-volume approach was chosen both for mass preservation and to ease parallelization. To overcome the pole restriction on time-steps for a regular latitude-longitude grid, Pangolin uses a quasi-area-preserving reduced latitude-longitude grid. The features of the regular grid are exploited to improve parallel performances and a custom domain decomposition algorithm is presented. To assess the validity of the transport scheme, its results are compared with state-of-the-art models on analytical test cases. Finally, parallel performances are shown in terms of strong scaling and confirm the efficient scalability up to a few hundred of cores.
NASA Astrophysics Data System (ADS)
El Kadi Abderrezzak, Kamal; Die Moran, Andrés; Tassi, Pablo; Ata, Riadh; Hervouet, Jean-Michel
2016-07-01
Bank erosion can be an important form of morphological adjustment in rivers. With the advances made in computational techniques, two-dimensional (2D) depth-averaged numerical models have become valuable tools for resolving many engineering problems dealing with sediment transport. The objective of this research work is to present a simple, new, bank-erosion operator that is integrated into a 2D Saint-Venant-Exner morphodynamic model. The numerical code is based on an unstructured grid of triangular elements and finite-element algorithms. The slope of each element in the grid is compared to the angle of repose of the bank material. Elements for which the slope is too steep are tilted to bring them to the angle of repose along a horizontal axis defined such that the volume loss above the axis is equal to the volume gain below, thus ensuring mass balance. The model performance is assessed using data from laboratory flume experiments and a scale model of the Old Rhine. For the flume experiment case with uniform bank material, relevant results are obtained for bank geometry changes. For the more challenging case (i.e. scale model of the Old Rhine with non-uniform bank material), the numerical model is capable of reproducing the main features of the bank failure, induced by the newly designed groynes, as well as the transport of the mobilized sediment material downstream. Some deviations between the computed results and measured data are, however, observed. They are ascribed to the effects of three-dimensional (3D) flow structures, pore pressure and cohesion, which are not considered in the present 2D model.
Kraloua, B.; Hennad, A.
2008-09-23
The aim of this paper is to determine electric and physical properties by 2D modelling of glow discharge low pressure in continuous regime maintained by term constant source. This electric discharge is confined in reactor plan-parallel geometry. This reactor is filled by Argon monatomic gas. Our continuum model the order two is composed the first three moments the Boltzmann's equations coupled with Poisson's equation by self consistent method. These transport equations are discretized by the finite volumes method. The equations system is resolved by a new technique, it is about the N-BEE explicit scheme using the time splitting method.
Voss, Clifford I.; Provost, A.M.
2002-01-01
SUTRA (Saturated-Unsaturated Transport) is a computer program that simulates fluid movement and the transport of either energy or dissolved substances in a subsurface environment. This upgraded version of SUTRA adds the capability for three-dimensional simulation to the former code (Voss, 1984), which allowed only two-dimensional simulation. The code employs a two- or three-dimensional finite-element and finite-difference method to approximate the governing equations that describe the two interdependent processes that are simulated: 1) fluid density-dependent saturated or unsaturated ground-water flow; and 2) either (a) transport of a solute in the ground water, in which the solute may be subject to: equilibrium adsorption on the porous matrix, and both first-order and zero-order production or decay; or (b) transport of thermal energy in the ground water and solid matrix of the aquifer. SUTRA may also be used to simulate simpler subsets of the above processes. A flow-direction-dependent dispersion process for anisotropic media is also provided by the code and is introduced in this report. As the primary calculated result, SUTRA provides fluid pressures and either solute concentrations or temperatures, as they vary with time, everywhere in the simulated subsurface system. SUTRA flow simulation may be employed for two-dimensional (2D) areal, cross sectional and three-dimensional (3D) modeling of saturated ground-water flow systems, and for cross sectional and 3D modeling of unsaturated zone flow. Solute-transport simulation using SUTRA may be employed to model natural or man-induced chemical-species transport including processes of solute sorption, production, and decay. For example, it may be applied to analyze ground-water contaminant transport problems and aquifer restoration designs. In addition, solute-transport simulation with SUTRA may be used for modeling of variable-density leachate movement, and for cross sectional modeling of saltwater intrusion in
A Deterministic Transport Code for Space Environment Electrons
NASA Technical Reports Server (NTRS)
Nealy, John E.; Chang, C. K.; Norman, Ryan B.; Blattnig, Steve R.; Badavi, Francis F.; Adamczyk, Anne M.
2010-01-01
A deterministic computational procedure has been developed to describe transport of space environment electrons in various shield media. This code is an upgrade and extension of an earlier electron code. Whereas the former code was formulated on the basis of parametric functions derived from limited laboratory data, the present code utilizes well established theoretical representations to describe the relevant interactions and transport processes. The shield material specification has been made more general, as have the pertinent cross sections. A combined mean free path and average trajectory approach has been used in the transport formalism. Comparisons with Monte Carlo calculations are presented.
Recent developments in the Los Alamos radiation transport code system
Forster, R.A.; Parsons, K.
1997-06-01
A brief progress report on updates to the Los Alamos Radiation Transport Code System (LARTCS) for solving criticality and fixed-source problems is provided. LARTCS integrates the Diffusion Accelerated Neutral Transport (DANT) discrete ordinates codes with the Monte Carlo N-Particle (MCNP) code. The LARCTS code is being developed with a graphical user interface for problem setup and analysis. Progress in the DANT system for criticality applications include a two-dimensional module which can be linked to a mesh-generation code and a faster iteration scheme. Updates to MCNP Version 4A allow statistical checks of calculated Monte Carlo results.
NASA Astrophysics Data System (ADS)
Ivy, D. J.; Rigby, M. L.; Prinn, R. G.; Muhle, J.; Weiss, R. F.
2009-12-01
We present optimized annual global emissions from 1973-2008 of nitrogen trifluoride (NF3), a powerful greenhouse gas which is not currently regulated by the Kyoto Protocol. In the past few decades, NF3 production has dramatically increased due to its usage in the semiconductor industry. Emissions were estimated through the 'pulse-method' discrete Kalman filter using both a simple, flexible 2-D 12-box model used in the Advanced Global Atmospheric Gases Experiment (AGAGE) network and the Model for Ozone and Related Tracers (MOZART v4.5), a full 3-D atmospheric chemistry model. No official audited reports of industrial NF3 emissions are available, and with limited information on production, a priori emissions were estimated using both a bottom-up and top-down approach with two different spatial patterns based on semiconductor perfluorocarbon (PFC) emissions from the Emission Database for Global Atmospheric Research (EDGAR v3.2) and Semiconductor Industry Association sales information. Both spatial patterns used in the models gave consistent results, showing the robustness of the estimated global emissions. Differences between estimates using the 2-D and 3-D models can be attributed to transport rates and resolution differences. Additionally, new NF3 industry production and market information is presented. Emission estimates from both the 2-D and 3-D models suggest that either the assumed industry release rate of NF3 or industry production information is still underestimated.
Advanced Nodal P_{3}/SP_{3} Axial Transport Solvers for the MPACT 2D/1D Scheme
Stimpson, Shane G; Collins, Benjamin S
2015-01-01
As part of its initiative to provide multiphysics simulations of nuclear reactor cores, the Consortium for Advanced Simulation of Light Water Reactors (CASL) is developing the Virtual Environment for Reactor Applications Core Simulator (VERA-CS). The MPACT code, which is the primary neutron transport solver of VERA-CS, employs the two-dimensional/one-dimensional (2D/1D) method to solve 3-dimensional neutron transport problems and provide sub-pin-level resolution of the power distribution. While 2D method of characteristics is used to solve for the transport effects within each plane, 1D-nodal methods are used axially. There have been extensive studies of the 2D/1D method with a variety nodal methods, and the P_{3}/SP_{3} solver has proved to be an effective method of providing higher-fidelity solutions while maintaining a low computational burden.The current implementation in MPACT wraps a one-node nodal expansion method (NEM) kernel for each moment, iterating between them and performing multiple sweeps to resolve flux distributions. However, it has been observed that this approach is more sensitive to convergence problems. This paper documents the theory and application two new nodal P_{3}/SP_{3} approaches to be used within the 2D/1D method in MPACT. These two approaches aim to provide enhanced stability compared with the pre-existing one-node approach. Results from the HY-NEM-SP_{3} solver show that the accuracy is consistent with the one-node formulations and provides improved convergence for some problems; but the solver has issues with cases in thin planes. Although the 2N-SENM-SP_{3} solver is still under development, it is intended to resolve the issues with HY-NEM-SP_{3} but it will incur some additional computational burden by necessitating an additional 1D-CMFD-P_{3} solver to generate the second moment cell-averaged scalar flux.
2D coherent charge transport in highly ordered conducting polymers doped by solid state diffusion
NASA Astrophysics Data System (ADS)
Kang, Keehoon; Watanabe, Shun; Broch, Katharina; Sepe, Alessandro; Brown, Adam; Nasrallah, Iyad; Nikolka, Mark; Fei, Zhuping; Heeney, Martin; Matsumoto, Daisuke; Marumoto, Kazuhiro; Tanaka, Hisaaki; Kuroda, Shin-Ichi; Sirringhaus, Henning
2016-08-01
Doping is one of the most important methods to control charge carrier concentration in semiconductors. Ideally, the introduction of dopants should not perturb the ordered microstructure of the semiconducting host. In some systems, such as modulation-doped inorganic semiconductors or molecular charge transfer crystals, this can be achieved by spatially separating the dopants from the charge transport pathways. However, in conducting polymers, dopants tend to be randomly distributed within the conjugated polymer, and as a result the transport properties are strongly affected by the resulting structural and electronic disorder. Here, we show that in the highly ordered lamellar microstructure of a regioregular thiophene-based conjugated polymer, a small-molecule p-type dopant can be incorporated by solid state diffusion into the layers of solubilizing side chains without disrupting the conjugated layers. In contrast to more disordered systems, this allows us to observe coherent, free-electron-like charge transport properties, including a nearly ideal Hall effect in a wide temperature range, a positive magnetoconductance due to weak localization and the Pauli paramagnetic spin susceptibility.
NASA Astrophysics Data System (ADS)
Zhao, Dongmiao; Tang, Jun; Wu, Xiuguang; Lin, Changning; Liu, Lijun; Chen, Jian
2016-05-01
A 2D vertical (2DV) numerical model, without σ-coordinate transformation in the vertical direction, is developed for the simulation of fl ow and sediment transport in open channels. In the model, time-averaged Reynolds equations are closed by the k-ɛ nonlinear turbulence model. The modifi ed Youngs-VOF method is introduced to capture free surface dynamics, and the free surface slope is simulated using the ELVIRA method. Based on the power-law scheme, the k-ɛ model and the suspended-load transport model are solved numerically with an implicit scheme applied in the vertical plane and an explicit scheme applied in the horizontal plane. Bedload transport is modeled using the Euler-WENO scheme, and the grid-closing skill is adopted to deal with the moving channel bed boundary. Verifi cation of the model using laboratory data shows that the model is able to adequately simulate fl ow and sediment transport in open channels, and is a good starting point for the study of sediment transport dynamics in strong nonlinear fl ow scenarios.
Tupala, E; Hall, H; Bergström, K; Särkioja, T; Räsänen, P; Mantere, T; Callaway, J; Hiltunen, J; Tiihonen, J
2001-05-01
Alcohol acts through mechanisms involving the brain neurotransmitter dopamine (DA) with the nucleus accumbens as the key zone for mediating these effects. We evaluated the densities of DA D(2)/D(3) receptors and transporters in the nucleus accumbens and amygdala of post-mortem human brains by using [(125)l]epidepride and [(125)I]PE2I as radioligands in whole hemispheric autoradiography of Cloninger type 1 and 2 alcoholics and healthy controls. When compared with controls, the mean binding of [(125)I]epidepride to DA D(2)/D(3) receptors was 20% lower in the nucleus accumbens and 41% lower in the amygdala, and [(125)I]PE2I binding to DA transporters in the nucleus accumbens was 39% lower in type 1 alcoholics. These data indicate that dopaminergic functions in these limbic areas may be impaired among type 1 alcoholics, due to the substantially lower number of receptor sites. Our results suggest that such a reduction may result in the chronic overuse of alcohol as an attempt to stimulate DA function. PMID:11326293
Turcksin, Bruno Ragusa, Jean C.
2014-10-01
In this paper, a Diffusion Synthetic Acceleration (DSA) technique applied to the S{sub n} radiation transport equation is developed using Piece-Wise Linear Discontinuous (PWLD) finite elements on arbitrary polygonal grids. The discretization of the DSA equations employs an Interior Penalty technique, as is classically done for the stabilization of the diffusion equation using discontinuous finite element approximations. The penalty method yields a system of linear equations that is Symmetric Positive Definite (SPD). Thus, solution techniques such as Preconditioned Conjugate Gradient (PCG) can be effectively employed. Algebraic MultiGrid (AMG) and Symmetric Gauss–Seidel (SGS) are employed as conjugate gradient preconditioners for the DSA system. AMG is shown to be significantly more efficient than SGS. Fourier analyses are carried out and we show that this discontinuous finite element DSA scheme is always stable and effective at reducing the spectral radius for iterative transport solves, even for grids with high-aspect ratio cells. Numerical results are presented for different grid types: quadrilateral, hexagonal, and polygonal grids as well as grids with local mesh adaptivity.
Coherent heat transport in 2D phononic crystals with acoustic impedance mismatch
NASA Astrophysics Data System (ADS)
Arantes, A.; Anjos, V.
2016-03-01
In this work we have calculated the cumulative thermal conductivities of micro-phononic crystals formed by different combinations of inclusions and matrices at a sub-Kelvin temperature regime. The low-frequency phonon spectra (up to tens of GHz) were obtained by solving the generalized wave equation for inhomogeneous media with the plane wave expansion method. The thermal conductivity was calculated from Boltzmann transport theory highlighting the role of the low-frequency thermal phonons and neglecting phonon-phonon scattering. A purely coherent thermal transport regime was assumed throughout the structures. Our findings show that the cumulative thermal conductivity drops dramatically when compared with their bulk counterpart. Depending on the structural composition this reduction may be attributed to the phonon group velocity due to a flattening of the phonon dispersion relation, the extinction of phonon modes in the density of states or due to the presence of complete band gaps. According to the contrast between the inclusions and the matrices, three types of two dimensional phononic crystals were considered: carbon/epoxy, carbon/polyethylene and tungsten/silicon, which correspond respectively to a moderate, strong and very strong mismatch in the mechanical properties of these materials.
High conductance 2D transport around the Hall mobility peak in electrolyte-gated rubrene crystals.
Xie, Wei; Wang, Shun; Zhang, Xin; Leighton, C; Frisbie, C Daniel
2014-12-12
We report the observation of the Hall effect at hole densities up to 6×10¹³ cm⁻² (0.3 holes/molecule) on the surface of electrolyte-gated rubrene crystals. The perplexing peak in the conductance as a function of gate voltage is confirmed to result from a maximum in mobility, which reaches 4 cm² V⁻¹ s⁻¹ at 2.5×10¹³ cm⁻². Measurements to liquid helium temperatures reveal that this peak is markedly asymmetric, with bandlike and hopping-type transport occurring on the low density side, while unconventional, likely electrostatic-disorder-affected transport dominates the high density side. Most significantly, near the mobility peak the temperature coefficient of the resistance remains positive to as low as 120 K, the low temperature resistance becomes weakly temperature dependent, and the conductance reaches within a factor of 2 of e²/h, revealing conduction unprecedentedly close to a two-dimensional metallic state. PMID:25541790
NASA Technical Reports Server (NTRS)
Hanson, Donald B.
1994-01-01
A two dimensional linear aeroacoustic theory for rotor/stator interaction with unsteady coupling was derived and explored in Volume 1 of this report. Computer program CUP2D has been written in FORTRAN embodying the theoretical equations. This volume (Volume 2) describes the structure of the code, installation and running, preparation of the input file, and interpretation of the output. A sample case is provided with printouts of the input and output. The source code is included with comments linking it closely to the theoretical equations in Volume 1.
Transport of aurorally produced N/2D/ by winds in the high latitude thermosphere
NASA Technical Reports Server (NTRS)
Gerard, J.-C.; Roble, R. G.
1982-01-01
A time-dependent, two-dimensional model is developed for describing the meridional circulation of thermospheric odd nitrogen species produced in the auroral zone. The model is based on a previous model by Roble and Gary (1979) extended to upper altitude transport of the nitrogen species. Assumptions made include the existence of a steady neutral wind flowing from low to high latitudes, and an initial background due to scattered Lyman-beta and nightglow emissions. The aurora is also assumed as steady, along with a constant ion production. Predictions made using the model are compared with observations with the Atmosphere Explorer C spacecraft and rocket sounding measurements of the 5200 A distribution near the day-side polar cusp. The model requires thermospheric winds of 100-200 m/sec, flowing from day to nightside. Convective velocities near 1000 m/sec were detected by the Explorer spacecraft, as well as a day-to-nightside flow at the cusp.
Parallel FE Electron-Photon Transport Analysis on 2-D Unstructured Mesh
Drumm, C.R.; Lorenz, J.
1999-03-02
A novel solution method has been developed to solve the coupled electron-photon transport problem on an unstructured triangular mesh. Instead of tackling the first-order form of the linear Boltzmann equation, this approach is based on the second-order form in conjunction with the conventional multi-group discrete-ordinates approximation. The highly forward-peaked electron scattering is modeled with a multigroup Legendre expansion derived from the Goudsmit-Saunderson theory. The finite element method is used to treat the spatial dependence. The solution method is unique in that the space-direction dependence is solved simultaneously, eliminating the need for the conventional inner iterations, a method that is well suited for massively parallel computers.
TORUS: Radiation transport and hydrodynamics code
NASA Astrophysics Data System (ADS)
Harries, Tim
2014-04-01
TORUS is a flexible radiation transfer and radiation-hydrodynamics code. The code has a basic infrastructure that includes the AMR mesh scheme that is used by several physics modules including atomic line transfer in a moving medium, molecular line transfer, photoionization, radiation hydrodynamics and radiative equilibrium. TORUS is useful for a variety of problems, including magnetospheric accretion onto T Tauri stars, spiral nebulae around Wolf-Rayet stars, discs around Herbig AeBe stars, structured winds of O supergiants and Raman-scattered line formation in symbiotic binaries, and dust emission and molecular line formation in star forming clusters. The code is written in Fortran 2003 and is compiled using a standard Gnu makefile. The code is parallelized using both MPI and OMP, and can use these parallel sections either separately or in a hybrid mode.
NASA Astrophysics Data System (ADS)
Moustafa, Salli; Févotte, François; Lathuilière, Bruno; Plagne, Laurent
2014-06-01
The past few years have been marked by a noticeable increase in the interest in 3D whole-core heterogeneous deterministic neutron transport solvers for reference calculations. Due to the extremely large problem sizes tackled by such solvers, they need to use adapted numerical methods and need to be efficiently implemented to take advantage of the full computing power of modern systems. As for numerical methods, one possible approach consists in iterating over resolutions of 2D and 1D MOC problems by taking advantage of prismatic geometries. The MICADO solver, developed at EDF R&D, is a parallel implementation of such a method in distributed and shared memory systems. However it is currently unable to use SIMD vectorization to leverage the full computing power of modern CPUs. In this paper, we describe our first effort to support vectorization in MICADO, typically targeting Intel© SSE CPUs. Both the 2D and 1D algorithms are vectorized, allowing for high expected speedups for the whole spatial solver. We present benchmark computations, which show nearly optimal speedups for our vectorized implementation on the TAKEDA case.
NASA Astrophysics Data System (ADS)
Croissant, T.; Lague, D.; Davy, P.
2014-12-01
Numerical models of floodplain dynamics often use a simplified 1D description of flow hydraulics and sediment transport that cannot fully account for differential friction between vegetated banks and low friction in the main channel. Key parameters of such models are the friction coefficient and the description of the channel bathymetry which strongly influence predicted water depth and velocity, and therefore sediment transport capacity. In this study, we use a newly developed 2D hydrodynamic model, Floodos, whose efficiency is a major advantage for exploring channel morphodynamics from a flood event to millennial time scales. We evaluate the quality of Floodos predictions in the Whataroa river, New Zealand and assess the effect of a spatially distributed friction coefficient (SDFC) on long term sediment transport. Predictions from the model are compared to water depth data from a gauging station located on the Whataroa River in Southern Alps, New Zealand. The Digital Elevation Model (DEM) of the 2.5 km long studied reach is derived from a 2010 LiDAR acquisition with 2 m resolution and an interpolated bathymetry. The several large floods experienced by this river during 2010 allow us to access water depth for a wide range of possible river discharges and to retrieve the scaling between these two parameters. The high resolution DEM used has a non-negligible part of submerged bathymetry that airborne LiDAR was not able to capture. Bathymetry can be reconstructed by interpolation methods that introduce several uncertainties concerning water depth predictions. We address these uncertainties inherent to the interpolation using a simplified channel with a geometry (slope and width) similar to the Whataroa river. We then explore the effect of a SDFC on velocity pattern, water depth and sediment transport capacity and discuss its relevance on long term predictions of sediment transport and channel morphodynamics.
Kim, Young-Keun; Kim, Kyung-Soo
2014-10-15
Maritime transportation demands an accurate measurement system to track the motion of oscillating container boxes in real time. However, it is a challenge to design a sensor system that can provide both reliable and non-contact methods of 6-DOF motion measurements of a remote object for outdoor applications. In the paper, a sensor system based on two 2D laser scanners is proposed for detecting the relative 6-DOF motion of a crane load in real time. Even without implementing a camera, the proposed system can detect the motion of a remote object using four laser beam points. Because it is a laser-based sensor, the system is expected to be highly robust to sea weather conditions.
NASA Astrophysics Data System (ADS)
Kim, Young-Keun; Kim, Kyung-Soo
2014-10-01
Maritime transportation demands an accurate measurement system to track the motion of oscillating container boxes in real time. However, it is a challenge to design a sensor system that can provide both reliable and non-contact methods of 6-DOF motion measurements of a remote object for outdoor applications. In the paper, a sensor system based on two 2D laser scanners is proposed for detecting the relative 6-DOF motion of a crane load in real time. Even without implementing a camera, the proposed system can detect the motion of a remote object using four laser beam points. Because it is a laser-based sensor, the system is expected to be highly robust to sea weather conditions.
Shapiro, A.B.
1983-08-01
The computer code FACET calculates the radiation geometric view factor (alternatively called shape factor, angle factor, or configuration factor) between surfaces for axisymmetric, two-dimensional planar and three-dimensional geometries with interposed third surface obstructions. FACET was developed to calculate view factors for input to finite-element heat-transfer analysis codes. The first section of this report is a brief review of previous radiation-view-factor computer codes. The second section presents the defining integral equation for the geometric view factor between two surfaces and the assumptions made in its derivation. Also in this section are the numerical algorithms used to integrate this equation for the various geometries. The third section presents the algorithms used to detect self-shadowing and third-surface shadowing between the two surfaces for which a view factor is being calculated. The fourth section provides a user's input guide followed by several example problems.
Recent advances in the Mercury Monte Carlo particle transport code
Brantley, P. S.; Dawson, S. A.; McKinley, M. S.; O'Brien, M. J.; Stevens, D. E.; Beck, B. R.; Jurgenson, E. D.; Ebbers, C. A.; Hall, J. M.
2013-07-01
We review recent physics and computational science advances in the Mercury Monte Carlo particle transport code under development at Lawrence Livermore National Laboratory. We describe recent efforts to enable a nuclear resonance fluorescence capability in the Mercury photon transport. We also describe recent work to implement a probability of extinction capability into Mercury. We review the results of current parallel scaling and threading efforts that enable the code to run on millions of MPI processes. (authors)
Benchmarking NNWSI flow and transport codes: COVE 1 results
Hayden, N.K.
1985-06-01
The code verification (COVE) activity of the Nevada Nuclear Waste Storage Investigations (NNWSI) Project is the first step in certification of flow and transport codes used for NNWSI performance assessments of a geologic repository for disposing of high-level radioactive wastes. The goals of the COVE activity are (1) to demonstrate and compare the numerical accuracy and sensitivity of certain codes, (2) to identify and resolve problems in running typical NNWSI performance assessment calculations, and (3) to evaluate computer requirements for running the codes. This report describes the work done for COVE 1, the first step in benchmarking some of the codes. Isothermal calculations for the COVE 1 benchmarking have been completed using the hydrologic flow codes SAGUARO, TRUST, and GWVIP; the radionuclide transport codes FEMTRAN and TRUMP; and the coupled flow and transport code TRACR3D. This report presents the results of three cases of the benchmarking problem solved for COVE 1, a comparison of the results, questions raised regarding sensitivities to modeling techniques, and conclusions drawn regarding the status and numerical sensitivities of the codes. 30 refs.
Kasinathan, N.; Rajakumar, A.; Vaidyanathan, G.; Chetal, S.C.
1995-09-01
Post shutdown decay heat removal is an important safety requirement in any nuclear system. In order to improve the reliability of this function, Liquid metal (sodium) cooled fast breeder reactors (LMFBR) are equipped with redundant hot pool dipped immersion coolers connected to natural draught air cooled heat exchangers through intermediate sodium circuits. During decay heat removal, flow through the core, immersion cooler primary side and in the intermediate sodium circuits are also through natural convection. In order to establish the viability and validate computer codes used in making predictions, a 1:20 scale experimental model called RAMONA with water as coolant has been built and experimental simulation of decay heat removal situation has been performed at KfK Karlsruhe. Results of two such experiments have been compiled and published as benchmarks. This paper brings out the results of the numerical simulation of one of the benchmark case through a 1D/2D coupled code system, DHDYN-1D/THYC-2D and the salient features of the comparisons. Brief description of the formulations of the codes are also included.
Benchmarking of Neutron Production of Heavy-Ion Transport Codes
Remec, Igor; Ronningen, Reginald M.; Heilbronn, Lawrence
2012-01-01
Accurate prediction of radiation fields generated by heavy ion interactions is important in medical applications, space missions, and in design and operation of rare isotope research facilities. In recent years, several well-established computer codes in widespread use for particle and radiation transport calculations have been equipped with the capability to simulate heavy ion transport and interactions. To assess and validate these capabilities, we performed simulations of a series of benchmark-quality heavy ion experiments with the computer codes FLUKA, MARS15, MCNPX, and PHITS. We focus on the comparisons of secondary neutron production. Results are encouraging; however, further improvements in models and codes and additional benchmarking are required.
Bounce-averaged Fokker-Planck code for stellarator transport
Mynick, H.E.; Hitchon, W.N.G.
1985-07-01
A computer code for solving the bounce-averaged Fokker-Planck equation appropriate to stellarator transport has been developed, and its first applications made. The code is much faster than the bounce-averaged Monte-Carlo codes, which up to now have provided the most efficient numerical means for studying stellarator transport. Moreover, because the connection to analytic kinetic theory of the Fokker-Planck approach is more direct than for the Monte-Carlo approach, a comparison of theory and numerical experiment is now possible at a considerably more detailed level than previously.
NASA Astrophysics Data System (ADS)
Mani, Ramesh; Kriisa, A.
2015-03-01
Negative diagonal magneto-conductivity/resistivity is a spectacular- and thought provoking- property of driven, far-from-equilibrium, low dimensional electronic systems. The physical response of this exotic electronic state is not yet fully understood since it is rarely encountered in experiment. The microwave-radiation-induced zero-resistance state in the high mobility GaAs/AlGaAs 2D electron system is believed to be an example where negative magneto-conductivity/resistivity is responsible for the observed phenomena. Here, we examine the magneto-transport characteristics of this negative conductivity/resistivity state in the microwave photo-excited two-dimensional electron system (2DES) through a numerical solution of the associated boundary value problem. The results suggest, surprisingly, that a bare negative diagonal conductivity/resistivity state in the 2DES under photo-excitation should yield a positive diagonal resistance with a concomitant sign reversal in the Hall voltage. Transport measurements are supported by the DOE, Office of Basic Energy Sciences, Material Sciences and Engineering Division under DE-SC0001762. Additional support by the ARO under W911NF-07-01-015.
NASA Astrophysics Data System (ADS)
Ji, Jianhua; Wang, Yanfen; Wang, Ke; Xu, Ming; Zhang, Zhipeng; Yang, Shuwen
2013-09-01
A new two-dimensional OOC (optical orthogonal codes) named PC/MOCS is constructed, using PC (prime code) for time spreading and MOCS (modified one-coincidence sequence) for wavelength hopping. Compared with PC/PC, the number of wavelengths for PC/MOCS is not limited to a prime number. Compared with PC/OCS, the length of MOCS need not be expanded to the same length of PC. PC/MOCS can be constructed flexibly, and also can use available wavelengths effectively. Theoretical analysis shows that PC/MOCS can reduce the bit error rate (BER) of OCDMA system, and can support more users than PC/PC and PC/OCS.
TRAC code assessment using data from SCTF Core-III, a large-scale 2D/3D facility
Boyack, B.E.; Shire, P.R.; Harmony, S.C.; Rhee, G.
1988-01-01
Nine tests from the SCTF Core-III configuration have been analyzed using TRAC-PF1/MOD1. The objectives of these assessment activities were to obtain a better understanding of the phenomena occurring during the refill and reflood phases of a large-break loss-of-coolant accident, to determine the accuracy to which key parameters are calculated, and to identify deficiencies in key code correlations and models that provide closure for the differential equations defining thermal-hydraulic phenomena in pressurized water reactors. Overall, the agreement between calculated and measured values of peak cladding temperature is reasonable. In addition, TRAC adequately predicts many of the trends observed in both the integral effect and separate effect tests conducted in SCTF Core-III. The importance of assessment activities that consider potential contributors to discrepancies between the measured and calculated results arising from three sources are described as those related to (1) knowledge about the facility configuration and operation, (2) facility modeling for code input, and (3) deficiencies in code correlations and models. An example is provided. 8 refs., 7 figs., 2 tabs.
ZORNOC: a 1 1/2-D tokamak data analysis code for studying noncircular high beta plasmas
Zurro, B.; Wieland, R.M.; Murakami, M.; Swain, D.W.
1980-03-01
A new tokamak data analysis code, ZORNOC, was developed to study noncircular, high beta plasmas in the Impurity Study Experiment (ISX-B). These plasmas exhibit significant flux surface shifts and elongation in both ohmically heated and beam-heated discharges. The MHD equilibrium flux surface geometry is determined by solving the Grad-Shafranov equation based on: (1) the shape of the outermost flux surface, deduced from the magnetic loop probes; (2) a pressure profile, deduced by means of Thomson scattering data (electrons), charge exchange data (ions), and a Fokker-Planck model (fast ions); and (3) a safety factor profile, determined from the experimental data using a simple model (Z/sub eff/ = const) that is self-consistently altered while the plasma equilibrium is iterated. For beam-heated discharches the beam deposition profile is determined by means of a Monte Carlo scheme and the slowing down of the fast ions by means of an analytical solution of the Fokker-Planck equation. The code also carries out an electron power balance and calculates various confinement parameters. The code is described and examples of its operation are given.
Litvinenko, I. A.; Lykov, V. A.
1997-04-15
The results of numerical simulation of fast electrons motion and generated electro-magnetic fields at the picosecond pulse laser interaction with flat target are presented. The calculations were performed with PM2D code, where relativistic equation of electron motion joint with Maxwell equations is solved by particle method in cells. The efficiency of fast electrons energy conversion to the transverse electromagnetic wave of picosecond duration can reach the value 10{sup -4} for the intensity of ultrashort laser pulse at the target 10{sup 16}-10{sup 17} W/cm{sup 2}.
Jones, Edward H; Smith, Colin C
2005-12-01
This paper describes an investigation into non-equilibrium partitioning tracer transport and interaction with non-aqueous-phase liquid (NAPL) contaminated water-saturated porous media using a two-dimensional (2-D) physical modelling methodology. A fluorescent partitioning tracer is employed within a transparent porous model which when imaged by a CCD digital camera can provide full spatial tracer concentrations and tracer breakthrough curves. Quasi one-dimensional (1-D) benchmarking tests in models packed with various combinations of clean quartz sand and NAPL are described. These modelled residual NAPL saturations, S(n), of 0-15%. Results demonstrated that the fluorescent partitioning tracer was able to detect and quantify the presence of NAPL at low flow rates. At larger flow rates and/or higher NAPL saturations, the tracer increasingly underpredicted the NAPL volume as expected and this is attributed primarily to non-equilibrium partitioning. Despite little change in permeability, change in NAPL saturations from 4% to 8% resulted in significant NAPL saturation underestimates at the same flow rates implying coalescence of NAPL into wider separated but larger ganglia. A 2-D investigation of an idealised heterogeneous residual NAPL contaminated flow field indicated little permeability change in the NAPL contaminated zone and thus little flow bypassing, leading to reduced underpredictions of NAPL saturations than for equivalent quasi 1-D cases. This was attributed to increased 'sampling' of the NAPL by the tracer. The process is clearly visually identifiable from the experimental images. This rapid and relatively inexpensive experimental method is of value in laboratory studies of partitioning tracer behaviour in porous media; in particular, the ability to observe full field concentrations makes it valuable for the study of complex heterogeneous systems. PMID:16298415
Current status of the PSG Monte Carlo neutron transport code
Leppaenen, J.
2006-07-01
PSG is a new Monte Carlo neutron transport code, developed at the Technical Research Centre of Finland (VTT). The code is mainly intended for fuel assembly-level reactor physics calculations, such as group constant generation for deterministic reactor simulator codes. This paper presents the current status of the project and the essential capabilities of the code. Although the main application of PSG is in lattice calculations, the geometry is not restricted in two dimensions. This paper presents the validation of PSG against the experimental results of the three-dimensional MOX fuelled VENUS-2 reactor dosimetry benchmark. (authors)
Transport Code for Regular Triangular Geometry
Energy Science and Technology Software Center (ESTSC)
1993-06-09
DIAMANT2 solves the two-dimensional static multigroup neutron transport equation in planar regular triangular geometry. Both regular and adjoint, inhomogeneous and homogeneous problems subject to vacuum, reflective or input specified boundary flux conditions are solved. Anisotropy is allowed for the scattering source. Volume and surface sources are allowed for inhomogeneous problems.
Galián, Carmen; Manon, Florence; Dezi, Manuela; Torres, Cristina; Ebel, Christine; Lévy, Daniel; Jault, Jean-Michel
2011-01-01
Optimized protocols for achieving high-yield expression, purification and reconstitution of membrane proteins are required to study their structure and function. We previously reported high-level expression in Escherichia coli of active BmrC and BmrD proteins from Bacillus subtilis, previously named YheI and YheH. These proteins are half-transporters which belong to the ABC (ATP-Binding Cassette) superfamily and associate in vivo to form a functional transporter able to efflux drugs. In this report, high-yield purification and functional reconstitution were achieved for the heterodimer BmrC/BmrD. In contrast to other detergents more efficient for solubilizing the transporter, dodecyl-ß-D-maltoside (DDM) maintained it in a drug-sensitive and vanadate-sensitive ATPase-competent state after purification by affinity chromatography. High amounts of pure proteins were obtained which were shown either by analytical ultracentrifugation or gel filtration to form a monodisperse heterodimer in solution, which was notably stable for more than one month at 4°C. Functional reconstitution using different lipid compositions induced an 8-fold increase of the ATPase activity (kcat∼5 s−1). We further validated that the quality of the purified BmrC/BmrD heterodimer is suitable for structural analyses, as its reconstitution at high protein densities led to the formation of 2-D crystals. Electron microscopy of negatively stained crystals allowed the calculation of a projection map at 20 Å resolution revealing that BmrC/BmrD might assemble into oligomers in a lipidic environment. PMID:21602923
Galián, Carmen; Manon, Florence; Dezi, Manuela; Torres, Cristina; Ebel, Christine; Lévy, Daniel; Jault, Jean-Michel
2011-01-01
Optimized protocols for achieving high-yield expression, purification and reconstitution of membrane proteins are required to study their structure and function. We previously reported high-level expression in Escherichia coli of active BmrC and BmrD proteins from Bacillus subtilis, previously named YheI and YheH. These proteins are half-transporters which belong to the ABC (ATP-Binding Cassette) superfamily and associate in vivo to form a functional transporter able to efflux drugs. In this report, high-yield purification and functional reconstitution were achieved for the heterodimer BmrC/BmrD. In contrast to other detergents more efficient for solubilizing the transporter, dodecyl-ß-D-maltoside (DDM) maintained it in a drug-sensitive and vanadate-sensitive ATPase-competent state after purification by affinity chromatography. High amounts of pure proteins were obtained which were shown either by analytical ultracentrifugation or gel filtration to form a monodisperse heterodimer in solution, which was notably stable for more than one month at 4°C. Functional reconstitution using different lipid compositions induced an 8-fold increase of the ATPase activity (k(cat)∼5 s(-1)). We further validated that the quality of the purified BmrC/BmrD heterodimer is suitable for structural analyses, as its reconstitution at high protein densities led to the formation of 2-D crystals. Electron microscopy of negatively stained crystals allowed the calculation of a projection map at 20 Å resolution revealing that BmrC/BmrD might assemble into oligomers in a lipidic environment. PMID:21602923
Flow and Containment Transport Code for Modeling Variably Saturated Porous Media
Energy Science and Technology Software Center (ESTSC)
1998-05-14
FACT is a finite element based code designed to model subsurface flow and contaminant transport. It was designed to perform transient three-dimensional calculations that simulate isothermal groundwater flow, moisture movement, and solute transport in variably saturated and fully saturated subsurface porous media. The code is designed specifically to handle complex multi-layer and/or heterogenous aquifer systems in an efficient manner and accommodates a wide range of boundary conditions. Additionally 1-D and 2-D (in Cartesian coordinates) problemsmore » are handled in FACT by simply limiting the number of elements in a particular direction(s) to one. The governing equations in FACT are formulated only in Cartesian coordinates. FACT writes out both ascii and graphical binary files that are TECPLOT-ready. Special features are also available within FACT for handling the typical groundwater modeling needs for remediation efforts at the Savannah River Site.« less
Douard, Veronique; Sabbagh, Yves; Lee, Jacklyn; Patel, Chirag; Kemp, Francis W.; Bogden, John D.; Lin, Sheldon
2013-01-01
We recently discovered that chronic high fructose intake by lactating rats prevented adaptive increases in rates of active intestinal Ca2+ transport and in levels of 1,25-(OH)2D3, the active form of vitamin D. Since sufficient Ca2+ absorption is essential for skeletal growth, our discovery may explain findings that excessive consumption of sweeteners compromises bone integrity in children. We tested the hypothesis that 1,25-(OH)2D3 mediates the inhibitory effect of excessive fructose intake on active Ca2+ transport. First, compared with those fed glucose or starch, growing rats fed fructose for 4 wk had a marked reduction in intestinal Ca2+ transport rate as well as in expression of intestinal and renal Ca2+ transporters that was tightly associated with decreases in circulating levels of 1,25-(OH)2D3, bone length, and total bone ash weight but not with serum parathyroid hormone (PTH). Dietary fructose increased the expression of 24-hydroxylase (CYP24A1) and decreased that of 1α-hydroxylase (CYP27B1), suggesting that fructose might enhance the renal catabolism and impair the synthesis, respectively, of 1,25-(OH)2D3. Serum FGF23, which is secreted by osteocytes and inhibits CYP27B1 expression, was upregulated, suggesting a potential role of bone in mediating the fructose effects on 1,25-(OH)2D3 synthesis. Second, 1,25-(OH)2D3 treatment rescued the fructose effect and normalized intestinal and renal Ca2+ transporter expression. The mechanism underlying the deleterious effect of excessive fructose intake on intestinal and renal Ca2+ transporters is a reduction in serum levels of 1,25-(OH)2D3. This finding is significant because of the large amounts of fructose now consumed by Americans increasingly vulnerable to Ca2+ and vitamin D deficiency. PMID:23571713
DETERMINISTIC TRANSPORT METHODS AND CODES AT LOS ALAMOS
J. E. MOREL
1999-06-01
The purposes of this paper are to: Present a brief history of deterministic transport methods development at Los Alamos National Laboratory from the 1950's to the present; Discuss the current status and capabilities of deterministic transport codes at Los Alamos; and Discuss future transport needs and possible future research directions. Our discussion of methods research necessarily includes only a small fraction of the total research actually done. The works that have been included represent a very subjective choice on the part of the author that was strongly influenced by his personal knowledge and experience. The remainder of this paper is organized in four sections: the first relates to deterministic methods research performed at Los Alamos, the second relates to production codes developed at Los Alamos, the third relates to the current status of transport codes at Los Alamos, and the fourth relates to future research directions at Los Alamos.
HT2DINV: A 2D forward and inverse code for steady-state and transient hydraulic tomography problems
NASA Astrophysics Data System (ADS)
Soueid Ahmed, A.; Jardani, A.; Revil, A.; Dupont, J. P.
2015-12-01
Hydraulic tomography is a technique used to characterize the spatial heterogeneities of storativity and transmissivity fields. The responses of an aquifer to a source of hydraulic stimulations are used to recover the features of the estimated fields using inverse techniques. We developed a 2D free source Matlab package for performing hydraulic tomography analysis in steady state and transient regimes. The package uses the finite elements method to solve the ground water flow equation for simple or complex geometries accounting for the anisotropy of the material properties. The inverse problem is based on implementing the geostatistical quasi-linear approach of Kitanidis combined with the adjoint-state method to compute the required sensitivity matrices. For undetermined inverse problems, the adjoint-state method provides a faster and more accurate approach for the evaluation of sensitivity matrices compared with the finite differences method. Our methodology is organized in a way that permits the end-user to activate parallel computing in order to reduce the computational burden. Three case studies are investigated demonstrating the robustness and efficiency of our approach for inverting hydraulic parameters.
Modeling of Anomalous Transport in Tokamaks with FACETS code
NASA Astrophysics Data System (ADS)
Pankin, A. Y.; Batemann, G.; Kritz, A.; Rafiq, T.; Vadlamani, S.; Hakim, A.; Kruger, S.; Miah, M.; Rognlien, T.
2009-05-01
The FACETS code, a whole-device integrated modeling code that self-consistently computes plasma profiles for the plasma core and edge in tokamaks, has been recently developed as a part of the SciDAC project for core-edge simulations. A choice of transport models is available in FACETS through the FMCFM interface [1]. Transport models included in FMCFM have specific ranges of applicability, which can limit their use to parts of the plasma. In particular, the GLF23 transport model does not include the resistive ballooning effects that can be important in the tokamak pedestal region and GLF23 typically under-predicts the anomalous fluxes near the magnetic axis [2]. The TGLF and GYRO transport models have similar limitations [3]. A combination of transport models that covers the entire discharge domain is studied using FACETS in a realistic tokamak geometry. Effective diffusivities computed with the FMCFM transport models are extended to the region near the separatrix to be used in the UEDGE code within FACETS. 1. S. Vadlamani et al. (2009) %First time-dependent transport simulations using GYRO and NCLASS within FACETS (this meeting).2. T. Rafiq et al. (2009) %Simulation of electron thermal transport in H-mode discharges Submitted to Phys. Plasmas.3. C. Holland et al. (2008) %Validation of gyrokinetic transport simulations using %DIII-D core turbulence measurements Proc. of IAEA FEC (Switzerland, 2008)
A New Cell-Centered Implicit Numerical Scheme for Ions in the 2-D Axisymmetric Code Hall2de
NASA Technical Reports Server (NTRS)
Lopez Ortega, Alejandro; Mikellides, Ioannis G.
2014-01-01
We present a new algorithm in the Hall2De code to simulate the ion hydrodynamics in the acceleration channel and near plume regions of Hall-effect thrusters. This implementation constitutes an upgrade of the capabilities built in the Hall2De code. The equations of mass conservation and momentum for unmagnetized ions are solved using a conservative, finite-volume, cell-centered scheme on a magnetic-field-aligned grid. Major computational savings are achieved by making use of an implicit predictor/multi-corrector algorithm for time evolution. Inaccuracies in the prediction of the motion of low-energy ions in the near plume in hydrodynamics approaches are addressed by implementing a multi-fluid algorithm that tracks ions of different energies separately. A wide range of comparisons with measurements are performed to validate the new ion algorithms. Several numerical experiments with the location and value of the anomalous collision frequency are also presented. Differences in the plasma properties in the near-plume between the single fluid and multi-fluid approaches are discussed. We complete our validation by comparing predicted erosion rates at the channel walls of the thruster with measurements. Erosion rates predicted by the plasma properties obtained from simulations replicate accurately measured rates of erosion within the uncertainty range of the sputtering models employed.
Radiation transport phenomena and modeling - part A: Codes
Lorence, L.J.
1997-06-01
The need to understand how particle radiation (high-energy photons and electrons) from a variety of sources affects materials and electronics has motivated the development of sophisticated computer codes that describe how radiation with energies from 1.0 keV to 100.0 GeV propagates through matter. Predicting radiation transport is the necessary first step in predicting radiation effects. The radiation transport codes that are described here are general-purpose codes capable of analyzing a variety of radiation environments including those produced by nuclear weapons (x-rays, gamma rays, and neutrons), by sources in space (electrons and ions) and by accelerators (x-rays, gamma rays, and electrons). Applications of these codes include the study of radiation effects on electronics, nuclear medicine (imaging and cancer treatment), and industrial processes (food disinfestation, waste sterilization, manufacturing.) The primary focus will be on coupled electron-photon transport codes, with some brief discussion of proton transport. These codes model a radiation cascade in which electrons produce photons and vice versa. This coupling between particles of different types is important for radiation effects. For instance, in an x-ray environment, electrons are produced that drive the response in electronics. In an electron environment, dose due to bremsstrahlung photons can be significant once the source electrons have been stopped.
MORSE Monte Carlo radiation transport code system
Emmett, M.B.
1983-02-01
This report is an addendum to the MORSE report, ORNL-4972, originally published in 1975. This addendum contains descriptions of several modifications to the MORSE Monte Carlo Code, replacement pages containing corrections, Part II of the report which was previously unpublished, and a new Table of Contents. The modifications include a Klein Nishina estimator for gamma rays. Use of such an estimator required changing the cross section routines to process pair production and Compton scattering cross sections directly from ENDF tapes and writing a new version of subroutine RELCOL. Another modification is the use of free form input for the SAMBO analysis data. This required changing subroutines SCORIN and adding new subroutine RFRE. References are updated, and errors in the original report have been corrected. (WHK)
Monte-Carlo Impurity transport simulations in the edge of the DIII-D tokamak using the MCI code
Evans, T.E.; Mahdavi, M.A.; Sager, G.T.; West, W.P.; Fenstermacher, M.E.; Meyer, W.H.; Porter, G.D.
1995-07-01
A Monte-Carlo Impurity (MCI) transport code is used to follow trace impurities through multiple ionization states in realistic 2-D tokamak geometries. The MCI code is used to study impurity transport along the open magnetic field lines of the Scrape-off Layer (SOL) and to understand how impurities get into the core from the SOL. An MCI study concentrating on the entrainment of carbon impurities ions by deuterium background plasma into the DIII-D divertor is discussed. MCI simulation results are compared to experimental DIII-D carbon measurements.
SCTP as scalable video coding transport
NASA Astrophysics Data System (ADS)
Ortiz, Jordi; Graciá, Eduardo Martínez; Skarmeta, Antonio F.
2013-12-01
This study presents an evaluation of the Stream Transmission Control Protocol (SCTP) for the transport of the scalable video codec (SVC), proposed by MPEG as an extension to H.264/AVC. Both technologies fit together properly. On the one hand, SVC permits to split easily the bitstream into substreams carrying different video layers, each with different importance for the reconstruction of the complete video sequence at the receiver end. On the other hand, SCTP includes features, such as the multi-streaming and multi-homing capabilities, that permit to transport robustly and efficiently the SVC layers. Several transmission strategies supported on baseline SCTP and its concurrent multipath transfer (CMT) extension are compared with the classical solutions based on the Transmission Control Protocol (TCP) and the Realtime Transmission Protocol (RTP). Using ns-2 simulations, it is shown that CMT-SCTP outperforms TCP and RTP in error-prone networking environments. The comparison is established according to several performance measurements, including delay, throughput, packet loss, and peak signal-to-noise ratio of the received video.
Overview of development and design of MPACT: Michigan parallel characteristics transport code
Kochunas, B.; Collins, B.; Jabaay, D.; Downar, T. J.; Martin, W. R.
2013-07-01
MPACT (Michigan Parallel Characteristics Transport Code) is a new reactor analysis tool. It is being developed by students and research staff at the University of Michigan to be used for an advanced pin-resolved transport capability within VERA (Virtual Environment for Reactor Analysis). VERA is the end-user reactor simulation tool being produced by the Consortium for the Advanced Simulation of Light Water Reactors (CASL). The MPACT development project is itself unique for the way it is changing how students do research to achieve the instructional and research goals of an academic institution, while providing immediate value to industry. The MPACT code makes use of modern lean/agile software processes and extensive testing to maintain a level of productivity and quality required by CASL. MPACT's design relies heavily on object-oriented programming concepts and design patterns and is programmed in Fortran 2003. These designs are explained and illustrated as to how they can be readily extended to incorporate new capabilities and research ideas in support of academic research objectives. The transport methods currently implemented in MPACT include the 2-D and 3-D method of characteristics (MOC) and 2-D and 3-D method of collision direction probabilities (CDP). For the cross section resonance treatment, presently the subgroup method and the new embedded self-shielding method (ESSM) are implemented within MPACT. (authors)
A MARS-based method for estimating regional 2-D ionospheric VTEC and receiver differential code bias
NASA Astrophysics Data System (ADS)
Kao, Szu-Pyng; Chen, Yao-Chung; Ning, Fang-Shii
2014-01-01
The geometry-free linear combination of dual-frequency GNSS reference station ground observations are currently used to build the Vertical Total Electron Content (VTEC) model of the ionosphere. As it is known, besides ionospheric delays, there are differential code bias (DCB) of satellite (SDCB) and receiver (RDCB) in the geometry-free observation equation. The SDCB can be obtained using the International GNSS Service (IGS) analysis centers, but the RDCB for regional and local network receivers are not provided. Therefore, estimating the RDCB and VTEC model accurately and simultaneously is a critical factor investigated by researchers. This study uses Multivariate Adaptive Regression Splines (MARS) to estimate the VTEC approximate model and then substitutes this model in the observation equation to form the normal equation. The least squares method is used to solve the RDCB and VTEC model together. The research findings show that this method has good modeling effectiveness and the estimated RDCB has good reliability. The estimated VTEC model applied to GPS single-frequency precise point positioning has better positioning accuracy in comparison to the IGS global ionosphere map (GIM).
Update on the Development and Validation of MERCURY: A Modern, Monte Carlo Particle Transport Code
Procassini, R J; Taylor, J M; McKinley, M S; Greenman, G M; Cullen, D E; O'Brien, M J; Beck, B R; Hagmann, C A
2005-06-06
An update on the development and validation of the MERCURY Monte Carlo particle transport code is presented. MERCURY is a modern, parallel, general-purpose Monte Carlo code being developed at the Lawrence Livermore National Laboratory. During the past year, several major algorithm enhancements have been completed. These include the addition of particle trackers for 3-D combinatorial geometry (CG), 1-D radial meshes, 2-D quadrilateral unstructured meshes, as well as a feature known as templates for defining recursive, repeated structures in CG. New physics capabilities include an elastic-scattering neutron thermalization model, support for continuous energy cross sections and S ({alpha}, {beta}) molecular bound scattering. Each of these new physics features has been validated through code-to-code comparisons with another Monte Carlo transport code. Several important computer science features have been developed, including an extensible input-parameter parser based upon the XML data description language, and a dynamic load-balance methodology for efficient parallel calculations. This paper discusses the recent work in each of these areas, and describes a plan for future extensions that are required to meet the needs of our ever expanding user base.
Multidimensional electron-photon transport with standard discrete ordinates codes
Drumm, C.R.
1995-12-31
A method is described for generating electron cross sections that are compatible with standard discrete ordinates codes without modification. There are many advantages of using an established discrete ordinates solver, e.g. immediately available adjoint capability. Coupled electron-photon transport capability is needed for many applications, including the modeling of the response of electronics components to space and man-made radiation environments. The cross sections have been successfully used in the DORT, TWODANT and TORT discrete ordinates codes. The cross sections are shown to provide accurate and efficient solutions to certain multidimensional electronphoton transport problems.
RADTRAN 5: A computer code for transportation risk analysis
Neuhauser, K. S.; Kanipe, F. L.
1991-01-01
RADTRAN 5 is a computer code developed at Sandia National Laboratories (SNL) in Albuquerque, NM, to estimate radiological and nonradiological risks of radioactive materials transportation. RADTRAN 5 is written in ANSI Standard FORTRAN 77 and contains significant advances in the methodology for route-specific analysis first developed by SNL for RADTRAN 4 (Neuhauser and Kanipe, 1992). Like the previous RADTRAN codes, RADTRAN 5 contains two major modules for incident-free and accident risk amlysis, respectively. All commercially important transportation modes may be analyzed with RADTRAN 5: highway by combination truck; highway by light-duty vehicle; rail; barge; ocean-going ship; cargo air; and passenger air.
Validation of a coupled reactive transport code in porous media
NASA Astrophysics Data System (ADS)
Mugler, C.; Montarnal, P.; Dimier, A.
2003-04-01
The safety assessment of nuclear waste disposals needs to predict the migration of radionuclides and chemical species through a geological medium. It is therefore necessary to develop and assess qualified and validated tools which integrate both the transport mechanisms through the geological media and the chemical mechanisms governing the mobility of radionuclides. In this problem, both geochemical and hydrodynamic phenomena are tightly linked together. That is the reason why the French Nuclear Energy Agency (CEA) and the French Agency for the Management of Radioactive Wastes (ANDRA) are conjointly developping a coupled reactive transport code that solves simultaneously a geochemical model and a transport model. This code, which is part of the software project ALLIANCES, leans on the libraries of two geochemical codes solving the complex ensemble of reacting chemical species: CHESS and PHREEQC. Geochemical processes considered here include ion exchange, redox reactions, acid-base reactions, surface complexation and mineral dissolution and/or precipitation. Transport is simulated using the mixed-hybrid finite element scheme CAST3M or the finite volume scheme MT3D. All together solve Darcy's law and simulate several hydrological processes such as advection, diffusion and dispersion. The coupling algorithm is an iterative sequential algorithm. Several analytical test cases have been defined and used to validate the reactive transport code. Numerical results can be compared to analytical solutions.
NASA Astrophysics Data System (ADS)
Berger, Richard; Chapman, T.; Banks, J. W.; Brunner, S.
2015-11-01
We present 2D+2V Vlasov simulations of Ion Acoustic waves (IAWs) driven by an external traveling-wave potential, ϕ0 (x , t) , with frequency, ω, and wavenumber, k, obeying the kinetic dispersion relation. Both electrons and ions are treated kinetically. Simulations with ϕ0 (x , t) , localized transverse to the propagation direction, model IAWs driven in a laser speckle. The waves bow with a positive or negative curvature of the wave fronts that depends on the sign of the nonlinear frequency shift ΔωNL , which is in turn determined by the magnitude of ZTe /Ti where Z is the charge state and Te , i is the electron, ion temperature. These kinetic effects result can cause modulational and self-focusing instabilities that transfer wave energy to kinetic energy. Linear dispersion properties of IAWs are used in laser propagation codes that predict the amount of light reflected by stimulated Brillouin scattering. At high enough amplitudes, the linear dispersion is invalid and these kinetic effects should be incorporated. Including the spatial and time scales of these instabilities is computationally prohibitive. We report progress including kinetic models in laser propagation codes. This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract No. DE-AC52-07NA27344 and funded by the Laboratory Research and Development Program at LLNL under project tracking code 15.
Kraidith, Kamonshanok; Svasti, Saovaros; Teerapornpuntakit, Jarinthorn; Vadolas, Jim; Chaimana, Rattana; Lapmanee, Sarawut; Suntornsaratoon, Panan; Krishnamra, Nateetip; Fucharoen, Suthat; Charoenphandhu, Narattaphol
2016-07-01
Previously, β-thalassemia, an inherited anemic disorder with iron overload caused by loss-of-function mutation of β-globin gene, has been reported to induce osteopenia and impaired whole body calcium metabolism, but the pathogenesis of aberrant calcium homeostasis remains elusive. Herein, we investigated how β-thalassemia impaired intestinal calcium absorption and whether it could be restored by administration of 1,25-dihydroxyvitamin D3 [1,25(OH)2D3] or hepcidin, the latter of which was the liver-derived antagonist of intestinal iron absorption. The results showed that, in hemizygous β-globin knockout (BKO) mice, the duodenal calcium transport was lower than that in wild-type littermates, and severity was especially pronounced in female mice. Both active and passive duodenal calcium fluxes in BKO mice were found to be less than those in normal mice. This impaired calcium transport could be restored by 7-day 1,25(OH)2D3 treatment. The 1,25(OH)2D3-induced calcium transport was diminished by inhibitors of calcium transporters, e.g., L-type calcium channel, NCX1, and PMCA1b, as well as vesicular transport inhibitors. Interestingly, the duodenal calcium transport exhibited an inverse correlation with transepithelial iron transport, which was markedly enhanced in thalassemic mice. Thus, 3-day subcutaneous hepcidin injection and acute direct hepcidin exposure in the Ussing chamber were capable of restoring the thalassemia-associated impairment of calcium transport; however, the positive effect of hepcidin on calcium transport was completely blocked by proteasome inhibitors MG132 and bortezomib. In conclusion, both 1,25(OH)2D3 and hepcidin could be used to alleviate the β-thalassemia-associated impairment of calcium absorption. Therefore, our study has shed light on the development of a treatment strategy to rescue calcium dysregulation in β-thalassemia. PMID:27245334
NASA Astrophysics Data System (ADS)
Atzeni, Stefano; Marocchino, Alberto; Schiavi, Angelo
2016-03-01
Accurate descriptions of laser power coupling to the plasma and electron energy transport are crucial for designing shock-ignition targets and assessing their robustness (in particular with regard to laser and positioning errors). To this purpose, the 2D DUED laser fusion code has been improved with the inclusion of a 3D laser ray-tracing scheme and a model for non-local electron transport. 2D simulations with the upgraded code are presented; the dependence of the fusion yield vs target displacement is studied. Two different irradiation configurations are considered.
Differential Cross Section Kinematics for 3-dimensional Transport Codes
NASA Technical Reports Server (NTRS)
Norbury, John W.; Dick, Frank
2008-01-01
In support of the development of 3-dimensional transport codes, this paper derives the relevant relativistic particle kinematic theory. Formulas are given for invariant, spectral and angular distributions in both the lab (spacecraft) and center of momentum frames, for collisions involving 2, 3 and n - body final states.
The EGS4 Code System: Solution of Gamma-ray and Electron Transport Problems
DOE R&D Accomplishments Database
Nelson, W. R.; Namito, Yoshihito
1990-03-01
In this paper we present an overview of the EGS4 Code System -- a general purpose package for the Monte Carlo simulation of the transport of electrons and photons. During the last 10-15 years EGS has been widely used to design accelerators and detectors for high-energy physics. More recently the code has been found to be of tremendous use in medical radiation physics and dosimetry. The problem-solving capabilities of EGS4 will be demonstrated by means of a variety of practical examples. To facilitate this review, we will take advantage of a new add-on package, called SHOWGRAF, to display particle trajectories in complicated geometries. These are shown as 2-D laser pictures in the written paper and as photographic slides of a 3-D high-resolution color monitor during the oral presentation. 11 refs., 15 figs.
The EGS4 Code System: Solution of gamma-ray and electron transport problems
Nelson, W.R.; Namito, Yoshihito.
1990-02-09
In this paper we present an overview of the EGS4 Code System -- a general purpose package for the Monte Carlo simulation of the transport of electrons and photons. During the last 10-15 years EGS has been widely used to design accelerators and detectors for high-energy physics. More recently the code has been found to be of tremendous use in medical radiation physics and dosimetry. The problem-solving capabilities of EGS4 will be demonstrated by means of a variety of practical examples. To facilitate this review, we will take advantage of a new add-on package, called SHOWGRAF, to display particle trajectories in complicated geometries. These are shown as 2-D laser pictures in the written paper and as photographic slides of a 3-D high-resolution color monitor during the oral presentation. 11 refs., 15 figs.
Description of transport codes for space radiation shielding.
Kim, Myung-Hee Y; Wilson, John W; Cucinotta, Francis A
2012-11-01
Exposure to ionizing radiation in the space environment is one of the hazards faced by crews in space missions. As space radiations traverse spacecraft, habitat shielding, or tissues, their energies and compositions are altered by interactions with the shielding. Modifications to the radiation fields arise from atomic interactions of charged particles with orbital electrons and nuclear interactions leading to projectile and target fragmentation, including secondary particles such as neutrons, protons, mesons, and nuclear recoils. The transport of space radiation through shielding can be simulated using Monte Carlo techniques or deterministic solutions of the Boltzmann equation. To determine shielding requirements and to resolve radiation constraints for future human missions, the shielding evaluation of a spacecraft concept is required as an early step in the design process. To do this requires (1) accurate knowledge of space environmental models to define the boundary condition for transport calculations, (2) transport codes with detailed shielding and body geometry models to determine particle transmission into areas of internal shielding and at each critical body organ, and (3) the assessment of organ dosimetric quantities and biological risks by applying the corresponding response models for space radiation against the particle spectra that have been accurately determined from the transport code. This paper reviews current transport codes and analyzes their accuracy through comparison to laboratory and spaceflight data. This paper also introduces a probabilistic risk assessment approach for the evaluation of radiation shielding. PMID:23032892
Energy Science and Technology Software Center (ESTSC)
2005-07-01
Aniso2d is a two-dimensional seismic forward modeling code. The earth is parameterized by an X-Z plane in which the seismic properties Can have monoclinic with x-z plane symmetry. The program uses a user define time-domain wavelet to produce synthetic seismograms anrwhere within the two-dimensional media.
Boltzmann Transport Code Update: Parallelization and Integrated Design Updates
NASA Technical Reports Server (NTRS)
Heinbockel, J. H.; Nealy, J. E.; DeAngelis, G.; Feldman, G. A.; Chokshi, S.
2003-01-01
The on going efforts at developing a web site for radiation analysis is expected to result in an increased usage of the High Charge and Energy Transport Code HZETRN. It would be nice to be able to do the requested calculations quickly and efficiently. Therefore the question arose, "Could the implementation of parallel processing speed up the calculations required?" To answer this question two modifications of the HZETRN computer code were created. The first modification selected the shield material of Al(2219) , then polyethylene and then Al(2219). The modified Fortran code was labeled 1SSTRN.F. The second modification considered the shield material of CO2 and Martian regolith. This modified Fortran code was labeled MARSTRN.F.
Colloid transport code-nuclear user`s manual
Jain, R.
1992-04-03
This report describes the CTCN computer code, designed to solve the equations of transient colloidal transport of radionuclides in porous and fractured media. This Fortran 77 package solves systems of coupled nonlinear differential equations with a wide range of boundary conditions. The package uses the Method of Lines technique with a special section which forms finite-difference discretizations in up to four spatial dimensions to automatically convert the system into a set of ordinary differential equations. The CTCN code then solves these equations using a robust, efficient ODE solver. Thus CTCN can be used to solve population balance equations along with the usual transport equations to model colloid transport processes or as a general problem solver to treat up to four-dimensional differential systems.
NERO- a post-maximum supernova radiation transport code
NASA Astrophysics Data System (ADS)
Maurer, I.; Jerkstrand, A.; Mazzali, P. A.; Taubenberger, S.; Hachinger, S.; Kromer, M.; Sim, S.; Hillebrandt, W.
2011-12-01
The interpretation of supernova (SN) spectra is essential for deriving SN ejecta properties such as density and composition, which in turn can tell us about their progenitors and the explosion mechanism. A very large number of atomic processes are important for spectrum formation. Several tools for calculating SN spectra exist, but they mainly focus on the very early or late epochs. The intermediate phase, which requires a non-local thermodynamic equilibrium (NLTE) treatment of radiation transport has rarely been studied. In this paper, we present a new SN radiation transport code, NERO, which can look at those epochs. All the atomic processes are treated in full NLTE, under a steady-state assumption. This is a valid approach between roughly 50 and 500 days after the explosion depending on SN type. This covers the post-maximum photospheric and the early and the intermediate nebular phase. As a test, we compare NERO to the radiation transport code of Jerkstrand, Fransson & Kozma and to the nebular code of Mazzali et al. All three codes have been developed independently and a comparison provides a valuable opportunity to investigate their reliability. Currently, NERO is one-dimensional and can be used for predicting spectra of synthetic explosion models or for deriving SN properties by spectral modelling. To demonstrate this, we study the spectra of the 'normal' Type Ia supernova (SN Ia) 2005cf between 50 and 350 days after the explosion and identify most of the common SN Ia line features at post-maximum epochs.
Code System to Calculate Waste-Isolation Flow and Transport.
Energy Science and Technology Software Center (ESTSC)
2001-01-26
Version 00 Distribution is restricted to the United States Only. SWIFT2 (Sandia Waste Isolation Flow and Transport) is a fully transient, three-dimensional code that solves the coupled equations for transport in geologic media. The processes considered are fluid flow, heat transport, brine migration, and radionuclide-chain transport. Flow, heat and brine transport are coupled via fluid density, fluid viscosity, and porosity. Together they provide the velocity field on which the radionuclide transport depends. Both porous andmore » fractured media are considered. SWIFT2 was developed for use in the analysis of deep geologic nuclear waste-disposal facilities. However, it may be used in other areas such as waste injection into saline aquifers and heat storage in aquifers. Both dual-porosity and discrete-fracture conceptualizations may be considered for the fractured zones. A variable density is included throughout, and a variety of options are available to facilitate the various uses of the code.« less
3D unstructured-mesh radiation transport codes
Morel, J.
1997-12-31
Three unstructured-mesh radiation transport codes are currently being developed at Los Alamos National Laboratory. The first code is ATTILA, which uses an unstructured tetrahedral mesh in conjunction with standard Sn (discrete-ordinates) angular discretization, standard multigroup energy discretization, and linear-discontinuous spatial differencing. ATTILA solves the standard first-order form of the transport equation using source iteration in conjunction with diffusion-synthetic acceleration of the within-group source iterations. DANTE is designed to run primarily on workstations. The second code is DANTE, which uses a hybrid finite-element mesh consisting of arbitrary combinations of hexahedra, wedges, pyramids, and tetrahedra. DANTE solves several second-order self-adjoint forms of the transport equation including the even-parity equation, the odd-parity equation, and a new equation called the self-adjoint angular flux equation. DANTE also offers three angular discretization options: $S{_}n$ (discrete-ordinates), $P{_}n$ (spherical harmonics), and $SP{_}n$ (simplified spherical harmonics). DANTE is designed to run primarily on massively parallel message-passing machines, such as the ASCI-Blue machines at LANL and LLNL. The third code is PERICLES, which uses the same hybrid finite-element mesh as DANTE, but solves the standard first-order form of the transport equation rather than a second-order self-adjoint form. DANTE uses a standard $S{_}n$ discretization in angle in conjunction with trilinear-discontinuous spatial differencing, and diffusion-synthetic acceleration of the within-group source iterations. PERICLES was initially designed to run on workstations, but a version for massively parallel message-passing machines will be built. The three codes will be described in detail and computational results will be presented.
O'Dell, R.D.; Alcouffe, R.E.
1987-09-01
This report is for the serious user of discrete ordinates transport computer codes for performing nuclear analysis calculations. The first section after the introduction provides a reasonably thorough mathematical description of the analytic Boltzmann transport equation. Next is a section on the numerical discretization of the energy, angle, and space variables in the transport equation, along with an introduction to the source iteration method. The fourth section provides numerical details and features pertinent to discrete ordinates codes. That section details angular quadrature, spatial discretization methods, iteration acceleration methods, and search capabilities. The fifth section presents considerations in choosing a discrete ordinates code for use, and this is followed by a section on typical discrete ordinates codes available throughout the world. The report ends with some guidance for the user. 73 refs., 18 figs., 13 tabs.
Zhang, Yequn; Arabaci, Murat; Djordjevic, Ivan B
2012-04-01
Leveraging the advanced coherent optical communication technologies, this paper explores the feasibility of using four-dimensional (4D) nonbinary LDPC-coded modulation (4D-NB-LDPC-CM) schemes for long-haul transmission in future optical transport networks. In contrast to our previous works on 4D-NB-LDPC-CM which considered amplified spontaneous emission (ASE) noise as the dominant impairment, this paper undertakes transmission in a more realistic optical fiber transmission environment, taking into account impairments due to dispersion effects, nonlinear phase noise, Kerr nonlinearities, and stimulated Raman scattering in addition to ASE noise. We first reveal the advantages of using 4D modulation formats in LDPC-coded modulation instead of conventional two-dimensional (2D) modulation formats used with polarization-division multiplexing (PDM). Then we demonstrate that 4D LDPC-coded modulation schemes with nonbinary LDPC component codes significantly outperform not only their conventional PDM-2D counterparts but also the corresponding 4D bit-interleaved LDPC-coded modulation (4D-BI-LDPC-CM) schemes, which employ binary LDPC codes as component codes. We also show that the transmission reach improvement offered by the 4D-NB-LDPC-CM over 4D-BI-LDPC-CM increases as the underlying constellation size and hence the spectral efficiency of transmission increases. Our results suggest that 4D-NB-LDPC-CM can be an excellent candidate for long-haul transmission in next-generation optical networks. PMID:22513641
Livnat-Levanon, Nurit; I Gilson, Amy; Ben-Tal, Nir; Lewinson, Oded
2016-01-01
ABC transporters comprise a large and ubiquitous family of proteins. From bacteria to man they translocate solutes at the expense of ATP hydrolysis. Unlike other enzymes that use ATP as an energy source, ABC transporters are notorious for having high levels of basal ATPase activity: they hydrolyze ATP also in the absence of their substrate. It is unknown what are the effects of such prolonged and constant activity on the stability and function of ABC transporters or any other enzyme. Here we report that prolonged ATP hydrolysis is beneficial to the ABC transporter BtuC2D2. Using ATPase assays, surface plasmon resonance interaction experiments, and transport assays we observe that the constantly active transporter remains stable and functional for much longer than the idle one. Remarkably, during extended activity the transporter undergoes a slow conformational change (hysteresis) and gradually attains a hyperactive state in which it is more active than it was to begin with. This phenomenon is different from stabilization of enzymes by ligand binding: the hyperactive state is only reached through ATP hydrolysis, and not ATP binding. BtuC2D2 displays a strong conformational memory for this excited state, and takes hours to return to its basal state after catalysis terminates. PMID:26905293
Livnat-Levanon, Nurit; I. Gilson, Amy; Ben-Tal, Nir; Lewinson, Oded
2016-01-01
ABC transporters comprise a large and ubiquitous family of proteins. From bacteria to man they translocate solutes at the expense of ATP hydrolysis. Unlike other enzymes that use ATP as an energy source, ABC transporters are notorious for having high levels of basal ATPase activity: they hydrolyze ATP also in the absence of their substrate. It is unknown what are the effects of such prolonged and constant activity on the stability and function of ABC transporters or any other enzyme. Here we report that prolonged ATP hydrolysis is beneficial to the ABC transporter BtuC2D2. Using ATPase assays, surface plasmon resonance interaction experiments, and transport assays we observe that the constantly active transporter remains stable and functional for much longer than the idle one. Remarkably, during extended activity the transporter undergoes a slow conformational change (hysteresis) and gradually attains a hyperactive state in which it is more active than it was to begin with. This phenomenon is different from stabilization of enzymes by ligand binding: the hyperactive state is only reached through ATP hydrolysis, and not ATP binding. BtuC2D2 displays a strong conformational memory for this excited state, and takes hours to return to its basal state after catalysis terminates. PMID:26905293
NASA Astrophysics Data System (ADS)
Innocenti, M.; Beck, A.; Lapenta, G.; Markidis, S.
2012-12-01
The kinetic simulation of intrinsically multi scale processes such as magnetic reconnection events with realistic mass ratios is a daunting task for explicit Particle In Cell (PIC) codes, which require to use resolutions of the order of the electron Debye length even when simulating dramatically bigger domains. As an example, a simulation of reconnection in the magnetotail, with domain sizes of the order of 20 di x 10 di (˜ 7.2 106 m x 3.6 106 m, with di being the ion skin depth) and a resolution of λD,e= 687 m, with λD,e the electron Debye length, requires the astounding number of 10500 x 5240 cells. Higher grid spacings can be used if the simulation is performed with an implicit PIC code, which substitutes a much less strict accuracy constraint to the stability constraint of explicit PIC codes. The same reconnection problem as before can be simulated, with an implicit PIC code resolving the scale of interest of de /2 instead of the electron Debye length (de is the electron skin depth), with the much more manageable number of 1920 x 958 cells. However, an even smaller number of cells can be used if, instead of using the same, high resolution on the entire domain, the domain to simulate is divided into subdomains each resolved with a grid spacing related to the physical scale of interest in the specific subdomain. In the case of reconnection, the division which immediately springs to mind is between electron diffusion region, ion diffusion region and outer region, where resolutions respectively of the order of fractions of the electron skin depth, of the ion skin depth and bigger can be used. We present here a new Multi Level Multi Domain (MLMD) Implicit Moment Method (IMM) Particle In Cell (PIC) code, Parsek2D-MLMD, able to perform simulations of magnetic reconnection where the expensive high resolutions are used only when needed, while the rest of the domain is simulated with grid spacings chosen according to the local scales of interest. The major difference
Documentation for TRACE: an interactive beam-transport code
Crandall, K.R.; Rusthoi, D.P.
1985-01-01
TRACE is an interactive, first-order, beam-dynamics computer program. TRACE includes space-charge forces and mathematical models for a number of beamline elements not commonly found in beam-transport codes, such as permanent-magnet quadrupoles, rf quadrupoles, rf gaps, accelerator columns, and accelerator tanks. TRACE provides an immediate graphic display of calculative results, has a powerful and easy-to-use command procedure, includes eight different types of beam-matching or -fitting capabilities, and contains its own internal HELP package. This report describes the models and equations used for each of the transport elements, the fitting procedures, and the space-charge/emittance calculations, and provides detailed instruction for using the code.
Reduced Fast Ion Transport Model For The Tokamak Transport Code TRANSP
Podesta,, Mario; Gorelenkova, Marina; White, Roscoe
2014-02-28
Fast ion transport models presently implemented in the tokamak transport code TRANSP [R. J. Hawryluk, in Physics of Plasmas Close to Thermonuclear Conditions, CEC Brussels, 1 , 19 (1980)] are not capturing important aspects of the physics associated with resonant transport caused by instabilities such as Toroidal Alfv en Eigenmodes (TAEs). This work describes the implementation of a fast ion transport model consistent with the basic mechanisms of resonant mode-particle interaction. The model is formulated in terms of a probability distribution function for the particle's steps in phase space, which is consistent with the MonteCarlo approach used in TRANSP. The proposed model is based on the analysis of fast ion response to TAE modes through the ORBIT code [R. B. White et al., Phys. Fluids 27 , 2455 (1984)], but it can be generalized to higher frequency modes (e.g. Compressional and Global Alfv en Eigenmodes) and to other numerical codes or theories.
Simulation of neoclassical transport with the continuum gyrokinetic code COGENT
Dorf, M. A.; Cohen, R. H.; Dorr, M.; Rognlien, T.; Hittinger, J.; Compton, J.; Colella, P.; Martin, D.; McCorquodale, P.
2013-01-15
The development of the continuum gyrokinetic code COGENT for edge plasma simulations is reported. The present version of the code models a nonlinear axisymmetric 4D (R, v{sub Parallel-To }, {mu}) gyrokinetic equation coupled to the long-wavelength limit of the gyro-Poisson equation. Here, R is the particle gyrocenter coordinate in the poloidal plane, and v{sub Parallel-To} and {mu} are the guiding center velocity parallel to the magnetic field and the magnetic moment, respectively. The COGENT code utilizes a fourth-order finite-volume (conservative) discretization combined with arbitrary mapped multiblock grid technology (nearly field-aligned on blocks) to handle the complexity of tokamak divertor geometry with high accuracy. Topics presented are the implementation of increasingly detailed model collision operators, and the results of neoclassical transport simulations including the effects of a strong radial electric field characteristic of a tokamak pedestal under H-mode conditions.
Validation of a comprehensive space radiation transport code.
Shinn, J L; Cucinotta, F A; Simonsen, L C; Wilson, J W; Badavi, F F; Badhwar, G D; Miller, J; Zeitlin, C; Heilbronn, L; Tripathi, R K; Clowdsley, M S; Heinbockel, J H; Xapsos, M A
1998-12-01
The HZETRN code has been developed over the past decade to evaluate the local radiation fields within sensitive materials on spacecraft in the space environment. Most of the more important nuclear and atomic processes are now modeled and evaluation within a complex spacecraft geometry with differing material components, including transition effects across boundaries of dissimilar materials, are included. The atomic/nuclear database and transport procedures have received limited validation in laboratory testing with high energy ion beams. The codes have been applied in design of the SAGE-III instrument resulting in material changes to control injurious neutron production, in the study of the Space Shuttle single event upsets, and in validation with space measurements (particle telescopes, tissue equivalent proportional counters, CR-39) on Shuttle and Mir. The present paper reviews the code development and presents recent results in laboratory and space flight validation. PMID:11542474
Simulation of neoclassical transport with the continuum gyrokinetic code COGENT
Dorf, M. A.; Cohen, R. H.; Dorr, M.; Rognlien, T.; Hittinger, J.; Compton, J.; Colella, P.; Martin, D.; McCorquodale, P.
2013-01-25
The development of the continuum gyrokinetic code COGENT for edge plasma simulations is reported. The present version of the code models a nonlinear axisymmetric 4D (R, v∥, μ) gyrokinetic equation coupled to the long-wavelength limit of the gyro-Poisson equation. Here, R is the particle gyrocenter coordinate in the poloidal plane, and v∥ and μ are the guiding center velocity parallel to the magnetic field and the magnetic moment, respectively. The COGENT code utilizes a fourth-order finite-volume (conservative) discretization combined with arbitrary mapped multiblock grid technology (nearly field-aligned on blocks) to handle the complexity of tokamak divertor geometry with high accuracy.more » Furthermore, topics presented are the implementation of increasingly detailed model collision operators, and the results of neoclassical transport simulations including the effects of a strong radial electric field characteristic of a tokamak pedestal under H-mode conditions.« less
Simulation of neoclassical transport with the continuum gyrokinetic code COGENT
Dorf, M. A.; Cohen, R. H.; Dorr, M.; Rognlien, T.; Hittinger, J.; Compton, J.; Colella, P.; Martin, D.; McCorquodale, P.
2013-01-25
The development of the continuum gyrokinetic code COGENT for edge plasma simulations is reported. The present version of the code models a nonlinear axisymmetric 4D (R, v∥, μ) gyrokinetic equation coupled to the long-wavelength limit of the gyro-Poisson equation. Here, R is the particle gyrocenter coordinate in the poloidal plane, and v∥ and μ are the guiding center velocity parallel to the magnetic field and the magnetic moment, respectively. The COGENT code utilizes a fourth-order finite-volume (conservative) discretization combined with arbitrary mapped multiblock grid technology (nearly field-aligned on blocks) to handle the complexity of tokamak divertor geometry with high accuracy. Furthermore, topics presented are the implementation of increasingly detailed model collision operators, and the results of neoclassical transport simulations including the effects of a strong radial electric field characteristic of a tokamak pedestal under H-mode conditions.
New Parallel computing framework for radiation transport codes
Kostin, M.A.; Mokhov, N.V.; Niita, K.; /JAERI, Tokai
2010-09-01
A new parallel computing framework has been developed to use with general-purpose radiation transport codes. The framework was implemented as a C++ module that uses MPI for message passing. The module is significantly independent of radiation transport codes it can be used with, and is connected to the codes by means of a number of interface functions. The framework was integrated with the MARS15 code, and an effort is under way to deploy it in PHITS. Besides the parallel computing functionality, the framework offers a checkpoint facility that allows restarting calculations with a saved checkpoint file. The checkpoint facility can be used in single process calculations as well as in the parallel regime. Several checkpoint files can be merged into one thus combining results of several calculations. The framework also corrects some of the known problems with the scheduling and load balancing found in the original implementations of the parallel computing functionality in MARS15 and PHITS. The framework can be used efficiently on homogeneous systems and networks of workstations, where the interference from the other users is possible.
Parallelization of a Monte Carlo particle transport simulation code
NASA Astrophysics Data System (ADS)
Hadjidoukas, P.; Bousis, C.; Emfietzoglou, D.
2010-05-01
We have developed a high performance version of the Monte Carlo particle transport simulation code MC4. The original application code, developed in Visual Basic for Applications (VBA) for Microsoft Excel, was first rewritten in the C programming language for improving code portability. Several pseudo-random number generators have been also integrated and studied. The new MC4 version was then parallelized for shared and distributed-memory multiprocessor systems using the Message Passing Interface. Two parallel pseudo-random number generator libraries (SPRNG and DCMT) have been seamlessly integrated. The performance speedup of parallel MC4 has been studied on a variety of parallel computing architectures including an Intel Xeon server with 4 dual-core processors, a Sun cluster consisting of 16 nodes of 2 dual-core AMD Opteron processors and a 200 dual-processor HP cluster. For large problem size, which is limited only by the physical memory of the multiprocessor server, the speedup results are almost linear on all systems. We have validated the parallel implementation against the serial VBA and C implementations using the same random number generator. Our experimental results on the transport and energy loss of electrons in a water medium show that the serial and parallel codes are equivalent in accuracy. The present improvements allow for studying of higher particle energies with the use of more accurate physical models, and improve statistics as more particles tracks can be simulated in low response time.
NASA Astrophysics Data System (ADS)
Barreteau, C.; Michon, B.; Besnard, C.; Giannini, E.
2016-06-01
Silicon and Germanium monopnictides SiP, SiAs, GeP and GeAs form a family of 2D layered semiconductors. We have succeeded in growing bulk single crystals of these compounds by melt-growth under high pressure (0.5-1 GPa) in a cubic anvil hot press. Large (mm-size), shiny, micaceous crystals of GeP, GeAs and SiAs were obtained, and could be exfoliated into 2D flakes. Small and brittle crystals of SiP were yielded by this method. High-pressure sintered polycrystalline SiP and GeAs have also been successfully used as a precursor in the Chemical Vapor Transport growth of these crystals in the presence of I2 as a transport agent. All compounds are found to crystallize in the expected layered structure and do not undergo any structural transition at low temperature, as shown by Raman spectroscopy down to T=5 K. All materials exhibit a semiconducting behavior. The electrical resistivity of GeP, GeAs and SiAs is found to depend on temperature following a 2D-Variable Range Hopping conduction mechanism. The availability of bulk crystals of these compounds opens new perspectives in the field of 2D semiconducting materials for device applications.
Sun, Jianjun; Xu, Jinbin; Cairns, Nigel J.; Perlmutter, Joel S.; Mach, Robert H.
2012-01-01
The dopamine D1, D2, D3 receptors, vesicular monoamine transporter type-2 (VMAT2), and dopamine transporter (DAT) densities were measured in 11 aged human brains (aged 77–107.8, mean: 91 years) by quantitative autoradiography. The density of D1 receptors, VMAT2, and DAT was measured using [3H]SCH23390, [3H]dihydrotetrabenazine, and [3H]WIN35428, respectively. The density of D2 and D3 receptors was calculated using the D3-preferring radioligand, [3H]WC-10 and the D2-preferring radioligand [3H]raclopride using a mathematical model developed previously by our group. Dopamine D1, D2, and D3 receptors are extensively distributed throughout striatum; the highest density of D3 receptors occurred in the nucleus accumbens (NAc). The density of the DAT is 10–20-fold lower than that of VMAT2 in striatal regions. Dopamine D3 receptor density exceeded D2 receptor densities in extrastriatal regions, and thalamus contained a high level of D3 receptors with negligible D2 receptors. The density of dopamine D1 linearly correlated with D3 receptor density in the thalamus. The density of the DAT was negligible in the extrastriatal regions whereas the VMAT2 was expressed in moderate density. D3 receptor and VMAT2 densities were in similar level between the aged human and aged rhesus brain samples, whereas aged human brain samples had lower range of densities of D1 and D2 receptors and DAT compared with the aged rhesus monkey brain. The differential density of D3 and D2 receptors in human brain will be useful in the interpretation of PET imaging studies in human subjects with existing radiotracers, and assist in the validation of newer PET radiotracers having a higher selectivity for dopamine D2 or D3 receptors. PMID:23185343
77 FR 18716 - Transportation Security Administration Postal Zip Code Change; Technical Amendment
Federal Register 2010, 2011, 2012, 2013, 2014
2012-03-28
... SECURITY Transportation Security Administration 49 CFR Part 1572 Transportation Security Administration Postal Zip Code Change; Technical Amendment AGENCY: Transportation Security Administration, DHS. ACTION... Achuko, Office of the Chief Counsel, TSA-2, Transportation Security Administration, 601 South 12th...
Daniels, J.; Williams, J.; Asherson, P.; McGuffin, P.; Owen, M.
1995-02-27
It has been suggested that the cytochrome P450 mono-oxygenase, debrisoquine 4-hydroxylase, is involved in the catabolism and processing of neurotransmitters subsequent to their reuptake into target cells. It is also thought to be related to the dopamine transporter that acts to take released dopamine back up into presynaptic terminals. The present study used the association approach to test the hypothesis that mutations in the genes for debrisoquine 4-hydroxylase (CYP2D6) and the dopamine transporter (DAT) confer susceptibility to schizophrenia. There were no differences in allele or genotype frequencies between patients and controls in the mutations causing the poor metaboliser phenotype in CYP2D6. In addition there was no association found between schizophrenia and a 48 bp repeat within the 3{prime} untranslated region of DAT. 18 refs., 2 tabs.
Verification of ARES transport code system with TAKEDA benchmarks
NASA Astrophysics Data System (ADS)
Zhang, Liang; Zhang, Bin; Zhang, Penghe; Chen, Mengteng; Zhao, Jingchang; Zhang, Shun; Chen, Yixue
2015-10-01
Neutron transport modeling and simulation are central to many areas of nuclear technology, including reactor core analysis, radiation shielding and radiation detection. In this paper the series of TAKEDA benchmarks are modeled to verify the critical calculation capability of ARES, a discrete ordinates neutral particle transport code system. SALOME platform is coupled with ARES to provide geometry modeling and mesh generation function. The Koch-Baker-Alcouffe parallel sweep algorithm is applied to accelerate the traditional transport calculation process. The results show that the eigenvalues calculated by ARES are in excellent agreement with the reference values presented in NEACRP-L-330, with a difference less than 30 pcm except for the first case of model 3. Additionally, ARES provides accurate fluxes distribution compared to reference values, with a deviation less than 2% for region-averaged fluxes in all cases. All of these confirms the feasibility of ARES-SALOME coupling and demonstrate that ARES has a good performance in critical calculation.
NASA Astrophysics Data System (ADS)
Key, K.
2013-12-01
This work announces the public release of an open-source inversion code named MARE2DEM (Modeling with Adaptively Refined Elements for 2D Electromagnetics). Although initially designed for the rapid inversion of marine electromagnetic data, MARE2DEM now supports a wide variety of acquisition configurations for both offshore and onshore surveys that utilize electric and magnetic dipole transmitters or magnetotelluric plane waves. The model domain is flexibly parameterized using a grid of arbitrarily shaped polygonal regions, allowing for complicated structures such as topography or seismically imaged horizons to be easily assimilated. MARE2DEM efficiently solves the forward problem in parallel by dividing the input data parameters into smaller subsets using a parallel data decomposition algorithm. The data subsets are then solved in parallel using an automatic adaptive finite element method that iterative solves the forward problem on successively refined finite element meshes until a specified accuracy tolerance is met, thus freeing the end user from the burden of designing an accurate numerical modeling grid. Regularized non-linear inversion for isotropic or anisotropic conductivity is accomplished with a new implementation of Occam's method referred to as fast-Occam, which is able to minimize the objective function in much fewer forward evaluations than the required by the original method. This presentation will review the theoretical considerations behind MARE2DEM and use a few recent offshore EM data sets to demonstrate its capabilities and to showcase the software interface tools that streamline model building and data inversion.
Comparison of Space Radiation Calculations from Deterministic and Monte Carlo Transport Codes
NASA Technical Reports Server (NTRS)
Adams, J. H.; Lin, Z. W.; Nasser, A. F.; Randeniya, S.; Tripathi, r. K.; Watts, J. W.; Yepes, P.
2010-01-01
The presentation outline includes motivation, radiation transport codes being considered, space radiation cases being considered, results for slab geometry, results from spherical geometry, and summary. ///////// main physics in radiation transport codes hzetrn uprop fluka geant4, slab geometry, spe, gcr,
C5 Benchmark Problem with Discrete Ordinate Radiation Transport Code DENOVO
Yesilyurt, Gokhan; Clarno, Kevin T; Evans, Thomas M; Davidson, Gregory G; Fox, Patricia B
2011-01-01
The C5 benchmark problem proposed by the Organisation for Economic Co-operation and Development/Nuclear Energy Agency was modeled to examine the capabilities of Denovo, a three-dimensional (3-D) parallel discrete ordinates (S{sub N}) radiation transport code, for problems with no spatial homogenization. Denovo uses state-of-the-art numerical methods to obtain accurate solutions to the Boltzmann transport equation. Problems were run in parallel on Jaguar, a high-performance supercomputer located at Oak Ridge National Laboratory. Both the two-dimensional (2-D) and 3-D configurations were analyzed, and the results were compared with the reference MCNP Monte Carlo calculations. For an additional comparison, SCALE/KENO-V.a Monte Carlo solutions were also included. In addition, a sensitivity analysis was performed for the optimal angular quadrature and mesh resolution for both the 2-D and 3-D infinite lattices of UO{sub 2} fuel pin cells. Denovo was verified with the C5 problem. The effective multiplication factors, pin powers, and assembly powers were found to be in good agreement with the reference MCNP and SCALE/KENO-V.a Monte Carlo calculations.
VERIFICATION OF TRANSPORT CODES BY THE METHOD OF MANUFACTURED SOLUTIONS: THE ATTILA EXPERIENCE
S. D. PAUTZ
2001-03-19
We extend the Method of Manufactured Solutions (MMS) to the verification of transport codes. We derive analytic fixed sources required by the MMS procedure for several types of transport problems and apply the method to the Attila transport code. By means of this method we discover and correct several coding mistakes in Attila and ultimately verify its correct implementation for the problems studied. Our studies reveal that the MMS procedure is a useful tool for transport code development.
Application of the INSTANT-HPS PN Transport Code to the C5G7 Benchmark Problem
Y. Wang; H. Zhang; R. H. Szilard; R. C. Martineau
2011-06-01
INSTANT is the INL's next generation neutron transport solver to support high-fidelity multi-physics reactor simulation INSTANT is in continuous development to extend its capability Code is designed to take full advantage of middle to large cluster (10-1000 processors) Code is designed to focus on method adaptation while also mesh adaptation will be possible. It utilizes the most modern computing techniques to generate a neutronics tool of full-core transport calculations for reactor analysis and design. It can perform calculations on unstructured 2D/3D triangular, hexagonal and Cartesian geometries. Calculations can be easily extended to more geometries because of the independent mesh framework coded with the model Fortran. This code has a multigroup solver with thermal rebalance and Chebyshev acceleration. It employs second-order PN and Hybrid Finite Element method (PNHFEM) discretization scheme. Three different in-group solvers - preconditioned Conjugate Gradient (CG) method, preconditioned Generalized Minimal Residual Method (GMRES) and Red-Black iteration - have been implemented and parallelized with the spatial domain decomposition in the code. The input is managed with extensible markup language (XML) format. 3D variables including the flux distributions are outputted into VTK files, which can be visualized by tools such as VisIt and ParaView. An extension of the code named INSTANTHPS provides the capability to perform 3D heterogeneous transport calculations within fuel pins. C5G7 is an OECD/NEA benchmark problem created to test the ability of modern deterministic transport methods and codes to treat reactor core problems without spatial homogenization. This benchmark problem had been widely analyzed with various code packages. In this transaction, results of the applying the INSTANT-HPS code to the C5G7 problem are summarized.
NASA Astrophysics Data System (ADS)
Sawada, Y.; Uozumi, Y.; Nogamine, S.; Yamada, T.; Iwamoto, Y.; Sato, T.; Niita, K.
2012-11-01
The Intranuclear Cascade with Emission of Light Fragment (INC-ELF) code has been developed and implemented in the Particle and Heavy Ion Transport code System (PHITS). The INC-ELF code explicitly includes nucleon correlations within the framework of the INC model to describe light fragment emissions from nuclear spallation reactions by using the model in Phys. Rev. C 84, (2011) 064617. In addition to the degrees of freedom of nucleons, the developed code also accounts for pions, Δs, and N∗s, and can cover energy ranges up to 3 GeV. The predictive capabilities of the ELF/PHITS system have been verified through comparison with a diverse set of experimental observations. In particular, the verification was conducted with abundant double-differential cross-section data covering a wide range of reactions (e.g., (p, p'x), (p, nx), (p, dx), (p, 3Hex), (p, αx) and (p, πx) reactions) over a wide energy range (between 400 MeV and 1.5 GeV). As a result, our ELF/PHITS code has demonstrated strong predictive capability for all of these data, although areas requiring future study remain due to the lack of experimental data on high-energy cluster production.
Overview of Particle and Heavy Ion Transport Code System PHITS
NASA Astrophysics Data System (ADS)
Sato, Tatsuhiko; Niita, Koji; Matsuda, Norihiro; Hashimoto, Shintaro; Iwamoto, Yosuke; Furuta, Takuya; Noda, Shusaku; Ogawa, Tatsuhiko; Iwase, Hiroshi; Nakashima, Hiroshi; Fukahori, Tokio; Okumura, Keisuke; Kai, Tetsuya; Chiba, Satoshi; Sihver, Lembit
2014-06-01
A general purpose Monte Carlo Particle and Heavy Ion Transport code System, PHITS, is being developed through the collaboration of several institutes in Japan and Europe. The Japan Atomic Energy Agency is responsible for managing the entire project. PHITS can deal with the transport of nearly all particles, including neutrons, protons, heavy ions, photons, and electrons, over wide energy ranges using various nuclear reaction models and data libraries. It is written in Fortran language and can be executed on almost all computers. All components of PHITS such as its source, executable and data-library files are assembled in one package and then distributed to many countries via the Research organization for Information Science and Technology, the Data Bank of the Organization for Economic Co-operation and Development's Nuclear Energy Agency, and the Radiation Safety Information Computational Center. More than 1,000 researchers have been registered as PHITS users, and they apply the code to various research and development fields such as nuclear technology, accelerator design, medical physics, and cosmic-ray research. This paper briefly summarizes the physics models implemented in PHITS, and introduces some important functions useful for specific applications, such as an event generator mode and beam transport functions.
An object-oriented implementation of a parallel Monte Carlo code for radiation transport
NASA Astrophysics Data System (ADS)
Santos, Pedro Duarte; Lani, Andrea
2016-05-01
This paper describes the main features of a state-of-the-art Monte Carlo solver for radiation transport which has been implemented within COOLFluiD, a world-class open source object-oriented platform for scientific simulations. The Monte Carlo code makes use of efficient ray tracing algorithms (for 2D, axisymmetric and 3D arbitrary unstructured meshes) which are described in detail. The solver accuracy is first verified in testcases for which analytical solutions are available, then validated for a space re-entry flight experiment (i.e. FIRE II) for which comparisons against both experiments and reference numerical solutions are provided. Through the flexible design of the physical models, ray tracing and parallelization strategy (fully reusing the mesh decomposition inherited by the fluid simulator), the implementation was made efficient and reusable.
Energy Science and Technology Software Center (ESTSC)
2004-08-01
AnisWave2D is a 2D finite-difference code for a simulating seismic wave propagation in fully anisotropic materials. The code is implemented to run in parallel over multiple processors and is fully portable. A mesh refinement algorithm has been utilized to allow the grid-spacing to be tailored to the velocity model, avoiding the over-sampling of high-velocity materials that usually occurs in fixed-grid schemes.
Grant, K.E.; Taylor, K.E.; Ellis, J.S.; Wuebbles, D.J.
1987-07-01
The authors have implemented a series of state of the art radiation transport submodels in previously developed one dimensional and two dimensional chemical transport models of the troposphere and stratosphere. These submodels provide the capability of calculating accurate solar and infrared heating rates. They are a firm basis for further radiation submodel development as well as for studying interactions between radiation and model dynamics under varying conditions of clear sky, clouds, and aerosols. 37 refs., 3 figs.
Multidimensional electron-photon transport with standard discrete ordinates codes
Drumm, C.R.
1997-04-01
A method is described for generating electron cross sections that are comparable with standard discrete ordinates codes without modification. There are many advantages of using an established discrete ordinates solver, e.g. immediately available adjoint capability. Coupled electron-photon transport capability is needed for many applications, including the modeling of the response of electronics components to space and man-made radiation environments. The cross sections have been successfully used in the DORT, TWODANT and TORT discrete ordinates codes. The cross sections are shown to provide accurate and efficient solutions to certain multidimensional electron-photon transport problems. The key to the method is a simultaneous solution of the continuous-slowing-down (CSD) portion and elastic-scattering portion of the scattering source by the Goudsmit-Saunderson theory. The resulting multigroup-Legendre cross sections are much smaller than the true scattering cross sections that they represent. Under certain conditions, the cross sections are guaranteed positive and converge with a low-order Legendre expansion.
Time-dependent recycling modeling with edge plasma transport codes
NASA Astrophysics Data System (ADS)
Pigarov, A.; Krasheninnikov, S.; Rognlien, T.; Taverniers, S.; Hollmann, E.
2013-10-01
First,we discuss extensions to Macroblob approach which allow to simulate more accurately dynamics of ELMs, pedestal and edge transport with UEDGE code. Second,we present UEDGE modeling results for H mode discharge with infrequent ELMs and large pedestal losses on DIII-D. In modeled sequence of ELMs this discharge attains a dynamic equilibrium. Temporal evolution of pedestal plasma profiles, spectral line emission, and surface temperature matching experimental data over ELM cycle is discussed. Analysis of dynamic gas balance highlights important role of material surfaces. We quantified the wall outgassing between ELMs as 3X the NBI fueling and the recycling coefficient as 0.8 for wall pumping via macroblob-wall interactions. Third,we also present results from multiphysics version of UEDGE with built-in, reduced, 1-D wall models and analyze the role of various PMI processes. Progress in framework-coupled UEDGE/WALLPSI code is discussed. Finally, implicit coupling schemes are important feature of multiphysics codes and we report on the results of parametric analysis of convergence and performance for Picard and Newton iterations in a system of coupled deterministic-stochastic ODE and proposed modifications enhancing convergence.
NASA Astrophysics Data System (ADS)
Cappelli, Mark; Young, Chris; Cha, Eusnun; Fernandez, Eduardo; Stanford Plasma Physics Laboratory Collaboration; Eckerd College Collaboration
2015-09-01
We present a simple, zero-equation turbulence model for electron transport across the magnetic field of a plasma Hall thruster and integrate this model into 2-D hybrid particle-in-cell simulations of a 72 mm diameter laboratory thruster operating at 400 W. The turbulent transport model is based on the assumption that the primary means of electron energy dissipation is the turbulent eddy cascade in the electron fluid to smaller scales. Implementing the model into 2-D hybrid simulations is relatively straightforward and leverages the existing framework for solving the electron fluid equations. We find that the model captures the strong axial variation in the mobility seen in experiments. In particular, it predicts the existence of a strong transport barrier which anchors the region of plasma acceleration. The model also captures the time-averaged experimental discharge current and its fluctuations due to ionization instabilities. We observe quantitative agreement with recent laser induced fluorescence measurements of time-averaged xenon ion and neutral velocities along the channel centerline. This work was supported by the Air Force Office of Scientific Research.
Dimensionality-dependent charge transport in close-packed nanoparticle arrays: from 2D to 3D
Wang, Ying; Duan, Chao; Peng, Lianmao; Liao, Jianhui
2014-01-01
Charge transport properties in close-packed nanoparticle arrays with thickness crossing over from two dimensions to three dimensions have been studied. The dimensionality transition of nanoparticle arrays was realized by continually printing spatially well-defined nanoparticle monolayers on top of the device in situ. The evolution of charge transport properties depending on the dimensionality has been investigated in both the Efros-Shaklovskii variable-range-hopping (ES-VRH) (low temperature) regime and the sequential hopping (SH) (medium temperature) regime. We find that the energy barriers to transport decrease when the thickness of nanoparticle arrays increases from monolayer to multilayers, but start to level off at the thickness of 4–5 monolayers. The energy barriers are characterized by the coefficient βD at ES-VRH regime and the activation energy Ea at SH regime. Moreover, a turning point for the temperature coefficient of conductance was observed in multilayer nanoparticle arrays at high temperature, which is attributed to the increasing mobility with decreasing temperature of hopping transport in three dimensions. PMID:25523836
Dimensionality-dependent charge transport in close-packed nanoparticle arrays: from 2D to 3D.
Wang, Ying; Duan, Chao; Peng, Lianmao; Liao, Jianhui
2014-01-01
Charge transport properties in close-packed nanoparticle arrays with thickness crossing over from two dimensions to three dimensions have been studied. The dimensionality transition of nanoparticle arrays was realized by continually printing spatially well-defined nanoparticle monolayers on top of the device in situ. The evolution of charge transport properties depending on the dimensionality has been investigated in both the Efros-Shaklovskii variable-range-hopping (ES-VRH) (low temperature) regime and the sequential hopping (SH) (medium temperature) regime. We find that the energy barriers to transport decrease when the thickness of nanoparticle arrays increases from monolayer to multilayers, but start to level off at the thickness of 4-5 monolayers. The energy barriers are characterized by the coefficient βD at ES-VRH regime and the activation energy Ea at SH regime. Moreover, a turning point for the temperature coefficient of conductance was observed in multilayer nanoparticle arrays at high temperature, which is attributed to the increasing mobility with decreasing temperature of hopping transport in three dimensions. PMID:25523836
Dimensionality-dependent charge transport in close-packed nanoparticle arrays: from 2D to 3D
NASA Astrophysics Data System (ADS)
Wang, Ying; Duan, Chao; Peng, Lianmao; Liao, Jianhui
2014-12-01
Charge transport properties in close-packed nanoparticle arrays with thickness crossing over from two dimensions to three dimensions have been studied. The dimensionality transition of nanoparticle arrays was realized by continually printing spatially well-defined nanoparticle monolayers on top of the device in situ. The evolution of charge transport properties depending on the dimensionality has been investigated in both the Efros-Shaklovskii variable-range-hopping (ES-VRH) (low temperature) regime and the sequential hopping (SH) (medium temperature) regime. We find that the energy barriers to transport decrease when the thickness of nanoparticle arrays increases from monolayer to multilayers, but start to level off at the thickness of 4-5 monolayers. The energy barriers are characterized by the coefficient βD at ES-VRH regime and the activation energy Ea at SH regime. Moreover, a turning point for the temperature coefficient of conductance was observed in multilayer nanoparticle arrays at high temperature, which is attributed to the increasing mobility with decreasing temperature of hopping transport in three dimensions.
A massively parallel method of characteristic neutral particle transport code for GPUs
Boyd, W. R.; Smith, K.; Forget, B.
2013-07-01
Over the past 20 years, parallel computing has enabled computers to grow ever larger and more powerful while scientific applications have advanced in sophistication and resolution. This trend is being challenged, however, as the power consumption for conventional parallel computing architectures has risen to unsustainable levels and memory limitations have come to dominate compute performance. Heterogeneous computing platforms, such as Graphics Processing Units (GPUs), are an increasingly popular paradigm for solving these issues. This paper explores the applicability of GPUs for deterministic neutron transport. A 2D method of characteristics (MOC) code - OpenMOC - has been developed with solvers for both shared memory multi-core platforms as well as GPUs. The multi-threading and memory locality methodologies for the GPU solver are presented. Performance results for the 2D C5G7 benchmark demonstrate 25-35 x speedup for MOC on the GPU. The lessons learned from this case study will provide the basis for further exploration of MOC on GPUs as well as design decisions for hardware vendors exploring technologies for the next generation of machines for scientific computing. (authors)
Entin, M. V.; Magarill, L. I.; Olshanetsky, E. B. Kvon, Z. D.; Mikhailov, N. N.; Dvoretsky, S. A.
2013-11-15
The influence of e-h scattering on the conductivity and magnetotransport of 2D semimetallic HgTe is studied both theoretically and experimentally. The presence of e-h scattering leads to the friction between electrons and holes resulting in a large temperature-dependent contribution to the transport coefficients. The coefficient of friction between electrons and holes is determined. The comparison of experimental data with the theory shows that the interaction between electrons and holes based on the long-range Coulomb potential strongly underestimates the e-h friction. The experimental results are in agreement with the model of strong short-range e-h interaction.
NASA Technical Reports Server (NTRS)
Chang, Chau-Lyan
2003-01-01
During the past two decades, our understanding of laminar-turbulent transition flow physics has advanced significantly owing to, in a large part, the NASA program support such as the National Aerospace Plane (NASP), High-speed Civil Transport (HSCT), and Advanced Subsonic Technology (AST). Experimental, theoretical, as well as computational efforts on various issues such as receptivity and linear and nonlinear evolution of instability waves take part in broadening our knowledge base for this intricate flow phenomenon. Despite all these advances, transition prediction remains a nontrivial task for engineers due to the lack of a widely available, robust, and efficient prediction tool. The design and development of the LASTRAC code is aimed at providing one such engineering tool that is easy to use and yet capable of dealing with a broad range of transition related issues. LASTRAC was written from scratch based on the state-of-the-art numerical methods for stability analysis and modem software technologies. At low fidelity, it allows users to perform linear stability analysis and N-factor transition correlation for a broad range of flow regimes and configurations by using either the linear stability theory (LST) or linear parabolized stability equations (LPSE) method. At high fidelity, users may use nonlinear PSE to track finite-amplitude disturbances until the skin friction rise. Coupled with the built-in receptivity model that is currently under development, the nonlinear PSE method offers a synergistic approach to predict transition onset for a given disturbance environment based on first principles. This paper describes the governing equations, numerical methods, code development, and case studies for the current release of LASTRAC. Practical applications of LASTRAC are demonstrated for linear stability calculations, N-factor transition correlation, non-linear breakdown simulations, and controls of stationary crossflow instability in supersonic swept wing boundary
Physics models in the toroidal transport code PROCTR
Howe, H.C.
1990-08-01
The physics models that are contained in the toroidal transport code PROCTR are described in detail. Time- and space-dependent models are included for the plasma hydrogenic-ion, helium, and impurity densities, the electron and ion temperatures, the toroidal rotation velocity, and the toroidal current profile. Time- and depth-dependent models for the trapped and mobile hydrogenic particle concentrations in the wall and a time-dependent point model for the number of particles in the limiter are also included. Time-dependent models for neutral particle transport, neutral beam deposition and thermalization, fusion heating, impurity radiation, pellet injection, and the radial electric potential are included and recalculated periodically as the time-dependent models evolve. The plasma solution is obtained either in simple flux coordinates, where the radial shift of each elliptical, toroidal flux surface is included to maintain an approximate pressure equilibrium, or in general three-dimensional torsatron coordinates represented by series of helical harmonics. The detailed coupling of the plasma, scrape-off layer, limiter, and wall models through the neutral transport model makes PROCTR especially suited for modeling of recycling and particle control in toroidal plasmas. The model may also be used in a steady-state profile analysis mode for studying energy and particle balances starting with measured plasma profiles.
PRESTO low-level waste transport and risk assessment code
Little, C.A.; Fields, D.E.; McDowell-Boyer, L.M.; Emerson, C.J.
1981-01-01
PRESTO (Prediction of Radiation Effects from Shallow Trench Operations) is a computer code developed under US Environmental Protection Agency (EPA) funding to evaluate possible health effects from shallow land burial trenches. The model is intended to be generic and to assess radionuclide transport, ensuing exposure, and health impact to a static local population for a 1000-y period following the end of burial operations. Human exposure scenarios considered by the model include normal releases (including leaching and operational spillage), human intrusion, and site farming or reclamation. Pathways and processes of transit from the trench to an individual or population inlude: groundwater transport, overland flow, erosion, surface water dilution, resuspension, atmospheric transport, deposition, inhalation, and ingestion of contaminated beef, milk, crops, and water. Both population doses and individual doses are calculated as well as doses to the intruder and farmer. Cumulative health effects in terms of deaths from cancer are calculated for the population over the thousand-year period using a life-table approach. Data bases are being developed for three extant shallow land burial sites: Barnwell, South Carolina; Beatty, Nevada; and West Valley, New York.
Fevotte, F.; Lathuiliere, B.
2013-07-01
The large increase in computing power over the past few years now makes it possible to consider developing 3D full-core heterogeneous deterministic neutron transport solvers for reference calculations. Among all approaches presented in the literature, the method first introduced in [1] seems very promising. It consists in iterating over resolutions of 2D and ID MOC problems by taking advantage of prismatic geometries without introducing approximations of a low order operator such as diffusion. However, before developing a solver with all industrial options at EDF, several points needed to be clarified. In this work, we first prove the convergence of this iterative process, under some assumptions. We then present our high-performance, parallel implementation of this algorithm in the MICADO solver. Benchmarking the solver against the Takeda case shows that the 2D-1D coupling algorithm does not seem to affect the spatial convergence order of the MOC solver. As for performance issues, our study shows that even though the data distribution is suited to the 2D solver part, the efficiency of the ID part is sufficient to ensure a good parallel efficiency of the global algorithm. After this study, the main remaining difficulty implementation-wise is about the memory requirement of a vector used for initialization. An efficient acceleration operator will also need to be developed. (authors)
Schaffranek, Raymond W.
2004-01-01
A numerical model for simulation of surface-water integrated flow and transport in two (horizontal-space) dimensions is documented. The model solves vertically integrated forms of the equations of mass and momentum conservation and solute transport equations for heat, salt, and constituent fluxes. An equation of state for salt balance directly couples solution of the hydrodynamic and transport equations to account for the horizontal density gradient effects of salt concentrations on flow. The model can be used to simulate the hydrodynamics, transport, and water quality of well-mixed bodies of water, such as estuaries, coastal seas, harbors, lakes, rivers, and inland waterways. The finite-difference model can be applied to geographical areas bounded by any combination of closed land or open water boundaries. The simulation program accounts for sources of internal discharges (such as tributary rivers or hydraulic outfalls), tidal flats, islands, dams, and movable flow barriers or sluices. Water-quality computations can treat reactive and (or) conservative constituents simultaneously. Input requirements include bathymetric and topographic data defining land-surface elevations, time-varying water level or flow conditions at open boundaries, and hydraulic coefficients. Optional input includes the geometry of hydraulic barriers and constituent concentrations at open boundaries. Time-dependent water level, flow, and constituent-concentration data are required for model calibration and verification. Model output consists of printed reports and digital files of numerical results in forms suitable for postprocessing by graphical software programs and (or) scientific visualization packages. The model is compatible with most mainframe, workstation, mini- and micro-computer operating systems and FORTRAN compilers. This report defines the mathematical formulation and computational features of the model, explains the solution technique and related model constraints, describes the
2D coherent charge transport in highly ordered conducting polymers doped by solid state diffusion.
Kang, Keehoon; Watanabe, Shun; Broch, Katharina; Sepe, Alessandro; Brown, Adam; Nasrallah, Iyad; Nikolka, Mark; Fei, Zhuping; Heeney, Martin; Matsumoto, Daisuke; Marumoto, Kazuhiro; Tanaka, Hisaaki; Kuroda, Shin-Ichi; Sirringhaus, Henning
2016-08-01
Doping is one of the most important methods to control charge carrier concentration in semiconductors. Ideally, the introduction of dopants should not perturb the ordered microstructure of the semiconducting host. In some systems, such as modulation-doped inorganic semiconductors or molecular charge transfer crystals, this can be achieved by spatially separating the dopants from the charge transport pathways. However, in conducting polymers, dopants tend to be randomly distributed within the conjugated polymer, and as a result the transport properties are strongly affected by the resulting structural and electronic disorder. Here, we show that in the highly ordered lamellar microstructure of a regioregular thiophene-based conjugated polymer, a small-molecule p-type dopant can be incorporated by solid state diffusion into the layers of solubilizing side chains without disrupting the conjugated layers. In contrast to more disordered systems, this allows us to observe coherent, free-electron-like charge transport properties, including a nearly ideal Hall effect in a wide temperature range, a positive magnetoconductance due to weak localization and the Pauli paramagnetic spin susceptibility. PMID:27159015
NASA Astrophysics Data System (ADS)
Eckert, Dominik; Kürzinger, Petra; Bauer, Robert; Griebler, Christian; Cirpka, Olaf A.
2015-01-01
Biodegradation in contaminated aquifers has been shown to be most pronounced at the fringe of contaminant plumes, where mixing of contaminated water and ambient groundwater, containing dissolved electron acceptors, stimulates microbial activity. While physical mixing of contaminant and electron acceptor by transverse dispersion has been shown to be the major bottleneck for biodegradation in steady-state plumes, so far little is known on the effect of flow and transport dynamics (caused, e.g., by a seasonally fluctuating groundwater table) on biodegradation in these systems. Towards this end we performed experiments in quasi-two-dimensional flow-through microcosms on aerobic toluene degradation by Pseudomonas putida F1. Plume dynamics were simulated by vertical alteration of the toluene plume position and experimental results were analyzed by reactive-transport modeling. We found that, even after disappearance of the toluene plume for two weeks, the majority of microorganisms stayed attached to the sediment and regained their full biodegradation potential within two days after reappearance of the toluene plume. Our results underline that besides microbial growth, also maintenance and dormancy are important processes that affect biodegradation performance under transient environmental conditions and therefore deserve increased consideration in future reactive-transport modeling.
NASA Astrophysics Data System (ADS)
Hue, V.; Cavalié, T.; Dobrijevic, M.; Hersant, F.; Greathouse, T. K.
2015-09-01
Saturn's axial tilt of 26.7° produces seasons in a similar way as on Earth. Both the stratospheric temperature and composition are affected by this latitudinally varying insolation along Saturn's orbital path. A new time-dependent 2D photochemical model is presented to study the seasonal evolution of Saturn's stratospheric composition. This study focuses on the impact of the seasonally variable thermal field on the main stratospheric C2-hydrocarbon chemistry (C2H2 and C2H6) using a realistic radiative climate model. Meridional mixing and advective processes are implemented in the model but turned off in the present study for the sake of simplicity. The results are compared to a simple study case where a latitudinally and temporally steady thermal field is assumed. Our simulations suggest that, when the seasonally variable thermal field is accounted for, the downward diffusion of the seasonally produced hydrocarbons is faster due to the seasonal compression of the atmospheric column during winter. This effect increases with increasing latitudes which experience the most important thermal changes in the course of the seasons. The seasonal variability of C2H2 and C2H6 therefore persists at higher-pressure levels with a seasonally-variable thermal field. Cassini limb-observations of C2H2 and C2H6 (Guerlet, S. et al. [2009]. Icarus 203, 214-232) are reasonably well-reproduced from the equator to 40° in both hemispheres in the 0.1-1 mbar pressure range. At lower pressure levels, the models only fit the Cassini observations in the northern hemisphere, from the equator to 40°N. Beyond 40° in both hemispheres, deviations from the pure photochemical predictions, mostly in the southern hemisphere, suggest the presence of large-scale stratospheric dynamics.
Bandy, P.J.; Hall, L.F.
1993-03-01
This report presents information on computer codes for numerical and analytical models that have been used at the Idaho National Engineering Laboratory (INEL) to model ground water and surface water flow and contaminant transport. Organizations conducting modeling at the INEL include: EG G Idaho, Inc., US Geological Survey, and Westinghouse Idaho Nuclear Company. Information concerning computer codes included in this report are: agency responsible for the modeling effort, name of the computer code, proprietor of the code (copyright holder or original author), validation and verification studies, applications of the model at INEL, the prime user of the model, computer code description, computing environment requirements, and documentation and references for the computer code.
NASA Astrophysics Data System (ADS)
Baitsch-Ghirardello, Bettina; Stracke, Andreas; Connolly, James A. D.; Nikolaeva, Ksenia M.; Gerya, Taras V.
2014-11-01
Understanding the physical-chemical mechanisms and pathways of geochemical transport in subduction zones remains a long-standing goal of subduction-related research. In this study, we perform fully coupled geochemical-thermo-mechanical (GcTM) numerical simulations to investigate Pb isotopic signatures of the two key "outputs" of subduction zones: (A) serpentinite mélanges and (B) arc basalts. With this approach we analyze three different geodynamic regimes of intra-oceanic subduction systems: (1) retreating subduction with backarc spreading, (2) stable subduction with high fluid-related weakening, and (3) stable subduction with low fluid-related weakening. Numerical results suggest a three-stage Pb geochemical transport in subduction zones: (I) from subducting sediments and oceanic crust to serpentinite mélanges, (II) from subducting serpentinite mélanges to subarc asthenospheric wedge and (III) from the mantle wedge to arc volcanics. Mechanical mixing and fluid-assisted geochemical transport above slabs result in spatially and temporarily variable Pb concentrations in the serpentinized forearc mantle as well as in arc basalts. The Pb isotopic ratios are strongly heterogeneous and show five types of geochemical mixing trends: (i) binary mantle-MORB, (i) binary MORB-sediments, (iii) double binary MORB-mantle and MORB-sediments, (iv) double binary MORB-mantle and mantle-sediments and (v) triple MORB-sediment-mantle. Double binary and triple mixing trends are transient and characterize relatively early stages of subduction. In contrast, steady-state binary mantle-MORB and MORB-sediments trends are typical for mature subduction zones with respectively low and high intensity of sedimentary melange subduction. Predictions from our GcTM models are in agreement with Pb isotopic data from some natural subduction zones.
Bailey, T S; Adams, M L; Chang, J H
2008-10-01
We present a new spatial discretization of the discrete-ordinates transport equation in two-dimensional cylindrical (RZ) geometry for arbitrary polygonal meshes. This discretization is a discontinuous finite element method that utilizes the piecewise linear basis functions developed by Stone and Adams. We describe an asymptotic analysis that shows this method to be accurate for many problems in the thick diffusion limit on arbitrary polygons, allowing this method to be applied to radiative transfer problems with these types of meshes. We also present numerical results for multiple problems on quadrilateral grids and compare these results to the well-known bi-linear discontinuous finite element method.
NASA Astrophysics Data System (ADS)
Derakhshan, V.; Ketabi, S. A.; Moghaddam, A. G.
2016-09-01
We employed the formalism of bond currents, expressed in terms of non-equilibrium Green’s function to obtain the local currents and transport features of zigzag silicene ribbon in the presence of magnetic impurity. When only intrinsic and Rashba spin–orbit interactions are present, silicene behaves as a two-dimensional topological insulator with gapless edge states. But in the presence of finite intrinsic spin–orbit interaction, the edge states start to penetrate into the bulk of the sample by increasing Rashba interaction strength. The exchange interaction induced by local impurities breaks the time-reversal symmetry of the gapless edge states and influences the topological properties strongly. Subsequently, the singularity of partial Berry curvature disappears and the silicene nanoribbon becomes a trivial insulator. On the other hand, when the concentration of the magnetic impurities is low, the edge currents are not affected significantly. In this case, when the exchange field lies in the x–y plane, the spin mixing around magnetic impurity is more profound rather than the case in which the exchange field is directed along the z-axis. Nevertheless, when the exchange field of magnetic impurities is placed in the x–y plane, a spin-polarized conductance is observed. The resulting conductance polarization can be tuned by the concentration of the impurities and even completely polarized spin transport is achievable.
Status of the MORSE multigroup Monte Carlo radiation transport code
Emmett, M.B.
1993-06-01
There are two versions of the MORSE multigroup Monte Carlo radiation transport computer code system at Oak Ridge National Laboratory. MORSE-CGA is the most well-known and has undergone extensive use for many years. MORSE-SGC was originally developed in about 1980 in order to restructure the cross-section handling and thereby save storage. However, with the advent of new computer systems having much larger storage capacity, that aspect of SGC has become unnecessary. Both versions use data from multigroup cross-section libraries, although in somewhat different formats. MORSE-SGC is the version of MORSE that is part of the SCALE system, but it can also be run stand-alone. Both CGA and SGC use the Multiple Array System (MARS) geometry package. In the last six months the main focus of the work on these two versions has been on making them operational on workstations, in particular, the IBM RISC 6000 family. A new version of SCALE for workstations is being released to the Radiation Shielding Information Center (RSIC). MORSE-CGA, Version 2.0, is also being released to RSIC. Both SGC and CGA have undergone other revisions recently. This paper reports on the current status of the MORSE code system.
NASA Astrophysics Data System (ADS)
Olson, Gordon L.
2012-04-01
When using polynomial expansions for the angular variables in the radiation transport equation, the usual procedure is to truncate the series by setting all higher order terms to zero. At low order, such simple closures may not give the optimum solution. This work tests alternate closures that scale either the time- or spatial-derivatives in the highest order equation. These scale factors can be chosen such that waves propagate at exactly the speed of light in optically thin media. Alternatively, they may be chosen to significantly improve the accuracy of low-order solutions with no additional computational cost. The same scaling procedure and scale factors work in one- and multi-dimensions. In multidimensions, reducing the order of a solution can save significant amounts of computer time.
NASA Astrophysics Data System (ADS)
Walker, D. I.; Wolfand, J.; Wang, Y.; Bai, C.; Li, Y.; Abriola, L. M.; Pennell, K. D.
2011-12-01
Understanding the processes governing nanomaterial fate in subsurface environments is important due to their widespread use in commercial and manufacturing settings. To date, the majority of studies on nanoparticle transport in porous media have been conducted under idealized conditions in columns packed with uniform, clean sands. To evaluate the transport of nanoparticles in more representative subsurface systems, column and aquifer cell experiments were completed to investigate the transport of fullerene aggregates (nC60) in unwashed Ottawa sands and in the presence of low permeability lenses. To obtain independent measurements of attachment parameters, nC60 breakthrough and retention profiles were measured in 1-D columns (2.5 dia x 10.8 cm length) with three different size fractions (40-50 mesh, d50=335 μm; 80-100 mesh, d50=165 μm; and 100-140 mesh, d50=125 μm) of unwashed Ottawa sand. A three pore volume pulse of nC60 (ca. 2.5 mg/L) suspension containing 6 mM NaCl resulted in 71, 5 and 3 percent breakthrough in the 40-50, 80-100 and 100-140 mesh unwashed Ottawa sand, respectively. Solid phase retention profiles for all three size fractions exhibited hyper-exponential decline with distance from the column inlet, consistent with a physical straining process. The aquifer cells (63 x 36 x 1.4 cm) were packed with unwashed 40-50 mesh Ottawa sand as the background medium, in which rectangular lenses (2.4 cm x 9.6 cm) of lower permeability (80-100 and 100-140 mesh) were emplaced. Following completion of nonreactive tracer tests, 200-400 mL pulses of nC60 (ca. 2.5 mg/L) suspension were injected through side-ports up-gradient of the lenses. Only trace amounts of nC60 were detected in aqueous samples collected from down-gradient sampling ports. Post-experiment dissection of the aquifer cell demonstrated that solid phase concentrations of nC60 were largest surrounding the injection ports and decreased monotonically with no preferential accumulation observed at the
Edge Transport Modeling using the 3D EMC3-Eirene code on Tokamaks and Stellarators
NASA Astrophysics Data System (ADS)
Lore, J. D.; Ahn, J. W.; Briesemeister, A.; Ferraro, N.; Labombard, B.; McLean, A.; Reinke, M.; Shafer, M.; Terry, J.
2015-11-01
The fluid plasma edge transport code EMC3-Eirene has been applied to aid data interpretation and understanding the results of experiments with 3D effects on several tokamaks. These include applied and intrinsic 3D magnetic fields, 3D plasma facing components, and toroidally and poloidally localized heat and particle sources. On Alcator C-Mod, a series of experiments explored the impact of toroidally and poloidally localized impurity gas injection on core confinement and asymmetries in the divertor fluxes, with the differences between the asymmetry in L-mode and H-mode qualitatively reproduced in the simulations due to changes in the impurity ionization in the private flux region. Modeling of NSTX experiments on the effect of 3D fields on detachment matched the trend of a higher density at which the detachment occurs when 3D fields are applied. On DIII-D, different magnetic field models were used in the simulation and compared against the 2D Thomson scattering diagnostic. In simulating each device different aspects of the code model are tested pointing to areas where the model must be further developed. The application to stellarator experiments will also be discussed. Work supported by U.S. DOE: DE-AC05-00OR22725, DE AC02-09CH11466, DE-FC02-99ER54512, and DE-FC02-04ER54698.
Final Report for National Transport Code Collaboration PTRANSP
Arnold H. Kritz
2012-06-14
PTRANSP, which is the predictive version of the TRANSP code, was developed in a collaborative effort involving the Princeton Plasma Physics Laboratory, General Atomics Corporation, Lawrence Livermore National Laboratory, and Lehigh University. The PTRANSP/TRANSP suite of codes is the premier integrated tokamak modeling software in the United States. A production service for PTRANSP/TRANSP simulations is maintained at the Princeton Plasma Physics Laboratory; the server has a simple command line client interface and is subscribed to by about 100 researchers from tokamak projects in the US, Europe, and Asia. This service produced nearly 13000 PTRANSP/TRANSP simulations in the four year period FY 2005 through FY 2008. Major archives of TRANSP results are maintained at PPPL, MIT, General Atomics, and JET. Recent utilization, counting experimental analysis simulations as well as predictive simulations, more than doubled from slightly over 2000 simulations per year in FY 2005 and FY 2006 to over 4300 simulations per year in FY 2007 and FY 2008. PTRANSP predictive simulations applied to ITER increased eight fold from 30 simulations per year in FY 2005 and FY 2006 to 240 simulations per year in FY 2007 and FY 2008, accounting for more than half of combined PTRANSP/TRANSP service CPU resource utilization in FY 2008. PTRANSP studies focused on ITER played a key role in journal articles. Examples of validation studies carried out for momentum transport in PTRANSP simulations were presented at the 2008 IAEA conference. The increase in number of PTRANSP simulations has continued (more than 7000 TRANSP/PTRANSP simulations in 2010) and results of PTRANSP simulations appear in conference proceedings, for example the 2010 IAEA conference, and in peer reviewed papers. PTRANSP provides a bridge to the Fusion Simulation Program (FSP) and to the future of integrated modeling. Through years of widespread usage, each of the many parts of the PTRANSP suite of codes has been thoroughly