SEARCH · Search NASA

Results for “Particle simulations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Novel data interpretation method for DIII-D divertor retarding field energy analyzer with 3-D particle-in-cell simulations

A novel data interpretation process that utilizes comprehensive particle-in-cell (PIC) simulations is developed for the new retarding field energy analyzer (RFEA) currently being constructed at DIII-D for the lower divertor using the Divertor Material Evaluation System. Furthermore, this probe is expected to survive a heat load of up to 100 MW/m 2 for up to 5 s and reliably measure the main ion temperature (T i ) on the divertor target ranging from 10 to 200 eV. These extreme conditions posed significant engineering limitations on the probe geometry, thus extensive validation work has been performed. The conventional fitting method for the RFEA I–V characteristics is based on a simplified 1-D model without considering the ion space charge inside the probe cavity and may not be sufficient for probes designed for the DIII-D divertor environment. In this article, a more realistic description of the particle propagation process within the RFEA cavity is achieved by including both 3-D geometric effects and ion space charge in the PIC simulations, and the capability to reconstruct the ion energy distribution functions is demonstrated with reasonable consistency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1990

Simulation studies of plasma waves in the electron foreshock - The generation of Langmuir waves by a gentle bump-on-tail electron distribution

Particle simulation experiments were used to study the basic physical ingredients needed for building a global model of foreshock wave phenomena. In particular, the generation of Langmuir waves by a gentle bump-on-tail electron distribution is analyzed. It is shown that, with appropriately designed simulations experiments, quasi-linear theory can be quantitatively verified for parameters corresponding to the electron foreshock.

Dum, C. T.↗

NASA NTRS1991

Simulation of whistler waves excited in the presence of a cold plasma cloud - Implications for the CRRES mission

A one-dimensional electromagnetic particle simulation model is constructed to study the excitation of whistler waves in the presence of a cold plasma cloud for conditions representative of those after the release of lithium in the inner plasma sheet during the Combined Release and Radiation Effect Satellite mission. The results indicate that a standing-wave pattern with discrete wave frequencies is formed within the cloud. The magnetic wave amplitude inside the cloud, which is limited by quasi-linear diffusion, is of the order of several nanoteslas. Assuming a magnetospheric loss cone of 5 deg, the observed pitch angle diffusion produced by the whistler waves is sufficient to put the electrons on strong diffusion.

Pritchett, P. L.↗

NASA NTRSDate not supplied

Comparing Particle Flow Regimes in the L2K Arcjet with Martian Entry Conditions

In this study, the US3D-DUST code will be used to simulate particles in the L2K dusty-flow test environments and compare them to conditions in the shock layer of the Schiaparelli capsule during its entry into the Martian atmosphere in 2016. The particle Reynolds, Mach, and Knudsen numbers, key non-dimensional numbers that characterize the dusty-flow environments, will be compared between the L2K and Martian shock layers. This study will provide a preliminary assessment of how well-suited L2K dusty-flow conditions are for particle drag model validation and the development of particle impact surface damage correlations. The results of this study may help to guide the design of future dusty flow environments in the L2K facility

Dusty flow↗

NASA NTRS1985

A simulation study of kilometric radiation generation along an auroral field line

A model for the high-altitude auroral zone is combined with local particle simulations in order to assess the quantitative implications of relativistic dispersion for the generation of auroral kilometric radiation. The auroral zone model is described, including the distribution functions for the primary and secondary electron populations and the variations of the total number density with altitude. The analysis shows that the most intense cyclotron maser emission should occur in the 1.75-2.0 R(E) altitude range. At lower altitudes, the secondary electrons are dominant, and the instability saturates at a low level. At higher altitudes, the energy of the primary electron is lower, and again the radiation level drops. The simulations indicate that linear growth rates of the order of 0.002 Omega(e) occur in this region and that the conversion of primary electron energy into AKR is about one percent.

Pritchett, P. L.↗

NASA NTRS1986

A simulation study of the critical ionization velocity process

The critical ionization velocity process is studied by first investigating a coupled system of equations describing the production of several ion species and electrons by impact ionization, their collisions with neutrals, and the heating of electrons. Analytic relations derived from this were tested with the help of a particle simulation, including collisional processes between neutrals and plasma particles. It was found that resistive heating of electrons plays an important role when the density of the neutrals is high, and that electron heating due to lower hybrid waves is significant when the neutral density is low. In both cases, the control of the plasma production rate by the ratio of the beam velocity to the critical velocity was verified.

Machida, S.↗

NASA NTRS1993

