Search NASA⌕ Search

SEARCH · Search NASA

Results for “kinetic 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Kinetic simulation of magnetic field generation and collisionless shock formation in expanding laboratory plasmas

Recent laboratory experiments with laser-produced plasmas have observed and studied a number of fundamental physical processes relevant to magnetized astrophysical plasmas, including magnetic reconnection, collisionless shocks, and magnetic field generation by Weibel instability, opening up new experimental platforms for laboratory astrophysics. We develop a fully kinetic simulation model for first-principles simulation of these systems including the dynamics of magnetic fields—magnetic field generation by the Biermann battery effect or Weibel instability; advection by the ion flow, Hall effect, and Nernst effect; and destruction of the field by dissipative mechanisms. Key dimensionless parameters describing the system are derived for scaling between kinetic simulation, recent experiments, and astrophysical plasmas. First, simulations are presented which model Biermann battery magnetic field generation in plasmas expanding from a thin target. Ablation of two neighboring plumes leads to the formation of a current sheet as the opposing Biermann-generated fields collide, modeling recent laser-driven magnetic reconnection experiments. Second, we simulate recent experiments on collisionless magnetized shock generation, by expanding a piston plasma into a pre-magnetized ambient plasma. For parameters considered, the Biermann effect generates additional magnetic fields in the curved shock front and thereby increases shock particle reflection. In conclusion, both cases show the importance of kinetic processes in the interaction of plasmas with magnetic fields and open opportunities to benchmark these important processes through comparison of theory and experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Experimental data based combinatorial kinetic simulations for predictions of synergistic catalyst mixtures

Physical mixtures of catalysts can sometimes show synergistic activity which exceeds the sum of the activities of the same amount of catalysts employed separately. However, identifying such synergistic mixtures is non-trivial. Not all mixtures display synergy, and the number of combinations that are possible (even for binary mixtures of equal portions) scales very rapidly, and thus they would be costly to screen experimentally. In this work, we show that it is possible to predict synergistic mixtures using combinatorial kinetic simulations based on experimental data collected on individual catalysts. The data was first collected for conditions-of-interest under low conversions (also called near differential conditions, such that each condition approximates a small volume in a non-differential conditions reactor) to build a library of kinetic models (one model for each catalyst). This data was then used for combinatorial kinetic simulations of non-differential reactor scale conversions, and successfully predicted the qualitative behavior of two synergistic physical mixtures. The capability was utilized in the present work in the context of converting CO, C 3 H 6 , and NO species in simulated car exhaust, but is a general approach. We provide equations to estimate the costs of the screening method as well as the combinatorial kinetics method and show that the costs of identifying synergistic physical mixtures become much lower with the combinatorial kinetics method even before ten catalysts of interest. The cost of using prediction by the combinatorial kinetics method continue to become lower (relative to screening) when more catalysts are added -- and even more so if mixtures beyond two components, or consideration of multiple catalyst ratios, is of interest.

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

A review of low-rank methods for time-dependent kinetic simulations

Time-dependent kinetic models are ubiquitous in computational science and engineering. The underlying integro-differential equations in these models are high-dimensional, comprised of a six–dimensional phase space, making simulations of such phenomena extremely expensive. In this article we demonstrate that in many situations, the solution to kinetics problems lives on a low dimensional manifold that can be described by a low-rank matrix or tensor approximation. We then review the recent development of so-called low-rank methods that evolve the solution on this manifold. The two classes of methods we review are the dynamical low-rank (DLR) method, which derives differential equations for the low-rank factors, and a Step-and-Truncate (SAT) approach, which projects the solution onto the low-rank representation after each time step. Thorough discussions of time integrators, tensor decompositions, and method properties such as structure preservation and computational efficiency are included. We further show examples of low-rank methods as applied to particle transport and plasma dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Pilgrim : A thermal rate constant calculator and a chemical kinetics simulator

Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Explicit structure-preserving geometric particle-in-cell algorithm in curvilinear orthogonal coordinate systems and its applications to whole-device 6D kinetic simulations of tokamak physics

