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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The NASA ACTIVATE Mission

The NASA Aerosol Cloud Meteorology Interactions over the Western Atlantic Experiment (ACTIVATE) conducted 162 joint flights with two aircraft over the northwest Atlantic to study aerosol–cloud interactions (ACIs), which represent the largest uncertainty in estimating total anthropogenic radiative forcing. The combination of a high-flying King Air and low-flying HU-25 Falcon, equipped with remote sensing and in situ instruments, characterized trace gases, aerosol particles, clouds, and meteorological variables with data collected nearly simultaneously below, within, and above marine boundary layer (MBL) clouds. Flights spanning warm and cold seasons across 3 years (2020–22) provided a broad range of conditions associated with aerosol particles, cloud properties (including particle size and phase), and meteorology, ideally suited for robust ACI calculations and assessing how well models simulate a wide range of MBL clouds from stratiform to cumulus. ACTIVATE data suggest that drivers of cloud droplet number concentration N d , including aerosol particles and MBL dynamics, vary between winter and summer months with a stronger potential to convert aerosol particles into cloud droplets in winter. Models of varying complexity not only highlight some skills in simulating winter and summer cloud types but also identify challenges that still need to be addressed such as treatment of turbulence, wet scavenging, and mesoscale organization. Remote sensing advances range from new retrieval methods for N d , cloud phase classification, vertically resolved aerosol and cloud condensation nuclei number concentration, and ocean surface wind speed. This work describes these scientific and technological advances along with efforts in outreach and open data science.

aerosol indirect effect↗
DOE OSTI2024

Improving the reliability of machine learned potentials for modeling inhomogeneous liquids

The atomic-scale response of inhomogeneous fluids at interfaces and surrounding solute particles plays a critical role in governing chemical, electrochemical, and biological processes. Classical molecular dynamics simulations have been applied extensively to simulate the response of fluids to inhomogeneities directly, but are limited by the accuracy of the underlying interatomic potentials. Here, we use neural network potentials (NNPs) trained to ab initio simulations to accurately predict the inhomogeneous responses of two distinct fluids: liquid water and molten NaCl. Although NNPs can be readily trained to model complex bulk systems across a range of state points, we show that to appropriately model a fluid's response at an interface, relevant inhomogeneous configurations must be included in the training data. In order to sufficiently sample appropriate configurations of such inhomogeneous fluids, we develop protocols based on molecular dynamics simulations in the presence of external potentials. We demonstrate that NNPs trained on inhomogeneous fluid configurations can more accurately predict several key properties of fluids—including the density response, surface tension and size-dependent cavitation free energies—for liquid water and molten NaCl, compared to both empirical interatomic potentials and NNPs that are not trained on such inhomogeneous configurations. This work therefore provides a first demonstration and framework to extract the response of inhomogeneous fluids from first principles for classical density-functional treatment of fluids free from empirical potentials.

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

Ionization-Induced Subcycle Metallization of Nanoparticles in Few-Cycle Pulses

Strong-field laser-matter interactions in nanoscale targets offer unique avenues for the generation and detailed characterization of matter under extreme conditions. Field-driven, sub-cycle ionization-induced metallization of nanoscale solids in intense laser fields has been predicted (Peltz et al. Time-Resolved X-ray Imaging of Anisotropic Nano-plasma Expansion. Phys. Rev. Lett. 2014, 113, 133401), but its observation was hampered by a lack of a smoking gun. Here, we report the ultrafast metallization of isolated dielectric and semiconducting nanoparticles under intense few-cycle laser pulses. The highest-energy electron emission is found to be a decisive proof that shows a characteristic cutoff modification to a metallic limit for intensities high enough to ignite carrier avalanching in the volume of the particles. Furthermore, semiclassical Mean-field Mie Monte-Carlo transport simulations reveal the underlying dynamics and explain the observed evolution by near-field driven electron backscattering from the metallizing target.

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

Creation and characterization of warm dense matter isochorically heated by an intense laser-driven proton beam to temperatures exceeding 100 eV