Low-frequency mobility response functions for the central plasma sheet with application to tearing modes

Consideration is given to the effect of constant cross-tail magnetic field By on the collisionless conductivity produced by chaotic scattering and stochastic diffusion of particles in the current sheet for a parabolic geometry. It is shown that the correlation time scales as (By/Bz)-squared, and from this strong By scaling a strong tendency toward stabilization of the linear tearing modes with increasing values of By is inferred. This effect of increased dawn-dusk mobility is particularly dramatic when electrons are introduced in the calculation, and is in agreement with the results of kinetic particle simulations. The collisionless conductivity is expressed in terms of the ensemble-averaged power spectrum of the single particle trajectories, which makes it possible to calculate directly the linear conductivity instead of deriving it from the calculation of the irreversible heating rates.

Hernandez, J.↗

NASA NTRS1996

Excitation of Plasma Waves in Aurora by Electron Beams

In this paper, we study numerically the excitation of plasma waves by electron beams, in the auroral region above 2000 km of altitude. We have solved the fully kinetic dispersion relation, using numerical method and found the real frequency and the growth rate of the plasma wave modes. We have examined the instability properties of low-frequency waves such as the Electromagnetic Ion Cyclotron (EMIC) wave as well as Lower-Hybrid (LH) wave in the range of high-frequency. In all cases, the source of free energy are electron beams propagating parallel to the geomagnetic field. We present some features of the growth rate modes, when the cold plasma parameters are changed, such as background electrons and ions species (H(+) and O(+)) temperature, density or the electron beam density and/or drift velocity. These results can be used in a test-particle simulation code, to investigate the ion acceleration and their implication in the auroral acceleration processes, by wave-particle interaction.

daSilva, C. E.↗

DOE OSTI2025

SciDAC ISEP: Integrated Simulation of Energetic Particles in Burning Plasmas

The objective of the SciDAC Center for Integrated Simulation of Energetic Particles in Burning Plasmas (ISEP) is to improve physics understanding of energetic particle (EP) confinement and EP interactions with burning thermal plasmas through large-scale simulations. The ISEP center will develop a multiscale and multiphysics ISEP framework for a predictive capability of EP physics and deliver an EP module incorporating both first-principles simulations and high fidelity reduced transport models to the fusion whole device modeling (WDM) project. The ISEP framework will enable us to perform long time, global kinetic simulations of EP physics in burning plasmas, by utilizing the full power of the next generation supercomputers. Our research and development activities will build on fruitful collaborations with computer scientists and applied mathematicians to offer enabling technologies for performance scalability, portability, solvers, coupling for integration with the fusion WDM project, and long-term preservation of data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Simulations of neutron unbound physics for Geant4

The study of neutron unbound systems via the invariant mass technique is the primary focus of the MoNA Collaboration, which built and operates the Modular Neutron Array (MoNA) and the Large multi-Institutional Scintillator Array (LISA) at FRIB. Advancements in nuclear structure from theory and experiment along the neutron dripline have presented opportunities to understand the nature of unbound systems in higher mass nuclei. The GEometry ANd Tracking (Geant4) platform has been used in high-energy and nuclear physics to simulate particle interactions with as much detail as the user desires. Geant4 currently does not have a physics class to simulate neutron unbound systems. Given the advancement of accelerator facilities and active searches along the neutron dripline, detailed simulations to study the breakup of neutron unbound systems, are necessary. Furthermore, the implementation of the breakup of neutron unbound systems in Geant4 will be presented.

Neutron dripline↗

DOE OSTI2019

High-Frequency Plasma Waves and Pitch Angle Scattering Induced by Pulsed Electron Beams

Cherenkov radiation from a pulse of charge propagating along the magnetic field in a magnetized plasma is analyzed using theory and fluid–kinetic simulations. Besides radiation into whistler modes, the subject of many previous investigations in laboratory and space, radiation can occur through extraordinary (X) modes. Theory and simulations demonstrate that X mode radiation efficiencies can be orders of magnitude higher than those into whistler modes. Test particle simulations of the dynamics of energetic electrons in the beam–generated wavefield show that X modes can also induce pitch angle scattering much more efficiently than whistlers. While coherence effects associated with spreading of realistic beam pulses may limit the size of the X mode source region, a simple model of beam dynamics suggests that the size of this region could be substantial (hundreds of meters for ionospheric conditions). Furthermore, these results have potentially important implications for many problems, including understanding losses in the near–Earth environment and radiation belt remediation.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1983

Recent advances in numerical simulation of space-plasma-physics problems