Explicit structure-preserving geometric particle-in-cell (PIC) algorithm in curvilinear orthogonal coordinate systems is developed. The work reported represents a further development of the structure-preserving geometric PIC algorithm achieving the goal of practical applications in magnetic fusion research. The algorithm is constructed by discretizing the field theory for the system of charged particles and electromagnetic field using Whitney forms, discrete exterior calculus, and explicit non-canonical symplectic integration. In addition to the truncated infinitely dimensional symplectic structure, the algorithm preserves exactly many important physical symmetries and conservation laws, such as local energy conservation, gauge symmetry and the corresponding local charge conservation. As a result, the algorithm possesses the long-term accuracy and fidelity required for first-principles-based simulations of the multiscale tokamak physics. The algorithm has been implemented in the SymPIC code, which is designed for high-efficiency massively-parallel PIC simulations in modern clusters. The code has been applied to carry out whole-device 6D kinetic simulation studies of tokamak physics. A self-consistent kinetic steady state for fusion plasma in the tokamak geometry is numerically found with a predominately diagonal and anisotropic pressure tensor. The state also admits a steady-state sub-sonic ion flow in the range of 10 km s -1 , agreeing with experimental observations and analytical calculations Kinetic ballooning instability in the self-consistent kinetic steady state is simulated. It is shown that high-n ballooning modes have larger growth rates than low-n global modes, and in the nonlinear phase the modes saturate approximately in 5 ion transit times at the 2% level by the E × B flow generated by the instability. These results are consistent with early and recent electromagnetic gyrokinetic simulations.

43 PARTICLE ACCELERATORS↗
NASA NTRS2010

Jovian Plasma Torus Interaction with Europa: 3D Hybrid Kinetic Simulation. First results

The hybrid kinetic model supports comprehensive simulation of the interaction between different spatial and energetic elements of the Europa-moon-magnetosphere system with respect to variable upstream magnetic field and flux or density distributions of plasma and energetic ions, electrons, and neutral atoms. This capability is critical for improving the interpretation of the existing Europa flyby measurements from the Galileo orbiter mission, and for planning flyby and orbital measurements, (including the surface and atmospheric compositions) for future missions. The simulations are based on recent models of the atmosphere of Europa (Cassidy etal.,2007;Shematovichetal.,2005). In contrast to previous approaches with MHD simulations, the hybrid model allows us to fully take into account the finite gyro radius effect and electron pressure, and to correctly estimate the ion velocity distribution and the fluxes along the magnetic field (assuming an initial Maxwellian velocity distribution for upstream background ions).Non-thermal distributions of upstream plasma will be addressed in future work. Photoionization,electron-impact ionization, charge exchange and collisions between the ions and neutrals are also included in our model. We consider two models for background plasma:(a) with O(++) ions; (b) with O(++) and S(++) ions. The majority of O2 atmosphere is thermal with an extended cold population (Cassidyetal.,2007). A few first simulations already include an induced magnetic dipole; however, several important effects of induced magnetic fields arising from oceanic shell conductivity will be addressed in later work.

Lipatov, A. S.↗
DOE OSTI2015

Kinetic simulations of X-B and O-X-B mode conversion

We have performed fully-kinetic simulations of X-B and O-X-B mode conversion in one and two dimensional setups using the PIC code EPOCH. We have recovered the linear dispersion relation for electron Bernstein waves by employing relatively low amplitude incoming waves. The setups presented here can be used to study non-linear regimes of X-B and O-X-B mode conversion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Kinetic simulations underestimate the effects of waves during magnetic reconnection

Collisionless plasma systems are often studied using fully kinetic simulations, where protons and electrons are treated as particles. Due to their computational expense, it is necessary to reduce the ion-to-electron mass ratio $m_i/m_e$ or the ratio between plasma and cyclotron frequencies in simulations of large systems. In this Letter we show that when electron-scale waves are present in larger-scale systems, numerical parameters affect their amplitudes and effects on the larger system. Using lower-hybrid drift waves during magnetic reconnection as an example, we find that the ratio between the wave electric field and the reconnection electric field scales as $\sqrt{m_i/m_e}$, while the phase relationship is also affected. The combination of these effects means that the anomalous drag that contributes to momentum balance in the reconnection region can be underestimated by an order of magnitude. The results are relevant to the coupling of electron-scale waves to ion-scale reconnection regions, and other systems such as collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Multispecies Ion Acceleration in 3D Magnetic Reconnection with Hybrid-Kinetic Simulations