The warm dense matter (WDM) is an exotic state of matter encountered in inertial confinement implosions for fusion energy, as well as the interiors of giant planets like Jupiter, brown dwarfs, the atmospheres of white dwarfs, neutron star crusts, and newly discovered exo-planets. One efficient way to create WDM is to use protons accelerated by a high-intensity short-pulse laser to isochorically heat dense samples to WDM states. Despite its importance, direct temperature measurements within WDM targets are scarce. This study utilizes an intense proton beam generated by the kilojoule EP laser further focused and guided by a curved foil and cone structure to efficiently heat a thin copper sample. A high-resolution streaked spectrometer tuned to copper K α fluorescence lines provided bulk temperature measurements every ~2 ps, revealing temperatures exceeding 100 eV in under 50 ps. Particle-in-cell simulations of proton transport and energy deposition closely matched the observed heating dynamics, including transverse temperature gradients revealed by the broadening of K α lines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Adaptive Spatially Aware I/O for Multiresolution Particle Data Layouts

Large-scale simulations on nonuniform particle distributions that evolve over time are widely used in cosmology, molecular dynamics, and engineering. Such data are often saved in an unstructured format that neither preserves spatial locality nor provides metadata for accelerating spatial or attribute subset queries, leading to poor performance of visualization tasks. Furthermore, the parallel I/O strategy used typically writes a file per process or a single shared file, neither of which is portable or scalable across different HPC systems. We present a portable technique for scalable, spatially aware adaptive aggregation that preserves spatial locality in the output. We evaluate our approach on two supercomputers, Stampede2 and Summit, and demonstrate that it outperforms prior approaches at scale, achieving up to 2.5× faster writes and reads for nonuniform distributions. Furthermore, the layout written by our method is directly suitable for visual analytics, supporting low-latency reads and attribute-based filtering with little overhead.

Usher, Will↗
NASA NTRS1991

Interrelation of soft and hard X-ray emissions during solar flares. II - Simulation model

Two-dimensional electrostatic particle simulations are presented which incorporate the effect of quasi-static electric fields on particle dynamics as well as effects associated with wave-particle interactions induced by the accelerated particles. The properties of the soft and hard X-ray and microwave emissions from such systems are examined. In particular, it is shown that acceleration by quasi-static electric fields and heating via wave-particle interactions produces electron distributions with a broken-power law, similar to those inferred from hard X-ray spectra. Also, heating of the ambient plasma gives rise to a region of hot plasma propagating down to the chromosphere at about the ion sound speed.

Winglee, R. M.↗
NASA NTRS2010

Simulation and 'TWINS Observations of the 22 July 2009 Storm

TWINS is the first mission to perform stereo imaging of the Earth's ring current. The magnetic storm on 22 July 2009 is the largest storm observed since TWINS began routine stereo imaging in June 2008. On 22 July 2009, the Dst dropped to nearly -80nT at 7:00 and 10:00 UT. During the main phase and at the peak of the storm, TWINS 1 and 2 were near apogee and moving from pre-dawn to post-dawn local time. The energetic neutral atom (ENA) imagers on the 2 spacecraft captured the storm intensification and the formation of the partial ring current. The peak of the ENA emissions was seen in the midnight-to-dawn local-time sector. The development of this storm has been simulated using the Comprehensive Ring Current Model (CRCM) to understand and interpret the observed signatures. We perform CRCM runs with constant and time-varying magnetic field. The model calculations are validated by comparing the simulated ENA and ion flux intensities with TWINS ENA images and in-situ ion data from THEMIS satellites. Simulation with static magnetic field produces a strong shielding electric field that skews the ion drift trajectories toward dawn. The model's corresponding peak ENA emissions are always eastward than those in the observed TWINS images. On the other hand, simulation with a dynamic magnetic field gives better spatial agreements with both ENA and insitu particle data, suggesting that temporal variations of the geomagnetic field exert a significant influence upon global ring current ion dynamics.

Fok, Mei-Ching↗
NASA NTRSDate not supplied