Computer simulations have become an increasingly popular, important and insightful tool for studying space plasmas. This review describes MHD and particle simulations, both of which treat the plasma and the electromagnetic field in which it moves in a self consistent fashion but on drastically different spatial and temporal scales. The complementary roles of simulation, observations and theory are stressed. Several examples of simulations being carried out in the area of magnetospheric plasma physics are described to illustrate the power, potential and limitations of the approach.

Birmingham, T. J.↗

DOE OSTI2025

Generating multi-scale Li-ion battery cathode particles with radial grain architectures using stereological generative adversarial networks

Abstract Understanding structure-property relationships of Li-ion battery cathodes is crucial for optimizing rate-performance and cycle-life resilience. However, correlating the morphology of cathode particles, such as in LiNi0.8Mn0.1Co0.1O2 (NMC811), and their inner grain architecture with electrode performance is challenging, particularly, due to the significant length-scale difference between grain and particle sizes. Experimentally, it is not feasible to image such a high number of particles with full granular detail. A second challenge is that sufficiently high-resolution 3D imaging techniques remain expensive and are sparsely available at research institutions. Here, we present a stereological generative adversarial network-based model fitting approach to tackle this, that generates representative 3D information from 2D data, enabling characterization of materials in 3D using cost-effective 2D data. Once calibrated, this multi-scale model can rapidly generate virtual cathode particles that are statistically similar to experimental data, and thus is suitable for virtual characterization and materials testing through numerical simulations. A large dataset of simulated particles with inner grain architecture has been made publicly available.

25 ENERGY STORAGE↗

NASA NTRS1987

On the generation of broadband electrostatic noise

Linear theory analysis and particle simulation studies were carried out to understand the mechanisms generating broadband electrostatic noise (BEN) and the effects of BEN on particles as well as the nonlinear saturation level in the geomagnetic tail region. Streaming ion beams as well as warm ions and electrons are used to model the plasma-sheet particle population. When the beam ion temperature is comparable to the temperature of the warm plasma-sheet particles, electrostatic ion cyclotron instabilities become unstable, giving rise to low-frequency noise at omega less than about n(Omega i), where Omega i is the ion gyrofrequency.

Ashour-Abdalla, M.↗

NASA NTRS1993

Current loop coalescence studied by 3-D electromagnetic particle code

Solar flare plasma data from the Yohkoh satellite is analyzed. The interactions of current loops were observed in the active regions on the Sun. This observation pointed out the importance of the idea that the solar flare is generated by the coalescence of current loops. The three dimensional electromagnetic particle simulations are to help in understanding the global interaction between two current loops including the evolution of the twist of loops due to instabilities. Associated rapid dynamics of current loop coalescence such as reconnection, shock waves and associated kinetic processes such as energy transfer, acceleration of particles, and electromagnetic emissions are to be studied by the code to complement analytical theories and magnetohydrodynamic simulations of the current loop coalescence. The simulation results show the strong interactions between two current loops, beam and whistler instabilities, and associated parallel and perpendicular particle heating.

Nishikawa, Ken-Ichi↗

DOE OSTI2022

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗

DOE OSTI2021

Computing the Paschen curve for argon with speed-limited particle-in-cell simulation

Upon inclusion of collisions, the speed-limited particle-in-cell (SLPIC) simulation method successfully computed the Paschen curve for argon. The 1D3V simulations modeled an electron cascade across an argon-filled capacitor, including electron-neutral ionization, electron-neutral elastic collisions, electron-neutral excitation, and ion-induced secondary electron emission. In electrical breakdown, the timescale difference between ion and electron motion makes traditional PIC methods computationally slow. To decrease this timescale difference and speed up computation, here we used SLPIC, a time-domain algorithm that limits the speed of the fastest electrons in the simulation. The SLPIC algorithm facilitates a straightforward, fully kinetic treatment of dynamics and collisions. SLPIC was as accurate as PIC, but ran up to 200 times faster. SLPIC accurately computed the Paschen curve for argon over three orders of magnitude in pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Moment-preserving Monte-Carlo Coulomb collision method for particle codes

Binary-pairing Monte-Carlo methods are widely used in particle-in-cell codes to capture effects of small angle Coulomb collisions. These methods preserve momentum and energy exactly when the simulation particles have equal weights. However, when the interacting particles are of varying weight, these physical conservation laws are only preserved on average. Here, we 1) extend these methods to weighted particles such that the scattering physics is correct on average, and 2) describe a new method for adjusting the particle velocities post scatter to restore exact conservation of momentum and energy. In conclusion, the efficacy of the model is illustrated with various test problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