Search NASA⌕ Search

SEARCH · Search NASA

Results for “particle dynamics simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Collaborative Research: Effects of the magnetic field shear and flow shear on the kinetic physics and particle acceleration in relativistic magnetic reconnection (Final Report)

Magnetic reconnection is a fundamental plasma process that allows rapid changes of magnetic field topology and the conversion of magnetic energy into plasma kinetic energy. There has been a strong surge of interest in relativistic reconnection over the past few years in plasma astrophysics, but the rich physics of collisionless reconnection and its associated particle acceleration in the relativistic regime remain less studied compared to the non-relativistic counterparts. Magnetic reconnection in a realistic system is often accompanied with the magnetic field shear and shear flows. This setting is more general and likely the more common situation where most reconnection takes place. Within the past few decades, researchers had made significant progress in understanding how magnetic reconnection works under these conditions in the non-relativistic regime. However, how the effect of special relativity, among other effects, influences the dynamics of reconnection in the strongly magnetized astrophysical plasmas under these conditions is largely unknown. Our proposed research is to use fully kinetic simulations to understand kinetic physics and particle acceleration in relativistic magnetic reconnection in more realistic configurations and physical conditions. The primary goal of this project is to identify the fundamental effects caused by special relativity under the presence of field shear and flow shear.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Machine-Learned Force Field for Molecular Dynamics Simulations of Nonequilibrium Ammonia Synthesis on Iron Catalysts

Ammonia (NH 3 ) is one of the most important industrial chemicals. The conventional NH 3 synthesis method-the Haber–Bosch process-converts atmospheric nitrogen (N 2 ) into NH 3 using H 2 with an iron (Fe) catalyst. However, this process requires high pressures (100–200 atm) and temperatures (700–800 K) near thermal equilibrium. Recently, Fe-based nanocatalysts have been reported to produce promising NH 3 yields under atmospheric pressures and temperature-modulated nonequilibrium conditions. Understanding the mechanism of nonequilibrium catalysis with programmed temperature variation could help to optimize this fully electrified and less energy-intensive process. Although reactive molecular dynamics (RMD) simulations can be a useful tool to model nonequilibrium catalytic processes, they require the development of accurate force fields (i.e., interatomic potentials). Here, we present a machine-learned (ML) force field within the Deep Potential MD (DPMD) framework, trained using periodic density functional theory (DFT) calculations, to model NH 3 synthesis on Fe catalysts with various surface adsorbates such as *N, *H, *N 2 , *H 2 , *NH, *NH 2 , and *NH 3 . Here, we generated the DFT data from static models of elementary reactions on the most stable (110) surface of body-centered cubic Fe, which then were augmented by data from constant number of particles–volume–temperature (NVT) DFT-MD trajectories at various temperatures. Finally, we utilized the fully optimized ML force field to investigate reaction dynamics at an Fe(110) surface at linearly increasing temperatures using NVT-DPMD simulations. Our simulations indicate that pulsed temperature ramping could prove favorable for NH3 synthesis. For example, we conducted ramping under multiple sets of conditions: (i) from 900 to 1200 K over periods of 0.1–0.3 ns for Fe surfaces precovered with N or NH along with H; and (ii) from 300 to 600 K over 0.1–0.3 ns for Fe surfaces precovered with NH 3 . While our simulations so far are limited to short time scales (very rapid heating), these observations shed light on the mechanism of the high NH 3 synthesis rate achieved in a novel temperature-modulated nonequilibrium catalytic reactor using pulsed heating and cooling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Topological phase transition without single particle gap closing in strongly correlated systems

Here, in this study, we show two models where changing topology does not necessarily close the bulk insulating charge gap as demanded in the standard noninteracting picture. From extensive determinantal and dynamical cluster quantum Monte Carlo simulations of the half-filled and quarter-filled Kane-Mele-Hubbard model, we show that, for sufficiently strong interactions at either half- or quarter-filling, a transition between topological and trivial insulators occurs without the closing of a charge gap. To shed light on this behavior, we illustrate that an exactly solvable model reveals that while the single-particle gap remains, the many-body gap does, in fact, close. These two gaps are the same in the noninteracting system but depart from each other as the interaction turns on. We purport that for interacting systems, the proper probe of topological phase transitions is the closing of the many-body rather than the single-particle gap.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Experimental and computational investigations of ethane and ethylene kinetics with copper oxide particles for Chemical Looping Combustion