Magnetic reconnection drives multispecies particle acceleration broadly in space and astrophysics. Here we perform the first 3D hybrid simulations (fluid electrons, kinetic ions) that contain sufficient scale separation to produce nonthermal heavy-ion acceleration, with fragmented flux ropes critical for accelerating all species. We demonstrate the acceleration of all ion species (up to Fe) into power-law spectra with similar indices, by a common Fermi acceleration mechanism. The upstream ion velocities influence the first Fermi reflection for injection. The subsequent onsets of Fermi acceleration are delayed for ions with lower charge-mass ratios (Q/M), until growing flux ropes magnetize them. This leads to a species-dependent maximum energy/nucleon ∝(Q/M) α . These findings are consistent with in situ observations in reconnection regions, suggesting Fermi acceleration as the dominant multispecies ion acceleration mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Particle Acceleration in Kinetic Simulations of Nonrelativistic Magnetic Reconnection with Different Ion–Electron Mass Ratios

By means of fully kinetic particle-in-cell simulations, we study whether the proton-to-electron mass ratio m i /m e influences the energy spectrum and underlying acceleration mechanism during magnetic reconnection. While kinetic simulations are critical to studying particle acceleration during magnetic reconnection, a reduced m i /m e is often used to alleviate the demanding computing resources, which leads to artificial scale separation between electron and proton scales. Recent kinetic simulations with high mass ratios have indicated new regimes of reconnection, as electron pressure anisotropy develops in the exhaust region and supports extended current layers. Here, we study whether different m i /m e changes the particle acceleration processes by performing a series of simulations with different mass ratio (m i /m e = 25–400) and guide field strength in a low-β plasma. We find that mass ratio does not strongly influence reconnection rate, magnetic energy conversion, ion internal energy gain, plasma energization processes, ion energy spectra, and the acceleration mechanisms for high-energy ions. Simulations with different mass ratios are different in electron acceleration processes, including electron internal energy gain, electron energy spectrum, and the acceleration efficiencies for high-energy electrons. We find that high-energy electron acceleration becomes less efficient when the mass ratio gets larger because the Fermi -like mechanism associated with particle curvature drift becomes less efficient. These results indicate that when particle curvature drift dominates high-energy particle acceleration, the further the particle kinetic scales are from the magnetic field curvature scales (~d i ), the weaker the acceleration will be.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Particle collisionality in scaled kinetic plasma simulations

Kinetic plasma processes, such as magnetic reconnection, collisionless shocks, and turbulence, are fundamental to the dynamics of astrophysical and laboratory plasmas. Simulating these processes often requires particle-in-cell (PIC) methods, but the computational cost of fully kinetic simulations can necessitate the use of artificial parameters, such as a reduced speed of light and ion-to-electron mass ratio, to decrease expense. While these approximations can preserve overall dynamics under specific conditions, they introduce nontrivial impacts on particle collisionality that are not yet well understood. In this work, we develop a method to scale particle collisionality in simulations employing an artificial speed of light and/or an artificial ion-to-electron mass ratio. By introducing species-dependent scaling factors, we independently adjust inter- and intra-species collision rates to better replicate the collisional properties of the physical system. Our approach maintains the fidelity of electron and ion transport properties while preserving critical relaxation rates, such as energy exchange timescales, within the limits of weakly collisional plasma theory. Furthermore, we demonstrate the accuracy of this scaling method through benchmarking tests against theoretical relaxation rates and connecting to fluid theory, highlighting its ability to retain key transport properties. Existing collisional PIC implementations can be easily modified to include this scaling, which will enable deeper insights into the behavior of marginally collisional plasmas across various contexts.