Sensitivity Study of Dust-Laden Flows in the DLR GBK Facility using the DUST-US3D Framework

An unstructured integrated simulation framework DUST-US3D has been used to study particle-laden flows in the German Aerospace Center (DLR) GBK facility. Air/particle inflow conditions and the overall simulation methodology has been adjusted to ensure consonance with recent experimental work. The resulting particle dynamics are found to be extremely sensitive to the configuration of the particle injection process. However, owing to relatively low particle mass flow rates incorporating 2-way coupling (particles-gas exchange momentum/energy) or particle-particle collisions has no significant impact on predictions. Numerical estimates for particle velocity magnitude are consistently lower than those observed during experiments which could be indicative of particles coalescing into larger, possibly non-spherical aggregates.

Dusty Flows↗
DOE OSTI2025

Nonlinear Alfvén instability simulation and EP transport for ITER reversed shear (steady-state) and monotonic q-profile regimes *

The nonlinear evolution and dynamics of energetic particle (EP) driven instabilities in both steady-state (reversed shear q-profile) and monotonic q-profile regimes of ITER are examined using the global gyro-Landau closure model FAR3d. Both neutral beam and alpha components are included, leading to synergistic effects between the two populations. In order to maintain computational feasibility, the present nonlinear simulation includes all toroidal mode numbers ranging from n = 0 to 15. It cannot be excluded that n’s above this range could introduce larger linear growth rates; however, as indicated from the results, the n = 0 to 15 range is sufficient to provide a strong linear drive and significant nonlinear transport effects. While the toroidal mode numbers n = 6,13 dominate the linear instability growth phase, nonlinear energy transfers in the saturated phase reverse this trend, leading to the dominance of lower n’s in the saturated phase of the simulation. Zonal flow structure generation and EP density flattening are further consequences of the nonlinear phase. The density profile flattening is caused by collectively driven EP radial transport fluxes, which have non-local characteristics. Instantaneous alpha particle transport fluxes are significant (Γ α ~ 2 x 10 20 m –2 s –1 ) for the reversed shear regime and lower (Γ α ~ 1.1 x 10 20 m –2 s –1 ) for the monotonic q-profile case.

AE↗
DOE OSTI2022

Parametric studies of stream instability-induced higher harmonics in plasma ionization breakdown near an emissive dielectric surface

Here, in this work, we comprehensively investigate the generation of higher harmonic (HH) electric fields normal to the applied rf electric field in multipactor-coexisting plasma breakdown by fully kinetic particle-in-cell (PIC) simulations and a theoretical model. Firstly, a base case at driving frequency 1 GHz, transverse rf electric field amplitude 3 MV m –1 , and background gas pressure 0.2 Torr, is studied in detail. The enhanced harmonic frequency observed is around ten times the fundamental rf frequency, significantly lower than the Langmuir frequency. A theoretical model reveals that the fundamental mechanism of HHs generation is stream–plasma instability, which originates from stream-like secondary electron emission interacting with plasma. The resulting HH frequency and the growth rate of its oscillating amplitude from the theoretical model, agree well with the PIC simulations. With increasing pressure, the HH oscillations are found to be significantly reduced. This is because at higher pressure the gas ionization rate is higher, which causes a more rapidly increasing plasma density, leaving less time for the growth of instability. Furthermore, the parameter space in terms of background gas pressure and rf field amplitude within which the HHs appear is revealed. Finally, the effect of the driving rf frequency on HHs is also investigated, and it is found that the instability-induced oscillating HHs field remains when the driving frequency is increased to 2 GHz, however, it is highly reduced at higher driving frequency of 5 GHz, as oscillations at the fundamental frequency start playing a more important role.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2025

Rheology of Lunar Regolith Simulant Under Varying Gravitational Conditions

