Kalman Filter Tracking on Parallel Architectures
NASA Astrophysics Data System (ADS)
Cerati, Giuseppe; Elmer, Peter; Lantz, Steven; McDermott, Kevin; Riley, Dan; Tadel, Matevž; Wittich, Peter; Würthwein, Frank; Yagil, Avi
2015-12-01
Power density constraints are limiting the performance improvements of modern CPUs. To address this we have seen the introduction of lower-power, multi-core processors, but the future will be even more exciting. In order to stay within the power density limits but still obtain Moore's Law performance/price gains, it will be necessary to parallelize algorithms to exploit larger numbers of lightweight cores and specialized functions like large vector units. Example technologies today include Intel's Xeon Phi and GPGPUs. Track finding and fitting is one of the most computationally challenging problems for event reconstruction in particle physics. At the High Luminosity LHC, for example, this will be by far the dominant problem. The need for greater parallelism has driven investigations of very different track finding techniques including Cellular Automata or returning to Hough Transform. The most common track finding techniques in use today are however those based on the Kalman Filter [2]. Significant experience has been accumulated with these techniques on real tracking detector systems, both in the trigger and offline. They are known to provide high physics performance, are robust and are exactly those being used today for the design of the tracking system for HL-LHC. Our previous investigations showed that, using optimized data structures, track fitting with Kalman Filter can achieve large speedup both with Intel Xeon and Xeon Phi. We report here our further progress towards an end-to-end track reconstruction algorithm fully exploiting vectorization and parallelization techniques in a realistic simulation setup.
Visualization and Tracking of Parallel CFD Simulations
NASA Technical Reports Server (NTRS)
Vaziri, Arsi; Kremenetsky, Mark
1995-01-01
We describe a system for interactive visualization and tracking of a 3-D unsteady computational fluid dynamics (CFD) simulation on a parallel computer. CM/AVS, a distributed, parallel implementation of a visualization environment (AVS) runs on the CM-5 parallel supercomputer. A CFD solver is run as a CM/AVS module on the CM-5. Data communication between the solver, other parallel visualization modules, and a graphics workstation, which is running AVS, are handled by CM/AVS. Partitioning of the visualization task, between CM-5 and the workstation, can be done interactively in the visual programming environment provided by AVS. Flow solver parameters can also be altered by programmable interactive widgets. This system partially removes the requirement of storing large solution files at frequent time steps, a characteristic of the traditional 'simulate (yields) store (yields) visualize' post-processing approach.
Parallel transport on principal bundles over stacks
NASA Astrophysics Data System (ADS)
Collier, Brian; Lerman, Eugene; Wolbert, Seth
2016-09-01
In this paper we introduce a notion of parallel transport for principal bundles with connections over differentiable stacks. We show that principal bundles with connections over stacks can be recovered from their parallel transport thereby extending the results of Barrett, Caetano and Picken, and Schreiber and Waldorf from manifolds to stacks. In the process of proving our main result we simplify Schreiber and Waldorf's original definition of a transport functor for principal bundles with connections over manifolds and provide a more direct proof of the correspondence between principal bundles with connections and transport functors.
An optimization-based parallel particle filter for multitarget tracking
NASA Astrophysics Data System (ADS)
Sutharsan, S.; Sinha, A.; Kirubarajan, T.; Farooq, M.
2005-09-01
Particle filter based estimation is becoming more popular because it has the capability to effectively solve nonlinear and non-Gaussian estimation problems. However, the particle filter has high computational requirements and the problem becomes even more challenging in the case of multitarget tracking. In order to perform data association and estimation jointly, typically an augmented state vector of target dynamics is used. As the number of targets increases, the computation required for each particle increases exponentially. Thus, parallelization is a possibility in order to achieve the real time feasibility in large-scale multitarget tracking applications. In this paper, we present a real-time feasible scheduling algorithm that minimizes the total computation time for the bus connected heterogeneous primary-secondary architecture. This scheduler is capable of selecting the optimal number of processors from a large pool of secondary processors and mapping the particles among the selected processors. Furthermore, we propose a less communication intensive parallel implementation of the particle filter without sacrificing tracking accuracy using an efficient load balancing technique, in which optimal particle migration is ensured. In this paper, we present the mathematical formulations for scheduling the particles as well as for particle migration via load balancing. Simulation results show the tracking performance of our parallel particle filter and the speedup achieved using parallelization.
Tracking Connections: An Exercise about Series and Parallel Resistances
ERIC Educational Resources Information Center
Jankovic, Srdjan
2010-01-01
Unlike many other topics in basic physics, series and parallel resistances are rarely noticed in the real life of an ordinary individual, making it difficult to design a laboratory activity that can simulate something familiar. The activities described here entail minimal costs and are based on a puzzle-like game of tracking wire connections. A…
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.
Tracking Connections: An Exercise about Series and Parallel Resistances
NASA Astrophysics Data System (ADS)
Jankovic, Srdjan
2010-09-01
Unlike many other topics in basic physics, series and parallel resistances are rarely noticed in the real life of an ordinary individual, making it difficult to design a laboratory activity that can simulate something familiar. The activities described here entail minimal costs and are based on a puzzle-like game of tracking wire connections. A simple resistor-based device is built by students, which enables them to use a common multimeter to track down wire connections in a set of unmarked wires. A similar approach is sometimes used by electricians to identify wire connections.
An experiment in hurricane track prediction using parallel computing methods
NASA Technical Reports Server (NTRS)
Song, Chang G.; Jwo, Jung-Sing; Lakshmivarahan, S.; Dhall, S. K.; Lewis, John M.; Velden, Christopher S.
1994-01-01
The barotropic model is used to explore the advantages of parallel processing in deterministic forecasting. We apply this model to the track forecasting of hurricane Elena (1985). In this particular application, solutions to systems of elliptic equations are the essence of the computational mechanics. One set of equations is associated with the decomposition of the wind into irrotational and nondivergent components - this determines the initial nondivergent state. Another set is associated with recovery of the streamfunction from the forecasted vorticity. We demonstrate that direct parallel methods based on accelerated block cyclic reduction (BCR) significantly reduce the computational time required to solve the elliptic equations germane to this decomposition and forecast problem. A 72-h track prediction was made using incremental time steps of 16 min on a network of 3000 grid points nominally separated by 100 km. The prediction took 30 sec on the 8-processor Alliant FX/8 computer. This was a speed-up of 3.7 when compared to the one-processor version. The 72-h prediction of Elena's track was made as the storm moved toward Florida's west coast. Approximately 200 km west of Tampa Bay, Elena executed a dramatic recurvature that ultimately changed its course toward the northwest. Although the barotropic track forecast was unable to capture the hurricane's tight cycloidal looping maneuver, the subsequent northwesterly movement was accurately forecasted as was the location and timing of landfall near Mobile Bay.
Parallel transport and defects on nematic shells
NASA Astrophysics Data System (ADS)
Rosso, Riccardo; Virga, Epifanio G.; Kralj, Samo
2012-11-01
Nematic shells are thin films of nematic liquid crystal deposited on the boundary of colloidal particles, where liquid crystal molecules may freely glide, while remaining tangent to the surface substrate. The surface nematic order is described here by an appropriate tensor field Q, which vanishes wherever a defect occurs in the molecular order. We show how the classical concept of parallel transport on a manifold introduced by Levi-Civita can be adapted to this setting to define the topological charge m of a defect. We arrive at a simple formula to compute m from a generic representation of Q. In a number of separate appendices, we revisit in a unified language several, apparently disparate applications of Levi-Civita's parallel transport.
Nonisothermal multiphase subsurface transport on parallel computers
Martinez, M.J.; Hopkins, P.L.; Shadid, J.N.
1997-10-01
We present a numerical method for nonisothermal, multiphase subsurface transport in heterogeneous porous media. The mathematical model considers nonisothermal two-phase (liquid/gas) flow, including capillary pressure effects, binary diffusion in the gas phase, conductive, latent, and sensible heat transport. The Galerkin finite element method is used for spatial discretization, and temporal integration is accomplished via a predictor/corrector scheme. Message-passing and domain decomposition techniques are used for implementing a scalable algorithm for distributed memory parallel computers. An illustrative application is shown to demonstrate capabilities and performance.
Traditional Tracking with Kalman Filter on Parallel Architectures
NASA Astrophysics Data System (ADS)
Cerati, Giuseppe; Elmer, Peter; Lantz, Steven; MacNeill, Ian; McDermott, Kevin; Riley, Dan; Tadel, Matevž; Wittich, Peter; Würthwein, Frank; Yagil, Avi
2015-05-01
Power density constraints are limiting the performance improvements of modern CPUs. To address this, we have seen the introduction of lower-power, multi-core processors, but the future will be even more exciting. In order to stay within the power density limits but still obtain Moore's Law performance/price gains, it will be necessary to parallelize algorithms to exploit larger numbers of lightweight cores and specialized functions like large vector units. Example technologies today include Intel's Xeon Phi and GPGPUs. Track finding and fitting is one of the most computationally challenging problems for event reconstruction in particle physics. At the High Luminosity LHC, for example, this will be by far the dominant problem. The most common track finding techniques in use today are however those based on the Kalman Filter. Significant experience has been accumulated with these techniques on real tracking detector systems, both in the trigger and offline. We report the results of our investigations into the potential and limitations of these algorithms on the new parallel hardware.
Monte Carlo radiation transport¶llelism
Cox, L. J.; Post, S. E.
2002-01-01
This talk summarizes the main aspects of the LANL ASCI Eolus project and its major unclassified code project, MCNP. The MCNP code provide a state-of-the-art Monte Carlo radiation transport to approximately 3000 users world-wide. Almost all hardware platforms are supported because we strictly adhere to the FORTRAN-90/95 standard. For parallel processing, MCNP uses a mixture of OpenMp combined with either MPI or PVM (shared and distributed memory). This talk summarizes our experiences on various platforms using MPI with and without OpenMP. These platforms include PC-Windows, Intel-LINUX, BlueMountain, Frost, ASCI-Q and others.
TRAMP. Transport of Metallic Fission Products Along Multiple Parallel Paths
Hudritsch, W.; Richards, M.
1991-11-01
TRAMP is used to calculate the transport of metallic fission products along multiple parallel paths; the primary application is transport in and release from nuclear-grade graphite. The transport mechanisms are concentration-driven diffusion, thermal diffusion, and convection.
Beam dynamics calculations and particle tracking using massively parallel processors
Ryne, R.D.; Habib, S.
1995-12-31
During the past decade massively parallel processors (MPPs) have slowly gained acceptance within the scientific community. At present these machines typically contain a few hundred to one thousand off-the-shelf microprocessors and a total memory of up to 32 GBytes. The potential performance of these machines is illustrated by the fact that a month long job on a high end workstation might require only a few hours on an MPP. The acceptance of MPPs has been slow for a variety of reasons. For example, some algorithms are not easily parallelizable. Also, in the past these machines were difficult to program. But in recent years the development of Fortran-like languages such as CM Fortran and High Performance Fortran have made MPPs much easier to use. In the following we will describe how MPPs can be used for beam dynamics calculations and long term particle tracking.
Tracking moving radar targets with parallel, velocity-tuned filters
Bickel, Douglas L.; Harmony, David W.; Bielek, Timothy P.; Hollowell, Jeff A.; Murray, Margaret S.; Martinez, Ana
2013-04-30
Radar data associated with radar illumination of a movable target is processed to monitor motion of the target. A plurality of filter operations are performed in parallel on the radar data so that each filter operation produces target image information. The filter operations are defined to have respectively corresponding velocity ranges that differ from one another. The target image information produced by one of the filter operations represents the target more accurately than the target image information produced by the remainder of the filter operations when a current velocity of the target is within the velocity range associated with the one filter operation. In response to the current velocity of the target being within the velocity range associated with the one filter operation, motion of the target is tracked based on the target image information produced by the one filter operation.
Particle Transport in Parallel-Plate Reactors
Rader, D.J.; Geller, A.S.
1999-08-01
A major cause of semiconductor yield degradation is contaminant particles that deposit on wafers while they reside in processing tools during integrated circuit manufacturing. This report presents numerical models for assessing particle transport and deposition in a parallel-plate geometry characteristic of a wide range of single-wafer processing tools: uniform downward flow exiting a perforated-plate showerhead separated by a gap from a circular wafer resting on a parallel susceptor. Particles are assumed to originate either upstream of the showerhead or from a specified position between the plates. The physical mechanisms controlling particle deposition and transport (inertia, diffusion, fluid drag, and external forces) are reviewed, with an emphasis on conditions encountered in semiconductor process tools (i.e., sub-atmospheric pressures and submicron particles). Isothermal flow is assumed, although small temperature differences are allowed to drive particle thermophoresis. Numerical solutions of the flow field are presented which agree with an analytic, creeping-flow expression for Re < 4. Deposition is quantified by use of a particle collection efficiency, which is defined as the fraction of particles in the reactor that deposit on the wafer. Analytic expressions for collection efficiency are presented for the limiting case where external forces control deposition (i.e., neglecting particle diffusion and inertia). Deposition from simultaneous particle diffusion and external forces is analyzed by an Eulerian formulation; for creeping flow and particles released from a planar trap, the analysis yields an analytic, integral expression for particle deposition based on process and particle properties. Deposition from simultaneous particle inertia and external forces is analyzed by a Lagrangian formulation, which can describe inertia-enhanced deposition resulting from particle acceleration in the showerhead. An approximate analytic expression is derived for particle
Radiation transport on unstructured mesh with parallel computers
Fan, W.C.; Drumm, C.R.
2000-07-01
This paper summarizes the developmental work on a deterministic transport code that provides multidimensional radiation transport capabilities on an unstructured mesh. The second-order form of the Boltzmann transport equation is solved utilizing the discrete ordinates angular differencing and the Galerkin finite element spatial differencing. The discretized system, which couples the spatial-angular dependence, is solved simultaneously using a parallel conjugate-gradient (CG) iterative solver. This approach eliminates the need for the conventional inner iterations over the discrete directions and is well-suited for massively parallel computers.
Parallel heat transport in integrable and chaotic magnetic fields
Castillo-Negrete, D. del; Chacon, L.
2012-05-15
The study of transport in magnetized plasmas is a problem of fundamental interest in controlled fusion, space plasmas, and astrophysics research. Three issues make this problem particularly challenging: (i) The extreme anisotropy between the parallel (i.e., along the magnetic field), {chi}{sub ||} , and the perpendicular, {chi}{sub Up-Tack }, conductivities ({chi}{sub ||} /{chi}{sub Up-Tack} may exceed 10{sup 10} in fusion plasmas); (ii) Nonlocal parallel transport in the limit of small collisionality; and (iii) Magnetic field lines chaos which in general complicates (and may preclude) the construction of magnetic field line coordinates. Motivated by these issues, we present a Lagrangian Green's function method to solve the local and non-local parallel transport equation applicable to integrable and chaotic magnetic fields in arbitrary geometry. The method avoids by construction the numerical pollution issues of grid-based algorithms. The potential of the approach is demonstrated with nontrivial applications to integrable (magnetic island), weakly chaotic (Devil's staircase), and fully chaotic magnetic field configurations. For the latter, numerical solutions of the parallel heat transport equation show that the effective radial transport, with local and non-local parallel closures, is non-diffusive, thus casting doubts on the applicability of quasilinear diffusion descriptions. General conditions for the existence of non-diffusive, multivalued flux-gradient relations in the temperature evolution are derived.
Shift: A Massively Parallel Monte Carlo Radiation Transport Package
Pandya, Tara M; Johnson, Seth R; Davidson, Gregory G; Evans, Thomas M; Hamilton, Steven P
2015-01-01
This paper discusses the massively-parallel Monte Carlo radiation transport package, Shift, developed at Oak Ridge National Laboratory. It reviews the capabilities, implementation, and parallel performance of this code package. Scaling results demonstrate very good strong and weak scaling behavior of the implemented algorithms. Benchmark results from various reactor problems show that Shift results compare well to other contemporary Monte Carlo codes and experimental results.
Local and nonlocal parallel heat transport in general magnetic fields
Del-Castillo-Negrete, Diego B; Chacon, Luis
2011-01-01
A novel approach for the study of parallel transport in magnetized plasmas is presented. The method avoids numerical pollution issues of grid-based formulations and applies to integrable and chaotic magnetic fields with local or nonlocal parallel closures. In weakly chaotic fields, the method gives the fractal structure of the devil's staircase radial temperature profile. In fully chaotic fields, the temperature exhibits self-similar spatiotemporal evolution with a stretched-exponential scaling function for local closures and an algebraically decaying one for nonlocal closures. It is shown that, for both closures, the effective radial heat transport is incompatible with the quasilinear diffusion model.
Parallel heat transport in integrable and chaotic magnetic fields
Del-Castillo-Negrete, Diego B; Chacon, Luis
2012-01-01
The study of transport in magnetized plasmas is a problem of fundamental interest in controlled fusion, space plasmas, and astrophysics research. Three issues make this problem particularly chal- lenging: (i) The extreme anisotropy between the parallel (i.e., along the magnetic field), , and the perpendicular, , conductivities ( / may exceed 1010 in fusion plasmas); (ii) Magnetic field lines chaos which in general complicates (and may preclude) the construction of magnetic field line coordinates; and (iii) Nonlocal parallel transport in the limit of small collisionality. Motivated by these issues, we present a Lagrangian Green s function method to solve the local and non-local parallel transport equation applicable to integrable and chaotic magnetic fields in arbitrary geom- etry. The method avoids by construction the numerical pollution issues of grid-based algorithms. The potential of the approach is demonstrated with nontrivial applications to integrable (magnetic island chain), weakly chaotic (devil s staircase), and fully chaotic magnetic field configurations. For the latter, numerical solutions of the parallel heat transport equation show that the effective radial transport, with local and non-local closures, is non-diffusive, thus casting doubts on the appropriateness of the applicability of quasilinear diffusion descriptions. General conditions for the existence of non-diffusive, multivalued flux-gradient relations in the temperature evolution are derived.
Parallel heat transport in reversed shear magnetic field configurations
NASA Astrophysics Data System (ADS)
Blazevski, D.; Del-Castillo-Negrete, D.
2012-03-01
Transport in magnetized plasmas is a key problem in controlled fusion, space plasmas, and astrophysics. Three issues make this problem particularly challenging: (i) The extreme anisotropy between the parallel (i.e., along the magnetic field), χ, and the perpendicular, χ, conductivities (χ/χ may exceed 10^10 in fusion plasmas); (ii) Magnetic field lines chaos; and (iii) Nonlocal parallel transport. We have recently developed a Lagrangian Green's function (LG) method to solve the local and non-local parallel (χ/χ->∞) transport equation applicable to integrable and chaotic magnetic fields. footnotetext D. del-Castillo-Negrete, L. Chac'on, PRL, 106, 195004 (2011); D. del-Castillo-Negrete, L. Chac'on, Phys. Plasmas, APS Invited paper, submitted (2011). The proposed method overcomes many of the difficulties faced by standard finite different methods related to the three issues mentioned above. Here we apply the LG method to study transport in reversed shear configurations. We focus on the following problems: (i) separatrix reconnection of magnetic islands and transport; (ii) robustness of shearless, q'=0, transport barriers; (iii) leaky barriers and shearless Cantori.
Measures of effectiveness for BMD mid-course tracking on MIMD massively parallel computers
VanDyke, J.P.; Tomkins, J.L.; Furnish, M.D.
1995-05-01
The TRC code, a mid-course tracking code for ballistic missiles, has previously been implemented on a 1024-processor MIMD (Multiple Instruction -- Multiple Data) massively parallel computer. Measures of Effectiveness (MOE) for this algorithm have been developed for this computing environment. The MOE code is run in parallel with the TRC code. Particularly useful MOEs include the number of missed objects (real objects for which the TRC algorithm did not construct a track); of ghost tracks (tracks not corresponding to a real object); of redundant tracks (multiple tracks corresponding to a single real object); and of unresolved objects (multiple objects corresponding to a single track). All of these are expressed as a function of time, and tend to maximize during the time in which real objects are spawned (multiple reentry vehicles per post-boost vehicle). As well, it is possible to measure the track-truth separation as a function of time. A set of calculations is presented illustrating these MOEs as a function of time for a case with 99 post-boost vehicles, each of which spawns 9 reentry vehicles.
CONTAMINANT TRANSPORT IN PARALLEL FRACTURED MEDIA: SUDICKY AND FRIND REVISITED
This paper is concerned with a modified, nondimensional form of the parallel fracture, contaminant transport model of Sudicky and Frind (1982). The modifications include the boundary condition at the fracture wall, expressed by a parameter, and the power-law relationship between...
CONTAMINANT TRANSPORT IN PARALLEL FRACTURED MEDIA: SUDICKY AND FRIND REVISITED
This paper is concerned with a modified, nondimensional form of the parallel fracture, contaminant transport model of Sudicky and Frind (1982). The modifications include the boundary condition at the fracture wall, expressed by a parameter , and the power-law relationship betwe...
A parallel, volume-tracking algorithm for unstructured meshes
Mosso, S.J.; Swartz, B.K.; Kothe, D.B.; Ferrell, R.C.
1996-10-01
Many diverse areas of industry benefit from the use of volume of fluid methods to predict the movement of materials. Casting is a common method of part fabrication. The accurate prediction of the casting process is pivotal to industry. Mold design and casting is currently considered an art by industry. It typically involves many trial mold designs, and the rejection of defective parts is costly. Failure of cast parts, because residual stresses reduce the part`s strength, can be catastrophic. Cast parts should have precise geometric details that reduce or eliminate the need for machining after casting. Volume of fluid codes will help designers predict how the molten metal fills a mold and where ay trapped voids remain. Prediction of defects due to thermal contraction or expansion will eliminate defective, trial mold designs and speed the parts to market with fewer rejections. Increasing the predictability and therefore the accuracy of the casting process will reduce the art that is involved in mold design and parts casting. Here, recent enhancements to multidimensional volume-tracking algorithms are presented. Illustrations in two dimensions are given. The improvements include new, local algorithms for interface normal constructions and a new full remapping algorithm for time integration. These methods are used on structured and unstructured grids.
Transport of parallel momentum by drift-Alfven turbulence
McDevitt, C. J.; Diamond, P. H.
2009-01-15
An electromagnetic gyrokinetic formulation is utilized to calculate the turbulent radial flux of parallel momentum for a strongly magnetized plasma in the large aspect ratio limit. For low-{beta} plasmas, excluding regions of steep density gradients, the level of momentum transport induced by microturbulence is found to be well described within the electrostatic approximation. However, near regions of steep equilibrium profile gradients, strong electromagnetic contributions to the momentum flux are predicted. In particular, for sufficiently steep density gradient, the magnitude of transport induced by the off-diagonal residual stress component of the momentum flux induced by drift wave turbulence can be quenched. This quenching mechanism, which results from shielding of the parallel electric field by the inductive term, is distinct from ExB shear decorrelation, since it allows for the level of off-diagonal turbulent transport to be strongly reduced without extinguishing the underlying microturbulence. In contrast, the level of transport induced by a given Alfvenic branch of the drift-Alfven dispersion relationship typically increases as the density gradient steepens, allowing an alternate channel for momentum transport. A calculation of the momentum transport induced by Alfvenic turbulence in a homogeneous medium suggests that an imbalance in Elsasser populations is required in order to introduce a finite level of off-diagonal momentum transport for the case of the simplified geometry considered.
High-resolution reactive transport: A coupled parallel hydrogeochemical model
NASA Astrophysics Data System (ADS)
Beisman, J. J.; Maxwell, R. M.; Steefel, C. I.; Sitchler, A.; Molins, S.
2013-12-01
Subsurface hydrogeochemical systems are an especially complex component of the terrestrial environment and play host to a multitude of interactions. Parameterizations of these interactions are perhaps the least understood component of terrestrial systems, presenting uncertainties in the predictive understanding of biogeochemical cycling and transport. Thorough knowledge of biogeochemical transport processes is critical to the quantification of carbon/nutrient fluxes in the subsurface, and to the development of effective contaminant remediation techniques. Here we present a coupled parallel hydrogeochemical model, ParCrunchFlow, as a tool to further our understanding of governing processes and interactions in natural hydrogeochemical systems. ParCrunchFlow is a coupling of the reactive transport simulator CrunchFlow with the hydrologic model ParFlow. CrunchFlow is a multicomponent reactive flow and transport code that can be used to simulate a range of important processes and environments, including reactive contaminant transport, chemical weathering, carbon sequestration, biogeochemical cycling, and water-rock interaction. ParFlow is a parallel, three-dimensional, variably-saturated, coupled surface-subsurface flow and transport code with the ability to simulate complex topography, geology, and heterogeneity. ParCrunchflow takes advantage of the efficient parallelism built into Parflow, allowing the numerical simulation of reactive transport processes in chemically and physically heterogeneous media at high spatial resolutions. This model provides an ability to further examine the interactions and feedbacks between biogeochemical systems and complex subsurface flow fields. In addition to the details of model construction, results will be presented that show floodplain nutrient cycling and the effects of heterogeneity on small-scale mixing reactions at the Department of Energy's Old Rifle Legacy site.
Fast parallel tracking algorithm for the muon detector of the CBM experiment at fair
NASA Astrophysics Data System (ADS)
Lebedev, A.; Höhne, C.; Kisel, I.; Ososkov, G.
2010-07-01
Particle trajectory recognition is an important and challenging task in the Compressed Baryonic Matter (CBM) experiment at the future FAIR accelerator at Darmstadt. The tracking algorithms have to process terabytes of input data produced in particle collisions. Therefore, the speed of the tracking software is extremely important for data analysis. In this contribution, a fast parallel track reconstruction algorithm which uses available features of modern processors is presented. These features comprise a SIMD instruction set (SSE) and multithreading. The first allows one to pack several data items into one register and to operate on all of them in parallel thus achieving more operations per cycle. The second feature enables the routines to exploit all available CPU cores and hardware threads. This parallel version of the tracking algorithm has been compared to the initial serial scalar version which uses a similar approach for tracking. A speed-up factor of 487 was achieved (from 730 to 1.5 ms/event) for a computer with 2 × Intel Core i7 processors at 2.66 GHz.
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.
Efficient parallel algorithm for statistical ion track simulations in crystalline materials
NASA Astrophysics Data System (ADS)
Jeon, Byoungseon; Grønbech-Jensen, Niels
2009-02-01
We present an efficient parallel algorithm for statistical Molecular Dynamics simulations of ion tracks in solids. The method is based on the Rare Event Enhanced Domain following Molecular Dynamics (REED-MD) algorithm, which has been successfully applied to studies of, e.g., ion implantation into crystalline semiconductor wafers. We discuss the strategies for parallelizing the method, and we settle on a host-client type polling scheme in which a multiple of asynchronous processors are continuously fed to the host, which, in turn, distributes the resulting feed-back information to the clients. This real-time feed-back consists of, e.g., cumulative damage information or statistics updates necessary for the cloning in the rare event algorithm. We finally demonstrate the algorithm for radiation effects in a nuclear oxide fuel, and we show the balanced parallel approach with high parallel efficiency in multiple processor configurations.
Three-degree-of-freedom parallel manipulator to track the sun for concentrated solar power systems
NASA Astrophysics Data System (ADS)
Ashith Shyam, R. B.; Ghosal, A.
2015-07-01
In concentrated solar power(CSP) generating stations, incident solar energy is reflected from a large number of mirrors or heliostats to a faraway receiver. In typical CSP installations, the mirror needs to be moved about two axes independently using two actuators in series with the mirror effectively mounted at a single point. A three degree-of-freedom parallel manipulator, namely the 3-RPS parallel manipulator, is proposed to track the sun. The proposed 3-RPS parallel manipulator supports the load of the mirror, structure and wind loading at three points resulting in less deflection, and thus a much larger mirror can be moved with the required tracking accuracy and without increasing the weight of the support structure. The kinematics equations to determine motion of the actuated prismatic joints in the 3-RPS parallel manipulator such that the sun's rays are reflected on to a stationary receiver are developed. Using finite element analysis, it is shown that for same sized mirror, wind loading and maximum deflection requirement, the weight of the support structure is between 15% and 60% less with the 3-RPS parallel manipulator when compared to azimuth-elevation or the target-aligned configurations.
Provably optimal parallel transport sweeps on regular grids
Adams, M. P.; Adams, M. L.; Hawkins, W. D.; Smith, T.; Rauchwerger, L.; Amato, N. M.; Bailey, T. S.; Falgout, R. D.
2013-07-01
We have found provably optimal algorithms for full-domain discrete-ordinate transport sweeps on regular grids in 3D Cartesian geometry. We describe these algorithms and sketch a 'proof that they always execute the full eight-octant sweep in the minimum possible number of stages for a given P{sub x} x P{sub y} x P{sub z} partitioning. Computational results demonstrate that our optimal scheduling algorithms execute sweeps in the minimum possible stage count. Observed parallel efficiencies agree well with our performance model. An older version of our PDT transport code achieves almost 80% parallel efficiency on 131,072 cores, on a weak-scaling problem with only one energy group, 80 directions, and 4096 cells/core. A newer version is less efficient at present-we are still improving its implementation - but achieves almost 60% parallel efficiency on 393,216 cores. These results conclusively demonstrate that sweeps can perform with high efficiency on core counts approaching 10{sup 6}. (authors)
Pelegant : a parallel accelerator simulation code for electron generation and tracking.
Wang, Y.; Borland, M. D.; Accelerator Systems Division
2006-01-01
elegant is a general-purpose code for electron accelerator simulation that has a worldwide user base. Recently, many of the time-intensive elements were parallelized using MPI. Development has used modest Linux clusters and the BlueGene/L supercomputer at Argonne National Laboratory. This has provided very good performance for some practical simulations, such as multiparticle tracking with synchrotron radiation and emittance blow-up in the vertical rf kick scheme. The effort began with development of a concept that allowed for gradual parallelization of the code, using the existing beamline-element classification table in elegant. This was crucial as it allowed parallelization without major changes in code structure and without major conflicts with the ongoing evolution of elegant. Because of rounding error and finite machine precision, validating a parallel program against a uniprocessor program with the requirement of bitwise identical results is notoriously difficult. We will report validating simulation results of parallel elegant against those of serial elegant by applying Kahan's algorithm to improve accuracy dramatically for both versions. The quality of random numbers in a parallel implementation is very important for some simulations. Some practical experience with generating parallel random numbers by offsetting the seed of each random sequence according to the processor ID will be reported.
Novel Parallel Numerical Methods for Radiation& Neutron Transport
Brown, P N
2001-03-06
In many of the multiphysics simulations performed at LLNL, transport calculations can take up 30 to 50% of the total run time. If Monte Carlo methods are used, the percentage can be as high as 80%. Thus, a significant core competence in the formulation, software implementation, and solution of the numerical problems arising in transport modeling is essential to Laboratory and DOE research. In this project, we worked on developing scalable solution methods for the equations that model the transport of photons and neutrons through materials. Our goal was to reduce the transport solve time in these simulations by means of more advanced numerical methods and their parallel implementations. These methods must be scalable, that is, the time to solution must remain constant as the problem size grows and additional computer resources are used. For iterative methods, scalability requires that (1) the number of iterations to reach convergence is independent of problem size, and (2) that the computational cost grows linearly with problem size. We focused on deterministic approaches to transport, building on our earlier work in which we performed a new, detailed analysis of some existing transport methods and developed new approaches. The Boltzmann equation (the underlying equation to be solved) and various solution methods have been developed over many years. Consequently, many laboratory codes are based on these methods, which are in some cases decades old. For the transport of x-rays through partially ionized plasmas in local thermodynamic equilibrium, the transport equation is coupled to nonlinear diffusion equations for the electron and ion temperatures via the highly nonlinear Planck function. We investigated the suitability of traditional-solution approaches to transport on terascale architectures and also designed new scalable algorithms; in some cases, we investigated hybrid approaches that combined both.
Nonlinear parallel momentum transport in strong electrostatic turbulence
NASA Astrophysics Data System (ADS)
Wang, Lu; Wen, Tiliang; Diamond, P. H.
2015-05-01
Most existing theoretical studies of momentum transport focus on calculating the Reynolds stress based on quasilinear theory, without considering the nonlinear momentum flux- ⟨ v ˜ r n ˜ u ˜ ∥ ⟩ . However, a recent experiment on TORPEX found that the nonlinear toroidal momentum flux induced by blobs makes a significant contribution as compared to the Reynolds stress [Labit et al., Phys. Plasmas 18, 032308 (2011)]. In this work, the nonlinear parallel momentum flux in strong electrostatic turbulence is calculated by using a three dimensional Hasegawa-Mima equation, which is relevant for tokamak edge turbulence. It is shown that the nonlinear diffusivity is smaller than the quasilinear diffusivity from Reynolds stress. However, the leading order nonlinear residual stress can be comparable to the quasilinear residual stress, and so may be important to intrinsic rotation in tokamak edge plasmas. A key difference from the quasilinear residual stress is that parallel fluctuation spectrum asymmetry is not required for nonlinear residual stress.
Parallel Monte Carlo Synthetic Acceleration methods for discrete transport problems
NASA Astrophysics Data System (ADS)
Slattery, Stuart R.
This work researches and develops Monte Carlo Synthetic Acceleration (MCSA) methods as a new class of solution techniques for discrete neutron transport and fluid flow problems. Monte Carlo Synthetic Acceleration methods use a traditional Monte Carlo process to approximate the solution to the discrete problem as a means of accelerating traditional fixed-point methods. To apply these methods to neutronics and fluid flow and determine the feasibility of these methods on modern hardware, three complementary research and development exercises are performed. First, solutions to the SPN discretization of the linear Boltzmann neutron transport equation are obtained using MCSA with a difficult criticality calculation for a light water reactor fuel assembly used as the driving problem. To enable MCSA as a solution technique a group of modern preconditioning strategies are researched. MCSA when compared to conventional Krylov methods demonstrated improved iterative performance over GMRES by converging in fewer iterations when using the same preconditioning. Second, solutions to the compressible Navier-Stokes equations were obtained by developing the Forward-Automated Newton-MCSA (FANM) method for nonlinear systems based on Newton's method. Three difficult fluid benchmark problems in both convective and driven flow regimes were used to drive the research and development of the method. For 8 out of 12 benchmark cases, it was found that FANM had better iterative performance than the Newton-Krylov method by converging the nonlinear residual in fewer linear solver iterations with the same preconditioning. Third, a new domain decomposed algorithm to parallelize MCSA aimed at leveraging leadership-class computing facilities was developed by utilizing parallel strategies from the radiation transport community. The new algorithm utilizes the Multiple-Set Overlapping-Domain strategy in an attempt to reduce parallel overhead and add a natural element of replication to the algorithm. It
Final Report - Ferrographic Tracking of Bacterial Transport
William P. Johnson
2002-10-10
The work performed during the past three years has been extremely productive. Ferrographic capture was utilized in analysis of several thousand field samples collected from arrays of multilevel samplers during three intensive field campaigns conducted at two shallow sandy aquifer sites in Oyster, VA. This work has shown resulted in three important conclusions: (1) Ferrographic capture provides unparalleled low quantitation limits for bacterial cell enumeration (Johnson et al., 2000). (2) The high-resolution analyses provided by ferrographic capture allowed observation of increased bacterial removal rates (from groundwater) that corresponded to increased populations of protozoa in the groundwater (Zhang et al., 2001). This novel data allowed determination of bacterial predation rates by protists in the field, a consideration that will be important for successful bioaugmentation strategies. (3) The high-resolution analyses provided by ferrographic capture allowed observation of detachment of indigenous cells in response to breakthrough of injected cells in groundwater (Johnson et al., 2001). The implication of this unique observation is that bacterial transport, specifically bacterial attachment and detachment, may be much more dynamic than has been indicated by short-term laboratory and field studies. Dynamic attachment and detachment of bacteria in groundwater may lead to greatly increased transport distances over long terms relative to what has been indicated by short-term laboratory and field studies.
Improved parallel solution techniques for the integral transport matrix method
Zerr, Robert J; Azmy, Yousry Y
2010-11-23
Alternative solution strategies to the parallel block Jacobi (PBJ) method for the solution of the global problem with the integral transport matrix method operators have been designed and tested. The most straightforward improvement to the Jacobi iterative method is the Gauss-Seidel alternative. The parallel red-black Gauss-Seidel (PGS) algorithm can improve on the number of iterations and reduce work per iteration by applying an alternating red-black color-set to the subdomains and assigning multiple sub-domains per processor. A parallel GMRES(m) method was implemented as an alternative to stationary iterations. Computational results show that the PGS method can improve on the PBJ method execution by up to {approx}50% when eight sub-domains per processor are used. However, compared to traditional source iterations with diffusion synthetic acceleration, it is still approximately an order of magnitude slower. The best-performing case are opticaUy thick because sub-domains decouple, yielding faster convergence. Further tests revealed that 64 sub-domains per processor was the best performing level of sub-domain division. An acceleration technique that improves the convergence rate would greatly improve the ITMM. The GMRES(m) method with a diagonal block preconditioner consumes approximately the same time as the PBJ solver but could be improved by an as yet undeveloped, more efficient preconditioner.
Particle Tracking Model and Abstraction of Transport Processes
B. Robinson
2000-04-07
The purpose of the transport methodology and component analysis is to provide the numerical methods for simulating radionuclide transport and model setup for transport in the unsaturated zone (UZ) site-scale model. The particle-tracking method of simulating radionuclide transport is incorporated into the FEHM computer code and the resulting changes in the FEHM code are to be submitted to the software configuration management system. This Analysis and Model Report (AMR) outlines the assumptions, design, and testing of a model for calculating radionuclide transport in the unsaturated zone at Yucca Mountain. In addition, methods for determining colloid-facilitated transport parameters are outlined for use in the Total System Performance Assessment (TSPA) analyses. Concurrently, process-level flow model calculations are being carrier out in a PMR for the unsaturated zone. The computer code TOUGH2 is being used to generate three-dimensional, dual-permeability flow fields, that are supplied to the Performance Assessment group for subsequent transport simulations. These flow fields are converted to input files compatible with the FEHM code, which for this application simulates radionuclide transport using the particle-tracking algorithm outlined in this AMR. Therefore, this AMR establishes the numerical method and demonstrates the use of the model, but the specific breakthrough curves presented do not necessarily represent the behavior of the Yucca Mountain unsaturated zone.
Ruh, Dominic; Tränkle, Benjamin; Rohrbach, Alexander
2011-10-24
Multi-dimensional, correlated particle tracking is a key technology to reveal dynamic processes in living and synthetic soft matter systems. In this paper we present a new method for tracking micron-sized beads in parallel and in all three dimensions - faster and more precise than existing techniques. Using an acousto-optic deflector and two quadrant-photo-diodes, we can track numerous optically trapped beads at up to tens of kHz with a precision of a few nanometers by back-focal plane interferometry. By time-multiplexing the laser focus, we can calibrate individually all traps and all tracking signals in a few seconds and in 3D. We show 3D histograms and calibration constants for nine beads in a quadratic arrangement, although trapping and tracking is easily possible for more beads also in arbitrary 2D arrangements. As an application, we investigate the hydrodynamic coupling and diffusion anomalies of spheres trapped in a 3 × 3 arrangement. PMID:22109012
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
Parallel Transport Quantum Logic Gates with Trapped Ions.
de Clercq, Ludwig E; Lo, Hsiang-Yu; Marinelli, Matteo; Nadlinger, David; Oswald, Robin; Negnevitsky, Vlad; Kienzler, Daniel; Keitch, Ben; Home, Jonathan P
2016-02-26
We demonstrate single-qubit operations by transporting a beryllium ion with a controlled velocity through a stationary laser beam. We use these to perform coherent sequences of quantum operations, and to perform parallel quantum logic gates on two ions in different processing zones of a multiplexed ion trap chip using a single recycled laser beam. For the latter, we demonstrate individually addressed single-qubit gates by local control of the speed of each ion. The fidelities we observe are consistent with operations performed using standard methods involving static ions and pulsed laser fields. This work therefore provides a path to scalable ion trap quantum computing with reduced requirements on the optical control complexity. PMID:26967401
XTP as a transport protocol for distributed parallel processing
Strayer, W.T.; Lewis, M.J.; Cline, R.E. Jr.
1994-12-31
The Xpress Transfer Protocol (XTP) is a flexible transport layer protocol designed to provide efficient service without dictating the communication paradigm or the delivery characteristics that quality the paradigm. XTP provides the tools to build communication services appropriate to the application. Current data delivery solutions for many popular cluster computing environments use TCP and UDP. We examine TCP, UDP, and XTP with respect to the communication characteristics typical of parallel applications. We perform measurements of end-to-end latency for several paradigms important to cluster computing. An implementation of XTP is shown to be comparable to TCP in end-to-end latency on preestablished connections, and does better for paradigms where connections must be constructed on the fly.
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.
Plane parallel radiance transport for global illumination in vegetation
Max, N.; Mobley, C.; Keating, B.; Wu, E.H.
1997-01-05
This paper applies plane parallel radiance transport techniques to scattering from vegetation. The leaves, stems, and branches are represented as a volume density of scattering surfaces, depending only on height and the vertical component of the surface normal. Ordinary differential equations are written for the multiply scattered radiance as a function of the height above the ground, with the sky radiance and ground reflectance as boundary conditions. They are solved using a two-pass integration scheme to unify the two-point boundary conditions, and Fourier series for the dependence on the azimuthal angle. The resulting radiance distribution is used to precompute diffuse and specular `ambient` shading tables, as a function of height and surface normal, to be used in rendering, together with a z-buffer shadow algorithm for direct solar illumination.
Parallel Transport Quantum Logic Gates with Trapped Ions
NASA Astrophysics Data System (ADS)
de Clercq, Ludwig E.; Lo, Hsiang-Yu; Marinelli, Matteo; Nadlinger, David; Oswald, Robin; Negnevitsky, Vlad; Kienzler, Daniel; Keitch, Ben; Home, Jonathan P.
2016-02-01
We demonstrate single-qubit operations by transporting a beryllium ion with a controlled velocity through a stationary laser beam. We use these to perform coherent sequences of quantum operations, and to perform parallel quantum logic gates on two ions in different processing zones of a multiplexed ion trap chip using a single recycled laser beam. For the latter, we demonstrate individually addressed single-qubit gates by local control of the speed of each ion. The fidelities we observe are consistent with operations performed using standard methods involving static ions and pulsed laser fields. This work therefore provides a path to scalable ion trap quantum computing with reduced requirements on the optical control complexity.
Nonlinear parallel momentum transport in strong electrostatic turbulence
Wang, Lu Wen, Tiliang; Diamond, P. H.
2015-05-15
Most existing theoretical studies of momentum transport focus on calculating the Reynolds stress based on quasilinear theory, without considering the nonlinear momentum flux-〈v{sup ~}{sub r}n{sup ~}u{sup ~}{sub ∥}〉. However, a recent experiment on TORPEX found that the nonlinear toroidal momentum flux induced by blobs makes a significant contribution as compared to the Reynolds stress [Labit et al., Phys. Plasmas 18, 032308 (2011)]. In this work, the nonlinear parallel momentum flux in strong electrostatic turbulence is calculated by using a three dimensional Hasegawa-Mima equation, which is relevant for tokamak edge turbulence. It is shown that the nonlinear diffusivity is smaller than the quasilinear diffusivity from Reynolds stress. However, the leading order nonlinear residual stress can be comparable to the quasilinear residual stress, and so may be important to intrinsic rotation in tokamak edge plasmas. A key difference from the quasilinear residual stress is that parallel fluctuation spectrum asymmetry is not required for nonlinear residual stress.
PFLOTRAN: Recent Developments Facilitating Massively-Parallel Reactive Biogeochemical Transport
NASA Astrophysics Data System (ADS)
Hammond, G. E.
2015-12-01
With the recent shift towards modeling carbon and nitrogen cycling in support of climate-related initiatives, emphasis has been placed on incorporating increasingly mechanistic biogeochemistry within Earth system models to more accurately predict the response of terrestrial processes to natural and anthropogenic climate cycles. PFLOTRAN is an open-source subsurface code that is specialized for simulating multiphase flow and multicomponent biogeochemical transport on supercomputers. The object-oriented code was designed with modularity in mind and has been coupled with several third-party simulators (e.g. CLM to simulate land surface processes and E4D for coupled hydrogeophysical inversion). Central to PFLOTRAN's capabilities is its ability to simulate tightly-coupled reactive transport processes. This presentation focuses on recent enhancements to the code that enable the solution of large parameterized biogeochemical reaction networks with numerous chemical species. PFLOTRAN's "reaction sandbox" is described, which facilitates the implementation of user-defined reaction networks without the need for a comprehensive understanding of PFLOTRAN software infrastructure. The reaction sandbox is written in modern Fortran (2003-2008) and leverages encapsulation, inheritance, and polymorphism to provide the researcher with a flexible workspace for prototyping reactions within a massively parallel flow and transport simulation framework. As these prototypical reactions mature into well-accepted implementations, they can be incorporated into PFLOTRAN as native biogeochemistry capability. Users of the reaction sandbox are encouraged to upload their source code to PFLOTRAN's main source code repository, including the addition of simple regression tests to better ensure the long-term code compatibility and validity of simulation results.
OPTIMIZING COLLAGEN TRANSPORT THROUGH TRACK-ETCHED NANOPORES
Bueno, Ericka M.; Ruberti, Jeffrey W.
2008-01-01
Polymer transport through nanopores is a potentially powerful tool for separation and organization of molecules in biotechnology applications. Our goal is to produce aligned collagen fibrils by mimicking cell-mediated collagen assembly: driving collagen monomers in solution through the aligned nanopores in track-etched membranes followed by fibrillogenesis at the pore exit. We examined type I atelo-collagen monomer transport in neutral, cold solution through polycarbonate track-etched membranes comprising 80-nm-diameter, 6-μm-long pores at 2% areal fraction. Source concentrations of 1.0, 2.8 and 7.0 mg/ml and pressure differentials of 0, 10 and 20 inH2O were used. Membrane surfaces were hydrophilized via covalent poly(ethylene-glycol) binding to limit solute-membrane interaction. Collagen transport through the nanopores was a non-intuitive process due to the complex behavior of this associating molecule in semi-dilute solution. Nonetheless, a modified open pore model provided reasonable predictions of transport parameters. Transport rates were concentration- and pressure-dependent, with diffusivities across the membrane in semi-dilute solution two-fold those in dilute solution, possibly via cooperative diffusion or polymer entrainment. The most significant enhancement of collagen transport was accomplished by membrane hydrophilization. The highest concentration transported (5.99±2.58 mg/ml) with the highest monomer flux (2.60±0.49 ×103 molecules s-1 pore-1) was observed using 2.8 mg collagen/ml, 10 inH2O and hydrophilic membranes. PMID:21394216
T839 fiber tracking transporter at New Muon Lab
Krider, J.
1991-10-23
A darkbox and its transporter have been designed for T839 fiber tracking tests. The darkbox is 3.35 m {times} 0.76 m {times} 0.25 m (1{center dot}w{center dot}h) and contains a scintillating fiber ribbon suspension system and mechanical hardware to support the readout electronics. The transporter provides 3.0 m of horizontal motion transverse to the beam for linear scans of fiber characteristics. In addition, 70 degrees of rotation about a vertical axis is provided to simulate tracking of particles emanating from a collision point at lab angles in the range 0{degrees}--70{degrees}. The transporter, which is located inside a radiation area, is remotely controlled to permit scanning the fiber array through the region defined by four small stationary triggering scintillators without disabling beam. The transporter rails extend 20 feet to the west beyond a gate in the radiation enclosure fencing. This provides a staging area to work on the apparatus, while the beam is on. 4 figs.
T839 fiber tracking transporter at New Muon Lab
Krider, J.
1991-10-23
A darkbox and its transporter have been designed for T839 fiber tracking tests. The darkbox is 3.35 m {times} 0.76 m {times} 0.25 m (1{center_dot}w{center_dot}h) and contains a scintillating fiber ribbon suspension system and mechanical hardware to support the readout electronics. The transporter provides 3.0 m of horizontal motion transverse to the beam for linear scans of fiber characteristics. In addition, 70 degrees of rotation about a vertical axis is provided to simulate tracking of particles emanating from a collision point at lab angles in the range 0{degrees}--70{degrees}. The transporter, which is located inside a radiation area, is remotely controlled to permit scanning the fiber array through the region defined by four small stationary triggering scintillators without disabling beam. The transporter rails extend 20 feet to the west beyond a gate in the radiation enclosure fencing. This provides a staging area to work on the apparatus, while the beam is on. 4 figs.
D`Azevedo, E.F.; Gwo, Jin-Ping
1997-02-01
The objectives of this research are (1) to parallelize a suite of multiregion groundwater flow and solute transport codes that use Galerkin and Lagrangian- Eulerian finite element methods, (2) to test the compatibility of a global shared memory emulation software with a High Performance FORTRAN (HPF) compiler, and (3) to obtain performance characteristics and scalability of the parallel codes. The suite of multiregion flow and transport codes, 3DMURF and 3DMURT, were parallelized using the DOLIB shared memory emulation, in conjunction with the PGI HPF compiler, to run on the Intel Paragons at the Oak Ridge National Laboratory (ORNL) and a network of workstations. The novelty of this effort is first in the use of HPF and global shared memory emulation concurrently to facilitate the conversion of a serial code to a parallel code, and secondly the shared memory library enables efficient implementation of Lagrangian particle tracking along flow characteristics. The latter allows long-time-step-size simulation with particle tracking and dynamic particle redistribution for load balancing, thereby reducing the number of time steps needed for most transient problems. The parallel codes were applied to a pumping well problem to test the efficiency of the domain decomposition and particle tracking algorithms. The full problem domain consists of over 200,000 degrees of freedom with highly nonlinear soil property functions. Relatively good scalability was obtained for a preliminary test run on the Intel Paragons at the Center for Computational Sciences (CCS), ORNL. However, due to the difficulties we encountered in the PGI HPF compiler, as of the writing of this manuscript we are able to report results from 3DMURF only.
Trajectory tracking control of parallel robots in the presence of joint drive flexibility
NASA Astrophysics Data System (ADS)
Ider, S. Kemal; Korkmaz, Ozan
2009-01-01
Trajectory tracking control of parallel manipulators is aimed in the presence of flexibility at the joint drives. Joint structural damping is also considered in the dynamic model. The system is first converted into an open-tree structure by disconnecting a sufficient number of unactuated joints. The closed loops are then expressed by constraint equations. It is shown that, in a parallel robot with flexible joint drives, the acceleration level inverse dynamics equations are singular because the control torques do not have an instantaneous effect on the end-effector accelerations due to the elastic media. Eliminating the Lagrange multipliers and the intermediate variables, a fourth-order input-output relation is obtained between the actuator torques and the end-effector position variables. The proposed control law decouples and linearizes the system and achieves asymptotic stability by feedback of positions and velocities of the actuated joints and rotors. As a case study, a three degree of freedom, two legged planar parallel manipulator is simulated to illustrate the performance of the method. The end-effector desired trajectory is chosen such that the kinematic and drive singular positions are avoided.
Active catheter tracking using parallel MRI and real-time image reconstruction.
Bock, Michael; Müller, Sven; Zuehlsdorff, Sven; Speier, Peter; Fink, Christian; Hallscheidt, Peter; Umathum, Reiner; Semmler, Wolfhard
2006-06-01
In this work active MR catheter tracking with automatic slice alignment was combined with an autocalibrated parallel imaging technique. Using an optimized generalized autocalibrating partially parallel acquisitions (GRAPPA) algorithm with an acceleration factor of 2, we were able to reduce the acquisition time per image by 34%. To accelerate real-time GRAPPA image reconstruction, the coil sensitivities were updated only after slice reorientation. For a 2D trueFISP acquisition (160 x 256 matrix, 80% phase matrix, half Fourier acquisition, TR = 3.7 ms, GRAPPA factor = 2) real-time image reconstruction was achieved with up to six imaging coils. In a single animal experiment the method was used to steer a catheter from the vena cava through the beating heart into the pulmonary vasculature at an image update rate of about five images per second. Under all slice orientations, parallel image reconstruction was accomplished with only minor image artifacts, and the increased temporal resolution provided a sharp delineation of intracardial structures, such as the papillary muscle. PMID:16683261
The optimization of improved mean shift object tracking in embedded multicore DSP parallel system
NASA Astrophysics Data System (ADS)
Tian, Li; Zhou, Fugen; Meng, Cai; Hu, Congliang
2014-11-01
This paper proposes a more robust and efficient Mean Shift object tracking algorithm which is optimized for embedded multicore DSP Parallel system. Firstly, the RGB image is transformed into HSV image which is robust in many aspects such as lighting changes. Then, the color histogram model is used in the back projection process to generate the color probability distribution. Secondly, the size and position of search window are initialized in the first frame, and Mean Shift algorithm calculates the center position of the target and adjusts the search window automatically both in size and location, according to the result of the previous frame. Finally, since the multicore DSP system is commonly adopted in the embedded application such as seeker and an optical scout system, we implement the proposed algorithm in the TI multicore DSP system to meet the need of large amount computation. For multicore parallel computing, the explicit IPC based multicore framework is designed which outperforms OpenMP standard. Moreover, the parallelisms of 8 functional units and cross path data fetch capability of C66 core are utilized to accelerate the computation of iteration in Mean Shift algorithm. The experimental results show that the algorithm has good performance in complex scenes such as deformation, scale change and occlusion, simultaneously the proposed optimization method can significantly reduce the computation time.
Geometry of matrix product states: Metric, parallel transport, and curvature
Haegeman, Jutho Verstraete, Frank; Mariën, Michaël; Osborne, Tobias J.
2014-02-15
We study the geometric properties of the manifold of states described as (uniform) matrix product states. Due to the parameter redundancy in the matrix product state representation, matrix product states have the mathematical structure of a (principal) fiber bundle. The total space or bundle space corresponds to the parameter space, i.e., the space of tensors associated to every physical site. The base manifold is embedded in Hilbert space and can be given the structure of a Kähler manifold by inducing the Hilbert space metric. Our main interest is in the states living in the tangent space to the base manifold, which have recently been shown to be interesting in relation to time dependence and elementary excitations. By lifting these tangent vectors to the (tangent space) of the bundle space using a well-chosen prescription (a principal bundle connection), we can define and efficiently compute an inverse metric, and introduce differential geometric concepts such as parallel transport (related to the Levi-Civita connection) and the Riemann curvature tensor.
Kinetic Modelling of Parallel Electron Transport in TdeV.
NASA Astrophysics Data System (ADS)
Shoucri, M.; Shkarofsky, I.; Stansfield, B.; Batishchev, O.; Batishcheva, A.; Krasheninnikov, S.; Sigmar
1996-11-01
The 1D2V Fokker-Planck code ALLA [1] is used to model parallel electron transport in the SOL of TdeV. Our model uses detached and attached experimental data [2] given by Langmuir probes, Li ablation and He spectroscopy. We obtain the electron distribution function on a precise 257x65x60 non-uniform grid. Strong deviations of hydrogen and carbon excitation rates, and heat conduction coefficient from their Maxwellian values are shown. We compare the calculated variation of the effective temperature at the reciprocating probe position with experimental measurements. We also explain by non-local effects why different experimental techniques show differences in the electron temperature. altaffiltext [1] A.A.Batishcheva et al., Physics of Plasmas 3 (1996) 1634 altaffiltext [2] B.L.Stansfield et al., Proc. 22 Eur.Conf., Bornemouth, 19C pIII-101. *Supported by Government of Canada, Hydro-Québec and INRS. Work performed under USDoE contracts DE-FG02-91-ER-54109 at MIT and DE-FG02-88-ER-53263 at Lodestar.
Particle Tracking Model and Abstraction of Transport Processes
B. Robinson
2004-10-21
The purpose of this report is to document the abstraction model being used in total system performance assessment (TSPA) model calculations for radionuclide transport in the unsaturated zone (UZ). The UZ transport abstraction model uses the particle-tracking method that is incorporated into the finite element heat and mass model (FEHM) computer code (Zyvoloski et al. 1997 [DIRS 100615]) to simulate radionuclide transport in the UZ. This report outlines the assumptions, design, and testing of a model for calculating radionuclide transport in the UZ at Yucca Mountain. In addition, methods for determining and inputting transport parameters are outlined for use in the TSPA for license application (LA) analyses. Process-level transport model calculations are documented in another report for the UZ (BSC 2004 [DIRS 164500]). Three-dimensional, dual-permeability flow fields generated to characterize UZ flow (documented by BSC 2004 [DIRS 169861]; DTN: LB03023DSSCP9I.001 [DIRS 163044]) are converted to make them compatible with the FEHM code for use in this abstraction model. This report establishes the numerical method and demonstrates the use of the model that is intended to represent UZ transport in the TSPA-LA. Capability of the UZ barrier for retarding the transport is demonstrated in this report, and by the underlying process model (BSC 2004 [DIRS 164500]). The technical scope, content, and management of this report are described in the planning document ''Technical Work Plan for: Unsaturated Zone Transport Model Report Integration'' (BSC 2004 [DIRS 171282]). Deviations from the technical work plan (TWP) are noted within the text of this report, as appropriate. The latest version of this document is being prepared principally to correct parameter values found to be in error due to transcription errors, changes in source data that were not captured in the report, calculation errors, and errors in interpretation of source data.
NASA Astrophysics Data System (ADS)
Ney, Michael; Safrani, Avner; Abdulhlaim, Ibrahim
2016-09-01
High resolution fast focus tracking and vibrometery system based on parallel phase shift polarization interferometry using three detectors is presented. The basic design and algorithm are described, followed by an experimental demonstration showing sub nm resolution of different controlled motion profiles instantaneously monitored at a feedback rate of 100 kHz. The fact that the method does not rely on active optical components, potentially allows extremely high vibration rates to be measured; limited only by the detector bandwidth and sampling rate. In addition, the relatively simple design relies only on standard optical equipment, combined with the simple algorithm, makes the task of setting up a high performance vibrometry system cheap and readily available.
NASA Astrophysics Data System (ADS)
Zhang, Dong Ping; Howes, Lee
2013-02-01
We present a parallel multi-hypothesis template tracking algorithm on heterogeneous platforms using a layered dispatch programming model. The contributions of this work are: an architecture-specific optimised solution for vasculature structure enhancement, an approach to segment the vascular lumen network from volumetric CTA images and a layered dispatch programming model to free the developers from hand-crafting mappings to particularly constrained execution domains on high throughput architecture. This abstraction is demonstrated through a vasculature segmentation application and can also be applied in other real-world applications. Current GPGPU programming models define a grouping concept which may lead to poorly scoped lo cal/ shared memory regions and an inconvenient approach to projecting complicated iterations spaces. To improve on this situation, we propose a simpler and more flexible programming model that leads to easier computation projections and hence a more convenient mapping of the same algorithm to a wide range of architectures. We first present an optimised image enhancement solution step- by-step, then solve a separable nonlinear least squares problem using a parallel Levenberg-Marquardt algorithm for template matching, and perform the energy efficiency analysis and performance comparison on a variety of platforms, including multi-core CPUs, discrete GPUs and APUs. We propose and discuss the efficiency of a layered-dispatch programming abstraction for mapping algorithms onto heterogeneous architectures.
Su, Hao; Dickstein-Fischer, Laurie; Harrington, Kevin; Fu, Qiushi; Lu, Weina; Huang, Haibo; Cole, Gregory; Fischer, Gregory S
2010-01-01
This paper presents the development of new prismatic actuation approach and its application in human-safe humanoid head design. To reduce actuator output impedance and mitigate unexpected external shock, the prismatic actuation method uses cables to drive a piston with preloaded spring. By leveraging the advantages of parallel manipulator and cable-driven mechanism, the developed neck has a parallel manipulator embodiment with two cable-driven limbs embedded with preloaded springs and one passive limb. The eye mechanism is adapted for low-cost webcam with succinct "ball-in-socket" structure. Based on human head anatomy and biomimetics, the neck has 3 degree of freedom (DOF) motion: pan, tilt and one decoupled roll while each eye has independent pan and synchronous tilt motion (3 DOF eyes). A Kalman filter based face tracking algorithm is implemented to interact with the human. This neck and eye structure is translatable to other human-safe humanoid robots. The robot's appearance reflects a non-threatening image of a penguin, which can be translated into a possible therapeutic intervention for children with Autism Spectrum Disorders. PMID:21095653
Tracking lithium transport and electrochemical reactions in nanoparticles.
Wang, Feng; Yu, Hui-Chia; Chen, Min-Hua; Wu, Lijun; Pereira, Nathalie; Thornton, Katsuyo; Van der Ven, Anton; Zhu, Yimei; Amatucci, Glenn G; Graetz, Jason
2012-01-01
Expectations for the next generation of lithium batteries include greater energy and power densities along with a substantial increase in both calendar and cycle life. Developing new materials to meet these goals requires a better understanding of how electrodes function by tracking physical and chemical changes of active components in a working electrode. Here we develop a new, simple in-situ electrochemical cell for the transmission electron microscope and use it to track lithium transport and conversion in FeF(2) nanoparticles by nanoscale imaging, diffraction and spectroscopy. In this system, lithium conversion is initiated at the surface, sweeping rapidly across the FeF(2) particles, followed by a gradual phase transformation in the bulk, resulting in 1-3 nm iron crystallites mixed with amorphous LiF. The real-time imaging reveals a surprisingly fast conversion process in individual particles (complete in a few minutes), with a morphological evolution resembling spinodal decomposition. This work provides new insights into the inter- and intra-particle lithium transport and kinetics of lithium conversion reactions, and may help to pave the way to develop high-energy conversion electrodes for lithium-ion batteries. PMID:23149745
Improvements in fast-response flood modeling: desktop parallel computing and domain tracking
Judi, David R; Mcpherson, Timothy N; Burian, Steven J
2009-01-01
It is becoming increasingly important to have the ability to accurately forecast flooding, as flooding accounts for the most losses due to natural disasters in the world and the United States. Flood inundation modeling has been dominated by one-dimensional approaches. These models are computationally efficient and are considered by many engineers to produce reasonably accurate water surface profiles. However, because the profiles estimated in these models must be superimposed on digital elevation data to create a two-dimensional map, the result may be sensitive to the ability of the elevation data to capture relevant features (e.g. dikes/levees, roads, walls, etc...). Moreover, one-dimensional models do not explicitly represent the complex flow processes present in floodplains and urban environments and because two-dimensional models based on the shallow water equations have significantly greater ability to determine flow velocity and direction, the National Research Council (NRC) has recommended that two-dimensional models be used over one-dimensional models for flood inundation studies. This paper has shown that two-dimensional flood modeling computational time can be greatly reduced through the use of Java multithreading on multi-core computers which effectively provides a means for parallel computing on a desktop computer. In addition, this paper has shown that when desktop parallel computing is coupled with a domain tracking algorithm, significant computation time can be eliminated when computations are completed only on inundated cells. The drastic reduction in computational time shown here enhances the ability of two-dimensional flood inundation models to be used as a near-real time flood forecasting tool, engineering, design tool, or planning tool. Perhaps even of greater significance, the reduction in computation time makes the incorporation of risk and uncertainty/ensemble forecasting more feasible for flood inundation modeling (NRC 2000; Sayers et al
Tracking Vehicle in GSM Network to Support Intelligent Transportation Systems
NASA Astrophysics Data System (ADS)
Koppanyi, Z.; Lovas, T.; Barsi, A.; Demeter, H.; Beeharee, A.; Berenyi, A.
2012-07-01
The penetration of GSM capable devices is very high, especially in Europe. To exploit the potential of turning these mobile devices into dynamic data acquisition nodes that provides valuable data for Intelligent Transportation Systems (ITS), position information is needed. The paper describes the basic operation principles of the GSM system and provides an overview on the existing methods for deriving location data in the network. A novel positioning solution is presented that rely on handover (HO) zone measurements; the zone geometry properties are also discussed. A new concept of HO zone sequence recognition is introduced that involves application of Probabilistic Deterministic Finite State Automata (PDFA). Both the potential commercial applications and the use of the derived position data in ITS is discussed for tracking vehicles and monitoring traffic flow. As a practical cutting edge example, the integration possibility of the technology in the SafeTRIP platform (developed in an EC FP7 project) is presented.
Effects of parallel electron dynamics on plasma blob transport
Angus, Justin R.; Krasheninnikov, Sergei I.; Umansky, Maxim V.
2012-08-15
The 3D effects on sheath connected plasma blobs that result from parallel electron dynamics are studied by allowing for the variation of blob density and potential along the magnetic field line and using collisional Ohm's law to model the parallel current density. The parallel current density from linear sheath theory, typically used in the 2D model, is implemented as parallel boundary conditions. This model includes electrostatic 3D effects, such as resistive drift waves and blob spinning, while retaining all of the fundamental 2D physics of sheath connected plasma blobs. If the growth time of unstable drift waves is comparable to the 2D advection time scale of the blob, then the blob's density gradient will be depleted resulting in a much more diffusive blob with little radial motion. Furthermore, blob profiles that are initially varying along the field line drive the potential to a Boltzmann relation that spins the blob and thereby acts as an addition sink of the 2D potential. Basic dimensionless parameters are presented to estimate the relative importance of these two 3D effects. The deviation of blob dynamics from that predicted by 2D theory in the appropriate limits of these parameters is demonstrated by a direct comparison of 2D and 3D seeded blob simulations.
A portable, parallel, object-oriented Monte Carlo neutron transport code in C++
Lee, S.R.; Cummings, J.C.; Nolen, S.D. |
1997-05-01
We have developed a multi-group Monte Carlo neutron transport code using C++ and the Parallel Object-Oriented Methods and Applications (POOMA) class library. This transport code, called MC++, currently computes k and {alpha}-eigenvalues and is portable to and runs parallel on a wide variety of platforms, including MPPs, clustered SMPs, and individual workstations. It contains appropriate classes and abstractions for particle transport and, through the use of POOMA, for portable parallelism. Current capabilities of MC++ are discussed, along with physics and performance results on a variety of hardware, including all Accelerated Strategic Computing Initiative (ASCI) hardware. Current parallel performance indicates the ability to compute {alpha}-eigenvalues in seconds to minutes rather than hours to days. Future plans and the implementation of a general transport physics framework are also discussed.
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)
Modeling reactive transport with particle tracking and kernel estimators
NASA Astrophysics Data System (ADS)
Rahbaralam, Maryam; Fernandez-Garcia, Daniel; Sanchez-Vila, Xavier
2015-04-01
Groundwater reactive transport models are useful to assess and quantify the fate and transport of contaminants in subsurface media and are an essential tool for the analysis of coupled physical, chemical, and biological processes in Earth Systems. Particle Tracking Method (PTM) provides a computationally efficient and adaptable approach to solve the solute transport partial differential equation. On a molecular level, chemical reactions are the result of collisions, combinations, and/or decay of different species. For a well-mixed system, the chem- ical reactions are controlled by the classical thermodynamic rate coefficient. Each of these actions occurs with some probability that is a function of solute concentrations. PTM is based on considering that each particle actually represents a group of molecules. To properly simulate this system, an infinite number of particles is required, which is computationally unfeasible. On the other hand, a finite number of particles lead to a poor-mixed system which is limited by diffusion. Recent works have used this effect to actually model incomplete mix- ing in naturally occurring porous media. In this work, we demonstrate that this effect in most cases should be attributed to a defficient estimation of the concentrations and not to the occurrence of true incomplete mixing processes in porous media. To illustrate this, we show that a Kernel Density Estimation (KDE) of the concentrations can approach the well-mixed solution with a limited number of particles. KDEs provide weighting functions of each particle mass that expands its region of influence, hence providing a wider region for chemical reactions with time. Simulation results show that KDEs are powerful tools to improve state-of-the-art simulations of chemical reactions and indicates that incomplete mixing in diluted systems should be modeled based on alternative conceptual models and not on a limited number of particles.
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.
Solving parallel transport equations in the higher-dimensional Kerr-NUT-(A)dS spacetimes
Connell, Patrick; Frolov, Valeri P.; Kubiznak, David
2008-07-15
We obtain and study the equations describing the parallel transport of orthonormal frames along geodesics in a spacetime admitting a nondegenerate, principal, conformal Killing-Yano tensor h. We demonstrate that the operator F, obtained by a projection of h to a subspace orthogonal to the velocity, has in a generic case eigenspaces of dimension not greater than 2. Each of these eigenspaces is independently parallel propagated. This allows one to reduce the parallel transport equations to a set of first order, ordinary, differential equations for the angles of rotation in the 2D eigenspaces. General analysis is illustrated by studying the equations of the parallel transport in the Kerr-NUT-(A)dS metrics. Examples of three-, four-, and five-dimensional Kerr-NUT-(A)dS are considered, and it is shown that the obtained first order equations can be solved by a separation of variables.
Solving parallel transport equations in the higher-dimensional Kerr-NUT-(A)dS spacetimes
NASA Astrophysics Data System (ADS)
Connell, Patrick; Frolov, Valeri P.; Kubizňák, David
2008-07-01
We obtain and study the equations describing the parallel transport of orthonormal frames along geodesics in a spacetime admitting a nondegenerate, principal, conformal Killing-Yano tensor h. We demonstrate that the operator F, obtained by a projection of h to a subspace orthogonal to the velocity, has in a generic case eigenspaces of dimension not greater than 2. Each of these eigenspaces is independently parallel propagated. This allows one to reduce the parallel transport equations to a set of first order, ordinary, differential equations for the angles of rotation in the 2D eigenspaces. General analysis is illustrated by studying the equations of the parallel transport in the Kerr-NUT-(A)dS metrics. Examples of three-, four-, and five-dimensional Kerr-NUT-(A)dS are considered, and it is shown that the obtained first order equations can be solved by a separation of variables.
3D multi-scale analysis of coupled heat and moisture transport and its parallel implementation
NASA Astrophysics Data System (ADS)
Kruis, Jaroslav
2016-06-01
Parallel implementation of two-scale model of coupled heat and moisture transport is described. The coupled heat and moisture transport is based on the Künzel model. Motivation for the two-scale analysis comes from the requirement to describe distribution of the relative humidity and temperature in historical masonry structures.
A massively parallel track-finding system for the LEVEL 2 trigger in the CLAS detector at CEBAF
Doughty, D.C. Jr.; Collins, P.; Lemon, S. ); Bonneau, P. )
1994-02-01
The track segment finding subsystem of the LEVEL 2 trigger in the CLAS detector has been designed and prototyped. Track segments will be found in the 35,076 wires of the drift chambers using a massively parallel array of 768 Xilinx XC-4005 FPGA's. These FPGA's are located on daughter cards attached to the front-end boards distributed around the detector. Each chip is responsible for finding tracks passing through a 4 x 6 slice of an axial superlayer, and reports two segment found bits, one for each pair of cells. The algorithm used finds segments even when one or two layers or cells along the track is missing (this number is programmable), while being highly resistant to false segments arising from noise hits. Adjacent chips share data to find tracks crossing cell and board boundaries. For maximum speed, fully combinatorial logic is used inside each chip, with the result that all segments in the detector are found within 150 ns. Segment collection boards gather track segments from each axial superlayer and pass them via a high speed link to the segment linking subsystem in an additional 400 ns for typical events. The Xilinx chips are ram-based and therefore reprogrammable, allowing for future upgrades and algorithm enhancements.
Le Crom, Stéphane; Schackwitz, Wendy; Pennacchio, Len; Magnuson, Jon K.; Culley, David E.; Collett, James R.; Martin, Joel; Druzhinina, Irina S.; Mathis, Hugues; Monot, Frédéric; Seiboth, Bernhard; Cherry, Barbara; Rey, Michael; Berka, Randy; Kubicek, Christian P.; Baker, Scott E.; Margeot, Antoine
2009-01-01
Trichoderma reesei (teleomorph Hypocrea jecorina) is the main industrial source of cellulases and hemicellulases harnessed for the hydrolysis of biomass to simple sugars, which can then be converted to biofuels such as ethanol and other chemicals. The highly productive strains in use today were generated by classical mutagenesis. To learn how cellulase production was improved by these techniques, we performed massively parallel sequencing to identify mutations in the genomes of two hyperproducing strains (NG14, and its direct improved descendant, RUT C30). We detected a surprisingly high number of mutagenic events: 223 single nucleotides variants, 15 small deletions or insertions, and 18 larger deletions, leading to the loss of more than 100 kb of genomic DNA. From these events, we report previously undocumented non-synonymous mutations in 43 genes that are mainly involved in nuclear transport, mRNA stability, transcription, secretion/vacuolar targeting, and metabolism. This homogeneity of functional categories suggests that multiple changes are necessary to improve cellulase production and not simply a few clear-cut mutagenic events. Phenotype microarrays show that some of these mutations result in strong changes in the carbon assimilation pattern of the two mutants with respect to the wild-type strain QM6a. Our analysis provides genome-wide insights into the changes induced by classical mutagenesis in a filamentous fungus and suggests areas for the generation of enhanced T. reesei strains for industrial applications such as biofuel production. PMID:19805272
Le Crom, Stéphane; Schackwitz, Wendy; Pennacchio, Len; Magnuson, Jon K; Culley, David E; Collett, James R; Martin, Joel; Druzhinina, Irina S; Mathis, Hugues; Monot, Frédéric; Seiboth, Bernhard; Cherry, Barbara; Rey, Michael; Berka, Randy; Kubicek, Christian P; Baker, Scott E; Margeot, Antoine
2009-09-22
Trichoderma reesei (teleomorph Hypocrea jecorina) is the main industrial source of cellulases and hemicellulases harnessed for the hydrolysis of biomass to simple sugars, which can then be converted to biofuels such as ethanol and other chemicals. The highly productive strains in use today were generated by classical mutagenesis. To learn how cellulase production was improved by these techniques, we performed massively parallel sequencing to identify mutations in the genomes of two hyperproducing strains (NG14, and its direct improved descendant, RUT C30). We detected a surprisingly high number of mutagenic events: 223 single nucleotides variants, 15 small deletions or insertions, and 18 larger deletions, leading to the loss of more than 100 kb of genomic DNA. From these events, we report previously undocumented non-synonymous mutations in 43 genes that are mainly involved in nuclear transport, mRNA stability, transcription, secretion/vacuolar targeting, and metabolism. This homogeneity of functional categories suggests that multiple changes are necessary to improve cellulase production and not simply a few clear-cut mutagenic events. Phenotype microarrays show that some of these mutations result in strong changes in the carbon assimilation pattern of the two mutants with respect to the wild-type strain QM6a. Our analysis provides genome-wide insights into the changes induced by classical mutagenesis in a filamentous fungus and suggests areas for the generation of enhanced T. reesei strains for industrial applications such as biofuel production. PMID:19805272
Le Crom, Stphane; Schackwitz, Wendy; Pennacchiod, Len; Magnuson, Jon K.; Culley, David E.; Collett, James R.; Martin, Joel X.; Druzhinina, Irina S.; Mathis, Hugues; Monot, Frdric; Seiboth, Bernhard; Cherry, Barbara; Rey, Michael; Berka, Randy; Kubicek, Christian P.; Baker, Scott E.; Margeot, Antoine
2009-09-22
Trichoderma reesei (teleomorph Hypocrea jecorina) is the main industrial source of cellulases and hemicellulases harnessed for the hydrolysis of biomass to simple sugars, which can then be converted to biofuels, such as ethanol, and other chemicals. The highly productive strains in use today were generated by classical mutagenesis. To learn how cellulase production was improved by these techniques, we performed massively parallel sequencing to identify mutations in the genomes of two hyperproducing strains (NG14, and its direct improved descendant, RUT C30). We detected a surprisingly high number of mutagenic events: 223 single nucleotides variants, 15 small deletions or insertions and 18 larger deletions leading to the loss of more than 100 kb of genomic DNA. From these events we report previously undocumented non-synonymous mutations in 43 genes that are mainly involved in nuclear transport, mRNA stability, transcription, secretion/vacuolar targeting, and metabolism. This homogeneity of functional categories suggests that multiple changes are necessary to improve cellulase production and not simply a few clear-cut mutagenic events. Phenotype microarrays show that some of these mutations result in strong changes in the carbon assimilation pattern of the two mutants with respect to the wild type strain QM6a. Our analysis provides the first genome-wide insights into the changes induced by classical mutagenesis in a filamentous fungus, and suggests new areas for the generation of enhanced T. reesei strains for industrial applications such as biofuel production.
Supersonic Transport Analysis on the IBM Parallel System SP2
NASA Technical Reports Server (NTRS)
Jameson, Antony; Cliff, Susan; Thomas, Scott; Baker, Timothy; Cheng, Wu-Sun
1995-01-01
Several studies of supersonic transport (SST) configurations have been undertaken by members of the High Speed Aerodynamics branch at NASA Ames. These computational investigations involved the analysis of shapes to study the sonic boom signatures, aerodynamic performance characteristics, as well as studies of nacelle/airframe integration. A variety of different computer codes were employed including both structured and unstructured codes. The AIRPLANE code has been used extensively in these investigations. This computer code solves the Euler equations for inviscid flow by exploiting an explicit finite volume method on a mesh of tetrahedral cells. AIRPLANE is capable of handling complete aircraft configurations including nacelles and diverters. An example of a generic SST configuration is shown and a comparison of computed and experimental force coefficients is presented. Most of the computations in support of the SST investigations have been run on the YMP and C-90 computers currently installed at NASA Ames. Additional information is contained in the original extended abstract.
Parallel processing implementation for the coupled transport of photons and electrons using OpenMP
NASA Astrophysics Data System (ADS)
Doerner, Edgardo
2016-05-01
In this work the use of OpenMP to implement the parallel processing of the Monte Carlo (MC) simulation of the coupled transport for photons and electrons is presented. This implementation was carried out using a modified EGSnrc platform which enables the use of the Microsoft Visual Studio 2013 (VS2013) environment, together with the developing tools available in the Intel Parallel Studio XE 2015 (XE2015). The performance study of this new implementation was carried out in a desktop PC with a multi-core CPU, taking as a reference the performance of the original platform. The results were satisfactory, both in terms of scalability as parallelization efficiency.
On the adequacy of message-passing parallel supercomputers for solving neutron transport problems
Azmy, Y.Y.
1990-01-01
A coarse-grained, static-scheduling parallelization of the standard iterative scheme used for solving the discrete-ordinates approximation of the neutron transport equation is described. The parallel algorithm is based on a decomposition of the angular domain along the discrete ordinates, thus naturally producing a set of completely uncoupled systems of equations in each iteration. Implementation of the parallel code on Intcl's iPSC/2 hypercube, and solutions to test problems are presented as evidence of the high speedup and efficiency of the parallel code. The performance of the parallel code on the iPSC/2 is analyzed, and a model for the CPU time as a function of the problem size (order of angular quadrature) and the number of participating processors is developed and validated against measured CPU times. The performance model is used to speculate on the potential of massively parallel computers for significantly speeding up real-life transport calculations at acceptable efficiencies. We conclude that parallel computers with a few hundred processors are capable of producing large speedups at very high efficiencies in very large three-dimensional problems. 10 refs., 8 figs.
Evaluation of Baltic Sea transport properties using particle tracking
NASA Astrophysics Data System (ADS)
Dargahi, Bijan; Cvetkovic, Vladimir
2014-05-01
Particle tracking model (PTM) is an effective tool for quantifying transport properties of large water bodies such as the Baltic Sea. We have applied PTM to our fully calibrated and validated Baltic Sea 3D hydrodynamic model for a 10-years period (2000-9). One hundred particles were released at a constant rate during an initial 10-days period from all the Baltic Sea sub-basins, the major rivers, and the open boundary in the Arkona Basin. In each basin, the particles were released at two different depths corresponding to the deep water and middle water layers. The objectives of the PTM simulations were to analyse the intra-exchange processes between the Baltic Sea basins and to estimate the arrival times and the paths of particles released from the rivers. The novel contribution of this study is determining the paths and arrival times of deeper water masses rather than the surface masses. Advective and diffusive transport processes in the Bornholm and Arkona basins are both driven by the interacting flows of the northern basins of the Baltic Sea and the North Sea. Particles released from Arkona basin flows northwards along the Stople Channel. The Gotland basins are the major contributors to the exchange process in the Baltic Sea. We find high values of the advection ratio, indicative of a forced advective transport process. The Bay of Gdansk is probably the most vulnerable region in the Baltic Sea. This is despite the fact that the main exchanging basins are the Bornholm Sea and the Easter Gotland Basin. The main reason is the intensive supply of the particles from the northern basins that normally take about 3000 days to reach the Bay of Gdansk. The process maintains a high level of particle concentration (90%) along its coastlines even after the 10-years period. Comparing the particle paths in the Western and Eastern Gotland basins two interesting features were found. Particles travelled in all four directions in the former basin and the middle layer particles
Xu, Linfeng; Lee, Hun; Panchapakesan, Rajagopal; Oh, Kwang W
2012-10-21
We propose a robust droplet fusion and sorting method for two parallel trains of droplets that is relatively insensitive to frequency and phase mismatch. Conventional methods of droplet fusion require an extremely precise control of aqueous/oil flows for perfect frequency matching between two trains of droplets. In this work, by combining our previous two methods (i.e., droplet synchronization using railroad-like channels and manipulation of shape-dependent droplets using guiding tracks), we realized an error-free droplet fusion/sorting device for the two parallel trains of droplets. If droplet pairs are synchronized through a railroad-like channel, they are electrically fused and the fused droplets transit to a middle guiding track to flow in a middle channel; otherwise non-synchronized non-fused droplets will be discarded into the side waste channels by flowing through their own guiding tracks. The simple droplet synchronization, fusion, and sorting technology will have widespread application in droplet-based chemical or biological experiments, where two trains of the chemically or biologically treated or pre-formed droplets yield a train of 100% one-to-one fused droplets at the desired outlet channel by sorting all the non-synchronized non-fused droplets into waste outlets. PMID:22814673
Modeling of tokamak divertor plasma for weakly collisional parallel electron transport
NASA Astrophysics Data System (ADS)
Umansky, M. V.; Dimits, A. M.; Joseph, I.; Omotani, J. T.; Rognlien, T. D.
2015-08-01
The parallel electron heat transport in a weakly collisional regime can be represented in the framework of the Landau-fluid model (Hammett et al., 1990). Practical implementation of Landau-fluid transport has become possible due to the recent invention of an efficient non-spectral method for the non-local closure operators (Dimits et al., 2014). Here the implementation of a Landau-fluid based model for the parallel plasma transport is described, and the model is tested for different collisionality regimes against Fokker-Planck simulations. The new method appears to represent the weakly collisional electron transport more accurately than the conventional flux-limiter based models, on the other hand it is computationally efficient enough to be incorporated in comprehensive edge plasma simulations.
Measurement of parallel transport of B+ and Al2+ impurities in a linear He+ plasma
NASA Astrophysics Data System (ADS)
Hollmann, E. M.; Umstadter, K. R.; Doerner, R. P.; Munoz, J.; Nishijima, D.; Pigarov, A. Yu.
2011-08-01
Laser ablation is used to inject trace quantities (˜1%) of boron or aluminum impurities into steady-state, cylindrical He+ plasmas. The He+ plasmas have an axial flow velocity of about 2 × 103 m/s, corresponding to Mach number M ˜ 0.2. The parallel flow of the resulting B+ or Al2+ impurity ions is measured directly using photomultiplier tubes with line filters. Perpendicular loss of the impurity ions is estimated indirectly from particle conservation. Rapid (<100 μs) entrainment of B+ in the background flow is observed; Al2+ ions appear to accelerate up to perhaps half the background flow velocity in the same time period. Preliminary comparisons with expected transport rates indicate that the observed parallel transport is ≈2-3× faster than expected from classical collisional transport, while perpendicular transport appears slightly slower than expected.
penORNL: a parallel Monte Carlo photon and electron transport package using PENELOPE
Bekar, Kursat B.; Miller, Thomas Martin; Patton, Bruce W.; Weber, Charles F.
2015-01-01
The parallel Monte Carlo photon and electron transport code package penORNL was developed at Oak Ridge National Laboratory to enable advanced scanning electron microscope (SEM) simulations on high-performance computing systems. This paper discusses the implementations, capabilities and parallel performance of the new code package. penORNL uses PENELOPE for its physics calculations and provides all available PENELOPE features to the users, as well as some new features including source definitions specifically developed for SEM simulations, a pulse-height tally capability for detailed simulations of gamma and x-ray detectors, and a modified interaction forcing mechanism to enable accurate energy deposition calculations. The parallel performance of penORNL was extensively tested with several model problems, and very good linear parallel scaling was observed with up to 512 processors. penORNL, along with its new features, will be available for SEM simulations upon completion of the new pulse-height tally implementation.
Lateral charge transport from heavy-ion tracks in integrated circuit chips
NASA Technical Reports Server (NTRS)
Zoutendyk, J. A.; Schwartz, H. R.; Nevill, L. R.
1988-01-01
A 256K DRAM has been used to study the lateral transport of charge (electron-hole pairs) induced by direct ionization from heavy-ion tracks in an IC. The qualitative charge transport has been simulated using a two-dimensional numerical code in cylindrical coordinates. The experimental bit-map data clearly show the manifestation of lateral charge transport in the creation of adjacent multiple-bit errors from a single heavy-ion track. The heavy-ion data further demonstrate the occurrence of multiple-bit errors from single ion tracks with sufficient stopping power. The qualitative numerical simulation results suggest that electric-field-funnel-aided (drift) collection accounts for single error generated by an ion passing through a charge-collecting junction, while multiple errors from a single ion track are due to lateral diffusion of ion-generated charge.
ERIC Educational Resources Information Center
Farmer, Thomas A.; Cargill, Sarah A.; Hindy, Nicholas C.; Dale, Rick; Spivey, Michael J.
2007-01-01
Although several theories of online syntactic processing assume the parallel activation of multiple syntactic representations, evidence supporting simultaneous activation has been inconclusive. Here, the continuous and non-ballistic properties of computer mouse movements are exploited, by recording their streaming x, y coordinates to procure…
Massively parallel simulation of flow and transport in variably saturated porous and fractured media
Wu, Yu-Shu; Zhang, Keni; Pruess, Karsten
2002-01-15
This paper describes a massively parallel simulation method and its application for modeling multiphase flow and multicomponent transport in porous and fractured reservoirs. The parallel-computing method has been implemented into the TOUGH2 code and its numerical performance is tested on a Cray T3E-900 and IBM SP. The efficiency and robustness of the parallel-computing algorithm are demonstrated by completing two simulations with more than one million gridblocks, using site-specific data obtained from a site-characterization study. The first application involves the development of a three-dimensional numerical model for flow in the unsaturated zone of Yucca Mountain, Nevada. The second application is the study of tracer/radionuclide transport through fracture-matrix rocks for the same site. The parallel-computing technique enhances modeling capabilities by achieving several-orders-of-magnitude speedup for large-scale and high resolution modeling studies. The resulting modeling results provide many new insights into flow and transport processes that could not be obtained from simulations using the single-CPU simulator.
Azmy, Yousry
2014-06-10
We employ the Integral Transport Matrix Method (ITMM) as the kernel of new parallel solution methods for the discrete ordinates approximation of the within-group neutron transport equation. The ITMM abandons the repetitive mesh sweeps of the traditional source iterations (SI) scheme in favor of constructing stored operators that account for the direct coupling factors among all the cells' fluxes and between the cells' and boundary surfaces' fluxes. The main goals of this work are to develop the algorithms that construct these operators and employ them in the solution process, determine the most suitable way to parallelize the entire procedure, and evaluate the behavior and parallel performance of the developed methods with increasing number of processes, P. The fastest observed parallel solution method, Parallel Gauss-Seidel (PGS), was used in a weak scaling comparison with the PARTISN transport code, which uses the source iteration (SI) scheme parallelized with the Koch-baker-Alcouffe (KBA) method. Compared to the state-of-the-art SI-KBA with diffusion synthetic acceleration (DSA), this new method- even without acceleration/preconditioning-is completitive for optically thick problems as P is increased to the tens of thousands range. For the most optically thick cells tested, PGS reduced execution time by an approximate factor of three for problems with more than 130 million computational cells on P = 32,768. Moreover, the SI-DSA execution times's trend rises generally more steeply with increasing P than the PGS trend. Furthermore, the PGS method outperforms SI for the periodic heterogeneous layers (PHL) configuration problems. The PGS method outperforms SI and SI-DSA on as few as P = 16 for PHL problems and reduces execution time by a factor of ten or more for all problems considered with more than 2 million computational cells on P = 4.096.
Makedonska, Nataliia; Painter, Scott L.; Bui, Quan M.; Gable, Carl W.; Karra, Satish
2015-09-16
The discrete fracture network (DFN) model is a method to mimic discrete pathways for fluid flow through a fractured low-permeable rock mass, and may be combined with particle tracking simulations to address solute transport. However, experience has shown that it is challenging to obtain accurate transport results in three-dimensional DFNs because of the high computational burden and difficulty in constructing a high-quality unstructured computational mesh on simulated fractures. We present a new particle tracking capability, which is adapted to control volume (Voronoi polygons) flow solutions on unstructured grids (Delaunay triangulations) on three-dimensional DFNs. The locally mass-conserving finite-volume approach eliminates mass balance-related problems during particle tracking. The scalar fluxes calculated for each control volume face by the flow solver are used to reconstruct a Darcy velocity at each control volume centroid. The groundwater velocities can then be continuously interpolated to any point in the domain of interest. The control volumes at fracture intersections are split into four pieces, and the velocity is reconstructed independently on each piece, which results in multiple groundwater velocities at the intersection, one for each fracture on each side of the intersection line. This technique enables detailed particle transport representation through a complex DFN structure. Verified for small DFNs, the new simulation capability enables numerical experiments on advective transport in large DFNs to be performed. As a result, we demonstrate this particle transport approach on a DFN model using parameters similar to those of crystalline rock at a proposed geologic repository for spent nuclear fuel in Forsmark, Sweden.
Makedonska, Nataliia; Painter, Scott L.; Bui, Quan M.; Gable, Carl W.; Karra, Satish
2015-09-16
The discrete fracture network (DFN) model is a method to mimic discrete pathways for fluid flow through a fractured low-permeable rock mass, and may be combined with particle tracking simulations to address solute transport. However, experience has shown that it is challenging to obtain accurate transport results in three-dimensional DFNs because of the high computational burden and difficulty in constructing a high-quality unstructured computational mesh on simulated fractures. We present a new particle tracking capability, which is adapted to control volume (Voronoi polygons) flow solutions on unstructured grids (Delaunay triangulations) on three-dimensional DFNs. The locally mass-conserving finite-volume approach eliminates massmore » balance-related problems during particle tracking. The scalar fluxes calculated for each control volume face by the flow solver are used to reconstruct a Darcy velocity at each control volume centroid. The groundwater velocities can then be continuously interpolated to any point in the domain of interest. The control volumes at fracture intersections are split into four pieces, and the velocity is reconstructed independently on each piece, which results in multiple groundwater velocities at the intersection, one for each fracture on each side of the intersection line. This technique enables detailed particle transport representation through a complex DFN structure. Verified for small DFNs, the new simulation capability enables numerical experiments on advective transport in large DFNs to be performed. As a result, we demonstrate this particle transport approach on a DFN model using parameters similar to those of crystalline rock at a proposed geologic repository for spent nuclear fuel in Forsmark, Sweden.« less
Massively Parallel, Three-Dimensional Transport Solutions for the k-Eigenvalue Problem
Davidson, Gregory G; Evans, Thomas M; Jarrell, Joshua J; Pandya, Tara M; Slaybaugh, R
2014-01-01
We have implemented a new multilevel parallel decomposition in the Denovo dis- crete ordinates radiation transport code. In concert with Krylov subspace iterative solvers, the multilevel decomposition allows concurrency over energy in addition to space-angle, enabling scalability beyond the limits imposed by the traditional KBA space-angle partitioning. Furthermore, a new Arnoldi-based k-eigenvalue solver has been implemented. The added phase-space concurrency combined with the high- performance Krylov and Arnoldi solvers has enabled weak scaling to O(100K) cores on the Jaguar XK6 supercomputer. The multilevel decomposition provides sucient parallelism to scale to exascale computing and beyond.
Track recognition in 4 [mu]s by a systolic trigger processor using a parallel Hough transform
Klefenz, F.; Noffz, K.H.; Conen, W.; Zoz, R.; Kugel, A. . Lehrstuhl fuer Informatik V); Maenner, R. . Lehrstuhl fuer Informatik V Univ. Heidelberg . Interdisziplinaeres Zentrum fuer Wissenschaftliches Rechnen)
1993-08-01
A parallel Hough transform processor has been developed that identifies circular particle tracks in a 2D projection of the OPAL jet chamber. The high-speed requirements imposed by the 8 bunch crossing mode of LEP could be fulfilled by computing the starting angle and the radius of curvature for each well defined track in less than 4 [mu]s. The system consists of a Hough transform processor that determines well defined tracks, and a Euler processor that counts their number by applying the Euler relation to the thresholded result of the Hough transform. A prototype of a systolic processor has been built that handles one sector of the jet chamber. It consists of 35 [times] 32 processing elements that were loaded into 21 programmable gate arrays (XILINX). This processor runs at a clock rate of 40 MHz. It has been tested offline with about 1,000 original OPAL events. No deviations from the off-line simulation have been found. A trigger efficiency of 93% has been obtained. The prototype together with the associated drift time measurement unit has been installed at the OPAL detector at LEP and 100k events have been sampled to evaluate the system under detector conditions.
NASA Astrophysics Data System (ADS)
Zerr, Robert Joseph
2011-12-01
The integral transport matrix method (ITMM) has been used as the kernel of new parallel solution methods for the discrete ordinates approximation of the within-group neutron transport equation. The ITMM abandons the repetitive mesh sweeps of the traditional source iterations (SI) scheme in favor of constructing stored operators that account for the direct coupling factors among all the cells and between the cells and boundary surfaces. The main goals of this work were to develop the algorithms that construct these operators and employ them in the solution process, determine the most suitable way to parallelize the entire procedure, and evaluate the behavior and performance of the developed methods for increasing number of processes. This project compares the effectiveness of the ITMM with the SI scheme parallelized with the Koch-Baker-Alcouffe (KBA) method. The primary parallel solution method involves a decomposition of the domain into smaller spatial sub-domains, each with their own transport matrices, and coupled together via interface boundary angular fluxes. Each sub-domain has its own set of ITMM operators and represents an independent transport problem. Multiple iterative parallel solution methods have investigated, including parallel block Jacobi (PBJ), parallel red/black Gauss-Seidel (PGS), and parallel GMRES (PGMRES). The fastest observed parallel solution method, PGS, was used in a weak scaling comparison with the PARTISN code. Compared to the state-of-the-art SI-KBA with diffusion synthetic acceleration (DSA), this new method without acceleration/preconditioning is not competitive for any problem parameters considered. The best comparisons occur for problems that are difficult for SI DSA, namely highly scattering and optically thick. SI DSA execution time curves are generally steeper than the PGS ones. However, until further testing is performed it cannot be concluded that SI DSA does not outperform the ITMM with PGS even on several thousand or tens of
A Novel Implementation of Massively Parallel Three Dimensional Monte Carlo Radiation Transport
NASA Astrophysics Data System (ADS)
Robinson, P. B.; Peterson, J. D. L.
2005-12-01
The goal of our summer project was to implement the difference formulation for radiation transport into Cosmos++, a multidimensional, massively parallel, magneto hydrodynamics code for astrophysical applications (Peter Anninos - AX). The difference formulation is a new method for Symbolic Implicit Monte Carlo thermal transport (Brooks and Szöke - PAT). Formerly, simultaneous implementation of fully implicit Monte Carlo radiation transport in multiple dimensions on multiple processors had not been convincingly demonstrated. We found that a combination of the difference formulation and the inherent structure of Cosmos++ makes such an implementation both accurate and straightforward. We developed a "nearly nearest neighbor physics" technique to allow each processor to work independently, even with a fully implicit code. This technique coupled with the increased accuracy of an implicit Monte Carlo solution and the efficiency of parallel computing systems allows us to demonstrate the possibility of massively parallel thermal transport. This work was performed under the auspices of the U.S. Department of Energy by University of California Lawrence Livermore National Laboratory under contract No. W-7405-Eng-48
Rykovanov, S. G.; Chen, M.; Geddes, C. G. R.; Schroeder, C. B.; Esarey, E.; Leemans, W. P.
2012-12-21
The Virtual Detector for Synchrotron Radiation (VDSR) is a parallel C++ code developed to calculate the incoherent radiation from a single charged particle or a beam moving in given external electro-magnetic fields. In this proceedings the code structure and features are introduced. An example of radiation generation from the betatron motion of a beam in the focusing fields of the wake in a laser-plasma accelerator is presented.
NASA Astrophysics Data System (ADS)
Coulette, David; Hirstoaga, Sever A.; Manfredi, Giovanni
2016-08-01
We develop a hybrid model to describe the parallel transport in a tokamak scrape-off layer following an edge-localized mode (ELM) event. The parallel dynamics is treated with a kinetic Vlasov–Poisson model, while the evolution of the perpendicular temperature {{T}\\bot} is governed by a fluid equation. The coupling is ensured by isotropising collisions. The model generalises an earlier approach where {{T}\\bot} was constant in space and time (Manfredi et al 2011 Plasma Phys. Control. Fusion 53 015012). Numerical results show that the main effect comes from electron–electron collisions, which limit the decrease of the parallel electron temperature and increase the potential drop in the Debye sheath in front of the surface. Ion–ion collisions have an almost negligible impact. The net effect is an increased peak power load on the target plates.
Johnson, William P.; Zhang, Pengfei; Fuller, Mark E.; Scheibe, Timothy D. ); Mailloux, Brian J.; Onstott, Tullis C.; Deflaun, Mary F.; Hubbard, Susan; Radtke, Jon; Kovacic, William P.; Holben, William
2001-01-01
The first results from an innovative bacterial tracking technique, ferrographic capture, applied to bacterial transport in groundwater are reported in this paper. Ferrographic capture was used to analyze samples during an October 1999 bacterial injection experiment at the Narrow Channel Focus Area of the South Oyster Site, VA. Data obtained using this method showed that the timing of bacterial breakthrough was controlled by physical (hydraulic conductivity) heterogeneity in the vertical dimension, as opposed to variation in sediment surface or aqueous chemical properties. Ferrographic tracking yielded results that compared well with results from other tracking techniques over a concentration range of eight orders of magnitude, and provided a low detection limit relative to most other bacterial tracking techniques. The low detection limit of this method allowed observation of transport of an adhesion-deficient bacterium over distances greater than 20 m in the fine sand aquifer under lying this site.
Pandya, Tara M.; Johnson, Seth R.; Evans, Thomas M.; Davidson, Gregory G.; Hamilton, Steven P.; Godfrey, Andrew T.
2015-12-21
This paper discusses the implementation, capabilities, and validation of Shift, a massively parallel Monte Carlo radiation transport package developed and maintained at Oak Ridge National Laboratory. It has been developed to scale well from laptop to small computing clusters to advanced supercomputers. Special features of Shift include hybrid capabilities for variance reduction such as CADIS and FW-CADIS, and advanced parallel decomposition and tally methods optimized for scalability on supercomputing architectures. Shift has been validated and verified against various reactor physics benchmarks and compares well to other state-of-the-art Monte Carlo radiation transport codes such as MCNP5, CE KENO-VI, and OpenMC. Somemore » specific benchmarks used for verification and validation include the CASL VERA criticality test suite and several Westinghouse AP1000® problems. These benchmark and scaling studies show promising results.« less
Pandya, Tara M.; Johnson, Seth R.; Evans, Thomas M.; Davidson, Gregory G.; Hamilton, Steven P.; Godfrey, Andrew T.
2015-12-21
This paper discusses the implementation, capabilities, and validation of Shift, a massively parallel Monte Carlo radiation transport package developed and maintained at Oak Ridge National Laboratory. It has been developed to scale well from laptop to small computing clusters to advanced supercomputers. Special features of Shift include hybrid capabilities for variance reduction such as CADIS and FW-CADIS, and advanced parallel decomposition and tally methods optimized for scalability on supercomputing architectures. Shift has been validated and verified against various reactor physics benchmarks and compares well to other state-of-the-art Monte Carlo radiation transport codes such as MCNP5, CE KENO-VI, and OpenMC. Some specific benchmarks used for verification and validation include the CASL VERA criticality test suite and several Westinghouse AP1000^{®} problems. These benchmark and scaling studies show promising results.
NASA Astrophysics Data System (ADS)
He, W.; Beyer, C.; Fleckenstein, J. H.; Jang, E.; Kolditz, O.; Naumov, D.; Kalbacher, T.
2015-03-01
This technical paper presents an efficient and performance-oriented method to model reactive mass transport processes in environmental and geotechnical subsurface systems. The open source scientific software packages OpenGeoSys and IPhreeqc have been coupled, to combine their individual strengths and features to simulate thermo-hydro-mechanical-chemical coupled processes in porous and fractured media with simultaneous consideration of aqueous geochemical reactions. Furthermore, a flexible parallelization scheme using MPI (Message Passing Interface) grouping techniques has been implemented, which allows an optimized allocation of computer resources for the node-wise calculation of chemical reactions on the one hand, and the underlying processes such as for groundwater flow or solute transport on the other hand. The coupling interface and parallelization scheme have been tested and verified in terms of precision and performance.
Parallel-burn options for dual-fuel single-stage orbital transports
NASA Technical Reports Server (NTRS)
Martin, J. A.
1978-01-01
A parallel-burn version of a single-stage vehicle for transport from the earth to low-earth orbit using two fuels and rocket propulsion is considered. New engine results were incorporated in vehicle performance and design studies. The results indicate that a hydrogen-cooled gas generator cycle engine provides attractive vehicle performance and that there is little incentive for increasing the chamber pressure beyond 27 MPa.
Goulet, T.; Keszei, E.; Jay-Gerin, J.
1988-03-15
We present a three-dimensional probabilistic model of particle transport in a medium where the particles suffer quasielastic collisions. The model accounts for bulk and surface scattering, as well as partial reflections at the boundaries of the medium. We give analytical and numerical methods for the evaluation of the particle transmission probability in the case of a medium with a plane-parallel geometry. The influence of the various parameters of the model on this probability is also discussed.
Stochastic simulation of charged particle transport on the massively parallel processor
NASA Technical Reports Server (NTRS)
Earl, James A.
1988-01-01
Computations of cosmic-ray transport based upon finite-difference methods are afflicted by instabilities, inaccuracies, and artifacts. To avoid these problems, researchers developed a Monte Carlo formulation which is closely related not only to the finite-difference formulation, but also to the underlying physics of transport phenomena. Implementations of this approach are currently running on the Massively Parallel Processor at Goddard Space Flight Center, whose enormous computing power overcomes the poor statistical accuracy that usually limits the use of stochastic methods. These simulations have progressed to a stage where they provide a useful and realistic picture of solar energetic particle propagation in interplanetary space.
A massively parallel semi-Lagrangian algorithm for solving the transport equation
Manson, Russell; Wang, Dali
2010-01-01
The scalar transport equation underpins many models employed in science, engineering, technology and business. Application areas include, but are not restricted to, pollution transport, weather forecasting, video analysis and encoding (the optical flow equation), options and stock pricing (the Black-Scholes equation) and spatially explicit ecological models. Unfortunately finding numerical solutions to this equation which are fast and accurate is not trivial. Moreover, finding such numerical algorithms that can be implemented on high performance computer architectures efficiently is challenging. In this paper the authors describe a massively parallel algorithm for solving the advection portion of the transport equation. We present an approach here which is different to that used in most transport models and which we have tried and tested for various scenarios. The approach employs an intelligent domain decomposition based on the vector field of the system equations and thus automatically partitions the computational domain into algorithmically autonomous regions. The solution of a classic pure advection transport problem is shown to be conservative, monotonic and highly accurate at large time steps. Additionally we demonstrate that the algorithm is highly efficient for high performance computer architectures and thus offers a route towards massively parallel application.
A parallel computational framework for integrated surface-subsurface flow and transport simulations
NASA Astrophysics Data System (ADS)
Park, Y.; Hwang, H.; Sudicky, E. A.
2010-12-01
HydroGeoSphere is a 3D control-volume finite element hydrologic model describing fully-integrated surface and subsurface water flow and solute and thermal energy transport. Because the model solves tighly-coupled highly-nonlinear partial differential equations, often applied at regional and continental scales (for example, to analyze the impact of climate change on water resources), high performance computing (HPC) is essential. The target parallelization includes the composition of the Jacobian matrix for the iterative linearization method and the sparse-matrix solver, a preconditioned Bi-CGSTAB. The matrix assembly is parallelized by using a coarse-grained scheme in that the local matrix compositions can be performed independently. The preconditioned Bi-CGSTAB algorithm performs a number of LU substitutions, matrix-vector multiplications, and inner products, where the parallelization of the LU substitution is not trivial. The parallelization of the solver is achieved by partitioning the domain into equal-size subdomains, with an efficient reordering scheme. The computational flow of the Bi-CGSTAB solver is also modified to reduce the parallelization overhead and to be suitable for parallel architectures. The parallelized model is tested on several benchmark simulations which include linear and nonlinear flow problems involving various domain sizes and degrees of hydrologic complexities. The performance is evaluated in terms of computational robustness and efficiency, using standard scaling performance measures. The results of simulation profiling indicate that the efficiency becomes higher with an increasing number of nodes/elements in the mesh, for increasingly nonlinear transient simulations, and with domains of irregular geometry. These characteristics are promising for the large-scale analysis water resources problems involved integrated surface/subsurface flow regimes.
Rajaram, Harihar; Brutz, Michael; Klein, Dylan R; Mallikamas, Wasin
2014-09-18
Matrix Diffusion and Adsorption within a rock matrix are important mechanisms for retarding transport of radionuclides in fractured rock. Due to computational limitations and difficulties in characterizing complex subsurface systems, diffusive exchange between a fracture network and surrounding rock matrix is often modeled using simplified conceptual representations. There is significant uncertainty in “effective” parameters used in these models, such as the “effective matrix diffusivity”. Often, these parameters are estimated by fitting sparse breakthrough data, and estimated values fall outside meaningful ranges, because simplified interpretive models do not consider complex three-dimensional flow. There is limited understanding of the relationship between the effective parameters and rock mass characteristics including network structure and matrix properties. There is also evidence for an apparent scale-dependence in “effective matrix diffusion” coefficients. These observations raise questions on whether fracture-matrix interaction parameters estimated from small-scale tracer tests can be used for predicting radionuclide fate and transport at the scale of DOE field sites. High-resolution three-dimensional Discrete-Fracture-Network-Matrix (DFNM) models based on well-defined local scale transport equations can help to address some of these questions. Due to tremendous advances in computational technology over the last 10 years, DFNM modeling in relatively large domains is now feasible. The overarching objective of our research is to use DFNM modeling to improve fundamental understanding of how effective parameters in conceptual models are related to fracture network structure and matrix properties. An advanced three-dimensional DFNM model is being developed, which combines upscaled particle-tracking algorithms for fracture-matrix interaction and a parallel fracture-network flow simulator. The particle-tracking algorithms allow complexity in flow fields
Recent Improvements to the IMPACT-T Parallel Particle TrackingCode
Qiang, J.; Pogorelov, I.V.; Ryne, R.
2006-11-16
The IMPACT-T code is a parallel three-dimensional quasi-static beam dynamics code for modeling high brightness beams in photoinjectors and RF linacs. Developed under the US DOE Scientific Discovery through Advanced Computing (SciDAC) program, it includes several key features including a self-consistent calculation of 3D space-charge forces using a shifted and integrated Green function method, multiple energy bins for beams with large energy spread, and models for treating RF standing wave and traveling wave structures. In this paper, we report on recent improvements to the IMPACT-T code including modeling traveling wave structures, short-range transverse and longitudinal wakefields, and longitudinal coherent synchrotron radiation through bending magnets.
Urban, Michael; Vor der Brüggen, Marc; Tampé, Robert
2016-01-01
Membrane protein transport on the single protein level still evades detailed analysis, if the substrate translocated is non-electrogenic. Considerable efforts have been made in this field, but techniques enabling automated high-throughput transport analysis in combination with solvent-free lipid bilayer techniques required for the analysis of membrane transporters are rare. This class of transporters however is crucial in cell homeostasis and therefore a key target in drug development and methodologies to gain new insights desperately needed. The here presented manuscript describes the establishment and handling of a novel biochip for the analysis of membrane protein mediated transport processes at single transporter resolution. The biochip is composed of microcavities enclosed by nanopores that is highly parallel in its design and can be produced in industrial grade and quantity. Protein-harboring liposomes can directly be applied to the chip surface forming self-assembled pore-spanning lipid bilayers using SSM-techniques (solid supported lipid membranes). Pore-spanning parts of the membrane are freestanding, providing the interface for substrate translocation into or out of the cavity space, which can be followed by multi-spectral fluorescent readout in real-time. The establishment of standard operating procedures (SOPs) allows the straightforward establishment of protein-harboring lipid bilayers on the chip surface of virtually every membrane protein that can be reconstituted functionally. The sole prerequisite is the establishment of a fluorescent read-out system for non-electrogenic transport substrates. High-content screening applications are accomplishable by the use of automated inverted fluorescent microscopes recording multiple chips in parallel. Large data sets can be analyzed using the freely available custom-designed analysis software. Three-color multi spectral fluorescent read-out furthermore allows for unbiased data discrimination into different
Quantum dots for tracking cellular transport of lectin-functionalized nanoparticles
Gao Xiaoling; Wang Tao; Wu Bingxian; Chen Jun; Chen Jiyao; Yue Yang; Dai Ning; Chen Hongzhuan Jiang Xinguo
2008-12-05
Successful drug delivery by functionalized nanocarriers largely depends on their efficient intracellular transport which has not yet been fully understood. We developed a new tracking technique by encapsulating quantum dots into the core of wheat germ agglutinin-conjugated nanoparticles (WGA-NP) to track cellular transport of functionalized nanocarriers. The resulting nanoparticles showed no changes in particle size, zeta potential or biobinding activity, and the loaded probe presented excellent photostability and tracking ability. Taking advantage of these properties, cellular transport profiles of WGA-NP in Caco-2 cells was demonstrated. The cellular uptake begins with binding of WGA to its receptor at the cell surface. The subsequent endocytosis happened in a cytoskeleton-dependent manner and by means of clathrin and caveolae-mediated mechanisms. After endosome creating, transport occurs to both trans-Golgi and lysosome. Our study provides new evidences for quantum dots as a cellular tracking probe of nanocarriers and helps understand intracellular transport profile of lectin-functionalized nanoparticles.
A two-phase thermal model for subsurface transport on massively parallel computers
Martinez, M.J.; Hopkins, P.L.
1997-12-01
Many research activities in subsurface transport require the numerical simulation of multiphase flow in porous media. This capability is critical to research in environmental remediation (e.g. contaminations with dense, non-aqueous-phase liquids), nuclear waste management, reservoir engineering, and to the assessment of the future availability of groundwater in many parts of the world. This paper presents an unstructured grid numerical algorithm for subsurface transport in heterogeneous porous media implemented for use on massively parallel (MP) computers. The mathematical model considers nonisothermal two-phase (liquid/gas) flow, including capillary pressure effects, binary diffusion in the gas phase, conductive, latent, and sensible heat transport. The Galerkin finite element method is used for spatial discretization, and temporal integration is accomplished via a predictor/corrector scheme. Message-passing and domain decomposition techniques are used for implementing a scalable algorithm for distributed memory parallel computers. Illustrative applications are shown to demonstrate capabilities and performance, one of which is modeling hydrothermal transport at the Yucca Mountain site for a radioactive waste facility.
Pumpless Transport of Low Surface Tension Liquids in Surface Tension Confined (STC) Tracks
NASA Astrophysics Data System (ADS)
Megaridis, Constantine; Schutzius, Thomas; Elsharkawy, Mohamed; Tiwari, Manish
2012-11-01
Surfaces with patterned wettability have potential applications in microfluidics, fog capture, pool boiling, etc. With recent fabrication advancements, surfaces with adjacent superhydrophobic and superhydrophilic regions are feasible at a reasonable cost; with properly designed patterns, one can produce microfluidic paths (a.k.a. surface tension confined or STC tracks) where a liquid is confined and transported by surface tension alone. The surface tension of water is relatively high (72 mN/m), as compared with oils (~25 mN/m) and organic solvents (~20 mN/m). This makes the design of STC channels for oils and organic solvents far more difficult. In this study, open STC tracks for pumpless transport of low-surface tension liquids (acetone, ethanol, and hexadecane) on microfluidic chips are fabricated using a large-area, wet-processing technique. Wettable, wax-based, submillimeter-wide tracks are applied by a fountain-pen procedure on superoleophobic, fluoroacrylic carbon nanofiber (CNF) composite coatings. The fabricated anisotropic wetting patterns confine the low-surface tension liquids onto the flow tracks, driving them with meniscus velocities exceeding 3 cm/s. Scaling arguments and Washburn's equation provide estimates of the liquid velocities measured in these tracks, which also act as rails for directional sliding control of mm-sized water droplets. The present facile patterned wettability approach can be extended to deposit micrometer-wide tracks.
Parallel Implementation and Scaling of an Adaptive Mesh Discrete Ordinates Algorithm for Transport
Howell, L H
2004-11-29
Block-structured adaptive mesh refinement (AMR) uses a mesh structure built up out of locally-uniform rectangular grids. In the BoxLib parallel framework used by the Raptor code, each processor operates on one or more of these grids at each refinement level. The decomposition of the mesh into grids and the distribution of these grids among processors may change every few timesteps as a calculation proceeds. Finer grids use smaller timesteps than coarser grids, requiring additional work to keep the system synchronized and ensure conservation between different refinement levels. In a paper for NECDC 2002 I presented preliminary results on implementation of parallel transport sweeps on the AMR mesh, conjugate gradient acceleration, accuracy of the AMR solution, and scalar speedup of the AMR algorithm compared to a uniform fully-refined mesh. This paper continues with a more in-depth examination of the parallel scaling properties of the scheme, both in single-level and multi-level calculations. Both sweeping and setup costs are considered. The algorithm scales with acceptable performance to several hundred processors. Trends suggest, however, that this is the limit for efficient calculations with traditional transport sweeps, and that modifications to the sweep algorithm will be increasingly needed as job sizes in the thousands of processors become common.
Monte Carlo photon transport on vector and parallel supercomputers: Final report
Martin, W.R.; Nowak, P.F.
1986-12-01
The University of Michigan has been investigating the implementation of vectorized and parallelized Monte Carlo algorithms for the analysis of photon transport in an inertially-confined fusion (ICF) plasma. The goal of this work is to develop and test Monte Carlo algorithms for vector/parallel supercomputers such as the Cray X-MP and Cray-2. Previous effort has resulted in the development of a vectorized photon transport code, named VPHOT, and a companion scalar code, named SPHOT, that performs the same analysis and is used for comparative purposes to assess the performance of the vectorized algorithm. A test problem, denoted the ICF test problem, has been created and tested with the VPHOT and SPHOT codes. By comparison with a reference LLNL calculation of the ICF test problem, the VPHOT/SPHOT codes have been verified to predict the correct results. Performance results with VPHOT versus SPHOT and the reference LLNL code have been reported previously and indicate that speedups in the range of 6 to 12 can be achieved with the vectorized algorithm versus the conventional scalar algorithm on the Cray X-MP. This report summarizes the progress made during the last year to continue the investigation of vectorized Monte Carlo (parameter studies, alternative vectorized algorithm, alternative target machines) and to extend the work into the area of parallel processing. 5 refs.
Quantum transport through the system of parallel quantum dots with Majorana bound states
Wang, Ning; Li, Yuxian; Lv, Shuhui
2014-02-28
We study the tunneling transport properties through a system of parallel quantum dots which are coupled to Majorana bound states (MBSs). The conductance and spectral function are computed using the retarded Green's function method based on the equation of motion. The conductance of the system is 2e{sup 2}/h at zero Fermi energy and is robust against the coupling between the MBSs and the quantum dots. The dependence of the Fermi energy on the spectral function is different for the first dot (dot1) than for the second dot (dot2) with fixed dot2-MBSs coupling. The influence of the Majorana bound states on the spectral function was studied for the series and parallel configurations of the system. It was found that when the configuration is in series, the Majorana bound states play an important role, resulting in a spectral function with three peaks. However, the spectral function shows two peaks when the system is in a parallel configuration. The zero Fermi energy spectral function is always 1/2 not only in series but also in the parallel configuration and robust against the coupling between the MBSs and the quantum dots. The phase diagram of the Fermi energy versus the quantum dot energy levels was also investigated.
Modeling nitrogen transport and transformation in aquifers using a particle-tracking approach
NASA Astrophysics Data System (ADS)
Cui, Zhengtao; Welty, Claire; Maxwell, Reed M.
2014-09-01
We have integrated multispecies biodegradation and geochemical reactions into an existing particle-tracking code to simulate reactive transport in three-dimensional variably saturated media, with a focus on nitrification and denitrification processes. This new numerical model includes reactive air-phase transport so that gases such as N2 and CO2 can be tracked. Although nitrogen biodegradation is the primary problem addressed here, the method presented is also applicable to other reactive multispecies transport problems. We verified the model by comparison with (1) analytical solutions for saturated one- and two-dimensional cases; (2) a finite element model for a one-dimensional unsaturated case; and (3) laboratory observations for a one-dimensional saturated case. Good agreement between the new code and the verification problems is demonstrated. The new model can simulate nitrogen transport and transformation in a heterogeneous permeability field where sharp concentration gradients are present. An example application to nitrogen species biodegradation and transport of a plume emanating from a leaking sewer in a heterogeneous, variably saturated aquifer is presented to illustrate this capability. This example is a novel application of coupling unsaturated/saturated zone transport with nitrogen species biodegradation. The code has the computational advantages of particle-tracking algorithms, including local and global mass conservation and minimal numerical dispersion. We also present new methods for improving particle code efficiency by implementing the concept of tracking surplus/deficit particles and particle recycling in order to control the growth of particle numbers. The new model retains the advantages of the particle tracking approach such as allowing relatively low spatial and temporal resolutions to be used, while incorporating the robustness of grid-based Monod kinetics to simulate biogeochemical reactions.
Parallel transport of long mean-free-path plasma along open magnetic field lines: Parallel heat flux
Guo Zehua; Tang Xianzhu
2012-06-15
In a long mean-free-path plasma where temperature anisotropy can be sustained, the parallel heat flux has two components with one associated with the parallel thermal energy and the other the perpendicular thermal energy. Due to the large deviation of the distribution function from local Maxwellian in an open field line plasma with low collisionality, the conventional perturbative calculation of the parallel heat flux closure in its local or non-local form is no longer applicable. Here, a non-perturbative calculation is presented for a collisionless plasma in a two-dimensional flux expander bounded by absorbing walls. Specifically, closures of previously unfamiliar form are obtained for ions and electrons, which relate two distinct components of the species parallel heat flux to the lower order fluid moments such as density, parallel flow, parallel and perpendicular temperatures, and the field quantities such as the magnetic field strength and the electrostatic potential. The plasma source and boundary condition at the absorbing wall enter explicitly in the closure calculation. Although the closure calculation does not take into account wave-particle interactions, the results based on passing orbits from steady-state collisionless drift-kinetic equation show remarkable agreement with fully kinetic-Maxwell simulations. As an example of the physical implications of the theory, the parallel heat flux closures are found to predict a surprising observation in the kinetic-Maxwell simulation of the 2D magnetic flux expander problem, where the parallel heat flux of the parallel thermal energy flows from low to high parallel temperature region.
4. Launch closure, concrete apron and tracks at left, transporter/erector ...
4. Launch closure, concrete apron and tracks at left, transporter/erector horizontal at right, view towards northwest - Ellsworth Air Force Base, Delta Flight, Launch Facility, On County Road T512, south of Exit 116 off I-90, Interior, Jackson County, SD
Making Tracks 1.0: Action Researching an Active Transportation Education Program
ERIC Educational Resources Information Center
Robinson, Daniel; Foran, Andrew; Robinson, Ingrid
2014-01-01
This paper reports on the results of the first cycle of an action research project. The objective of this action research was to examine the implementation of a school-based active transportation education program (Making Tracks). A two-cycle action research design was employed in which elementary school students' (ages 7-9), middle school…
NASA Astrophysics Data System (ADS)
Pandya, Tara M.; Johnson, Seth R.; Evans, Thomas M.; Davidson, Gregory G.; Hamilton, Steven P.; Godfrey, Andrew T.
2016-03-01
This work discusses the implementation, capabilities, and validation of Shift, a massively parallel Monte Carlo radiation transport package authored at Oak Ridge National Laboratory. Shift has been developed to scale well from laptops to small computing clusters to advanced supercomputers and includes features such as support for multiple geometry and physics engines, hybrid capabilities for variance reduction methods such as the Consistent Adjoint-Driven Importance Sampling methodology, advanced parallel decompositions, and tally methods optimized for scalability on supercomputing architectures. The scaling studies presented in this paper demonstrate good weak and strong scaling behavior for the implemented algorithms. Shift has also been validated and verified against various reactor physics benchmarks, including the Consortium for Advanced Simulation of Light Water Reactors' Virtual Environment for Reactor Analysis criticality test suite and several Westinghouse AP1000® problems presented in this paper. These benchmark results compare well to those from other contemporary Monte Carlo codes such as MCNP5 and KENO.
Parallel-Batch Scheduling and Transportation Coordination with Waiting Time Constraint
Gong, Hua; Chen, Daheng; Xu, Ke
2014-01-01
This paper addresses a parallel-batch scheduling problem that incorporates transportation of raw materials or semifinished products before processing with waiting time constraint. The orders located at the different suppliers are transported by some vehicles to a manufacturing facility for further processing. One vehicle can load only one order in one shipment. Each order arriving at the facility must be processed in the limited waiting time. The orders are processed in batches on a parallel-batch machine, where a batch contains several orders and the processing time of the batch is the largest processing time of the orders in it. The goal is to find a schedule to minimize the sum of the total flow time and the production cost. We prove that the general problem is NP-hard in the strong sense. We also demonstrate that the problem with equal processing times on the machine is NP-hard. Furthermore, a dynamic programming algorithm in pseudopolynomial time is provided to prove its ordinarily NP-hardness. An optimal algorithm in polynomial time is presented to solve a special case with equal processing times and equal transportation times for each order. PMID:24883385
Parallel transport gates in a mixed-species ion trap processor
NASA Astrophysics Data System (ADS)
Home, Jonathan
Scaled up quantum information processors will require large numbers of parallel gate operations. For ion trap quantum processing, a promising approach is to perform these operations in separated regions of a multi-zone processing chip between which quantum information is transported either by distributed photonic entanglement or by deterministic shuttling of the ions through the array. However scaling the technology for controlling pulsed laser beams which address each of multiple regions appears challenging. I will describe recent work on the control of both beryllium and calcium ions by transporting ions through static laser beams. We have demonstrated both parallel individually addressed operations as well as sequences of operations. Work is in progress towards multi-qubit gates, which requires good control of the ion transport velocity. We have developed a number of techniques for measuring and optimizing velocities in our trap, enabling significant improvements in performance. In addition to direct results, I will give an overview of our multi-species apparatus, including recent results on high fidelity multi-qubit gates. We are grateful for funding from the Swiss National Science Foundation and the ETH Zurich.
NASA Astrophysics Data System (ADS)
Zerr, Robert Joseph
2011-12-01
The integral transport matrix method (ITMM) has been used as the kernel of new parallel solution methods for the discrete ordinates approximation of the within-group neutron transport equation. The ITMM abandons the repetitive mesh sweeps of the traditional source iterations (SI) scheme in favor of constructing stored operators that account for the direct coupling factors among all the cells and between the cells and boundary surfaces. The main goals of this work were to develop the algorithms that construct these operators and employ them in the solution process, determine the most suitable way to parallelize the entire procedure, and evaluate the behavior and performance of the developed methods for increasing number of processes. This project compares the effectiveness of the ITMM with the SI scheme parallelized with the Koch-Baker-Alcouffe (KBA) method. The primary parallel solution method involves a decomposition of the domain into smaller spatial sub-domains, each with their own transport matrices, and coupled together via interface boundary angular fluxes. Each sub-domain has its own set of ITMM operators and represents an independent transport problem. Multiple iterative parallel solution methods have investigated, including parallel block Jacobi (PBJ), parallel red/black Gauss-Seidel (PGS), and parallel GMRES (PGMRES). The fastest observed parallel solution method, PGS, was used in a weak scaling comparison with the PARTISN code. Compared to the state-of-the-art SI-KBA with diffusion synthetic acceleration (DSA), this new method without acceleration/preconditioning is not competitive for any problem parameters considered. The best comparisons occur for problems that are difficult for SI DSA, namely highly scattering and optically thick. SI DSA execution time curves are generally steeper than the PGS ones. However, until further testing is performed it cannot be concluded that SI DSA does not outperform the ITMM with PGS even on several thousand or tens of
A parallel finite element simulator for ion transport through three-dimensional ion channel systems.
Tu, Bin; Chen, Minxin; Xie, Yan; Zhang, Linbo; Eisenberg, Bob; Lu, Benzhuo
2013-09-15
A parallel finite element simulator, ichannel, is developed for ion transport through three-dimensional ion channel systems that consist of protein and membrane. The coordinates of heavy atoms of the protein are taken from the Protein Data Bank and the membrane is represented as a slab. The simulator contains two components: a parallel adaptive finite element solver for a set of Poisson-Nernst-Planck (PNP) equations that describe the electrodiffusion process of ion transport, and a mesh generation tool chain for ion channel systems, which is an essential component for the finite element computations. The finite element method has advantages in modeling irregular geometries and complex boundary conditions. We have built a tool chain to get the surface and volume mesh for ion channel systems, which consists of a set of mesh generation tools. The adaptive finite element solver in our simulator is implemented using the parallel adaptive finite element package Parallel Hierarchical Grid (PHG) developed by one of the authors, which provides the capability of doing large scale parallel computations with high parallel efficiency and the flexibility of choosing high order elements to achieve high order accuracy. The simulator is applied to a real transmembrane protein, the gramicidin A (gA) channel protein, to calculate the electrostatic potential, ion concentrations and I - V curve, with which both primitive and transformed PNP equations are studied and their numerical performances are compared. To further validate the method, we also apply the simulator to two other ion channel systems, the voltage dependent anion channel (VDAC) and α-Hemolysin (α-HL). The simulation results agree well with Brownian dynamics (BD) simulation results and experimental results. Moreover, because ionic finite size effects can be included in PNP model now, we also perform simulations using a size-modified PNP (SMPNP) model on VDAC and α-HL. It is shown that the size effects in SMPNP can
Hybrid Parallel Programming Models for AMR Neutron Monte-Carlo Transport
NASA Astrophysics Data System (ADS)
Dureau, David; Poëtte, Gaël
2014-06-01
This paper deals with High Performance Computing (HPC) applied to neutron transport theory on complex geometries, thanks to both an Adaptive Mesh Refinement (AMR) algorithm and a Monte-Carlo (MC) solver. Several Parallelism models are presented and analyzed in this context, among them shared memory and distributed memory ones such as Domain Replication and Domain Decomposition, together with Hybrid strategies. The study is illustrated by weak and strong scalability tests on complex benchmarks on several thousands of cores thanks to the petaflopic supercomputer Tera100.
Alfven wave transport effects in the time evolution of parallel cosmic-ray-modified shocks
NASA Technical Reports Server (NTRS)
Jones, T. W.
1993-01-01
This paper presents a numerical study of the time evolution of plane, cosmic-ray modified shocks with magnetic field parallel to the shock normal, based on the diffusive shock acceleration formalism and including the effects from the finite propagation speed and energy of Alfven waves responsible for controlling the transport of the cosmic rays. The simulations discussed are based on a three-fluid model for the dynamics, but a more complete formalism is laid out for future work. The results of the simulations confirm earlier steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and the gas sound speed are comparable, i.e., when the plasma and magnetic pressures are similar. It is also clear, however, that the impact of Alfven transport effects, which tend to slow shock evolution and reduce the time asymptotic cosmic-ray pressure in the shock, is strongly dependent upon uncertain details in the transport models. Both cosmic-ray advection tied to streaming Alfven waves and dissipation of wave energy are important to include in the models. Further, Alfven transport properties on both sides of the shock are also influential.
Suk, Heejun
2012-01-01
Abstract In articles published in 2009 and 2010, Suk and Yeh reported the development of an accurate and efficient particle tracking algorithm for simulating a path line under complicated unsteady flow conditions, using a range of elements within finite elements in multidimensions. Here two examples, an aquifer storage and recovery (ASR) example and a landfill leachate migration example, are examined to enhance the practical implementation of the proposed particle tracking method, known as Suk's method, to a real field of groundwater flow and transport. Results obtained by Suk's method are compared with those obtained by Pollock's method. Suk's method produces superior tracking accuracy, which suggests that Suk's method can describe more accurately various advection-dominated transport problems in a real field than existing popular particle tracking methods, such as Pollock's method. To illustrate the wide and practical applicability of Suk's method to random-walk particle tracking (RWPT), the original RWPT has been modified to incorporate Suk's method. Performance of the modified RWPT using Suk's method is compared with the original RWPT scheme by examining the concentration distributions obtained by the modified RWPT and the original RWPT under complicated transient flow systems. PMID:22476629
Rosa, M.; Warsa, J. S.; Chang, J. H.
2006-07-01
A Fourier analysis is conducted for the discrete-ordinates (SN) approximation of the neutron transport problem solved with Richardson iteration (Source Iteration) and Richardson iteration preconditioned with Transport Synthetic Acceleration (TSA), using the Parallel Block-Jacobi (PBJ) algorithm. Both 'traditional' TSA (TTSA) and a 'modified' TSA (MTSA), in which only the scattering in the low order equations is reduced by some non-negative factor {beta} and < 1, are considered. The results for the un-accelerated algorithm show that convergence of the PBJ algorithm can degrade. The PBJ algorithm with TTSA can be effective provided the {beta} parameter is properly tuned for a given scattering ratio c, but is potentially unstable. Compared to TTSA, MTSA is less sensitive to the choice of {beta}, more effective for the same computational effort (c'), and it is unconditionally stable. (authors)
Investigation of current transport normal and parallel to the tape plane in BSCCO/Ag tapes
Maley, M.P.; Cho, J.H.; Willis, J.O.; Bulaevskii, L.N.
1995-07-01
We have performed transport, resistivity and critical current measurements on Bi-2223/Ag and Bi-2212/Ag tapes with current directions both parallel and perpendicular to the tape plane in magnetic fields up to 7 T and 50
Global restructuring of the CPM-2 transport algorithm for vector and parallel processing
Vujic, J.L.; Martin, W.R. )
1989-11-01
The CPM-2 code is an assembly transport code based on the collision probability (CP) method. It can in principle be applied to global reactor problems, but its excessive computational demands prevent this application. Therefore, a new transport algorithm for CPM-2 has been developed for vector-parallel architectures, which has resulted in an overall factor of 20 speedup (wall clock) on the IBM 3090-600E. This paper presents the detailed results of this effort as well as a brief description of ongoing effort to remove some of the modeling limitations in CPM-2 that inhibit its use for global applications, such as the use of the pure CP treatment and the assumption of isotropic scattering.
Adaptive finite element simulation of flow and transport applications on parallel computers
NASA Astrophysics Data System (ADS)
Kirk, Benjamin Shelton
The subject of this work is the adaptive finite element simulation of problems arising in flow and transport applications on parallel computers. Of particular interest are new contributions to adaptive mesh refinement (AMR) in this parallel high-performance context, including novel work on data structures, treatment of constraints in a parallel setting, generality and extensibility via object-oriented programming, and the design/implementation of a flexible software framework. This technology and software capability then enables more robust, reliable treatment of multiscale--multiphysics problems and specific studies of fine scale interaction such as those in biological chemotaxis (Chapter 4) and high-speed shock physics for compressible flows (Chapter 5). The work begins by presenting an overview of key concepts and data structures employed in AMR simulations. Of particular interest is how these concepts are applied in the physics-independent software framework which is developed here and is the basis for all the numerical simulations performed in this work. This open-source software framework has been adopted by a number of researchers in the U.S. and abroad for use in a wide range of applications. The dynamic nature of adaptive simulations pose particular issues for efficient implementation on distributed-memory parallel architectures. Communication cost, computational load balance, and memory requirements must all be considered when developing adaptive software for this class of machines. Specific extensions to the adaptive data structures to enable implementation on parallel computers is therefore considered in detail. The libMesh framework for performing adaptive finite element simulations on parallel computers is developed to provide a concrete implementation of the above ideas. This physics-independent framework is applied to two distinct flow and transport applications classes in the subsequent application studies to illustrate the flexibility of the
Parallelization and load balancing of a comprehensive atmospheric chemistry transport model
NASA Astrophysics Data System (ADS)
Elbern, Hendrik
Chemistry transport models are generally claimed to be well suited for massively parallel processing on distributed memory architectures since the arithmetic-to-communication ratio is usually high. However, this observation proves insufficient to account for an efficient parallel performance with increasing complexity of the model. The modeling of the local state of the atmosphere ensues very different branches of the modules' code and greater differences in the computational work load and, consequently, runtime of individual processors occur to a much larger extent during a time step than reported for meteorological models. Variable emissions, changes in actinic fluxes, and all processes associated with cloud modeling are highly variable in time and space and are identified to induce large load imbalances which severely affect the parallel efficiency. This is more so, when the model domain encompasses more heterogeneous meteorological or regional regimes, which impinge dissimilarly on simulations of atmospheric chemistry processes. These conditions hold for the EURAD model applied in this study, which covers the European continental scale as integration domain. Based on a master-worker configuration with a horizontal grid partitioning approach, a method is proposed where the integration domain of the individual processors is locally adjusted to accommodate for load imbalances. This ensures a minimal communication volume and data exchange only with the next neighbors. The interior boundary adjustments of the processors are combined with routine boundary exchange which is required each time step anyway. Two dynamic load balancing schemes were implemented and compared against a conventional equal area partition and a static load balancing scheme. The methods are devised for massively parallel distributed memory computers of both, Single and Multiple Instruction stream Multiple Data stream (SIMD, MIMD) types. A midsummer episode of highly elevated ozone concentrations
Random walk particle tracking simulations of non-Fickian transport in heterogeneous media
Srinivasan, G. Tartakovsky, D.M. Dentz, M. Viswanathan, H.; Berkowitz, B.; Robinson, B.A.
2010-06-01
Derivations of continuum nonlocal models of non-Fickian (anomalous) transport require assumptions that might limit their applicability. We present a particle-based algorithm, which obviates the need for many of these assumptions by allowing stochastic processes that represent spatial and temporal random increments to be correlated in space and time, be stationary or non-stationary, and to have arbitrary distributions. The approach treats a particle trajectory as a subordinated stochastic process that is described by a set of Langevin equations, which represent a continuous time random walk (CTRW). Convolution-based particle tracking (CBPT) is used to increase the computational efficiency and accuracy of these particle-based simulations. The combined CTRW-CBPT approach enables one to convert any particle tracking legacy code into a simulator capable of handling non-Fickian transport.
Procassini, R.J.
1997-12-31
The fine-scale, multi-space resolution that is envisioned for accurate simulations of complex weapons systems in three spatial dimensions implies flop-rate and memory-storage requirements that will only be obtained in the near future through the use of parallel computational techniques. Since the Monte Carlo transport models in these simulations usually stress both of these computational resources, they are prime candidates for parallelization. The MONACO Monte Carlo transport package, which is currently under development at LLNL, will utilize two types of parallelism within the context of a multi-physics design code: decomposition of the spatial domain across processors (spatial parallelism) and distribution of particles in a given spatial subdomain across additional processors (particle parallelism). This implementation of the package will utilize explicit data communication between domains (message passing). Such a parallel implementation of a Monte Carlo transport model will result in non-deterministic communication patterns. The communication of particles between subdomains during a Monte Carlo time step may require a significant level of effort to achieve a high parallel efficiency.
OpenGeoSys-GEMS: Hybrid parallelization of a reactive transport code with MPI and threads
NASA Astrophysics Data System (ADS)
Kosakowski, G.; Kulik, D. A.; Shao, H.
2012-04-01
OpenGeoSys-GEMS is a generic purpose reactive transport code based on the operator splitting approach. The code couples the Finite-Element groundwater flow and multi-species transport modules of the OpenGeoSys (OGS) project (http://www.ufz.de/index.php?en=18345) with the GEM-Selektor research package to model thermodynamic equilibrium of aquatic (geo)chemical systems utilizing the Gibbs Energy Minimization approach (http://gems.web.psi.ch/). The combination of OGS and the GEM-Selektor kernel (GEMS3K) is highly flexible due to the object-oriented modular code structures and the well defined (memory based) data exchange modules. Like other reactive transport codes, the practical applicability of OGS-GEMS is often hampered by the long calculation time and large memory requirements. • For realistic geochemical systems which might include dozens of mineral phases and several (non-ideal) solid solutions the time needed to solve the chemical system with GEMS3K may increase exceptionally. • The codes are coupled in a sequential non-iterative loop. In order to keep the accuracy, the time step size is restricted. In combination with a fine spatial discretization the time step size may become very small which increases calculation times drastically even for small 1D problems. • The current version of OGS is not optimized for memory use and the MPI version of OGS does not distribute data between nodes. Even for moderately small 2D problems the number of MPI processes that fit into memory of up-to-date workstations or HPC hardware is limited. One strategy to overcome the above mentioned restrictions of OGS-GEMS is to parallelize the coupled code. For OGS a parallelized version already exists. It is based on a domain decomposition method implemented with MPI and provides a parallel solver for fluid and mass transport processes. In the coupled code, after solving fluid flow and solute transport, geochemical calculations are done in form of a central loop over all finite
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)
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)
Regional Atmospheric Transport Code for Hanford Emission Tracking, Version 2(RATCHET2)
Ramsdell, James V.; Rishel, Jeremy P.
2006-07-01
This manual describes the atmospheric model and computer code for the Atmospheric Transport Module within SAC. The Atmospheric Transport Module, called RATCHET2, calculates the time-integrated air concentration and surface deposition of airborne contaminants to the soil. The RATCHET2 code is an adaptation of the Regional Atmospheric Transport Code for Hanford Emissions Tracking (RATCHET). The original RATCHET code was developed to perform the atmospheric transport for the Hanford Environmental Dose Reconstruction Project. Fundamentally, the two sets of codes are identical; no capabilities have been deleted from the original version of RATCHET. Most modifications are generally limited to revision of the run-specification file to streamline the simulation process for SAC.
Transport properties of track-etched membranes having variable effective pore-lengths
NASA Astrophysics Data System (ADS)
Nguyen, Quoc Hung; Ali, Mubarak; Nasir, Saima; Ensinger, Wolfgang
2015-12-01
The transport rate of molecules through polymeric membranes is normally limited because of their micrometer-scale thickness which restricts their suitability for more practical application. To study the effect of effective pore length on the transport behavior, polymer membranes containing cylindrical and asymmetric-shaped nanopores were prepared through a two-step ion track-etching technique. Permeation experiments were performed separately to investigate the transport properties (molecular flux and selectivity) of these track-etched membranes. The permeation data shows that the molecular flux across membranes containing asymmetric nanopores is higher compared to those having cylindrical pores. On the other hand, the cylindrical pore membranes exhibit higher selectivity than asymmetric pores for the permeation of charged molecules across the membrane. Current-voltage (I-V) measurements of single-pore membranes further verify that asymmetric pores exhibit lower resistance for the flow of ions and therefore show higher currents than cylindrical pores. Moreover, unmodified and polyethyleneimine (PEI) modified asymmetric-shaped pore membranes were successfully used for the separation of cationic and anionic analyte molecules from their mixture, respectively. In this study, two distinct effects (pore geometry and pore density, i.e. number of pores cm-2), which mainly influence membrane selectivity and molecular transport rates, were thoroughly investigated in order to optimize the membrane performance. In this context, we believe that membranes with high molecular transport rates could readily find their application in molecular separation and controlled drug delivery processes.
Transport properties of track-etched membranes having variable effective pore-lengths.
Nguyen, Quoc Hung; Ali, Mubarak; Nasir, Saima; Ensinger, Wolfgang
2015-12-01
The transport rate of molecules through polymeric membranes is normally limited because of their micrometer-scale thickness which restricts their suitability for more practical application. To study the effect of effective pore length on the transport behavior, polymer membranes containing cylindrical and asymmetric-shaped nanopores were prepared through a two-step ion track-etching technique. Permeation experiments were performed separately to investigate the transport properties (molecular flux and selectivity) of these track-etched membranes. The permeation data shows that the molecular flux across membranes containing asymmetric nanopores is higher compared to those having cylindrical pores. On the other hand, the cylindrical pore membranes exhibit higher selectivity than asymmetric pores for the permeation of charged molecules across the membrane. Current-voltage (I-V) measurements of single-pore membranes further verify that asymmetric pores exhibit lower resistance for the flow of ions and therefore show higher currents than cylindrical pores. Moreover, unmodified and polyethyleneimine (PEI) modified asymmetric-shaped pore membranes were successfully used for the separation of cationic and anionic analyte molecules from their mixture, respectively. In this study, two distinct effects (pore geometry and pore density, i.e. number of pores cm(-2)), which mainly influence membrane selectivity and molecular transport rates, were thoroughly investigated in order to optimize the membrane performance. In this context, we believe that membranes with high molecular transport rates could readily find their application in molecular separation and controlled drug delivery processes. PMID:26553245
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.
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.
Brogi, Bharat Bhushan Ahluwalia, P. K.; Chand, Shyam
2015-06-24
Theoretical study of the Coulomb blockade effect on transport properties (Transmission Probability and I-V characteristics) for varied configuration of coupled quantum dot system has been studied by using Non Equilibrium Green Function(NEGF) formalism and Equation of Motion(EOM) method in the presence of magnetic flux. The self consistent approach and intra-dot Coulomb interaction is being taken into account. As the key parameters of the coupled quantum dot system such as dot-lead coupling, inter-dot tunneling and magnetic flux threading through the system can be tuned, the effect of asymmetry parameter and magnetic flux on this tuning is being explored in Coulomb blockade regime. The presence of the Coulomb blockade due to on-dot Coulomb interaction decreases the width of transmission peak at energy level ε + U and by adjusting the magnetic flux the swapping effect in the Fano peaks in asymmetric and symmetric parallel configuration sustains despite strong Coulomb blockade effect.
Time-dependent 3-D dterministic transport on parallel architectures using Dantsys/MPI
Baker, R.S.; Alcouffe, R.E.
1996-12-31
In addition to the ability to solve the static transport equation, we have also incorporated time dependence into our parallel 3-D S{sub {ital N}} code DANTSYS/MPI. Using a semi-implicit scheme, DANTSYS/MPI is capable of performing time-dependent calculations for both fissioning and pure source driven problems. We have applied this to various types of problems such as nuclear well logging and prompt fission experiments. This paper describes the form of the time- dependent equations implemented, their solution strategies in DANTSYS/MPI including iteration acceleration, and the strategies used for time-step control. Results are presented for a model nuclear well logging calculation.
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.
BioFVM: an efficient, parallelized diffusive transport solver for 3-D biological simulations
Ghaffarizadeh, Ahmadreza; Friedman, Samuel H.; Macklin, Paul
2016-01-01
Motivation: Computational models of multicellular systems require solving systems of PDEs for release, uptake, decay and diffusion of multiple substrates in 3D, particularly when incorporating the impact of drugs, growth substrates and signaling factors on cell receptors and subcellular systems biology. Results: We introduce BioFVM, a diffusive transport solver tailored to biological problems. BioFVM can simulate release and uptake of many substrates by cell and bulk sources, diffusion and decay in large 3D domains. It has been parallelized with OpenMP, allowing efficient simulations on desktop workstations or single supercomputer nodes. The code is stable even for large time steps, with linear computational cost scalings. Solutions are first-order accurate in time and second-order accurate in space. The code can be run by itself or as part of a larger simulator. Availability and implementation: BioFVM is written in C ++ with parallelization in OpenMP. It is maintained and available for download at http://BioFVM.MathCancer.org and http://BioFVM.sf.net under the Apache License (v2.0). Contact: paul.macklin@usc.edu. Supplementary information: Supplementary data are available at Bioinformatics online. PMID:26656933
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)
Particle tracking approach for transport in three-dimensional discrete fracture networks
Makedonska, Nataliia; Painter, Scott L; Bui, Quan M; Gable, Carl; Karra, Satish
2015-01-01
The discrete fracture network (DFN) model is a method to mimic discrete pathways for fluid flow through a fractured low-permeable rock mass, and may be combined with particle tracking simulations to address solute transport. However, experience has shown that it is challenging to obtain accurate transport results in three-dimensional DFNs because of the high computational burden and difficulty in constructing a high-quality unstructured computational mesh on simulated fractures. We present a new particle tracking capability, which is adapted to control volume (Voronoi polygons) flow solutions on unstructured grids (Delaunay triangulations) on three-dimensional DFNs. The locally mass-conserving finite-volume approach eliminates massmore » balance related problems during particle tracking. The scalar fluxes calculated for each control volume face by the flow solver are used to reconstruct a Darcy velocity at each control volume centroid. The groundwater velocities can then be continuously interpolated to any point in the domain of interest. The control volumes at fracture intersections are split into four pieces and the velocity is reconstructed independently on each piece, which results in multiple groundwater velocities at the intersection, one for each fracture on each side of the intersection line. This technique enables detailed particle transport representation through a complex DFN structure. Verified for small DFNs, the new simulation capability enables numerical experiments on advective transport in large DFNs to be performed. We demonstrate this particle transport approach on a DFN model using parameters similar to those of crystalline rock at a proposed geologic repository for spent nuclear fuel in Forsmark, Sweden.« less
Particle tracking approach for transport in three-dimensional discrete fracture networks
Makedonska, Nataliia; Painter, Scott L; Bui, Quan M; Gable, Carl; Karra, Satish
2015-01-01
The discrete fracture network (DFN) model is a method to mimic discrete pathways for fluid flow through a fractured low-permeable rock mass, and may be combined with particle tracking simulations to address solute transport. However, experience has shown that it is challenging to obtain accurate transport results in three-dimensional DFNs because of the high computational burden and difficulty in constructing a high-quality unstructured computational mesh on simulated fractures. We present a new particle tracking capability, which is adapted to control volume (Voronoi polygons) flow solutions on unstructured grids (Delaunay triangulations) on three-dimensional DFNs. The locally mass-conserving finite-volume approach eliminates mass balance related problems during particle tracking. The scalar fluxes calculated for each control volume face by the flow solver are used to reconstruct a Darcy velocity at each control volume centroid. The groundwater velocities can then be continuously interpolated to any point in the domain of interest. The control volumes at fracture intersections are split into four pieces and the velocity is reconstructed independently on each piece, which results in multiple groundwater velocities at the intersection, one for each fracture on each side of the intersection line. This technique enables detailed particle transport representation through a complex DFN structure. Verified for small DFNs, the new simulation capability enables numerical experiments on advective transport in large DFNs to be performed. We demonstrate this particle transport approach on a DFN model using parameters similar to those of crystalline rock at a proposed geologic repository for spent nuclear fuel in Forsmark, Sweden.
Qiang, J.; Leitner, D.; Todd, D.
2005-05-16
The driver linac of the proposed Rare Isotope Accelerator (RIA) requires a great variety of high intensity, high charge state ion beams. In order to design and to optimize the low energy beamline optics of the RIA front end,we have developed a new parallel three-dimensional model to simulate the low energy, multi-species ion beam formation and transport from the ECR ion source extraction region to the focal plane of the analyzing magnet. A multisection overlapped computational domain has been used to break the original transport system into a number of each subsystem, macro-particle tracking is used to obtain the charge density distribution in this subdomain. The three-dimensional Poisson equation is solved within the subdomain and particle tracking is repeated until the solution converges. Two new Poisson solvers based on a combination of the spectral method and the multigrid method have been developed to solve the Poisson equation in cylindrical coordinates for the beam extraction region and in the Frenet-Serret coordinates for the bending magnet region. Some test examples and initial applications will also be presented.
Advanced quadratures and periodic boundary conditions in parallel 3D S{sub n} transport
Manalo, K.; Yi, C.; Huang, M.; Sjoden, G.
2013-07-01
Significant updates in numerical quadratures have warranted investigation with 3D Sn discrete ordinates transport. We show new applications of quadrature departing from level symmetric (S{sub 2}o). investigating 3 recently developed quadratures: Even-Odd (EO), Linear-Discontinuous Finite Element - Surface Area (LDFE-SA), and the non-symmetric Icosahedral Quadrature (IC). We discuss implementation changes to 3D Sn codes (applied to Hybrid MOC-Sn TITAN and 3D parallel PENTRAN) that can be performed to accommodate Icosahedral Quadrature, as this quadrature is not 90-degree rotation invariant. In particular, as demonstrated using PENTRAN, the properties of Icosahedral Quadrature are suitable for trivial application using periodic BCs versus that of reflective BCs. In addition to implementing periodic BCs for 3D Sn PENTRAN, we implemented a technique termed 'angular re-sweep' which properly conditions periodic BCs for outer eigenvalue iterative loop convergence. As demonstrated by two simple transport problems (3-group fixed source and 3-group reflected/periodic eigenvalue pin cell), we remark that all of the quadratures we investigated are generally superior to level symmetric quadrature, with Icosahedral Quadrature performing the most efficiently for problems tested. (authors)
NASA Astrophysics Data System (ADS)
Künzli, Pierre; Tsunematsu, Kae; Albuquerque, Paul; Falcone, Jean-Luc; Chopard, Bastien; Bonadonna, Costanza
2016-04-01
Volcanic ash transport and dispersal models typically describe particle motion via a turbulent velocity field. Particles are advected inside this field from the moment they leave the vent of the volcano until they deposit on the ground. Several techniques exist to simulate particles in an advection field such as finite difference Eulerian, Lagrangian-puff or pure Lagrangian techniques. In this paper, we present a new flexible simulation tool called TETRAS (TEphra TRAnsport Simulator) based on a hybrid Eulerian-Lagrangian model. This scheme offers the advantages of being numerically stable with no numerical diffusion and easily parallelizable. It also allows us to output particle atmospheric concentration or ground mass load at any given time. The model is validated using the advection-diffusion analytical equation. We also obtained a good agreement with field observations of the tephra deposit associated with the 2450 BP Pululagua (Ecuador) and the 1996 Ruapehu (New Zealand) eruptions. As this kind of model can lead to computationally intensive simulations, a parallelization on a distributed memory architecture was developed. A related performance model, taking into account load imbalance, is proposed and its accuracy tested.
Transport-parallel cross folds within a mid-crustal Caledonian thrust stack, northern Scotland
NASA Astrophysics Data System (ADS)
Alsop, G. I.; Holdsworth, R. E.; Strachan, R. A.
1996-06-01
Cross folds are typically associated with zones of anomalous foliation trend that lie at high angles to orogenic strike. This case study concentrates on a region of transport-parallel cross folding developed by buckling during Caledonian ductile thrusting within the Moine and Naver Nappes of northern Scotland. Detailed structural analysis reveals a systematic angular relationship between the trend of tectonic transport, fold axes, and the vergence of minor folds. Large-scale cross folding is considered here to be related to wrench-dominated differential shearing during thrust-sense displacements along an important, possibly out-of-sequence structure, here termed the Ben Blandy Shear Zone. This suggests that patterns of folding within the internal parts of the Caledonian orogen in Scotland are principally controlled by the kinematic constraints imposed by low-angle thrusting. Thus, early, 'main-phase' folding is associated with the initiation and propagation of ductile thrusts, whilst later, secondary structures, including cross folds can be related to the development of flow perturbations during displacement along well-established, regionally important detachments.
Pigarov, A Y; Krasheninnikov, S I; LaBombard, B; Rognlien, T D
2006-06-06
Large-Mach-number parallel plasma flows in the single-null SOL of different tokamaks are simulated with multi-fluid transport code UEDGE. The key role of poloidal asymmetry of cross-field plasma transport as the driving mechanism for such flows is discussed. The impact of ballooning-like diffusive and convective transport and plasma flows on divertor detachment, material migration, impurity flows, and erosion/deposition profiles is studied. The results on well-balanced double null plasma modeling that are indicative of strong asymmetry of cross-field transport are presented.
NASA Astrophysics Data System (ADS)
Masciopinto, Costantino; Volpe, Angela; Palmiotta, Domenico; Cherubini, Claudia
2010-09-01
A combination of a parallel fracture model with the PHREEQC-2 geochemical model was developed to simulate sequential flow and chemical transport with reactions in fractured media where both laminar and turbulent flows occur. The integration of non-laminar flow resistances in one model produced relevant effects on water flow velocities, thus improving model prediction capabilities on contaminant transport. The proposed conceptual model consists of 3D rock-blocks, separated by horizontal bedding plane fractures with variable apertures. Particle tracking solved the transport equations for conservative compounds and provided input for PHREEQC-2. For each cluster of contaminant pathways, PHREEQC-2 determined the concentration for mass-transfer, sorption/desorption, ion exchange, mineral dissolution/precipitation and biodegradation, under kinetically controlled reactive processes of equilibrated chemical species. Field tests have been performed for the code verification. As an example, the combined model has been applied to a contaminated fractured aquifer of southern Italy in order to simulate the phenol transport. The code correctly fitted the field available data and also predicted a possible rapid depletion of phenols as a result of an increased biodegradation rate induced by a simulated artificial injection of nitrates, upgradient to the sources.
Masciopinto, Costantino; Volpe, Angela; Palmiotta, Domenico; Cherubini, Claudia
2010-09-20
A combination of a parallel fracture model with the PHREEQC-2 geochemical model was developed to simulate sequential flow and chemical transport with reactions in fractured media where both laminar and turbulent flows occur. The integration of non-laminar flow resistances in one model produced relevant effects on water flow velocities, thus improving model prediction capabilities on contaminant transport. The proposed conceptual model consists of 3D rock-blocks, separated by horizontal bedding plane fractures with variable apertures. Particle tracking solved the transport equations for conservative compounds and provided input for PHREEQC-2. For each cluster of contaminant pathways, PHREEQC-2 determined the concentration for mass-transfer, sorption/desorption, ion exchange, mineral dissolution/precipitation and biodegradation, under kinetically controlled reactive processes of equilibrated chemical species. Field tests have been performed for the code verification. As an example, the combined model has been applied to a contaminated fractured aquifer of southern Italy in order to simulate the phenol transport. The code correctly fitted the field available data and also predicted a possible rapid depletion of phenols as a result of an increased biodegradation rate induced by a simulated artificial injection of nitrates, upgradient to the sources. PMID:20701994
A Many-Task Parallel Approach for Multiscale Simulations of Subsurface Flow and Reactive Transport
Scheibe, Timothy D.; Yang, Xiaofan; Schuchardt, Karen L.; Agarwal, Khushbu; Chase, Jared M.; Palmer, Bruce J.; Tartakovsky, Alexandre M.
2014-12-16
Continuum-scale models have long been used to study subsurface flow, transport, and reactions but lack the ability to resolve processes that are governed by pore-scale mixing. Recently, pore-scale models, which explicitly resolve individual pores and soil grains, have been developed to more accurately model pore-scale phenomena, particularly reaction processes that are controlled by local mixing. However, pore-scale models are prohibitively expensive for modeling application-scale domains. This motivates the use of a hybrid multiscale approach in which continuum- and pore-scale codes are coupled either hierarchically or concurrently within an overall simulation domain (time and space). This approach is naturally suited to an adaptive, loosely-coupled many-task methodology with three potential levels of concurrency. Each individual code (pore- and continuum-scale) can be implemented in parallel; multiple semi-independent instances of the pore-scale code are required at each time step providing a second level of concurrency; and Monte Carlo simulations of the overall system to represent uncertainty in material property distributions provide a third level of concurrency. We have developed a hybrid multiscale model of a mixing-controlled reaction in a porous medium wherein the reaction occurs only over a limited portion of the domain. Loose, minimally-invasive coupling of pre-existing parallel continuum- and pore-scale codes has been accomplished by an adaptive script-based workflow implemented in the Swift workflow system. We describe here the methods used to create the model system, adaptively control multiple coupled instances of pore- and continuum-scale simulations, and maximize the scalability of the overall system. We present results of numerical experiments conducted on NERSC supercomputing systems; our results demonstrate that loose many-task coupling provides a scalable solution for multiscale subsurface simulations with minimal overhead.
NASA Astrophysics Data System (ADS)
Patchimpattapong, Apisit
This dissertation develops an expert system for generating an effective spatial mesh distribution for the discrete ordinates particle transport method in a parallel environment. This expert system consists of two main parts: (1) an algorithm for generating an effective mesh distribution in a serial environment, and (2) an algorithm for inference of an effective domain decomposition strategy for parallel computing. The mesh generation algorithm consists of four steps: creation of a geometric model as partitioned into coarse meshes, determination of an approximate flux shape, selection of appropriate differencing schemes, and generation of an effective fine mesh distribution. A geometric model was created using AutoCAD. A parallel code PENFC (Parallel Environment Neutral-Particle First Collision) has been developed to calculate an uncollided flux in a 3-D Cartesian geometry. The appropriate differencing schemes were selected based on the uncollided flux distribution using a least squares methodology. A menu-driven serial code PENXMSH has been developed to generate an effective spatial mesh distribution that preserves problem geometry and physics. The domain decomposition selection process involves evaluation of the four factors that affect parallel performance, which include number of processors and memory available per processor, load balance, granularity, and degree-of-coupling among processors. These factors are used to derive a parallel-performance-index that provides expected performance of a parallel algorithm depending on computing environment and resources. A large index indicates a high granularity algorithm with relatively low coupling among processors. This expert system has been successfully tested within the PENTRAN (Parallel Environment Neutral-Particle Transport) code system for simulating real-life shielding problems: the VENUS-3 experimental facility and the BWR core shroud.
Microtubule doublets are double-track railways for intraflagellar transport trains.
Stepanek, Ludek; Pigino, Gaia
2016-05-01
The cilium is a large macromolecular machine that is vital for motility, signaling, and sensing in most eukaryotic cells. Its conserved core structure, the axoneme, contains nine microtubule doublets, each comprising a full A-microtubule and an incomplete B-microtubule. However, thus far, the function of this doublet geometry has not been understood. We developed a time-resolved correlative fluorescence and three-dimensional electron microscopy approach to investigate the dynamics of intraflagellar transport (IFT) trains, which carry ciliary building blocks along microtubules during the assembly and disassembly of the cilium. Using this method, we showed that each microtubule doublet is used as a bidirectional double-track railway: Anterograde IFT trains move along B-microtubules, and retrograde trains move along A-microtubules. Thus, the microtubule doublet geometry provides direction-specific rails to coordinate bidirectional transport of ciliary components. PMID:27151870
Electrical Transport Through Micro Porous Track Etch Membranes of same Porosity
NASA Astrophysics Data System (ADS)
Garg, Ravish; Kumar, Vijay; Kumar, Dinesh; Chakarvarti, S. K.
2012-12-01
Porosity, pore size and thickness of membrane are vital factors to influence the transport phenomena through micro porous track etch membranes (TEMs) and affect the various applications like separations, drug release, flow control, bio-sensing and cell size detection etc. based on transport process. Therefore, a better understanding of transport mechanism through TEMs is required for new applications in various thrust areas like biomedical devices and packaging of foods and drugs. Transport studies of electrolytic solutions of potassium chloride, through porous polycarbonate TEMS having cylindrical pores of size 0.2 μm and 0.4 μm with same porosity of 15%, have been carried out using an electrochemical cell. In this technique, the etched filter is sandwiched between two compartments of cell in such a way that the TEM acts as a membrane separating the cell into two chambers. The two chambers are then filled with electrolyte solution (KCl in distilled water). The current voltage characteristics have been drawn by stepping the voltage ranging 0 to 10 V using Keithley 2400 Series Source Measurement Unit. The results indicate that rate of ion transport through cylindrical pores although is independent of pore size of TEMs of same porosity but there seems to be effect of TEM aperture size exposed to the electrolyte used in conducting cell on ion transport magnitude. From the experimental studies, a large deviation in the conduction through TEMs was observed when compared with theoretical consideration which led to the need for modification in the applicability of simple Ohm's law to the conduction through TEMs. It is found that ion transport increases with increase in area of aperture of TEM but much lower than the expected theoretically value.
NASA Astrophysics Data System (ADS)
Henri, Christopher; Fernàndez-Garcia, Daniel
2015-04-01
Modeling multi-species reactive transport in natural systems with strong heterogeneities and complex biochemical reactions is a major challenge for assessing groundwater polluted sites with organic and inorganic contaminants. A large variety of these contaminants react according to serial-parallel reaction networks commonly simplified by a combination of first-order kinetic reactions. In this context, a random-walk particle tracking method is presented. This method is capable of efficiently simulating the motion of particles affected by first-order network reactions in three-dimensional systems, which are represented by spatially variable physical and biochemical coefficients described at high resolution. The approach is based on the development of transition probabilities that describe the likelihood that particles belonging to a given species and location at a given time will be transformed into and moved to another species and location afterwards. These probabilities are derived from the solution matrix of the spatial moments governing equations. The method is fully coupled with reactions, free of numerical dispersion and overcomes the inherent numerical problems stemming from the incorporation of heterogeneities to reactive transport codes. In doing this, we demonstrate that the motion of particles follows a standard random walk with time-dependent effective retardation and dispersion parameters that depend on the initial and final chemical state of the particle. The behavior of effective parameters develops as a result of differential retardation effects among species. Moreover, explicit analytic solutions of the transition probability matrix and related particle motions are provided for serial reactions. An example of the effect of heterogeneity on the dechlorination of organic solvents in a three-dimensional random porous media shows that the power-law behavior typically observed in conservative tracers breakthrough curves can be largely compromised by the
NASA Astrophysics Data System (ADS)
Henri, Christopher V.; Fernández-Garcia, Daniel
2014-09-01
Modeling multispecies reactive transport in natural systems with strong heterogeneities and complex biochemical reactions is a major challenge for assessing groundwater polluted sites with organic and inorganic contaminants. A large variety of these contaminants react according to serial-parallel reaction networks commonly simplified by a combination of first-order kinetic reactions. In this context, a random-walk particle tracking method is presented. This method is capable of efficiently simulating the motion of particles affected by first-order network reactions in three-dimensional systems, which are represented by spatially variable physical and biochemical coefficients described at high resolution. The approach is based on the development of transition probabilities that describe the likelihood that particles belonging to a given species and location at a given time will be transformed into and moved to another species and location afterward. These probabilities are derived from the solution matrix of the spatial moments governing equations. The method is fully coupled with reactions, free of numerical dispersion and overcomes the inherent numerical problems stemming from the incorporation of heterogeneities to reactive transport codes. In doing this, we demonstrate that the motion of particles follows a standard random walk with time-dependent effective retardation and dispersion parameters that depend on the initial and final chemical state of the particle. The behavior of effective parameters develops as a result of differential retardation effects among species. Moreover, explicit analytic solutions of the transition probability matrix and related particle motions are provided for serial reactions. An example of the effect of heterogeneity on the dechlorination of organic solvents in a three-dimensional random porous media shows that the power-law behavior typically observed in conservative tracers breakthrough curves can be largely compromised by the
Linear Quadratic Tracking Design for a Generic Transport Aircraft with Structural Load Constraints
NASA Technical Reports Server (NTRS)
Burken, John J.; Frost, Susan A.; Taylor, Brian R.
2011-01-01
When designing control laws for systems with constraints added to the tracking performance, control allocation methods can be utilized. Control allocations methods are used when there are more command inputs than controlled variables. Constraints that require allocators are such task as; surface saturation limits, structural load limits, drag reduction constraints or actuator failures. Most transport aircraft have many actuated surfaces compared to the three controlled variables (such as angle of attack, roll rate & angle of side slip). To distribute the control effort among the redundant set of actuators a fixed mixer approach can be utilized or online control allocation techniques. The benefit of an online allocator is that constraints can be considered in the design whereas the fixed mixer cannot. However, an online control allocator mixer has a disadvantage of not guaranteeing a surface schedule, which can then produce ill defined loads on the aircraft. The load uncertainty and complexity has prevented some controller designs from using advanced allocation techniques. This paper considers actuator redundancy management for a class of over actuated systems with real-time structural load limits using linear quadratic tracking applied to the generic transport model. A roll maneuver example of an artificial load limit constraint is shown and compared to the same no load limitation maneuver.
NASA Astrophysics Data System (ADS)
Romano, Paul Kollath
Monte Carlo particle transport methods are being considered as a viable option for high-fidelity simulation of nuclear reactors. While Monte Carlo methods offer several potential advantages over deterministic methods, there are a number of algorithmic shortcomings that would prevent their immediate adoption for full-core analyses. In this thesis, algorithms are proposed both to ameliorate the degradation in parallel efficiency typically observed for large numbers of processors and to offer a means of decomposing large tally data that will be needed for reactor analysis. A nearest-neighbor fission bank algorithm was proposed and subsequently implemented in the OpenMC Monte Carlo code. A theoretical analysis of the communication pattern shows that the expected cost is O( N ) whereas traditional fission bank algorithms are O(N) at best. The algorithm was tested on two supercomputers, the Intrepid Blue Gene/P and the Titan Cray XK7, and demonstrated nearly linear parallel scaling up to 163,840 processor cores on a full-core benchmark problem. An algorithm for reducing network communication arising from tally reduction was analyzed and implemented in OpenMC. The proposed algorithm groups only particle histories on a single processor into batches for tally purposes---in doing so it prevents all network communication for tallies until the very end of the simulation. The algorithm was tested, again on a full-core benchmark, and shown to reduce network communication substantially. A model was developed to predict the impact of load imbalances on the performance of domain decomposed simulations. The analysis demonstrated that load imbalances in domain decomposed simulations arise from two distinct phenomena: non-uniform particle densities and non-uniform spatial leakage. The dominant performance penalty for domain decomposition was shown to come from these physical effects rather than insufficient network bandwidth or high latency. The model predictions were verified with
Kostin, Mikhail; Mokhov, Nikolai; Niita, Koji
2013-09-25
A 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. It is intended to be used with older radiation transport codes implemented in Fortran77, Fortran 90 or C. 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 developed and tested in conjunction with the MARS15 code. It is possible to use it with other codes such as PHITS, FLUKA and MCNP after certain adjustments. 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. The framework 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.
Xie, Dexuan; Dash, Ranjan K.; Beard, Daniel A.
2009-01-01
Fast algorithms for simulating mathematical models of coupled blood-tissue transport and metabolism are critical for the analysis of data on transport and reaction in tissues. Here, by combining the method of characteristics with the standard grid discretization technique, a novel algorithm is introduced for solving a general blood-tissue transport and metabolism model governed by a large system of one-dimensional semilinear first order partial differential equations. The key part of the algorithm is to approximate the model as a group of independent ordinary differential equation (ODE) systems such that each ODE system has the same size as the model and can be integrated independently. Thus the method can be easily implemented in parallel on a large scale multiprocessor computer. The accuracy of the algorithm is demonstrated for solving a simple blood-tissue exchange model introduced by Sangren and Sheppard (Bull. Math. Biophys. 15:387–394, 1953), which has an analytical solution. Numerical experiments made on a distributed-memory parallel computer (an HP Linux cluster) and a shared-memory parallel computer (a SGI Origin 2000) demonstrate the parallel efficiency of the algorithm. PMID:20161089
Xie, Dexuan; Dash, Ranjan K; Beard, Daniel A
2009-11-01
Fast algorithms for simulating mathematical models of coupled blood-tissue transport and metabolism are critical for the analysis of data on transport and reaction in tissues. Here, by combining the method of characteristics with the standard grid discretization technique, a novel algorithm is introduced for solving a general blood-tissue transport and metabolism model governed by a large system of one-dimensional semilinear first order partial differential equations. The key part of the algorithm is to approximate the model as a group of independent ordinary differential equation (ODE) systems such that each ODE system has the same size as the model and can be integrated independently. Thus the method can be easily implemented in parallel on a large scale multiprocessor computer. The accuracy of the algorithm is demonstrated for solving a simple blood-tissue exchange model introduced by Sangren and Sheppard (Bull. Math. Biophys. 15:387-394, 1953), which has an analytical solution. Numerical experiments made on a distributed-memory parallel computer (an HP Linux cluster) and a shared-memory parallel computer (a SGI Origin 2000) demonstrate the parallel efficiency of the algorithm. PMID:20161089
Using CO observations from space to track long-range transport of pollution
NASA Astrophysics Data System (ADS)
Jacob, D. J.; Kim, S.; Wang, H.; Zoogman, P.
2012-12-01
The long and high-quality record of CO observations from satellite sensors provides a unique resource for tracking long-range transport from combustion sources, with many applications for our understanding of air quality and climate forcing. I will discuss some of these applications based on recent and ongoing work at Harvard. I will show how the high density of observations from AIRS and IASI has advanced our knowledge of intercontinental transport of pollution. The new TIR+NIR retrieval from MOPITT provides vertical profile information for CO over continents and I will examine its utility for testing boundary layer ventilation in global models. CO observations from space can be of great value in correlative analyses with other species to isolate the effect of combustion sources and synoptic transport. I will discuss how ozone-CO (OMI+AIRS) and aerosol-CO (MODIS+CALIOP+AIRS) correlations from space can thus provide constraints on ozone sources and on aerosol scavenging. Correlative information from CO observations is also useful in formal data assimilation by effectively reducing the forward model errors on transport and sources. I will discuss the application to CO2 source inversions and to future data assimilation of ozone air quality from geostationary orbit.
Clarno, Kevin T
2007-01-01
The NEWTRNX transport module solves the multigroup, discrete-ordinates source-driven or k-eigenvalue transport equation in parallel on a 3-D unstructured tetrahedral mesh using the extended step characteristics (ESC), also known as the slice-balance approach (SBA), spatial discretization. The spatial domains are decomposed using METIS. NEWTRNX is under development for nuclear reactor analysis on computer hardware ranging from clusters to massively parallel machines, like the Cray XT4. Transport methods that rely on full sweeps across the spatial domain have been shown to display poor scaling for thousands of processors. The Parallel Block-Jacobi (PBJ) algorithm allows each spatial partition to sweep over all discrete-ordinate directions and energies independently of all other domains, potentially allowing for much better scaling than possible with full sweeps. The PBJ algorithm has been implemented in NEWTRNX using a Gauss-Seidel iteration in energy and an asynchronous communication by an energy group, such that each partition utilizes the latest boundary solution available for each group before solving the withingroup scattering in a given group. For each energy group, the within-group scattering converges with a generalized minimum residual (GMRES) solver, preconditioned with beta transport synthetic acceleration ({beta}-TSA).
Le Roux, J. A.
2011-12-10
Earlier work based on nonlinear guiding center (NLGC) theory suggested that perpendicular cosmic-ray transport is diffusive when cosmic rays encounter random three-dimensional magnetohydrodynamic turbulence dominated by uniform two-dimensional (2D) turbulence with a minor uniform slab turbulence component. In this approach large-scale perpendicular cosmic-ray transport is due to cosmic rays microscopically diffusing along the meandering magnetic field dominated by 2D turbulence because of gyroresonant interactions with slab turbulence. However, turbulence in the solar wind is intermittent and it has been suggested that intermittent turbulence might be responsible for the observation of 'dropout' events in solar energetic particle fluxes on small scales. In a previous paper le Roux et al. suggested, using NLGC theory as a basis, that if gyro-scale slab turbulence is intermittent, large-scale perpendicular cosmic-ray transport in weak uniform 2D turbulence will be superdiffusive or subdiffusive depending on the statistical characteristics of the intermittent slab turbulence. In this paper we expand and refine our previous work further by investigating how both parallel and perpendicular transport are affected by intermittent slab turbulence for weak as well as strong uniform 2D turbulence. The main new finding is that both parallel and perpendicular transport are the net effect of an interplay between diffusive and nondiffusive (superdiffusive or subdiffusive) transport effects as a consequence of this intermittency.
NASA Technical Reports Server (NTRS)
Yu, Hongbin; Remer, Lorraine A.; Kahn, Ralph A.; Chin, Mian; Zhang, Yan
2012-01-01
Evidence of aerosol intercontinental transport (ICT) is both widespread and compelling. Model simulations suggest that ICT could significantly affect regional air quality and climate, but the broad inter-model spread of results underscores a need of constraining model simulations with measurements. Satellites have inherent advantages over in situ measurements to characterize aerosol ICT, because of their spatial and temporal coverage. Significant progress in satellite remote sensing of aerosol properties during the Earth Observing System (EOS) era offers opportunity to increase quantitative characterization and estimates of aerosol ICT, beyond the capability of pre-EOS era satellites that could only qualitatively track aerosol plumes. EOS satellites also observe emission strengths and injection heights of some aerosols, aerosol precursors, and aerosol-related gases, which can help characterize aerosol ICT. After an overview of these advances, we review how the current generation of satellite measurements have been used to (1) characterize the evolution of aerosol plumes (e.g., both horizontal and vertical transport, and properties) on an episodic basis, (2) understand the seasonal and inter-annual variations of aerosol ICT and their control factors, (3) estimate the export and import fluxes of aerosols, and (4) evaluate and constrain model simulations. Substantial effort is needed to further explore an integrated approach using measurements from on-orbit satellites (e.g., A-Train synergy) for observational characterization and model constraint of aerosol intercontinental transport and to develop advanced sensors for future missions.
NASA Astrophysics Data System (ADS)
Hansen, Scott K.; Scher, Harvey; Berkowitz, Brian
2014-07-01
Both Eulerian and Lagrangian reactive transport simulations in natural media require selection of a parameter that controls the “promiscuity” of the reacting particles. In Eulerian models, measurement of this parameter may be difficult because its value will generally differ between natural (diffusion-limited) systems and batch experiments, even though both are modeled by reaction terms of the same form. And in Lagrangian models, there previously has been no a priori way to compute this parameter. In both cases, then, selection is typically done by calibration, or ad hoc. This paper addresses the parameter selection problem for Fickian transport by deriving, from first principles and D (the diffusion constant) the reaction-rate-controlling parameters for particle tracking (PT) codes and for the diffusion-reaction equation (DRE). Using continuous time random walk analysis, exact reaction probabilities are derived for pairs of potentially reactive particles based on D and their probability of reaction provided that they collocate. Simultaneously, a second PT scheme directly employing collocation probabilities is derived. One-to-one correspondence between each of D, the reaction radius specified for a PT scheme, and the DRE decay constant are then developed. These results serve to ground reactive transport simulations in their underlying thermodynamics, and are confirmed by simulations.
McGhee, J.M.; Roberts, R.M.; Morel, J.E.
1997-06-01
A spherical harmonics research code (DANTE) has been developed which is compatible with parallel computer architectures. DANTE provides 3-D, multi-material, deterministic, transport capabilities using an arbitrary finite element mesh. The linearized Boltzmann transport equation is solved in a second order self-adjoint form utilizing a Galerkin finite element spatial differencing scheme. The core solver utilizes a preconditioned conjugate gradient algorithm. Other distinguishing features of the code include options for discrete-ordinates and simplified spherical harmonics angular differencing, an exact Marshak boundary treatment for arbitrarily oriented boundary faces, in-line matrix construction techniques to minimize memory consumption, and an effective diffusion based preconditioner for scattering dominated problems. Algorithm efficiency is demonstrated for a massively parallel SIMD architecture (CM-5), and compatibility with MPP multiprocessor platforms or workstation clusters is anticipated.
Modeling bimolecular reactions and transport in porous media via particle tracking
NASA Astrophysics Data System (ADS)
Ding, Dong; Benson, David A.; Paster, Amir; Bolster, Diogo
2013-03-01
We use a particle-tracking method to simulate several one-dimensional bimolecular reactive transport experiments. In our numerical scheme, the reactants are represented by particles: advection and dispersion dominate the flow, and molecular diffusion dictates, in large part, the reactions. The particle/particle reactions are determined by a combination of two probabilities dictated by the physics of transport and energetics of reaction. The first is that reactant particles occupy the same volume over a short time interval. The second is the conditional probability that two collocated particles favorably transform into a reaction. The first probability is a direct physical representation of the degree of mixing in an advancing interface between dissimilar waters, and as such lacks empirical parameters except for the user-defined number of particles. This number can be determined analytically from concentration autocovariance, if this type of data is available. The simulations compare favorably to two physical experiments. In one, the concentration of product, 1,2-naphthoquinoe-4-aminobenzene (NQAB) from reaction between 1,2-naphthoquinone-4-sulfonic acid (NQS) and aniline (AN), was measured at the outflow of a column filled with glass beads at different times. In the other, the concentration distribution of reactants (CuSO and EDTA) and product (CuEDTA) were quantified by snapshots of light transmitted through a column packed with cryolite sand. These snapshots allow us to estimate concentration statistics and calculate the required number of particles. The experiments differ significantly due to a ˜107 difference in thermodynamic rate coefficients, making the latter experiment effectively instantaneous. When compared to the solution of the advection-dispersion-reaction equation (ADRE) with the well-mixed reaction coefficient, the experiments and the particle-tracking simulations showed on the order of 20-40% less overall product, which is attributed to poor mixing
Modeling Bimolecular Reactions and Transport in Porous Media Via Particle Tracking
Dong Ding; David Benson; Amir Paster; Diogo Bolster
2012-01-01
We use a particle-tracking method to simulate several one-dimensional bimolecular reactive transport experiments. In this numerical method, the reactants are represented by particles: advection and dispersion dominate the flow, and molecular diffusion dictates, in large part, the reactions. The particle/particle reactions are determined by a combination of two probabilities dictated by the physics of transport and energetics of reaction. The first is that reactant particles occupy the same volume over a short time interval. The second is the conditional probability that two collocated particles favorably transform into a reaction. The first probability is a direct physical representation of the degree of mixing in an advancing displacement front, and as such lacks empirical parameters except for the user-defined number of particles. This number can be determined analytically from concentration autocovariance, if this type of data is available. The simulations compare favorably to two physical experiments. In one, the concentration of product, 1,2-naphthoquinoe-4-aminobenzene (NQAB) from reaction between 1,2-naphthoquinone-4-sulfonic acid (NQS) and aniline (AN), was measured at the outflow of a column filled with glass beads at different times. In the other, the concentration distribution of reactants (CuSO_4 and EDTA^{4-}) and products (CuEDTA^{4-}) were quantified by snapshots of transmitted light through a column packed with cryloite sand. The thermodynamic rate coefficient in the latter experiment was 10^7 times greater than the former experiment, making it essentially instantaneous. When compared to the solution of the advection-dispersion-reaction equation (ADRE) with the well-mixed reaction coefficient, the experiments and the particle-tracking simulations showed on the order of 20% to 40% less overall product, which is attributed to poor mixing. The poor mixing also leads to higher product concentrations on the edges of the mixing zones, which the particle
Reactor Dosimetry Applications Using RAPTOR-M3G:. a New Parallel 3-D Radiation Transport Code
NASA Astrophysics Data System (ADS)
Longoni, Gianluca; Anderson, Stanwood L.
2009-08-01
The numerical solution of the Linearized Boltzmann Equation (LBE) via the Discrete Ordinates method (SN) requires extensive computational resources for large 3-D neutron and gamma transport applications due to the concurrent discretization of the angular, spatial, and energy domains. This paper will discuss the development RAPTOR-M3G (RApid Parallel Transport Of Radiation - Multiple 3D Geometries), a new 3-D parallel radiation transport code, and its application to the calculation of ex-vessel neutron dosimetry responses in the cavity of a commercial 2-loop Pressurized Water Reactor (PWR). RAPTOR-M3G is based domain decomposition algorithms, where the spatial and angular domains are allocated and processed on multi-processor computer architectures. As compared to traditional single-processor applications, this approach reduces the computational load as well as the memory requirement per processor, yielding an efficient solution methodology for large 3-D problems. Measured neutron dosimetry responses in the reactor cavity air gap will be compared to the RAPTOR-M3G predictions. This paper is organized as follows: Section 1 discusses the RAPTOR-M3G methodology; Section 2 describes the 2-loop PWR model and the numerical results obtained. Section 3 addresses the parallel performance of the code, and Section 4 concludes this paper with final remarks and future work.
Crane, Jonathan M; Haggie, Peter M; Verkman, A S
2009-03-01
Single particle tracking (SPT) provides information about the microscopic motions of individual particles in live cells. We applied SPT to study the diffusion of membrane transport proteins in cell plasma membranes in which individual proteins are labeled with quantum dots at engineered extracellular epitopes. Software was created to deduce particle diffusive modes from quantum dot trajectories. SPT of aquaporin (AQP) water channels and cystic fibrosis transmembrane conductance regulator (CFTR) chloride channels revealed several types of diffusion. AQP1 was freely mobile in cell membranes, showing rapid, Brownian-type diffusion. The full-length (M1) isoform of AQP4 also diffused rapidly, though the diffusion of a shorter (M23) isoform of AQP4 was highly restricted due to its supermolecular assembly in raft-like orthogonal arrays. CFTR mobility was also highly restricted, in a spring-like potential, due to its tethering to the actin cytoskeleton through PDZ-domain C-terminus interactions. The biological significance of regulated diffusion of membrane transport proteins is a subject of active investigation. PMID:24386532
NASA Astrophysics Data System (ADS)
Crane, Jonathan M.; Haggie, Peter M.; Verkman, A. S.
2009-02-01
Single particle tracking (SPT) provides information about the microscopic motions of individual particles in live cells. We applied SPT to study the diffusion of membrane transport proteins in cell plasma membranes in which individual proteins are labeled with quantum dots at engineered extracellular epitopes. Software was created to deduce particle diffusive modes from quantum dot trajectories. SPT of aquaporin (AQP) water channels and cystic fibrosis transmembrane conductance regulator (CFTR) chloride channels revealed several types of diffusion. AQP1 was freely mobile in cell membranes, showing rapid, Brownian-type diffusion. The full-length (M1) isoform of AQP4 also diffused rapidly, though the diffusion of a shorter (M23) isoform of AQP4 was highly restricted due to its supermolecular assembly in raft-like orthogonal arrays. CFTR mobility was also highly restricted, in a spring-like potential, due to its tethering to the actin cytoskeleton through PDZ-domain C-terminus interactions. The biological significance of regulated diffusion of membrane transport proteins is a subject of active investigation.
NASA Astrophysics Data System (ADS)
Cadini, F.; Bertoli, I.; de Sanctis, J.; Zio, E.
2012-10-01
The estimation of the extent and timing of solute migration in a fractured medium is a fundamental task for verifying the level of protection against contaminant releases (e.g., toxic chemicals or radionuclides) offered by the engineered and natural barriers of a waste repository. In this paper we present a novel approach for modeling solute transport in a fractured medium, based on an extension of the Kolmogorov-Dmitriev theory of stochastic branching processes. The model equations for the expected values of the solute concentration take a form similar to that of classical dual-continua models. On the other hand, the stochastic nature of the modeling approach lends itself to a new particle tracking scheme of resolution, which allows accounting for realistic features of the transport process. The proposed stochastic modeling framework and simulation solution approach are illustrated with reference to the experimental results from a case study of literature. Some of the model parameters are optimally identified by means of a genetic algorithm search aimed at best fitting the experimental data.
NASA Technical Reports Server (NTRS)
Schwan, Karsten; Alyea, Fred; Ribarsky, M. William; Trauner, Mary; Eisenhauer, Greg; Jean, Yves; Gu, Weiming; Wang, Ray; Waldrop, Jeffrey; Schroeder, Beth; Heiner, Jeremy; Kindler, Thomas; Silva, Dilma; Vetter, Jeffrey; Martin, Vernard
1996-01-01
The three-dimensional, spectral transport model used in the current project was first successfully integrated over climatological time scales by Dr. Guang Ping Lou for the simulation of atmospheric N2O using the United Kingdom Meteorological Office (UKMO) 4-dimensional, assimilated wind and temperature data set. A non-parallel, FORTRAN version of this integration using a fairly simple N2O chemistry package containing only photo-chemical reactions was used to verify our initial parallel model results. The integrations reproduced the gross features of the observed stratospheric climatological N2O distributions but also simulated the structure of the stratospheric Antarctic vortex and its evolution. Subsequently, Dr. Thomas Kindler, who produced much of the parallel version of our model, enlarged the N2O model chemistry package to include N2O reactions involving O(D-1) and also introduced assimilated wind data from NASA as well as UKMO. Initially, transport calculations without chemistry were run using Carbon-14 as a non-reactive tracer gas with the result that large differences in the transport properties of the two assimilated wind data sets were apparent from the resultant Carbon-14 distributions. Subsequent calculations for N2O, including its chemistry, with the two input winds data sets with verification from UARS satellite observations have refined the transport differences between the two such that the model's steering capabilities could be used to infer the correct climatological vertical velocity fields required to support the N2O observations. During this process, it was also discovered that both the NASA and the UKMO data contained spurious values in some of the higher frequency wave components, leading to incorrect local transport calculations and ultimately affecting the large scale properties of the model's N2O distributions, particularly at tropical latitudes. Subsequent model runs with wind data that had been filtered to remove some of the high
NASA Astrophysics Data System (ADS)
Noury, M.; Bernet, M.; Jaillard, E.; Sempere, T.
2014-12-01
Proterozoic metamorphic rocks largely crop out between ~14 and ~17°S along the coast of southern Peru. Previous thermochronological studies in this belt focused on the Late Neogene cooling history and yielded Late Paleozoic to mid-Cretaceous zircon U-Th/He and fission-track ages, and Late Cretaceous to Early Paleocene apatite fission-track ages. Our study aims at extending the fission-track database in the ~300 km-long, high-grade metamorphic Atico-Mollendo Block (AMB; 15.8-17.1°S), in order to understand the lateral variations of thermochronologic ages along the Pacific coast of southern Peru. We present new apatite and zircon fission-track ages from the AMB and interpret them along with the previously published dates from this block and its surrounding areas. The apatite fission-track data show a longitudinal variation pattern with two localized shifts from younger (60 Ma) to older (~90 Ma) dates, corresponding to the northwestern and southeastern borders of the AMB. These shifts coincide in the field with two major normal fault zones that strike perpendicular to the coast and bound the AMB, and had been understated until now. Given the offsetted ages, the faults were active after 60 Ma, and tilted the AMB down to the NW. This interpretation is consistent with the fact that the metamorphic grade of the basement increases towards the SE, and with the deposition between ~55 and 30 Ma of forearc continental deposits (Lower Moquegua Group) that exhibit numerous synsedimentary extensional features. The zircon fission-track age pattern is more complicated but can be interpreted either as (i) reflecting a partial reset of the whole area due to burial during the Jurassic and Early Cretaceous, or (ii) as the reactivation of an important Carboniferous detachment fault during the ~100-80 Ma interval. We suggest that this coast-parallel Eocene-Oligocene extension accommodated the counterclockwise block rotation of the southern Peruvian forearc that accompanied the
Ramsdell, J.V. Jr.; Simonen, C.A.; Burk, K.W.
1994-02-01
The purpose of the Hanford Environmental Dose Reconstruction (HEDR) Project is to estimate radiation doses that individuals may have received from operations at the Hanford Site since 1944. This report deals specifically with the atmospheric transport model, Regional Atmospheric Transport Code for Hanford Emission Tracking (RATCHET). RATCHET is a major rework of the MESOILT2 model used in the first phase of the HEDR Project; only the bookkeeping framework escaped major changes. Changes to the code include (1) significant changes in the representation of atmospheric processes and (2) incorporation of Monte Carlo methods for representing uncertainty in input data, model parameters, and coefficients. To a large extent, the revisions to the model are based on recommendations of a peer working group that met in March 1991. Technical bases for other portions of the atmospheric transport model are addressed in two other documents. This report has three major sections: a description of the model, a user`s guide, and a programmer`s guide. These sections discuss RATCHET from three different perspectives. The first provides a technical description of the code with emphasis on details such as the representation of the model domain, the data required by the model, and the equations used to make the model calculations. The technical description is followed by a user`s guide to the model with emphasis on running the code. The user`s guide contains information about the model input and output. The third section is a programmer`s guide to the code. It discusses the hardware and software required to run the code. The programmer`s guide also discusses program structure and each of the program elements.
Bedload transport in a formerly glaciated mountain catchment constrained by particle tracking
NASA Astrophysics Data System (ADS)
Dell'Agnese, A.; Brardinoni, F.; Toro, M.; Mao, L.; Engel, M.; Comiti, F.
2015-11-01
In formerly glaciated mountain settings, Pleistocene glaciations are responsible for profound spatial reorganization of the landscape structure. By imposing local channel slope and the degree of hillslope-channel connectivity, glacial macro-forms can exert first-order controls on the downstream strength and continuity of the coarse sediment cascade. To estimate quantitatively these controls we trace bedload transport for 3 years along Strimm Creek, Eastern Italian Alps. Specifically, we monitor the travel distance of 490 PIT-tagged particles (b axis: 23-229 mm; weight: 83-6525 g) at two contrasting sites: Upper Strimm Creek (US; 4 km2), which flows through a fluvially dominated hanging valley, and Lower Strimm Creek (LS; 7.5 km2), located downstream, in a relict glacial trough where it experiences periodic colluvial sediment inputs from lateral tributaries. Tracer positioning within the streambed is periodically tracked in the field with a portable antenna in order to assess progressive travel distances, as well as the extent of the channel active layer, in relation to snowmelt and rainfall-driven peak flows. Interestingly, we show that tracer virtual velocities for selected inter-survey periods are independent of tracer weight at both study sites. Cumulatively, tracers in US have travelled across distances (i.e. inner quartiles) shorter than 2 m, which correspond to over 2 orders of magnitude less than what was observed in LS. These figures translate, after calculations of tracer inter-survey virtual velocities, into estimated bedload volumes equal to about 3 m3 in US and 600 m3 in LS, with most of the transport (75 % in US, and 93 % in LS) occurring during snowmelt. A similar contrast in bedload transport rates, even without considering the additional volumes of material mobilized by mass-wasting processes in LS, testifies the extent to which the glacial imprinting can still affect contemporary sediment transfer, and thus postglacial landscape evolution, in
Gao, Hao; Phan, Lan; Lin, Yuting
2012-09-01
A graphics processing unit-based parallel multigrid solver for a radiative transfer equation with vacuum boundary condition or reflection boundary condition is presented for heterogeneous media with complex geometry based on two-dimensional triangular meshes or three-dimensional tetrahedral meshes. The computational complexity of this parallel solver is linearly proportional to the degrees of freedom in both angular and spatial variables, while the full multigrid method is utilized to minimize the number of iterations. The overall gain of speed is roughly 30 to 300 fold with respect to our prior multigrid solver, which depends on the underlying regime and the parallelization. The numerical validations are presented with the MATLAB codes at https://sites.google.com/site/rtefastsolver/. PMID:23085905
NASA Astrophysics Data System (ADS)
Guo, L.; Huang, H.; Gaston, D.; Redden, G. D.; Fox, D. T.; Fujita, Y.
2010-12-01
Inducing mineral precipitation in the subsurface is one potential strategy for immobilizing trace metal and radionuclide contaminants. Generating mineral precipitates in situ can be achieved by manipulating chemical conditions, typically through injection or in situ generation of reactants. How these reactants transport, mix and react within the medium controls the spatial distribution and composition of the resulting mineral phases. Multiple processes, including fluid flow, dispersive/diffusive transport of reactants, biogeochemical reactions and changes in porosity-permeability, are tightly coupled over a number of scales. Numerical modeling can be used to investigate the nonlinear coupling effects of these processes which are quite challenging to explore experimentally. Many subsurface reactive transport simulators employ a de-coupled or operator-splitting approach where transport equations and batch chemistry reactions are solved sequentially. However, such an approach has limited applicability for biogeochemical systems with fast kinetics and strong coupling between chemical reactions and medium properties. A massively parallel, fully coupled, fully implicit Reactive Transport simulator (referred to as “RAT”) based on a parallel multi-physics object-oriented simulation framework (MOOSE) has been developed at the Idaho National Laboratory. Within this simulator, systems of transport and reaction equations can be solved simultaneously in a fully coupled, fully implicit manner using the Jacobian Free Newton-Krylov (JFNK) method with additional advanced computing capabilities such as (1) physics-based preconditioning for solution convergence acceleration, (2) massively parallel computing and scalability, and (3) adaptive mesh refinements for 2D and 3D structured and unstructured mesh. The simulator was first tested against analytical solutions, then applied to simulating induced calcium carbonate mineral precipitation in 1D columns and 2D flow cells as analogs
Philip, Bobby; Berrill, Mark A.; Allu, Srikanth; Hamilton, Steven P.; Sampath, Rahul S.; Clarno, Kevin T.; Dilts, Gary A.
2015-01-26
We describe an efficient and nonlinearly consistent parallel solution methodology for solving coupled nonlinear thermal transport problems that occur in nuclear reactor applications over hundreds of individual 3D physical subdomains. Efficiency is obtained by leveraging knowledge of the physical domains, the physics on individual domains, and the couplings between them for preconditioning within a Jacobian Free Newton Krylov method. Details of the computational infrastructure that enabled this work, namely the open source Advanced Multi-Physics (AMP) package developed by the authors are described. The details of verification and validation experiments, and parallel performance analysis in weak and strong scaling studies demonstrating the achieved efficiency of the algorithm are presented. Moreover, numerical experiments demonstrate that the preconditioner developed is independent of the number of fuel subdomains in a fuel rod, which is particularly important when simulating different types of fuel rods. Finally, we demonstrate the power of the coupling methodology by considering problems with couplings between surface and volume physics and coupling of nonlinear thermal transport in fuel rods to an external radiation transport code.
Philip, Bobby; Berrill, Mark A.; Allu, Srikanth; Hamilton, Steven P.; Sampath, Rahul S.; Clarno, Kevin T.; Dilts, Gary A.
2015-01-26
We describe an efficient and nonlinearly consistent parallel solution methodology for solving coupled nonlinear thermal transport problems that occur in nuclear reactor applications over hundreds of individual 3D physical subdomains. Efficiency is obtained by leveraging knowledge of the physical domains, the physics on individual domains, and the couplings between them for preconditioning within a Jacobian Free Newton Krylov method. Details of the computational infrastructure that enabled this work, namely the open source Advanced Multi-Physics (AMP) package developed by the authors are described. The details of verification and validation experiments, and parallel performance analysis in weak and strong scaling studies demonstratingmore » the achieved efficiency of the algorithm are presented. Moreover, numerical experiments demonstrate that the preconditioner developed is independent of the number of fuel subdomains in a fuel rod, which is particularly important when simulating different types of fuel rods. Finally, we demonstrate the power of the coupling methodology by considering problems with couplings between surface and volume physics and coupling of nonlinear thermal transport in fuel rods to an external radiation transport code.« less
Chen, Xueli; Gao, Xinbo; Qu, Xiaochao; Chen, Duofang; Ma, Bin; Wang, Lin; Peng, Kuan; Liang, Jimin; Tian, Jie
2010-01-01
During the past decade, Monte Carlo method has obtained wide applications in optical imaging to simulate photon transport process inside tissues. However, this method has not been effectively extended to the simulation of free-space photon transport at present. In this paper, a uniform framework for noncontact optical imaging is proposed based on Monte Carlo method, which consists of the simulation of photon transport both in tissues and in free space. Specifically, the simplification theory of lens system is utilized to model the camera lens equipped in the optical imaging system, and Monte Carlo method is employed to describe the energy transformation from the tissue surface to the CCD camera. Also, the focusing effect of camera lens is considered to establish the relationship of corresponding points between tissue surface and CCD camera. Furthermore, a parallel version of the framework is realized, making the simulation much more convenient and effective. The feasibility of the uniform framework and the effectiveness of the parallel version are demonstrated with a cylindrical phantom based on real experimental results. PMID:20689705
PARALLEL MEASUREMENT AND MODELING OF TRANSPORT IN THE DARHT II BEAMLINE ON ETA II
Chambers, F W; Raymond, B A; Falabella, S; Lee, B S; Richardson, R A; Weir, J T; Davis, H A; Schultze, M E
2005-05-31
To successfully tune the DARHT II transport beamline requires the close coupling of a model of the beam transport and the measurement of the beam observables as the beam conditions and magnet settings are varied. For the ETA II experiment using the DARHT II beamline components this was achieved using the SUICIDE (Simple User Interface Connecting to an Integrated Data Environment) data analysis environment and the FITS (Fully Integrated Transport Simulation) model. The SUICIDE environment has direct access to the experimental beam transport data at acquisition and the FITS predictions of the transport for immediate comparison. The FITS model is coupled into the control system where it can read magnet current settings for real time modeling. We find this integrated coupling is essential for model verification and the successful development of a tuning aid for the efficient convergence on a useable tune. We show the real time comparisons of simulation and experiment and explore the successes and limitations of this close coupled approach.
Transport of nonconductive and conductive droplets in a parallel plate array
NASA Astrophysics Data System (ADS)
Chatterjee, Debalina; Hetayothin, Boonta
2005-03-01
Electrowetting on dielectric technique is used to actuate conductive liquid droplets on electrodes patterned beneath a dielectric. Nonconductive liquids can be transported electrohydrodynamically inside channels. We show for the first time that it is possible to transport droplets of nonconductive liquids on dielectric surfaces, using modest voltages and frequencies (<100 V, <10 kHz). Ionic liquids, aqueous surfactants, buffers, and organic solutions can also be transported. Although conductive liquids show a significant change in liquid contact angle on application of potential, nonconductive liquids do not, suggesting a different mechanism of transport. The empirical criteria for moving droplets in a two-dimensional array are a liquid dielectric constant >= 4.3 and a molecular dipole moment >= 1.2 D. The transport mechanisms are discussed along with new microfluidic applications that these results suggest are now feasible.
Non-local approach to kinetic effects on parallel transport in fluid models of the scrape-off layer
NASA Astrophysics Data System (ADS)
Omotani, J. T.; Dudson, B. D.
2013-05-01
Using a non-local model, fluid simulations can capture kinetic effects in the parallel electron heat-flux better than is possible using flux limiters in the usual diffusive models. Non-local and diffusive models are compared using a test case representative of an edge-localized mode crash in the JET scrape-off layer (SOL), simulated in one dimension. The non-local model shows substantially enhanced electron temperature gradients, which cannot be achieved using a flux limiter. The performance of the implementation, in the BOUT++ framework, is also analysed to demonstrate its suitability for application in three-dimensional simulations of turbulent transport in the SOL.
Transport of Parallel Momentum Induced by Current-Symmetry Breaking in Toroidal Plasmas
Camenen, Y.; Peeters, A. G.; Casson, F. J.; Hornsby, W. A.; Snodin, A. P.; Angioni, C.; Strintzi, D.
2009-03-27
The symmetry of a physical system strongly impacts on its properties. In toroidal plasmas, the symmetry along a magnetic field line usually constrains the radial flux of parallel momentum to zero in the absence of background flows. By breaking the up-down symmetry of the toroidal currents, this constraint can be relaxed. The parallel asymmetry in the magnetic configuration then leads to an incomplete cancellation of the turbulent momentum flux across a flux surface. The magnitude of the subsequent toroidal rotation increases with the up-down asymmetry and its sign depends on the direction of the toroidal magnetic field and plasma current. Such a mechanism offers new insights in the interpretation and control of the intrinsic toroidal rotation in present day experiments.
Parallelization in SCALE continuous-energy resonance module GEMINEWTRN and transport module NEWT
Zhong, Z.; Downar, T. J.; DeHart, M. D.; Williams, M. L.
2006-07-01
A new resonance module, GEMINEWTRN, has been developed in SCALE, it can calculate the continuous-energy neutron flux within the whole two-dimensional geometry, providing us a rigorous solution. However, the new code needs tremendously amount of computation and memory for practical problem. To relieve the computational burden and memory requirement, parallelization has been implemented into GEMINEWTRN, both angular and spatial decomposition have been adopted so that both the computation and the memory requirement on each processor can be saved considerably, and this effort makes the new resonance method much feasible for practical use. Because the two-dimensional geometry capability and SN/ESC solver of GEMINEWTRN come from lattice physics code NEWT, the similar parallel technique has also been implemented into NEWT, which can also save the computation considerably. (authors)
Parallel domain decomposition methods in fluid models with Monte Carlo transport
Alme, H.J.; Rodrigues, G.H.; Zimmerman, G.B.
1996-12-01
To examine the domain decomposition code coupled Monte Carlo-finite element calculation, it is important to use a domain decomposition that is suitable for the individual models. We have developed a code that simulates a Monte Carlo calculation ( ) on a massively parallel processor. This code is used to examine the load balancing behavior of three domain decomposition ( ) for a Monte Carlo calculation. Results are presented.
Resolving magnetic field line stochasticity and parallel thermal transport in MHD simulations
Nishimura, Y.; Callen, J.D.; Hegna, C.C.
1998-12-31
Heat transport along braided, or chaotic magnetic field lines is a key to understand the disruptive phase of tokamak operations, both the major disruption and the internal disruption (sawtooth oscillation). Recent sawtooth experimental results in the Tokamak Fusion Test Reactor (TFTR) have inferred that magnetic field line stochasticity in the vicinity of the q = 1 inversion radius plays an important role in rapid changes in the magnetic field structures and resultant thermal transport. In this study, the characteristic Lyapunov exponents and spatial correlation of field line behaviors are calculated to extract the characteristic scale length of the microscopic magnetic field structure (which is important for net radial global transport). These statistical values are used to model the effect of finite thermal transport along magnetic field lines in a physically consistent manner.
Tracking log transport and deposition during a 20-year flood in a wide mountain river
NASA Astrophysics Data System (ADS)
Wyżga, Bartłomiej; Mikuś, Paweł; Zawiejska, Joanna; Ruiz-Villanueva, Virginia; Kaczka, Ryszard; Czech, Wiktoria
2016-04-01
Distance of large wood transport during floods and conditions for wood deposition in wide mountain rivers are still insufficiently recognised. Tracking logs tagged with radio transmitters was used to investigate differences in depositional conditions and the length of log displacement during a 20-year flood between channel reaches of different morphology in the Czarny Dunajec River, Polish Carpathians. During a rising limb of the flood, logs were placed into the river at the beginning of an incised reach, close to the beginning of a channelized reach, and 1 km upstream from the beginning of a wide, multi-thread reach. The incised, channelized, and multi-thread reaches retained 12.5%, 33%, and 94% of tagged logs introduced to these reaches, and all the logs retained in the multi-thread reach were deposited in its upstream half. Significant differences in the length of displacement existed between the logs delivered to the river at the three locations, with logs placed into the river at the beginning of the incised reach moved the longest distances and those delivered just upstream from the multi-thread reach the shortest ones. One-fourth of the logs were deposited in a low-flow channel or on channel margin, one-fifth on the floodplain and more than half on gravel bars. After the flood, river cross-sections with deposited logs and a set of cross-sections without wood deposits were surveyed to collect data for one-dimensional modelling of hydraulic conditions at the flood peak. The cross-sections with deposited logs were typified by significantly greater flow width and flow area, and significantly smaller mean flow depth, mean velocity, Froude number, mean bed shear stress and unit stream power. Principal component analysis of the hydraulic parameters in the analysed cross-sections grouped the two types of cross-sections in distinct clusters, indicating that multi-thread cross-sections differed in hydraulic parameters from all the other cross-sections. The experiment
Jolliet, S.; McMillan, B. F.; Vernay, T.; Villard, L.; Hatzky, R.; Bottino, A.; Angelino, P.
2009-07-15
In this paper, the influence of the parallel nonlinearity on zonal flows and heat transport in global particle-in-cell ion-temperature-gradient simulations is studied. Although this term is in theory orders of magnitude smaller than the others, several authors [L. Villard, P. Angelino, A. Bottino et al., Plasma Phys. Contr. Fusion 46, B51 (2004); L. Villard, S. J. Allfrey, A. Bottino et al., Nucl. Fusion 44, 172 (2004); J. C. Kniep, J. N. G. Leboeuf, and V. C. Decyck, Comput. Phys. Commun. 164, 98 (2004); J. Candy, R. E. Waltz, S. E. Parker et al., Phys. Plasmas 13, 074501 (2006)] found different results on its role. The study is performed using the global gyrokinetic particle-in-cell codes TORB (theta-pinch) [R. Hatzky, T. M. Tran, A. Koenies et al., Phys. Plasmas 9, 898 (2002)] and ORB5 (tokamak geometry) [S. Jolliet, A. Bottino, P. Angelino et al., Comput. Phys. Commun. 177, 409 (2007)]. In particular, it is demonstrated that the parallel nonlinearity, while important for energy conservation, affects the zonal electric field only if the simulation is noise dominated. When a proper convergence is reached, the influence of parallel nonlinearity on the zonal electric field, if any, is shown to be small for both the cases of decaying and driven turbulence.
Piras, P; Teresi, L; Traversetti, L; Varano, V; Gabriele, S; Kotsakis, T; Raia, P; Puddu, P E; Scalici, M
2016-05-01
Ontogeny is usually studied by analyzing a deformation series spanning over juvenile to adult shapes. In geometric morphometrics, this approach implies applying generalized Procrustes analysis coupled with principal component analysis on multiple individuals or multiple species datasets. The trouble with such a procedure is that it mixes intra- and inter-group variation. While MANCOVA models are relevant statistical/mathematical tools to draw inferences about the similarities of trajectories, if one wants to observe and interpret the morphological deformation alone by filtering inter-group variability, a particular tool, namely parallel transport, is necessary. In the context of ontogenetic trajectories, one should firstly perform separate multivariate regressions between shape and size, using regression predictions to estimate within-group deformations relative to the smallest individuals. These deformations are then applied to a common reference (the mean of per-group smallest individuals). The estimation of deformations can be performed on the Riemannian manifold by using sophisticated connection metrics. Nevertheless, parallel transport can be effectively achieved by estimating deformations in the Euclidean space via ordinary Procrustes analysis. This approach proved very useful in comparing ontogenetic trajectories of species presenting large morphological differences at early developmental stages. PMID:27161949
Scalable load balancing for massively parallel distributed Monte Carlo particle transport
O'Brien, M. J.; Brantley, P. S.; Joy, K. I.
2013-07-01
In order to run computer simulations efficiently on massively parallel computers with hundreds of thousands or millions of processors, care must be taken that the calculation is load balanced across the processors. Examining the workload of every processor leads to an unscalable algorithm, with run time at least as large as O(N), where N is the number of processors. We present a scalable load balancing algorithm, with run time 0(log(N)), that involves iterated processor-pair-wise balancing steps, ultimately leading to a globally balanced workload. We demonstrate scalability of the algorithm up to 2 million processors on the Sequoia supercomputer at Lawrence Livermore National Laboratory. (authors)
Generalized parallel heat transport equations in collisional to weakly collisional plasmas
NASA Astrophysics Data System (ADS)
Zawaideh, Emad; Kim, N. S.; Najmabadi, Farrokh
1988-11-01
A new set of two-fluid heat-transport equations for heat conduction in collisional to weakly collisional plasmas was derived on the basis of gyrokinetic equations in flux coordinates. In these equations, no restrictions on the anisotropy of the ion distribution function or the collisionality are imposed. In the highly collisional limit, these equations reduce to the classical heat conduction equation of Spitzer and Haerm (1953), while in the weakly collisional limit, they describe a saturated heat flux. Numerical examples comparing these equations with conventional heat transport equations are presented.
Modeling Bimolecular Reactions and Transport in Porous Media via Particle Tracking
NASA Astrophysics Data System (ADS)
Ding, D.; Benson, D. A.; Paster, A.; Bolster, D.
2012-12-01
We use a particle-tracking model to simulate several one-dimensional bi-molecular reactive transport experiments. In this numerical method, the reactants are represented by particles; advection and dispersion dominate the flow of particles, and molecular diffusion dictates, in large part, the reactions. The reactions are determined by a combination of two probabilities; one is that reactant particles occupy the same volume over a short time interval, which is dictated by diffusion; the other is the conditional probability that two collocated particles favorably transform into a reaction. The first probability is a direct physical representation of the degree of mixing in an advancing displacement front, and as such lacks empirical parameters except for the user-defined number of particles. This number is theoretically correlated to the concentration statistics, which can be estimated from concentration autocovariance in the experiment if concentration data is properly collected. The simulations compare favorably to two physical experiments. In one, Raje and Kapoor (2000) [1] measured the product concentrations at the end of a column at different times (the breakthrough curve). In the other, Gramling et al. [2] measured the distribution of reactants and products within a translucent column (snapshots). In addition, one experiment used reactants with a well-mixed thermodynamic rate coefficient 107 times greater than the other. The higher rate can be considered an essentially instantaneous reaction. When compared to the solution of the classical advection-dispersion-reaction equation with the well-mixed reaction coefficient, both experiments showed on the order of 20% to 40% slower reaction attributed to poor mixing. The Lagrangian model in this study accurately simulated the incomplete mixing for both the breakthrough curves [1] and product concentration profile [2]. In addition to model performance, the advantage is the lack of empirical parameters or assumptions
Tunable spin selective transport and quantum phase transition in parallel double dot system
NASA Astrophysics Data System (ADS)
Xiong, Yong-Chen; Wang, Wei-Zhong; Luo, Shi-Jun; Yang, Jun-Tao
2016-02-01
We study theoretically the spin selective transport and the quantum phase transition (QPT) in a double dot device by means of the numerical renormalization group technique. When the gate voltage ε is in the Kondo regime and the interdot hopping t is large enough, a first order QPT between local spin singlet and Sz=1 of the triplet is observed as the magnetic field B increases. Beyond the Kondo regime, the QPTs depend closely on ε and t, and perfect spin filter is found, where the effect of spin filtering could easily be manipulated by tuning external parameters. We show that the interplay between the Zeeman effect and the antiferromagnetic interdot hopping, and occupancy switching are responsible for the QPT and the spin selective transport.
Generalized parallel heat transport equations in collisional to weakly collisional plasmas
Zawaideh, E.; Kim, N.S.; Najmabadi, F.
1988-11-01
A new set of two-fluid heat transport equations that is valid from collisional to weakly collisional limits is derived. Starting from gyrokinetic equations in flux coordinates, a set of moment equations describing plasma energy transport along the field lines of a space- and time-dependent magnetic field is derived. No restrictions on the anisotropy of the ion distribution function or collisionality are imposed. In the highly collisional limit, these equations reduce to the classical heat conduction equation (e.g., Spitzer and Haerm or Braginskii), while in the weakly collisional limit, they describe a saturated heat flux (flux limited). Numerical examples comparing these equations with conventional heat transport equations show that in the limit where the ratio of the mean free path lambda to the scale length of the temperature gradient L/sub T/ approaches zero, there is no significant difference between the solutions of the new and conventional heat transport equations. As lambda/L/sub T/..-->..1, the conventional heat conduction equation contains a significantly larger error than (lambda/L/sub T/)/sup 2/. The error is found to be O(lambda/L)/sup 2/, where L is the smallest of the scale lengths of the gradient in the magnetic field, or the macroscopic plasma parameters (e.g., velocity scale length, temperature scale length, and density scale length). The accuracy of the flux-limited model depends significantly on the value of the flux limit parameter which, in general, is not known. The new set of equations shows that the flux-limited parameter is a function of the magnetic field and plasma parameter profiles.
libmpdata++ 0.1: a library of parallel MPDATA solvers for systems of generalised transport equations
NASA Astrophysics Data System (ADS)
Jaruga, A.; Arabas, S.; Jarecka, D.; Pawlowska, H.; Smolarkiewicz, P. K.; Waruszewski, M.
2014-11-01
This paper accompanies first release of libmpdata++, a C++ library implementing the Multidimensional Positive-Definite Advection Transport Algorithm (MPDATA). The library offers basic numerical solvers for systems of generalised transport equations. The solvers are forward-in-time, conservative and non-linearly stable. The libmpdata++ library covers the basic second-order-accurate formulation of MPDATA, its third-order variant, the infinite-gauge option for variable-sign fields and a flux-corrected transport extension to guarantee non-oscillatory solutions. The library is equipped with a non-symmetric variational elliptic solver for implicit evaluation of pressure gradient terms. All solvers offer parallelisation through domain decomposition using shared-memory parallelisation. The paper describes the library programming interface, and serves as a user guide. Supported options are illustrated with benchmarks discussed in the MPDATA literature. Benchmark descriptions include code snippets as well as quantitative representations of simulation results. Examples of applications include: homogeneous transport in one, two and three dimensions in Cartesian and spherical domains; shallow-water system compared with analytical solution (originally derived for a 2-D case); and a buoyant convection problem in an incompressible Boussinesq fluid with interfacial instability. All the examples are implemented out of the library tree. Regardless of the differences in the problem dimensionality, right-hand-side terms, boundary conditions and parallelisation approach, all the examples use the same unmodified library, which is a key goal of libmpdata++ design. The design, based on the principle of separation of concerns, prioritises the user and developer productivity. The libmpdata++ library is implemented in C++, making use of the Blitz++ multi-dimensional array containers, and is released as free/libre and open-source software.
libmpdata++ 1.0: a library of parallel MPDATA solvers for systems of generalised transport equations
NASA Astrophysics Data System (ADS)
Jaruga, A.; Arabas, S.; Jarecka, D.; Pawlowska, H.; Smolarkiewicz, P. K.; Waruszewski, M.
2015-04-01
This paper accompanies the first release of libmpdata++, a C++ library implementing the multi-dimensional positive-definite advection transport algorithm (MPDATA) on regular structured grid. The library offers basic numerical solvers for systems of generalised transport equations. The solvers are forward-in-time, conservative and non-linearly stable. The libmpdata++ library covers the basic second-order-accurate formulation of MPDATA, its third-order variant, the infinite-gauge option for variable-sign fields and a flux-corrected transport extension to guarantee non-oscillatory solutions. The library is equipped with a non-symmetric variational elliptic solver for implicit evaluation of pressure gradient terms. All solvers offer parallelisation through domain decomposition using shared-memory parallelisation. The paper describes the library programming interface, and serves as a user guide. Supported options are illustrated with benchmarks discussed in the MPDATA literature. Benchmark descriptions include code snippets as well as quantitative representations of simulation results. Examples of applications include homogeneous transport in one, two and three dimensions in Cartesian and spherical domains; a shallow-water system compared with analytical solution (originally derived for a 2-D case); and a buoyant convection problem in an incompressible Boussinesq fluid with interfacial instability. All the examples are implemented out of the library tree. Regardless of the differences in the problem dimensionality, right-hand-side terms, boundary conditions and parallelisation approach, all the examples use the same unmodified library, which is a key goal of libmpdata++ design. The design, based on the principle of separation of concerns, prioritises the user and developer productivity. The libmpdata++ library is implemented in C++, making use of the Blitz++ multi-dimensional array containers, and is released as free/libre and open-source software.
NASA Astrophysics Data System (ADS)
Hyang Park, Young; Vivier, Frédéric; Roquet, Fabien; Kestenare, Elodie; Sekma, Hela; Durand, Isabelle
2010-05-01
The Kerguelen Plateau is considered as a major barrier to the eastward flowing Antarctic Circumpolar Current (ACC), diverting about 2/3 of its transport to the north of the plateau and the remaining transport across the vast plateau developed between the Kerguelen Islands and Antarctica. However, due to the lack of systematic high-quality observations especially across the Fawn Trough, our knowledge of ACC branches and associated transports across the Kerguelen Plateau has been largely indirect and debated, with previous transport estimates ranging from 30 to 100 Sv. In order to fill this knowledge gap, we undertook full-depth CTD and LADCP measurements in the Fawn Trough area and moored three lines of a total of 9 current-meters across the Fawn Trough in February-March 2009 during the TRACK-1 cruise. The TRACK-1 data have permitted us to document for the first time a detailed vertical current structure especially across the Fawn Trough as well as the Deep Western Boundary Current (DWBC) on the eastern flank of the Kerguelen Plateau at 58°S. The Fawn Trough Current appears as a principal across-plateau jet associated with the Southern ACC Front, concentrating the majority (43 Sv) of the eastward transport passing to the south of the Kerguelen Islands (58 Sv), compared to ~92 Sv at the Subantarctic Front as estimated previously to the north of the islands. This yields a total transport of ~150 Sv for the ACC transiting through the Kerguelen longitude, consistent with other independent estimates in the other sectors of the Southern Ocean, but with an up-to-date partition of ~60% of the transport passing to the north of the Kerguelen Plateau and ~40% across the plateau. The DWBC transport at 58°S is estimated as 43 Sv, of which 36 Sv come from the northward turning of the westward flowing Antarctic Slope Current and the rest (7 Sv) originates from the southernmost branch of the ACC passing through the northern Princess Elizabeth Trough. Two other transport branches
Alfven wave transport effects in the time evolution of parallel cosmic-ray modified shocks
NASA Technical Reports Server (NTRS)
Jones, T. W.
1993-01-01
Some of the issues associated with a more complete treatment of Alfven transport in cosmic ray shocks are explored qualitatively. The treatment is simplified in some important respects, but some new issues are examined and for the first time a nonlinear, time dependent study of plane cosmic ray mediated shocks with both the entropy producing effects of wave dissipation and effects due to the Alfven wave advection of the cosmic ray relative to the gas is included. Examination of the direct consequences of including the pressure and energy of the Alfven waves in the formalism began.
Using Activated Transport in Parallel Nanowires for Energy Harvesting and Hot-Spot Cooling
NASA Astrophysics Data System (ADS)
Bosisio, Riccardo; Gorini, Cosimo; Fleury, Geneviève; Pichard, Jean-Louis
2015-05-01
We study arrays of parallel doped semiconductor nanowires in a temperature range where the electrons propagate through the nanowires by phonon-assisted hops between localized states. By solving the random-resistor-network problem, we compute the thermopower S , the electrical conductance G , and the electronic thermal conductance Ke of the device. We investigate how those quantities depend on the position—which can be tuned with a back gate—of the nanowire impurity band with respect to the equilibrium electrochemical potential. We show that large power factors can be reached near the band edges, when S self-averages to large values while G is small but scales with the number of wires. Calculating the amount of heat exchanged locally between the electrons inside the nanowires and the phonons of the environment, we show that phonons are mainly absorbed near one electrode and emitted near the other when a charge current is driven through the nanowires near their band edges. This phenomenon could be exploited for a field control of the heat exchange between the phonons and the electrons at submicron scales in electronic circuits. It could be also used for cooling hot spots.
NASA Astrophysics Data System (ADS)
Kihm, You Hong; Kim, Jeong Hwan
The rocks in the central Ogcheon Belt, South Korea, underwent three deformational phases (D 1, D 2 and D 3) during the Mesozoic Era. In the study area, NW- and NE-trending geological structures such as thrusts and folds are juxtaposed. The NW- and NE-trending folds deform isoclinal folds (F 1) with axial planar slaty cleavage (S 1), and have axial planar crenulation cleavage. All of these folds are overprinted by brittle faults (D 3). It is inferred from field data that the NW- and NE-trending structures are products of a single deformational phase (D 2). A very large-scale inclined NE-vergent isoclinal fold (F 1) produced an irregular boundary of Precambrian basement and a basement promontory, which acted as a structural obstacle against subsequent deformational phases. During the E-vergent phase (D 2), deformation partitioning occurred due to the irregular block boundary, and orogen-parallel and -orthogonal structures were produced. The D 3 phase is recognized as large-scale E-W trending structures including folds and faults. Thus, the structural evolution of the central Ogcheon Belt is related to the clockwise rotation of the maximum compressive stress axis from NE-SW to N-S during the Mesozoic Era. This study shows that the shape of colliding boundary is a very important factor in controlling the structural pattern and evolution in the study area.
Yaqi Wang; Cristian Rabiti; Giuseppe Palmiotti
2011-06-01
The Red-Black algorithm has been successfully applied on solving the second-order parity transport equation with the PN approximation in angle and the Hybrid Finite Element Method (HFEM) in space, i.e., the Variational Nodal Method (VNM) [1,2,3,4,5]. Any transport solving techniques, including the Red-Black algorithm, need to be parallelized in order to take the advantage of the development of supercomputers with multiple processors for the advanced modeling and simulation. To our knowledge, an attempt [6] was done to parallelize it, but it was devoted only to the z axis plans in three-dimensional calculations. General parallelization of the Red-Black algorithm with the spatial domain decomposition has not been reported in the literature. In this summary, we present our implementation of the parallelization of the Red-Black algorithm and its efficiency results.
Yoshida, Hiroyuki; Takatsuji; Wataru
2000-04-01
The parallel transport of an organic acid by solid-phase and macropore diffusion within a porous ion exchanger was studied by measuring equilibrium isotherms and uptake curves for adsorption of acetic acid and lactic acid on a weakly basic ion exchanger, DIAION WA30. Experimental adsorption isotherms were correlated by the Langmuir equation. The Langmuir equilibrium constant of acetic acid was close to that of lactic acid, and the saturation capacity of acetic acid was about 84% that of lactic acid. Intraparticle effective diffusivity D{sub eff} was determined using the homogeneous Fickian diffusion model. The value of D{sub eff} for acetic acid was about 1.5 times lactic acid. Because D{sub eff} increased with linearly increasing bulk phase concentration C{sub 0}, D{sub eff} was separated to the solid-phase diffusivity D{sub s} and the macropore diffusivity D{sub P} by applying the parallel diffusion model. The model agreed well with the experimental curves. The values of D{sub S} and D{sub P} for acetic acid were about 2 and 1.5 times those of lactic acid, respectively. The acetic acid and the lactic acid may be separated by the difference of the diffusion rates.
Toroidal momentum transport in a tokamak caused by symmetry breaking parallel derivatives
Sung, T.; Buchholz, R.; Grosshauser, S. R.; Hornsby, W. A.; Migliano, P.; Peeters, A. G.; Casson, F. J.; Fable, E.
2013-04-15
A new mechanism for toroidal momentum transport in a tokamak is investigated using the gyro-kinetic model. First, an analytic model is developed through the use of the ballooning transform. The terms that generate the momentum transport are then connected with the poloidal derivative of the ballooning envelope, which are one order smaller in the normalised Larmor radius, compared with the derivative of the eikonal. The mechanism, therefore, does not introduce an inhomogeneity in the radial direction, in contrast with the effect of profile shearing. Numerical simulations of the linear ion temperature gradient mode with adiabatic electrons, retaining the finite {rho}{sub *} effects in the E Multiplication-Sign B velocity, the drift, and the gyro-average, are presented. The momentum flux is found to be linear in the normalised Larmor radius ({rho}{sub *}) but is, nevertheless, generating a sizeable counter-current rotation. The total momentum flux scales linear with the aspect ratio of the considered magnetic surface, and increases with increasing magnetic shear, safety factor, and density and temperature gradients.
Parallel transport studies of high-Z impurities in the core of Alcator C-Mod plasmasa)
NASA Astrophysics Data System (ADS)
Reinke, M. L.; Hutchinson, I. H.; Rice, J. E.; Greenwald, M.; Howard, N. T.; Hubbard, A.; Hughes, J. W.; Terry, J. L.; Wolfe, S. M.
2013-05-01
Measurements of poloidal variation, ñz/⟨nz⟩, in high-Z impurity density have been made using photodiode arrays sensitive to vacuum ultraviolet and soft x-ray emission in Alcator C-Mod plasmas. In/out asymmetries in the range of -0.2
Li, Zhelong; Zhang, Dongxiao; Li, Xiqing
2010-02-15
Advances in pore structure characterization and lattice-Boltzmann (LB) simulations of flow fields in pore spaces are making mechanistic simulations of colloid transport in real porous media a realistic goal. The primary challenge to reach this goal may be the computational demand of LB flow simulations in discretized porous medium domains at an assemblage scale. In this work, flow fields in simple cubic and dense packing systems were simulated at different discretization resolutions using the LB method. The simulated flow fields were incorporated into to a three-dimensional particle tracking model to simulate colloid transport in the two systems. The simulated colloid deposition tended to become asymptotic at a critical discretization resolution (voxel-grain size ratio = 0.01) at groundwater flow regimes for colloids down to submicrometer level under favorable conditions and down to around 1 microm under unfavorable conditions. The average simulated fluid velocities near grain surfaces were extracted to explain the sensitivities of simulated depositions to space discretization under both conditions. At the critical discretization resolution, current computation capacity would allow flow simulations and particle tracking in assemblage porous medium domains. In addition, particle tracking simulations revealed that colloids may be retained in flow vortices under conditions both favorable and unfavorable for deposition. Colloid retention in flow vortices has been proposed only very recently. Here we provide a mechanistic confirmation to this novel retention process. PMID:20088544
Analysis of Massively Parallel Discrete-Ordinates Transport Sweep Algorithms with Collisions
Bailey, T S; Falgout, R D
2008-10-14
We present theoretical scaling models for a variety of discrete-ordinates sweep algorithms. In these models, we pay particular attention to the way each algorithm handles collisions. A collision is defined as a processor having multiple angles with ready to be swept during one stage of the sweep. The models also take into account how subdomains are assigned to processors and how angles are grouped during the sweep. We describe a data driven algorithm that resolves collisions efficiently during the sweep as well as other algorithms that have been designed to avoid collisions completely. Our models are validated using the ARGES and AMTRAN transport codes. We then use the models to study and predict scaling trends in all of the sweep algorithms.
Electron Beam-Blip Spectroscopic Diagnostics of the Scrape-off-Layer Parallel Transport in C-2
NASA Astrophysics Data System (ADS)
Osin, Dmitry; Thompson, Matthew; Garate, Eusebio; TAE Team
2015-11-01
C-2 is a microscopically stable, high-performance field-reversed configuration (FRC), where high plasma temperatures with significant fast ion population and record lifetimes were achieved by a combination of tangential neutral beam injection, electrically biased plasma guns at the ends and wall conditioning. FRC confinement depends on the properties of both the open and closed field lines, therefore, understanding the electron transport in the scrape-of-layer (SOL) is critical. To study parallel heat conduction in SOL, a high-energy pulsed electron beam (e-beam) was injected on-axis into C-2 to produce a heat pulse, which causes a fast rise and slower decay of the electron temperature, Te, in the SOL. The heat-blip was observed by means of He-jet spectroscopy. A small fraction of the total deposited e-beam energy is necessary to explain the measured Te increase. The electron thermal conductivity along the magnetic field lines can be inferred from the Te decay. Experiments suggest that a high energy e-beam pulse can serve as a direct diagnostic of heat transport in the SOL.
Archer, Charles Jens; Musselman, Roy Glenn; Peters, Amanda; Pinnow, Kurt Walter; Swartz, Brent Allen; Wallenfelt, Brian Paul
2010-11-16
A massively parallel computer system contains an inter-nodal communications network of node-to-node links. An automated routing strategy routes packets through one or more intermediate nodes of the network to reach a destination. Some packets are constrained to be routed through respective designated transporter nodes, the automated routing strategy determining a path from a respective source node to a respective transporter node, and from a respective transporter node to a respective destination node. Preferably, the source node chooses a routing policy from among multiple possible choices, and that policy is followed by all intermediate nodes. The use of transporter nodes allows greater flexibility in routing.
Parallel transport studies of high-Z impurities in the core of Alcator C-Mod plasmas
Reinke, M. L.; Hutchinson, I. H.; Rice, J. E.; Greenwald, M.; Howard, N. T.; Hubbard, A.; Hughes, J. W.; Terry, J. L.; Wolfe, S. M.
2013-05-15
Measurements of poloidal variation, ñ{sub z}/
NASA Astrophysics Data System (ADS)
Naeger, Aaron R.; Gupta, Pawan; Zavodsky, Bradley T.; McGrath, Kevin M.
2016-06-01
The primary goal of this study was to generate a near-real time (NRT) aerosol optical depth (AOD) product capable of providing a comprehensive understanding of the aerosol spatial distribution over the Pacific Ocean, in order to better monitor and track the trans-Pacific transport of aerosols. Therefore, we developed a NRT product that takes advantage of observations from both low-earth orbiting and geostationary satellites. In particular, we utilize AOD products from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Suomi National Polar-orbiting Partnership (NPP) Visible Infrared Imaging Radiometer Suite (VIIRS) satellites. Then, we combine these AOD products with our own retrieval algorithms developed for the NOAA Geostationary Operational Environmental Satellite (GOES-15) and Japan Meteorological Agency (JMA) Multi-functional Transport Satellite (MTSAT-2) to generate a NRT daily AOD composite product. We present examples of the daily AOD composite product for a case study of trans-Pacific transport of Asian pollution and dust aerosols in mid-March 2014. Overall, the new product successfully tracks this aerosol plume during its trans-Pacific transport to the west coast of North America as the frequent geostationary observations lead to a greater coverage of cloud-free AOD retrievals equatorward of about 35° N, while the polar-orbiting satellites provide a greater coverage of AOD poleward of 35° N. However, we note several areas across the domain of interest from Asia to North America where the GOES-15 and MTSAT-2 retrieval algorithms can introduce significant uncertainties into the new product.
NASA Astrophysics Data System (ADS)
Naeger, A. R.; Gupta, P.; Zavodsky, B.; McGrath, K. M.
2015-10-01
The primary goal of this study was to generate a near-real time (NRT) aerosol optical depth (AOD) product capable of providing a comprehensive understanding of the aerosol spatial distribution over the Pacific Ocean in order to better monitor and track the trans-Pacific transport of aerosols. Therefore, we developed a NRT product that takes advantage of observations from both low-earth orbiting and geostationary satellites. In particular, we utilize AOD products from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Suomi National Polar-orbiting Partnership (NPP) Visible Infrared Imaging Radiometer Suite (VIIRS) satellites. Then, we combine these AOD products with our own retrieval algorithms developed for the NOAA Geostationary Operational Environmental Satellite (GOES-15) and Japan Meteorological Agency (JMA) Multi-functional Transport Satellite (MTSAT-2) to generate a NRT daily AOD composite product. We present examples of the daily AOD composite product for a case study of trans-Pacific transport of Asian pollution and dust aerosols in mid-March 2014. Overall, the new product successfully tracks this aerosol plume during its trans-Pacific transport to the west coast of North America. However, we identify several areas across the domain of interest from Asia to North America where the new product can encounter significant uncertainties due to the inclusion of the geostationary AOD retrievals. The uncertainties associated with geostationary AOD retrievals are expected to be minimized after the successful launch of the next-generation advanced NOAA GOES-R and recently launched JMA Himawari satellites. Observations from these advanced satellites will ultimately provide an enhanced understanding of the spatial and temporal distribution of aerosols over the Pacific.
Rosa, Massimiliano; Warsa, James S; Perks, Michael
2010-12-14
We have implemented a cell-wise, block-Gauss-Seidel (bGS) iterative algorithm, for the solution of the S{sub n} transport equations on the Roadrunner hybrid, parallel computer architecture. A compute node of this massively parallel machine comprises AMD Opteron cores that are linked to a Cell Broadband Engine{trademark} (Cell/B.E.). LAPACK routines have been ported to the Cell/B.E. in order to make use of its parallel Synergistic Processing Elements (SPEs). The bGS algorithm is based on the LU factorization and solution of a linear system that couples the fluxes for all S{sub n} angles and energy groups on a mesh cell. For every cell of a mesh that has been parallel decomposed on the higher-level Opteron processors, a linear system is transferred to the Cell/B.E. and the parallel LAPACK routines are used to compute a solution, which is then transferred back to the Opteron, where the rest of the computations for the S{sub n} transport problem take place. Compared to standard parallel machines, a hundred-fold speedup of the bGS was observed on the hybrid Roadrunner architecture. Numerical experiments with strong and weak parallel scaling demonstrate the bGS method is viable and compares favorably to full parallel sweeps (FPS) on two-dimensional, unstructured meshes when it is applied to optically thick, multi-material problems. As expected, however, it is not as efficient as FPS in optically thin problems.
NASA Astrophysics Data System (ADS)
Kravets, L. I.; Dmitriev, S. N.; Drachev, A. I.; Gilman, A. B.; Lazea, A.; Dinescu, G.
2007-04-01
A process of plasma polymerization of dimethylaniline and acrylic acid vapours on the surface of poly(ethylene terephthalate) track membranes has been investigated. The surface and hydrodynamic properties of the composite membranes produced in this case have been studied. It is shown that the water permeability of the obtained polymeric membranes can be controlled by changing the filtrate pH. Membranes with such properties can be used for controllable drug delivery and in sensor control.
Joseph, D.D.; Bai, R.; Liao, T.Y.; Huang, A.; Hu, H.H.
1995-09-01
In this paper the authors introduce the idea of parallel pipelining for water lubricated transportation of oil (or other viscous material). A parallel system can have major advantages over a single pipe with respect to the cost of maintenance and continuous operation of the system, to the pressure gradients required to restart a stopped system and to the reduction and even elimination of the fouling of pipe walls in continuous operation. The authors show that the action of capillarity in small pipes is more favorable for restart than in large pipes. In a parallel pipeline system, they estimate the number of small pipes needed to deliver the same oil flux as in one larger pipe as N = (R/r){sup {alpha}}, where r and R are the radii of the small and large pipes, respectively, and {alpha} = 4 or 19/7 when the lubricating water flow is laminar or turbulent.
Supej, Matej; Cuk, Ivan
2014-01-01
Global Navigation Satellite Systems (GNSS) are, in addition to being most widely used vehicle navigation method, becoming popular in sport-related tests. There is a lack of knowledge regarding tracking speed using GNSS, therefore the aims of this study were to examine under dynamic conditions: (1) how accurate technologically different GNSS measure speed and (2) how large is latency in speed measurements in real time applications. Five GNSSs were tested. They were fixed to a car's roof-rack: a smart phone, a wrist watch, a handheld device, a professional system for testing vehicles and a high-end Real Time Kinematics (RTK) GNSS. The speed data were recorded and analyzed during rapid acceleration and deceleration as well as at steady speed. The study produced four main findings. Higher frequency and high quality GNSS receivers track speed at least at comparable accuracy to a vehicle speedometer. All GNSS systems measured maximum speed and movement at a constant speed well. Acceleration and deceleration have different level of error at different speeds. Low cost GNSS receivers operating at 1 Hz sampling rate had high latency (up to 2.16 s) and are not appropriate for tracking speed in real time, especially during dynamic movements. PMID:25494349
Supej, Matej; Čuk, Ivan
2014-01-01
Global Naavigation Satellite Systems (GNSS) are, in addition to being most widely used vehicle navigation method, becoming popular in sport-related tests. There is a lack of knowledge regarding tracking speed using GNSS, therefore the aims of this study were to examine under dynamic conditions: (1) how accurate technologically different GNSS measure speed and (2) how large is latency in speed measurements in real time applications. Five GNSSs were tested. They were fixed to a car's roof-rack: a smart phone, a wrist watch, a handheld device, a professional system for testing vehicles and a high-end Real Time Kinematics (RTK) GNSS. The speed data were recorded and analyzed during rapid acceleration and deceleration as well as at steady speed. The study produced four main findings. Higher frequency and high quality GNSS receivers track speed at least at comparable accuracy to a vehicle speedometer. All GNSS systems measured maximum speed and movement at a constant speed well. Acceleration and deceleration have different level of error at different speeds. Low cost GNSS receivers operating at 1 Hz sampling rate had high latency (up to 2.16 s) and are not appropriate for tracking speed in real time, especially during dynamic movements. PMID:25494349
Mair, L.O.; Superfine, R.
2014-01-01
Magnetic drug targeting has been proposed as a means of efficiently targeting drugs to tumors. However, the extracellular matrix (ECM) remains a significant barrier to long-range magnetophoretic transport through the tumor volume. While ensemble measurements of nanoparticle magnetophoresis have been reported, a single particle level understanding of magnetophoretic transport remains at large. We quantify nanorod magnetophoresis through ECM based on single particle observations. We find that smaller diameter particles achieve larger velocities through ECM despite experiencing smaller magnetic forces. Additionally, two interesting dynamics are elucidated. First, 18 nm diameter nanorods experience bimodal stick-slip motion through ECM during static field magnetophoresis, while similar bimodal transport is not observed for 55 nm nor 200 nm diameter nanorods. Second, smaller particles experience larger deviations in their orientation angle with respect to the magnetic field. This work elucidates important dynamics of nanoparticle transport through complex, porous biomaterials that may go unnoticed during ensemble measurements. PMID:24744160
Particle-tracking simulations of anomalous transport in hierarchically fractured rocks
NASA Astrophysics Data System (ADS)
Roubinet, Delphine; de Dreuzy, Jean-Raynald; Tartakovsky, Daniel M.
2013-01-01
Complex topology of fracture networks and interactions between transport processes in a fracture and the ambient un-fractured rock (matrix) combine to render modeling solute transport in fractured media a challenge. Classical approaches rely on both strong assumptions of either limited or full diffusion of solutes in the matrix and simplified fracture configurations. We analyze fracture-matrix transport in two-dimensional Sierpinski lattice structures, which display a wide range of matrix block sizes. The analysis is conducted in several transport regimes that are limited by either diffusion or block sizes. Our simulation results can be used to validate the simplifying assumptions that underpin classical analytical solutions and to benchmark other numerical methods. They also demonstrate that both hydraulic and structural properties of fractured rocks control the residence time distribution.
Well-to-Wheels Water Consumption: Tracking the Virtual Flow of Water into Transportation
NASA Astrophysics Data System (ADS)
Lampert, D. J.; Elgowainy, A.; Hao, C.
2015-12-01
Water and energy resources are fundamental to life on Earth and essential for the production of consumer goods and services in the economy. Energy and water resources are heavily interdependent—energy production consumes water, while water treatment and distribution consume energy. One example of this so-called energy-water nexus is the consumption of water associated with the production of transportation fuels. The Greenhouse Gases, Regulated Emissions, and Energy Use in Transportation (GREET) model is an analytical tool that can be used to compare the environmental impacts of different transportation fuels on a consistent basis. In this presentation, the expansion of GREET to perform life cycle water accounting or the "virtual flow" of water into transportation and other energy sectors and the associated implications will be discussed. The results indicate that increased usage of alternative fuels may increase freshwater resource consumption. The increased water consumption must be weighed against the benefits of decreased greenhouse gas and fossil energy consumption. Our analysis highlights the importance of regionality, co-product allocation, and consistent system boundaries when comparing the water intensity of alternative transportation fuel production pathways such as ethanol, biodiesel, compressed natural gas, hydrogen, and electricity with conventional petroleum-based fuels such as diesel and gasoline.
Qiang, J.; Leitner, D.; Todd, D.S.; Ryne, R.D.
2005-03-15
The superconducting ECR ion source VENUS serves as the prototype injector ion source for the Rare Isotope Accelerator (RIA) driver linac. The RIA driver linac requires a great variety of high charge state ion beams with up to an order of magnitude higher intensity than currently achievable with conventional ECR ion sources. In order to design the beam line optics of the low energy beam line for the RIA front end for the wide parameter range required for the RIA driver accelerator, reliable simulations of the ion beam extraction from the ECR ion source through the ion mass analyzing system are essential. The RIA low energy beam transport line must be able to transport intense beams (up to 10 mA) of light and heavy ions at 30 keV.For this purpose, LBNL is developing the parallel 3D particle-in-cell code IMPACT to simulate the ion beam transport from the ECR extraction aperture through the analyzing section of the low energy transport system. IMPACT, a parallel, particle-in-cell code, is currently used to model the superconducting RF linac section of RIA and is being modified in order to simulate DC beams from the ECR ion source extraction. By using the high performance of parallel supercomputing we will be able to account consistently for the changing space charge in the extraction region and the analyzing section. A progress report and early results in the modeling of the VENUS source will be presented.
NASA Astrophysics Data System (ADS)
Qiang, J.; Leitner, D.; Todd, D. S.; Ryne, R. D.
2005-03-01
The superconducting ECR ion source VENUS serves as the prototype injector ion source for the Rare Isotope Accelerator (RIA) driver linac. The RIA driver linac requires a great variety of high charge state ion beams with up to an order of magnitude higher intensity than currently achievable with conventional ECR ion sources. In order to design the beam line optics of the low energy beam line for the RIA front end for the wide parameter range required for the RIA driver accelerator, reliable simulations of the ion beam extraction from the ECR ion source through the ion mass analyzing system are essential. The RIA low energy beam transport line must be able to transport intense beams (up to 10 mA) of light and heavy ions at 30 keV. For this purpose, LBNL is developing the parallel 3D particle-in-cell code IMPACT to simulate the ion beam transport from the ECR extraction aperture through the analyzing section of the low energy transport system. IMPACT, a parallel, particle-in-cell code, is currently used to model the superconducting RF linac section of RIA and is being modified in order to simulate DC beams from the ECR ion source extraction. By using the high performance of parallel supercomputing we will be able to account consistently for the changing space charge in the extraction region and the analyzing section. A progress report and early results in the modeling of the VENUS source will be presented.
Wrong side of the tracks: the neglected human costs of transporting oil and gas.
Burton, Lloyd; Stretesky, Paul
2014-01-01
The connection between human rights and climate change is most evident when examining carbon dioxide emissions that result from burning fossil fuels (e.g., sea level rise and displaced coastal cultures). However, the transport of fossil fuels also has human rights implications for human rights and climate change. This research focuses on the health and safety risks inflicted on those residents who are adjacent to the railways that ship fossil fuels along the US-Canada transportation corridors. Applying sociological and jurisprudential perspectives, we review the environmental/climate justice literature as it pertains to industrial transport corridors, documenting the forms of heightened risk imposed on those living along these transportation paths. Next, we develop an illustrative case study of Canada's worst rail catastrophe. In so doing, we provide evidence of a decades-long failure of US and Canadian regulators to prevent such disasters. We interpret that disaster through a human rights case law suggesting that States have an affirmative duty to protect their citizens from foreseeable disasters. Based on this analysis, we propose specific rail safety regulatory reforms. We argue that if the US and Canadian governments implement these regulations as required under human rights law, they can more effectively honor their obligations to their citizens who are paying a high human cost for the material benefits associated with increased energy production. PMID:25474613
Lewis, Christina M.; Graves, Stephen A.; Hernandez, Reinier; Valdovinos, Hector F.; Barnhart, Todd E.; Cai, Weibo; Meyerand, Mary E.; Nickles, Robert J.; Suzuki, Masatoshi
2015-01-01
There is a growing demand for long-term in vivo stem cell imaging for assessing cell therapy techniques and guiding therapeutic decisions. This work develops the production of 52Mn and establishes proof of concept for the use of divalent metal transporter 1 (DMT1) as a positron emission tomography (PET) and magnetic resonance imaging (MRI) reporter gene for stem cell tracking in the rat brain. 52Mn was produced via proton irradiation of a natural chromium target. In a comparison of two 52Mn separation methods, solvent-solvent extraction was preferred over ion exchange chromatography because of reduced chromium impurities and higher 52Mn recovery. In vitro uptake of Mn-based PET and MRI contrast agents (52Mn2+ and Mn2+, respectively) was enhanced in DMT1 over-expressing human neural progenitor cells (hNPC-DMT1) compared to wild-type control cells (hNPC-WT). After cell transplantation in the rat striatum, increased uptake of Mn-based contrast agents in grafted hNPC-DMT1 was detected in in vivo manganese-enhanced MRI (MEMRI) and ex vivo PET and autoradiography. These initial studies indicate that this approach holds promise for dual-modality PET/MR tracking of transplanted stem cells in the central nervous system and prompt further investigation into the clinical applicability of this technique. PMID:25553111
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.
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.
Tracking aeolian transport patterns across a mega-nourishment using video imagery
NASA Astrophysics Data System (ADS)
Wijnberg, Kathelijne; van der Weerd, Lianne; Hulscher, Suzanne
2014-05-01
Coastal dune areas protect the hinterland from flooding. In order to maintain the safety level provided by the dunes, it may be necessary to artificially supply the beach-dune system with sand. How to best design these shore nourishments, amongst others with respect to optimal dune growth on the long-term (decadal scale), is not yet clear. One reason for this is that current models for aeolian transport on beaches appear to have limited predictive capabilities regarding annual onshore sediment supply. These limited capabilities may be attributed to the lack of appropriate input data, for instance on moisture content of the beach surface, or shortcomings in process understanding. However, it may also be argued that for the long-term prediction of onshore aeolian sand supply from the beach to the dunes, we may need to develop some aggregated-scale transport equations, because the detailed input data required for the application of process-scale transport equations may never be available in reality. A first step towards the development of such new concepts for aggregated-scale transport equations is to increase phenomenological insight into the characteristics and number of aeolian transport events that account for the annual volume changes of the foredunes. This requires high-frequency, long-term data sets to capture the only intermittently occurring aeolian transport events. Automated video image collection seems a promising way to collect such data. In the present study we describe the movement (direction and speed) of sand patches and aeolian bed forms across a nourished site, using video imagery, to characterize aeolian transport pathways and their variability in time. The study site is a mega-nourishment (21 Mm3 of sand) that was recently constructed at the Dutch coast. This mega-nourishment, also referred to as the Sand Motor, is a pilot project that may potentially replace current practice of more frequently applying small scale nourishments. The mega
Imaging and Quantitation Techniques for Tracking Cargo along Endosome-to-Golgi Transport Pathways
Chia, Pei Zhi Cheryl; Gleeson, Paul A.
2013-01-01
Recent improvements in the resolution of light microscopy, coupled with the development of a range of fluorescent-based probes, have provided new approaches to dissecting membrane domains and the regulation of membrane trafficking. Here, we review these advances, as well as highlight developments in quantitative image analysis and novel unbiased analytical approaches to quantitate protein localization. The application of these approaches to endosomal sorting and endosome-to-Golgi transport is discussed. PMID:24709647
Peeters, A. G.; Camenen, Y.; Casson, F. J.; Hornsby, W. A.; Snodin, A. P.; Strintzi, D.; Angioni, C.
2009-04-15
The paper derives the gyro-kinetic equation in the comoving frame of a toroidally rotating plasma, including both the Coriolis drift effect [A. G. Peeters et al., Phys. Rev. Lett. 98, 265003 (2007)] as well as the centrifugal force. The relation with the laboratory frame is discussed. A low field side gyro-fluid model is derived from the gyro-kinetic equation and applied to the description of parallel momentum transport. The model includes the effects of the Coriolis and centrifugal force as well as the parallel dynamics. The latter physics effect allows for a consistent description of both the Coriolis drift effect as well as the ExB shear effect [R. R. Dominguez and G. M. Staebler, Phys. Fluids B 5, 3876 (1993)] on the momentum transport. Strong plasma rotation as well as parallel dynamics reduce the Coriolis (inward) pinch of momentum and can lead to a sign reversal generating an outward pinch velocity. Also, the ExB shear effect is, in a similar manner, reduced by the parallel dynamics and stronger rotation.
NASA Astrophysics Data System (ADS)
Guo, L.; Huang, H.; Gaston, D.; Redden, G. D.
2009-12-01
One approach for immobilizing subsurface metal contaminants involves stimulating the in situ production of mineral phases that sequester or isolate contaminants. One example is using calcium carbonate to immobilize strontium. The success of such approaches depends on understanding how various processes of flow, transport, reaction and resulting porosity-permeability change couple in subsurface systems. Reactive transport models are often used for such purpose. Current subsurface reactive transport simulators typically involve a de-coupled solution approach, such as operator-splitting, that solves the transport equations for components and batch chemistry sequentially, which has limited applicability for many biogeochemical processes with fast kinetics and strong medium property-reaction interactions. A massively parallel, fully coupled, fully implicit reactive transport simulator has been developed based on a parallel multi-physics object oriented software environment computing framework (MOOSE) developed at the Idaho National Laboratory. Within this simulator, the system of transport and reaction equations is solved simultaneously in a fully coupled manner using the Jacobian Free Newton-Krylov (JFNK) method with preconditioning. The simulator was applied to model reactive transport in a one-dimensional column where conditions that favor calcium carbonate precipitation are generated by urea hydrolysis that is catalyzed by urease enzyme. Simulation results are compared to both laboratory column experiments and those obtained using the reactive transport simulator STOMP in terms of: the spatial and temporal distributions of precipitates and reaction rates and other major species in the reaction system; the changes in porosity and permeability; and the computing efficiency based on wall clock simulation time.
Baba, Shoji Sailer, Juergen; Deacon, Russell S.; Oiwa, Akira; Shibata, Kenji; Hirakawa, Kazuhiko; Tarucha, Seigo
2015-11-30
We report conductance and supercurrent measurements for InAs single and parallel double quantum dot Josephson junctions contacted with Nb or NbTiN superconducting electrodes. Large superconducting gap energy, high critical field, and large switching current are observed, all reflecting the features of Nb-based electrodes. For the parallel double dots, we observe an enhanced supercurrent when both dots are on resonance, which may reflect split Cooper pair tunneling.
NASA Astrophysics Data System (ADS)
Bian, Nicolas H.; Kontar, Eduard P.; Emslie, A. Gordon
2016-06-01
The transport of the energy contained in electrons, both thermal and suprathermal, in solar flares plays a key role in our understanding of many aspects of the flare phenomenon, from the spatial distribution of hard X-ray emission to global energetics. Motivated by recent RHESSI observations that point to the existence of a mechanism that confines electrons to the coronal parts of flare loops more effectively than Coulomb collisions, we here consider the impact of pitch-angle scattering off turbulent magnetic fluctuations on the parallel transport of electrons in flaring coronal loops. It is shown that the presence of such a scattering mechanism in addition to Coulomb collisional scattering can significantly reduce the parallel thermal and electrical conductivities relative to their collisional values. We provide illustrative expressions for the resulting thermoelectric coefficients that relate the thermal flux and electrical current density to the temperature gradient and the applied electric field. We then evaluate the effect of these modified transport coefficients on the flare coronal temperature that can be attained, on the post-impulsive-phase cooling of heated coronal plasma, and on the importance of the beam-neutralizing return current on both ambient heating and the energy loss rate of accelerated electrons. We also discuss the possible ways in which anomalous transport processes have an impact on the required overall energy associated with accelerated electrons in solar flares.
49 CFR 213.63 - Track surface.
Code of Federal Regulations, 2011 CFR
2011-10-01
... 49 Transportation 4 2011-10-01 2011-10-01 false Track surface. 213.63 Section 213.63..., DEPARTMENT OF TRANSPORTATION TRACK SAFETY STANDARDS Track Geometry § 213.63 Track surface. Each owner of the track to which this part applies shall maintain the surface of its track within the limits prescribed...
49 CFR 213.63 - Track surface.
Code of Federal Regulations, 2010 CFR
2010-10-01
... 49 Transportation 4 2010-10-01 2010-10-01 false Track surface. 213.63 Section 213.63..., DEPARTMENT OF TRANSPORTATION TRACK SAFETY STANDARDS Track Geometry § 213.63 Track surface. Each owner of the track to which this part applies shall maintain the surface of its track within the limits prescribed...
NASA Astrophysics Data System (ADS)
He, W.; Beyer, C.; Fleckenstein, J. H.; Jang, E.; Kolditz, O.; Naumov, D.; Kalbacher, T.
2015-10-01
The open-source scientific software packages OpenGeoSys and IPhreeqc have been coupled to set up and simulate thermo-hydro-mechanical-chemical coupled processes with simultaneous consideration of aqueous geochemical reactions faster and easier on high-performance computers. In combination with the elaborated and extendable chemical database of IPhreeqc, it will be possible to set up a wide range of multiphysics problems with numerous chemical reactions that are known to influence water quality in porous and fractured media. A flexible parallelization scheme using MPI (Message Passing Interface) grouping techniques has been implemented, which allows an optimized allocation of computer resources for the node-wise calculation of chemical reactions on the one hand and the underlying processes such as for groundwater flow or solute transport on the other. This technical paper presents the implementation, verification, and parallelization scheme of the coupling interface, and discusses its performance and precision.
NASA Technical Reports Server (NTRS)
Frank, David R.; Zolensky, M. E.; Le, L.; Weisberg, M. K.; Kimura, M.
2013-01-01
The Stardust Mission returned a large fraction of high-temperature, crystalline material that was radially transported from the inner solar system to the Kuiper Belt [1,2]. The mineralogical diversity found in this single cometary collection points to an even greater number of source materials than most primitive chondrites. In particular, the type II olivine found in Wild 2 includes the three distinct Fe/Mn ratios found in the matrix and chondrules of carbonaceous chondrites (CCs) and unequilibrated ordinary chondrites (UOCs) [3]. We also find that low-Ca pyroxene is quite variable (approximately Fs3-29) and is usually indistinguishable from CC, UOC, and EH3 pyroxene as well. However, occasional olivine and pyroxene compositions are found in Wild 2 that are inconsistent with chondrites. The Stardust track 61 terminal particle (TP) is one such example and is the focus of this study. It s highly reduced forsterite and enstatite is consistent only with that in Aubrites, in which FeO is essentially absent from these phases (less than approximately 0.1 wt.% FeO) [4].
Code of Federal Regulations, 2010 CFR
2010-10-01
... 49 Transportation 4 2010-10-01 2010-10-01 false Track gage. 213.323 Section 213.323 Transportation Other Regulations Relating to Transportation (Continued) FEDERAL RAILROAD ADMINISTRATION, DEPARTMENT OF TRANSPORTATION TRACK SAFETY STANDARDS Train Operations at Track Classes 6 and Higher § 213.323 Track gage....
Code of Federal Regulations, 2011 CFR
2011-10-01
... 49 Transportation 4 2011-10-01 2011-10-01 false Track gage. 213.323 Section 213.323 Transportation Other Regulations Relating to Transportation (Continued) FEDERAL RAILROAD ADMINISTRATION, DEPARTMENT OF TRANSPORTATION TRACK SAFETY STANDARDS Train Operations at Track Classes 6 and Higher § 213.323 Track gage....
NASA Astrophysics Data System (ADS)
Emslie, A. Gordon; Bian, Nicolas H.; Kontar, Eduard
2016-05-01
Motivated by recent RHESSI observations that point to the existence of a mechanism that confines electrons to the coronal parts of flare loops more effectively than Coulomb collisions, we consider the impact of pitch-angle scattering off turbulent magnetic fluctuations on the parallel transport of electrons in flaring coronal loops. It is shown that the presence of such a scattering mechanism in addition to Coulomb collisional scattering can significantly reduce the parallel thermal and electrical conductivities relative to their collisional values. We provide illustrative expressions for the resulting thermoelectric coefficients that relate the thermal flux and electrical current density to the temperature gradient and the applied electric field. We then evaluate the effect of these modified transport coefficients on several items of interest to the modeling of flares, including: the peak flare coronal temperature that can be attained, the post-impulsive-phase cooling time of heated coronal plasma, and the importance of the beam-neutralizing return current on both ambient heating and the energy loss rate of accelerated electrons. We also discuss the ways in which anomalous transport processes have an impact on the required overall energy content of accelerated electrons in solar flares.
NASA Astrophysics Data System (ADS)
le Roux, J. A.
2012-12-01
It has been argued that the acceleration of SEPS at a quasi-parallel CME-driven shock to GeV energies in the corona only occurs if strong wave-excitation by SEPs ahead of the shock reduces the parallel mean free path upstream, thus boosting the rate of diffusive shock acceleration. To investigate this issue, we modeled SEP acceleration at a fast parallel traveling shock in the corona with an existing time-dependent focused transport model. The model has been expanded recently to also feature time-dependent self-excitation and damping of Alfvén waves by SEP anisotropies ahead of the shock based on standard quasi-linear theory. Alfvén wave propagation near the traveling shock is modeled based on standard theory for wave transport in a slowly varying non-uniform plasma medium. Preliminary results will be shown to illustrate the increase in wave power driven by SEP anisotropies upstream, the effect of the shock wave in shortening the wave length and increasing the wave amplitude of Alfvén waves, and the associated acceleration of SEPs by 1st order Fermi acceleration to high energies. The role of the acceleration of the cross-shock solar wind flow, which was found to create a downstream population of shock pre-heated particles which forms an additional source for injection into 1st order Fermi acceleration, will be discussed in terms of how it affects self-excitation of Alfvén waves and the formation of high-energy SEPs by 1st order Fermi acceleration.
NASA Astrophysics Data System (ADS)
Ablyazimov, T. O.; Zyzak, M. V.; Ivanov, V. V.; Kisel, P. I.
2015-05-01
One of the main goals in the Compressed Baryonic Matter (CBM) experiment (GSI, Germany) is to find parameters of charged particle trajectories. An online full event reconstruction is planned to be carried out in this experiment, thus demanding fast algorithms be developed, which make the most of the capabilities of modern CPU and GPU architectures. This paper presents the results of an analysis of the Kalman filter-based track reconstruction for charged particles implemented by using various code parallelization methods. A multicore server located at the Laboratory of Information Technologies, Joint Institute for Nuclear Research (LIT JINR), with two CPU Intel Xeon X5660 processors and a GPU Nvidia GTX 480 video card is used.
Luanjing Guo; Chuan Lu; Hai Huang; Derek R. Gaston
2012-06-01
Systems of multicomponent reactive transport in porous media that are large, highly nonlinear, and tightly coupled due to complex nonlinear reactions and strong solution-media interactions are often described by a system of coupled nonlinear partial differential algebraic equations (PDAEs). A preconditioned Jacobian-Free Newton-Krylov (JFNK) solution approach is applied to solve the PDAEs in a fully coupled, fully implicit manner. The advantage of the JFNK method is that it avoids explicitly computing and storing the Jacobian matrix during Newton nonlinear iterations for computational efficiency considerations. This solution approach is also enhanced by physics-based blocking preconditioning and multigrid algorithm for efficient inversion of preconditioners. Based on the solution approach, we have developed a reactive transport simulator named RAT. Numerical results are presented to demonstrate the efficiency and massive scalability of the simulator for reactive transport problems involving strong solution-mineral interactions and fast kinetics. It has been applied to study the highly nonlinearly coupled reactive transport system of a promising in situ environmental remediation that involves urea hydrolysis and calcium carbonate precipitation.