Totorica, S. R. [Princeton Univ., NJ (United State↗
DOE OSTI2023

The application of parallel kinetic simulations to laser and electron transport through plasmas (Final technical report)

This is a final report for the grant entitled, “The application of parallel kinetic simulations to laser and electron transport through plasmas”. The objectives of this grant were to significantly advance the fundamental understanding of the nonlinear optics of plasmas and electron transport in high-energy-density laboratory plasmas (HEDLP), including conditions of relevance to Inertial Fusion Energy (IFE). The ultimate goal was to use the understanding to determine how to fully control laser plasma interactions. The primary research tools were our own kinetic particle-in-cell software, OSIRIS, that includes kinetic physics and can run effectively on leadership class computing facilities. Therefore, one objective was to ensure that OSIRIS in continually improved so that it was more accurate and could effectively utilize state-of-the-art computing facilities. Another objective was to attract and train young researchers into the field of high energy density plasma physics. To meet the research objectives, the funds from this proposal were used to conduct research on stimulated Raman scattering (SRS) and enhance our PIC software OSIRIS. It was found that small normalized magnetic fields can in some cases mitigate SRS and that speckles can mutually interact through SRS. It was also found that it is possible for instabilities drive near quarter critical (e.g., the high frequency hybrid instability-HFHI) can generate light waves that propagate back down a density gradient where they can rescatter into the HFHI at 1/16 of the original quarter critical density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Kinetic simulations of argon dusty plasma afterglow including metastable atom kinetics

The afterglow of a dusty plasma of rf discharge in argon is simulated by the particle-in-cell-Monte Carlo collision (PIC-MCC) method. The experimental observation that heavy dust contamination of plasma leads to an anomalous increase in the electron density at the beginning of afterglow is explained by release of electrons from the dust surface. Under the assumption that the floating potential of particles is in equilibrium with plasma conditions, the fast cooling of electrons in afterglow plasma due to a rapid escape of hot electrons from the volume leads to a decrease in the magnitude of the floating potential and hence to a loss of charge by dust. The intensive desorption of electrons from nanoparticles is the origin of anomalous behavior of the electron density. At the next stage of afterglow, when the electrons become cool, the plasma decay is defined by ambipolar diffusion. The effect of metastable argon atoms is also considered. Additional ionization due to metastable atom collisions affects the electron temperature but does not change the behavior of the electron density qualitatively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Kinetic Simulations of the Interruption of Large-Amplitude Shear-Alfvén Waves in a High- β Plasma

Using two-dimensional hybrid-kinetic simulations, we explore the nonlinear “interruption” of standing and traveling shear-Alfvén waves in collisionless plasmas. Interruption involves a self-generated pressure anisotropy removing the restoring force of a linearly polarized Alfvénic perturbation, and occurs for wave amplitudes δB ⊥ /B 0 ≳β –1/2 (where β is the ratio of thermal to magnetic pressure). We use highly elongated domains to obtain maximal scale separation between the wave and the ion gyroscale. For standing waves above the amplitude limit, we find that the large-scale magnetic field of the wave decays rapidly. The dynamics are strongly affected by the excitation of oblique firehose modes, which transition into long-lived parallel fluctuations at the ion gyroscale and cause significant particle scattering. Traveling waves are damped more slowly, but are also influenced by small-scale parallel fluctuations created by the decay of firehose modes. Our results demonstrate that collisionless plasmas cannot support linearly polarized Alfvén waves above δB ⊥ /B 0 ~β –1/2 . Here, they also provide a vivid illustration of two key aspects of low-collisionality plasma dynamics: (i) the importance of velocity-space instabilities in regulating plasma dynamics at high β, and (ii) how nonlinear collisionless processes can transfer mechanical energy directly from the largest scales into thermal energy and microscale fluctuations, without the need for a scale-by-scale turbulent cascade.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Kinetic Cascade in Solar-wind Turbulence: 3D3V Hybrid-kinetic Simulations with Electron Inertia

Understanding the nature of the turbulent fluctuations below the ion gyroradius in solar-wind (SW) turbulence is a great challenge. Recent studies have been mostly in favor of kinetic Alfvén wave (KAW)-type fluctuations, but other kinds of fluctuations with characteristics typical of magnetosonic, whistler, and ion-Bernstein modes could also play a role depending on the plasma parameters. Here, we investigate the properties of the subproton-scale cascade with high-resolution hybrid-kinetic simulations of freely decaying turbulence in 3D3V phase space, including electron inertia effects. Two proton plasma beta are explored: the “intermediate” β {sub p} = 1 and “low” β {sub p} = 0.2 regimes, both typically observed in the SW and corona. The magnetic energy spectum exhibits k{sub ⊥}{sup −8/3} and k{sub ∥}{sup −7/2} power laws at β {sub p} = 1, while they are slightly steeper at β {sub p} = 0.2. Nevertheless, both regimes develop a spectral anisotropy consistent with k{sub ∥}∼k{sub ⊥}{sup 2/3} at k{sub ⊥}ρ{sub p}>1 and pronounced small-scale intermittency. In this context, we find that the kinetic-scale cascade is dominated by KAW-like fluctuations at β {sub p} = 1, whereas the low-β case presents a more complex scenario suggesting the simultaneous presence of different types of fluctuations. In both regimes, however, a possible role of the ion-Bernstein-type fluctuations at the smallest scales cannot be excluded.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Particle heating and energy partition in low- β guide field reconnection with kinetic Riemann simulations

Kinetic Riemann simulations have been completed to investigate particle heating during guide field reconnection in the low-β environment of the inner heliosphere and the solar corona. The reconnection exhaust is bounded by two rotational discontinuities (RDs), and two slow shocks (SSs) form within the exhaust as in magnetohydrodynamic (MHD) models. At the RDs, ions are accelerated by the magnetic field tension to drive the reconnection outflow as well as flows in the out-of-plane direction. The out-of-plane flows stream toward the midplane and meet to drive the SSs. The SSs differ greatly from those in the MHD model. The turbulence at the shock fronts and both upstream and downstream is weak, and so the shocks are laminar and produce little dissipation. Downstream of the SSs, the counterstreaming ion beams lead to higher density, which leads to a positive potential between the SSs which acts to confine the downstream electrons to maintain charge neutrality. The potential accelerates electrons from upstream of the SSs to the downstream region and traps a small fraction but only modestly increases the downstream electron temperature above the upstream value. In the low-β limit, the released magnetic energy is split between bulk flow and ion heating with little energy going to electrons. That the model produces neither strong electron heating nor an energetic electron component implicates that other mechanisms, such as multiple x-line reconnection, are required to explain energetic electron production in large flares. The model can be tested with the expected data from the Parker Solar Probe.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Evaluating gyro-viscosity in the Kelvin-Helmholtz instability by kinetic simulations

In the present paper, the finite-Larmor-radius (gyro-viscous) term [K. V. Roberts and J. B. Taylor, Phys. Rev. Lett. 8, 197–198 (1962)] is evaluated by using a full kinetic Vlasov simulation result of the Kelvin-Helmholtz instability (KHI). The velocity field and the pressure tensor are calculated from the high-resolution data of the velocity distribution functions obtained by the Vlasov simulation, which are used to approximate the Finite-Larmor-Radius (FLR) term according to Roberts and Taylor [Phys. Rev. Lett. 8, 197–198 (1962)]. The direct comparison between the pressure tensor and the FLR term shows an agreement. It is also shown that the anisotropic pressure gradient enhanced the linear growth of the KHI when the inner product between the vorticity of the primary velocity shear layer and the magnetic field is negative, which is consistent with the previous FLR-magnetohydrodynamic simulation result. This result suggests that it is not sufficient for reproducing the kinetic simulation result by fluid simulations to include the FLR term (or the pressure tensor) only in the equation of motion for fluid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Kinetic simulations of stimulated Raman backscattering and related processes for the shock-ignition approach to inertial confinement fusion

A detailed description of stimulated Raman backscattering and related processes for the purpose of inertial confinement fusion requires multi-dimensional kinetic simulations of a full speckle in a high-temperature, large-scale, inhomogeneous plasma. In particular for the shock-ignition scheme operating at high laser intensities, kinetic aspects are predominant. High- (I{lambda}{sub o}{sup 2}{approx}5x10{sup 15}W{mu}m{sup 2}/cm{sup 2}) as well as low-intensity (I{lambda}{sub o}{sup 2}{approx}10{sup 15}W{mu}m{sup 2}/cm{sup 2}) cases show the predominance of collisionless, collective processes for the interaction. While the two-plasmon decay instability and the cavitation scenario are hardly affected by intensity variation, inflationary Raman backscattering proves to be very sensitive. Brillouin backscattering evolves on longer time scales and dominates the reflectivities, although it is sensitive to the intensity. Filamentation and self-focusing do occur for all cases but on time scales too long to affect Raman backscattering.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