Understanding structure-property-process relationship aka rheology of regolith in varying gravity conditions is critical for exploration and future missions. In this work, a framework for studying rheology of lunar regolith simulant in varying gravity conditions (Terrestrial, Lunar and Martian) is described theoretically/numerically, where the driving force of flow of simulants was gravitational acceleration. In the analysis particle-particle interaction forces were included. It was found that the dynamic behavior of granular material is extremely sensitive to external conditions by virtue of the myriad of forces present between particle grains. The results obtained were validated against results obtained in earth and Lunar gravity conditions. The theoretical/numerical and experimental results showed that the complex interaction of these forces can drastically change the dynamics of the material, which is not captured by standard design rules, generally used in industry, for variable gravity applications.

Simulation↗
DOE OSTI2021

Noise and error analysis and optimization in particle-based kinetic plasma simulations

In this paper we analyze the noise in macro-particle methods used in plasma physics and fluid dynamics, leading to approaches for minimizing the total error, focusing on electrostatic models in one dimension. We begin by describing kernel density estimation for continuous values of the spatial variable x, expressing the kernel in a form in which its shape and width are represented separately. The covariance matrix of the noise in the density is computed, first for uniform true density. The bandwidth of the covariance matrix C(x,y) is related to the width of the kernel. A feature that stands out is the presence of constant negative terms in the elements of the covariance matrix both on and off-diagonal. These negative correlations are related to the fact that the total number of particles is fixed at each time step; they also lead to the property ∫C(x,y)dy = 0. We investigate the effect of these negative correlations on the electric field computed by Gauss's law, finding that the noise in the electric field is related to a process called the Ornstein-Uhlenbeck bridge, leading to a covariance matrix of the electric field with variance significantly reduced relative to that of a Brownian process. For non-constant density, p(x), still with continuous x, we analyze the total error in the density estimation and discuss it in terms of bias-variance optimization (BVO). For some characteristic length l, determined by the density and its second derivative, and kernel width h, having too few particles within h leads to too much variance; for h that is large relative to l, there is too much smoothing of the density. The optimum between these two limits is found by BVO. For kernels of the same width, it is shown that this optimum (minimum) is weakly sensitive to the kernel shape. Next, we repeat the analysis for x discretized on a grid. In this case the charge deposition rule is determined by a particle shape. An important property to be respected in the discrete system is the exact preservation of total charge on the grid; this property is necessary to ensure that the electric field is equal at both ends, consistent with periodic boundary conditions. We find that if the particle shapes satisfy a partition of unity property, the particle charge deposited on the grid is conserved exactly. Further, if the particle shape is expressed as the convolution of a kernel with another kernel that satisfies the partition of unity, then the particle shape obeys the partition of unity. This property holds for kernels of arbitrary width, including widths that are not integer multiples of the grid spacing. Furthermore, we show results relaxing the approximations used to do BVO optimization analytically, by doing numerical computations of the total error as a function of the kernel width, on a grid in x. The comparison between numerical and analytical results shows good agreement over a range of particle shapes. We discuss the practical implications of our results, including the criteria for design and implementation of computationally efficient particle shapes that take advantage of the developed theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Quantum computing for fusion energy science applications

This is a review of recent research exploring and extending present-day quantum computing capabilities for fusion energy science applications. We begin with a brief tutorial on both ideal and open quantum dynamics, universal quantum computation, and quantum algorithms. Then, we explore the topic of using quantum computers to simulate both linear and nonlinear dynamics in greater detail. Because quantum computers can only efficiently perform linear operations on the quantum state, it is challenging to perform nonlinear operations that are generically required to describe the nonlinear differential equations of interest. In this work, we extend previous results on embedding nonlinear systems within linear systems by explicitly deriving the connection between the Koopman evolution operator, the Perron–Frobenius evolution operator, and the Koopman–von Neumann evolution (KvN) operator. We also explicitly derive the connection between the Koopman and Carleman approaches to embedding. Extension of the KvN framework to the complex-analytic setting relevant to Carleman embedding, and the proof that different choices of complex analytic reproducing kernel Hilbert spaces depend on the choice of Hilbert space metric are covered in the appendixes. Finally, we conclude with a review of recent quantum hardware implementations of algorithms on present-day quantum hardware platforms that may one day be accelerated through Hamiltonian simulation. We discuss the simulation of toy models of wave–particle interactions through the simulation of quantum maps and of wave–wave interactions important in nonlinear plasma dynamics.