In this work, reaction pathways for the oxidation of methane, ethane, and ethylene with CuO was obtained by ReaxFF Molecular Dynamics (MD) simulations between temperatures of 1000 K and 2000 K. Experiments in a fixed-bed flow reactor were preformed with methane, ethane, and ethylene at temperatures ranging from 500 K to 1000 K with time-dependent species measurements from an Electron-Ionization Molecular Beam Mass Spectrometer (MBMS), and species validation with Gas Chromatography (GC) for detection of complete and intermediate combustion products. The MBMS and GC allow for the detection of oxygenated species and larger species produced from radical reformation. The simulation and experiment agree on the production of such species as CH 3 CHO, CH 2 O, CO, and H 2 O, which allow for the creation of simple C1 and C2 reaction pathways, which can be used in kinetic models of C2 species and larger fuels such as biofuels, which inherently depend on C1 and C2 kinetics and reaction pathway. The simulation and experiment disagree on the formation of C 2 H 2 , CH 3 OH, and CO 2 with large amounts of C 2 H 2 being measured in the ethylene oxidation simulations and CH 3 OH being formed in methane oxidation simulations, while neither species were experimentally found. In the case of CO 2 large amounts of CO 2 are rapidly produced in experiments with C2 fuels at 800 K, while little-to-no CO 2 was observed in simulations. This is believed to be resulting from the extremely short timescale of the simulations, preventing total oxidation of the fuel. Here, the differences in products produced between simulation and experiment allow for the potential to modify the ReaxFF potential functions to more accurately model the experimental products of Cu–H–O–C reaction kinetics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Influence of defects on the femtosecond laser damage resistance of multilayer dielectric gratings

Multilayer dielectric (MLD) gratings with high diffraction efficiency and a high laser-induced damage (LID) threshold for pulse compressors are key to scaling the peak and average power of chirped pulse amplification lasers. However, surface defects introduced by manufacturing, storage, and handling processes can reduce the LID resistance of MLD gratings and impact the laser output. The underlying mechanisms of such defect-initiated LID remain unclear, especially in the femtosecond regime. In this Letter, we model dynamic processes in interactions of a 20-fs near-infrared (NIR) laser pulse and a MLD grating design in the presence of cylindrically symmetrical nodules and particle contaminants and cracks at the surface. Utilizing a dynamic model based on a 2D finite difference in time domain (FDTD) field solver coupled with photoionization, electron collision, and refractive index modification, we study the simulation results for the damage site distribution initiated by defects of various types and sizes and its impact on the LID threshold of the grating design.

Zhang, Simin (ORCID:0000000174948230)↗
DOE OSTI2025

Cosmological Hydrodynamics at Exascale: A Trillion-Particle Leap in Capability

Resolving the most fundamental questions in cosmology requires simulations that match the scale, fidelity, and physical complexity demanded by next-generation sky surveys. To achieve the realism needed for this critical scientific partnership, detailed gas dynamics must be treated self-consistently with gravity for end-to-end modeling of structure formation. Exascale computing enables simulations that span survey-scale volumes while incorporating key astrophysical processes that shape complex cosmic structures. We present results from CRK-HACC, a cosmological hydrodynamics code built for extreme scalability. Using separation-of-scale techniques, GPU-resident tree solvers, in situ analysis pipelines, and multi-tiered I/O, CRK-HACCexecuted Frontier-E: a four trillion particle full-sky simulation, over an order of magnitude larger than previous efforts. The run achieved 513.1 PFLOPs peak performance, processing 46.6 billion particles per second and writing more than 100 PB of data in just over one week of runtime. Frontier-E marks a significant advance in predictive modeling for next-generation cosmological science.