Joseph, I. (ORCID:0000000255400840)↗
DOE OSTI2021

A statistical mechanics approach to macroscopic limits of car-following traffic dynamics

Herein we study the derivation of macroscopic traffic models from car-following vehicle dynamics by means of hydrodynamic limits of an Enskog-type kinetic description. We consider the superposition of Follow-the-Leader (FTL) interactions and relaxation towards a traffic-dependent Optimal Velocity (OV) and we show that the resulting macroscopic models depend on the relative frequency between these two microscopic processes. If FTL interactions dominate then one gets an inhomogeneous Aw–Rascle–Zhang model, whose (pseudo) pressure and stability of the uniform flow are precisely defined by some features of the microscopic FTL and OV dynamics. Conversely, if the rate of OV relaxation is comparable to that of FTL interactions then one gets a Lighthill–Whitham–Richards model ruled only by the OV function. We further confirm these findings by means of numerical simulations of the particle system and the macroscopic models. Unlike other formally analogous results, our approach builds the macroscopic models as physical limits of particle dynamics rather than assessing the convergence of microscopic to macroscopic solutions under suitable numerical discretisations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1984

Computer simulation of surface and film processes

All the investigations which were performed employed in one way or another a computer simulation technique based on atomistic level considerations. In general, three types of simulation methods were used for modeling systems with discrete particles that interact via well defined potential functions: molecular dynamics (a general method for solving the classical equations of motion of a model system); Monte Carlo (the use of Markov chain ensemble averaging technique to model equilibrium properties of a system); and molecular statics (provides properties of a system at T = 0 K). The effects of three-body forces on the vibrational frequencies of triatomic cluster were investigated. The multilayer relaxation phenomena for low index planes of an fcc crystal was analyzed also as a function of the three-body interactions. Various surface properties for Si and SiC system were calculated. Results obtained from static simulation calculations for slip formation were presented. The more elaborate molecular dynamics calculations on the propagation of cracks in two-dimensional systems were outlined.

Tiller, W. A.↗
DOE OSTI2021

Asymmetry between deuterium and tritium turbulent particle flows

The asymmetry between deuterium (D) and tritium (T) turbulent particle fluxes in mixed D–T plasmas is studied with numerical simulations of nonlinear gyrokinetic turbulence in ion temperature gradient-dominated and trapped electron mode-dominated regimes. At 50–50 D–T concentration, the asymmetry, or flow separation, between D and T fluxes is such that the tritium is better confined than the deuterium in both regimes. To supplement the nonlinear simulations, an analytic quasilinear theory of the particle flux symmetry breaking is developed and is valid for general electron dynamics. Furthermore, this theory correctly predicts the numerically computed deviation in the ion density fraction (from 50–50), or ion density gradient (from the electron gradient), required to restore equal deuterium and tritium fluxes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2008

Relativistic Particle-In-Cell Simulations of Particle Accleration in Relativistic Jets

Highly accelerated particles are observed in astrophysical systems containing relativistic jets and shocks, e.g., active galactic nuclei (AGNs), microquasars, and Gamma-Ray Bursts (GRBs). Particle-In-Cell (PIC) simulations of relativistic electron-ion and electron-positron jets injected into a stationary medium show that efficient acceleration occurs downstream in the jet. In collisionless relativistic shocks particle acceleration is due to plasma waves and their associated instabilities, e.g., the Buneman instability, other two-stream instabilities, and the Weibel (filamentation) instability. Simulations show that the Weibel instability is responsible for generating and amplifying highly non-uniform, small-scale magnetic fields. The instability depends on strength and direction of the magnetic field. Particles in relativistic jets may be accelerated in a complicated dynamics of relativistic jets with magnetic field. We present results of our recent PIC simulations.

Nishikawa, K.-I.↗