Frontiere, Nicholas [Argonne National Laboratory (↗
DOE OSTI2020

Fission Product Diffusion in Silicon Carbide: As Revealed by Computer Simulations and Experimental Measurements

The diffusion and release of fission products through silicon carbide in the tri-structural isotropic (TRISO) fuel particles remain unsolved for decades. The underlying mechanism is quite challenging to be determined. To help unveil the mysterious story, the current work applies molecular dynamics method to show the stability of silver, palladium, ruthenium and iodine as an interstitial and their atomic diffusion along coincident site lattice (CSL) boundary, especially ?3 grain boundary (GB). The major finding presents a much faster diffusion along GB than in bulk for all elements considered. The reasonably close estimate to experiments and simulations where available has confirmed the important role of grain boundary diffusion of Ag and Pd in SiC. However, the discrepancy addressed in Ag with measurements from fuel studies suggest a more complicated mechanism, which might be in correlation with high energy grain boundaries or the presence of crack. The subsequent characterization of Ru and I distribution in SiC-PyC-SiC diffusion couples, which have been ion irradiated at 900?C to 10 dpa and 20 dpa, has performed by secondary ion mass spectrometry (SIMS) analysis. The experiment measurements correlate well with the grain boundary diffusion by simulation, which provide further evidences that the grain boundary diffusion cannot be neglected once the fission products are accessible at the grain boundary.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Cluster-dynamics-based parameterization for sulfuric acid–dimethylamine nucleation: comparison and selection through box and three-dimensional modeling

Clustering of gaseous sulfuric acid (SA) enhanced by dimethylamine (DMA) is a major mechanism for new particle formation (NPF) in polluted atmospheres. However, uncertainty remains regarding the SA–DMA nucleation parameterization that reasonably represents cluster dynamics and is applicable across various atmospheric conditions. This uncertainty hinders accurate three-dimensional (3-D) modeling of NPF and the subsequent assessment of its environmental and climatic impacts. Here we extensively compare different cluster-dynamics-based parameterizations for SA–DMA nucleation and identify the most reliable one through a combination of box model simulations, 3-D modeling, and in situ observations. Results show that the parameterization derived from Atmospheric Cluster Dynamic Code (ACDC) simulations, incorporating the latest theoretical insights (DLPNO-CCSD(T)/aug-cc-pVTZ//ωB97X-D/6-311++G(3df,3pd) level of theory) and adequate representation of cluster dynamics, exhibits dependable performance in 3-D NPF simulation for both winter and summer conditions in Beijing and shows promise for application in diverse atmospheric conditions. Another ACDC-derived parameterization, replacing the level of theory with RI-CC2/aug-cc-pV(T+d)Z//M06-2X/6–311++G(3df,3pd), also performs well in NPF modeling at relatively low temperatures around 280 K but exhibits limitations at higher temperatures due to inappropriate representation of SA–DMA cluster thermodynamics. Additionally, a previously reported parameterization incorporating simplifications is applicable for simulating NPF in polluted atmospheres but tends to overestimate particle formation rates under conditions of elevated temperature (>∼300 K) and low-condensation sink (<∼3×10 -3 s −1 ). Our findings highlight the applicability of the new ACDC-derived parameterization, which couples the latest SA–DMA nucleation theory and holistic cluster dynamics, in 3-D NPF modeling. The ACDC-derived parameterization framework provides a valuable reference for developing parameterizations for other nucleation systems.

Shen, Jiewen↗
NASA NTRS1990

Simulations of gas clouds in interacting galaxies

A companion can induce a variety of morphological changes in a galaxy. The author uses N-body simulations to study the effects of different kinds of perturbations on the dynamics of a disk galaxy. The model is two-dimensional, with a disk consisting of about 60,000 particles. Most of the particles (80%) represent the old stellar population with a high velocity dispersion, while the rest (20%) represent gas clouds with a low velocity dispersion. Initially, the velocity dispersion corresponds to Q = 1 for the star particles, and Q = O for the gas particles, where Q is Toomre's (1964) stability parameter. The gas clouds can collide inelastically. The disk is stabilized by a rigid halo potential, and by the random motions of the old star particles. To simulate the effect of an encounter on the disk, a companion galaxy, modelled as a point mass, can move in a co-planar orbit around the disk. A complete description of the N-body code is found in Thomasson (1989). The spiral structures caused by a companion in first a direct and then a retrograde (with respect to the rotation of the disk) parabolic orbit are presented. The associated velocity fields suggest a way to observationally distinguish between leading and trailing spiral arms. The stability of the gas component in a disk in which tidally triggered infall of gas to the center occurs is studied. Finally, the author shows how a ring of gas can form in a disk as a result of a co-planar encounter with another galaxy.

Thomasson, Magnus↗
DOE OSTI2025

Nyström type exponential integrators for strongly magnetized charged particle dynamics

Solving for charged particle motion in electromagnetic fields (i.e. the particle pushing problem) is a computationally intensive component of particle-in-cell (PIC) methods for plasma physics simulations. This task is especially challenging when the plasma is strongly magnetized due numerical stiffness arising from the wide range of time scales between highly oscillatory gyromotion and long term macroscopic behavior. A promising approach to solve these problems is by a class of methods known as exponential integrators that can solve linear problems exactly and are A-stable. This work extends the standard exponential integration framework to derive Nyström-type exponential integrators that integrates the Newtonian equations of motion as a second-order differential equation directly. In particular, we derive second-order and third-order Nyström-type exponential integrators for strongly magnetized particle pushing problems. Numerical experiments show that the Nyström-type exponential integrators exhibit significant improvement in computation speed over the standard exponential integrators.

general physics↗
DOE OSTI2025

Impact of the delayed energy deposition in TREAT on transient dynamics

Estimating the energy released during a reactivity insertion accident (RIA) is of prime importance. TREAT is an experimental reactor designed to simulate these transients. Typically, models of these transients assume instantaneous energy deposition. However, a delayed energy deposition component exists. This work leverages the advanced capabilities of ALFRED, a high-fidelity Monte Carlo code for particle transport, to study the impact of this delayed component on RIA transient dynamics in TREAT. The study reveals that the delayed component can affect the energy released by up to $3\%$ for natural pulse transients. Understanding the impact of the delayed energy deposition component will enable the development of more accurate models, which will assist in designing experiments in TREAT and evaluating their safety.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Comparison of multiobjective optimization methods for the $\mathrm{LCLS-II}$ photoinjector

Particle accelerators are among some of the largest science experiments in the world and can consist of thousands of components with a wide variety of input ranges. These systems can easily become unwieldy optimization problems during design and operations studies. Starting in the early 2000s, searching for better beam dynamics configurations became synonymous with heuristic optimization methods in the accelerator physics community. Genetic algorithms and particle swarm optimization are currently the most widely used. These algorithms can take thousands of simulation evaluations to find optimal solutions for one machine prototype. For large facilities such as the Linac Coherent Light Source (LCLS) and others, this equates to a limited exploration of many possible design configurations. In this paper, the LCLS-II photoinjector is optimized with three optimization algorithms. All optimizations were started from both a uniform random and Latin hypercube sample. In all cases, the optimizations started from Latin hypercube samples outperformed optimizations started from uniform samples. All three algorithms were able to optimize the photoinjector, with the model-based methods approximating the Pareto front in fewer simulation evaluations. This work, in combination with previous optimization observations, indicates objective penalties have a strong impact on the efficiency of such methods. In general, we recommend heuristic methods for initial optimizations and model-based methods when information about the objective space is available.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Structures of the mycobacterial membrane protein MmpL3 reveal its mechanism of lipid transport

The mycobacterial membrane protein large 3 (MmpL3) transporter is essential and required for shuttling the lipid trehalose monomycolate (TMM), a precursor of mycolic acid (MA)-containing trehalose dimycolate (TDM) and mycolyl arabinogalactan peptidoglycan (mAGP), in Mycobacterium species, including Mycobacterium tuberculosis and Mycobacterium smegmatis. However, the mechanism that MmpL3 uses to facilitate the transport of fatty acids and lipidic elements to the mycobacterial cell wall remains elusive. Here, we report 7 structures of the M. smegmatis MmpL3 transporter in its unbound state and in complex with trehalose 6-decanoate (T6D) or TMM using single-particle cryo-electron microscopy (cryo-EM) and X-ray crystallography. Combined with calculated results from molecular dynamics (MD) and target MD simulations, we reveal a lipid transport mechanism that involves a coupled movement of the periplasmic domain and transmembrane helices of the MmpL3 transporter that facilitates the shuttling of lipids to the mycobacterial cell wall.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Cesium Removal and Hydraulic Performance Comparisons of CST Ion Exchange Media Batches in Support of the Tank Closure Cesium Removal Project - 20522

Three batches of Crystalline Silicotitanate (CST) ion exchange media have been evaluated for cesium removal efficiency from a Savannah River Site (SRS) High Level Waste (HLW) supernate simulant in support of the Tank Closure Cesium Removal (TCCR) Project, which involves at-tank column waste treatment. Results will be utilized to select the CST media batch to be added to the next TCCR unit columns. Two of the CST media samples were recent production batches and the third was an archived batch studied extensively at SRS nearly two decades ago. The archived CST batch had a significantly smaller average particle diameter than more recent batches. Small-scale (∼23 mL), side-by-side column tests were conducted at 35 deg. C with each of the CST batches using the same simulant at flow rates near 1.2 mL/min (3.0-3.4 CST bed volumes/hr). The cesium ion exchange column performance results were lower than expected for all three CST batches, but the data indicated that 30% more simulant volume and 40% more equivalent CST bed volumes of solution can be processed with the older CST batch prior to reaching the 50% cesium breakthrough point than can be processed with the newer CST batches. Hydraulic evaluations of the CST columns were also conducted in simulant at 24 deg. C to determine the impact of the different particle size distributions and associated bed porosities on frictional pressure drop under dynamic flow conditions. Higher differential pressure drops were observed with the archived CST batch relative to recent production batches, due to the smaller particle size distribution and the packing characteristics of the archived batch. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Direct numerical simulation of turbulent dispersion of evaporative aerosol clouds produced by an intense expiratory event

Airborne particles are a major route for transmission of COVID-19 and many other infectious diseases. When a person talks, sings, coughs, or sneezes, nasal and throat secretions are spewed into the air. After a short initial fragmentation stage, the expelled material is mostly composed of spherical particles of different sizes. While the dynamics of the largest droplets are dominated by gravitational effects, the smaller aerosol particles, mostly transported by means of hydrodynamic drag, form clouds that can remain afloat for long times. In subsaturated air environments, the dependence of pathogen-laden particle dispersion on their size is complicated due to evaporation of the aqueous fraction. Particle dynamics can significantly change when ambient conditions favor rapid evaporation rates that result in a transition from buoyancyto-drag dominated dispersion regimes. To investigate the effect of particle size and evaporation on pathogen-laden cloud evolution, a direct numerical simulation of a mild cough was coupled with an evaporative Lagrangian particle advection model. The results suggest that while the dispersion of cough particles in the tails of the size distribution are unlikely to be disrupted by evaporative effects, preferential aerosol diameters (30-40 μm) may exhibit significant increases in the residence time and horizontal range under typical ambient conditions. Using estimations of the viral concentration in the spewed fluid and the number of ejected particles in a typical respiratory event, we obtained a map of viral load per volume of air at the end of the cough and the number of virus copies per inhalation in the emitter vicinity.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Saturation of fishbone instability through zonal flows driven by energetic particle transport in tokamak plasmas

Gyrokinetic and kinetic-MHD simulations are performed for the fishbone instability in the DIII-D discharge #178631, chosen for validation of first-principles simulations to predict the energetic particle (EP) transport in an ITER prefusion baseline scenario. Fishbone modes are found to generate zonal flows, which dominate the fishbone saturation. The underlying mechanisms of the two-way fishbone-zonal flows nonlinear interplay are discussed in details. Numerical and analytical analyses identify the fishbone-induced EP redistribution as the dominant generation mechanism for zonal flows. The zonal flows modify the nonlinear dynamics of phase space zonal structures, which reduces the amount of EPs able to resonate with the mode, leading to an early fishbone saturation. Simulation results including zonal flows agree quantitatively with DIII-D experimental measurements of the fishbone saturation amplitude and EP transport, supporting this novel saturation mechanism by self- generated zonal flows. Moreover, the wave-particle mode-locking mechanism is shown to determine quantitatively the fishbone frequency down-chirping, as evident in GTC simulation results in agreement with predictions from analytical theory. Finally, the&#xD;fishbone-induced zonal flows are possibly responsible for the formation of an ion-ITB in the DIII-D discharge. Based on the low EP transport and the large zonal flow shearing rates associated with the fishbone instability in gyrokinetic simulations of the ITER scenario, it is conjectured that high performance scenarios could be designed in ITER burning plasmas through fishbone-induced ITBs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2013

A Particle-in-Cell Simulation for the Traveling Wave Direct Energy Converter (TWDEC) for Fusion Propulsion

A Particle-in-cell simulation model has been developed to study the physics of the Traveling Wave Direct Energy Converter (TWDEC) applied to the conversion of charged fusion products into electricity. In this model the availability of a beam of collimated fusion products is assumed; the simulation is focused on the conversion of the beam kinetic energy into alternating current (AC) electric power. The model is electrostatic, as the electro-dynamics of the relatively slow ions can be treated in the quasistatic approximation. A two-dimensional, axisymmetric (radial-axial coordinates) geometry is considered. Ion beam particles are injected on one end and travel along the axis through ring-shaped electrodes with externally applied time-varying voltages, thus modulating the beam by forming a sinusoidal pattern in the beam density. Further downstream, the modulated beam passes through another set of ring electrodes, now electrically oating. The modulated beam induces a time alternating potential di erence between adjacent electrodes. Power can be drawn from the electrodes by connecting a resistive load. As energy is dissipated in the load, a corresponding drop in beam energy is measured. The simulation encapsulates the TWDEC process by reproducing the time-dependent transfer of energy and the particle deceleration due to the electric eld phase time variations.

Chap, Andrew↗
DOE OSTI2023

An implicit particle code with exact energy and charge conservation for electromagnetic studies of dense plasmas

A collisional particle code based on implicit energy- and charge-conserving methods is presented. A modified version of the particle-suppressed Jacobian-Free Newton-Krylov method that can enhance the solver efficiency is introduced. Mathematically, it is shown that this new approach can be viewed as a fixed-point iteration method for the particle positions. The model can exactly conserve global energy and local charge and can efficiently use time steps larger than the plasma period. In conclusion, the algorithm's ability to simulate dense plasmas accurately and efficiently is quantified by simulating the dynamic compression of a plasma slab via a magnetic piston in 1D planar geometry.

97 MATHEMATICS AND COMPUTING↗