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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

2D kinetic-ion simulations of inverted corona fusion targets

Laser-driven “inverted corona” fusion targets have attracted interest as a low-convergence neutron source and platform for studying kinetic physics. The scheme consists of a hollow or gas-filled spherical shell made of deuterated plastic. The shell has one or more laser entrance holes (LEH), resembling a spherical hohlraum. The laser passes through the LEH’s and illuminates the interior surface of the shell, ablating a plasma that travels inward towards the target center. Long ion mean free paths in the converging plasma can lead to significant interpenetration, atomic mix, and other kinetic effects. Here, in this work we report on numerical simulations of inverted corona targets using the kinetic-ion, fluid–electron hybrid particle-in-cell (PIC) approach in 2D RZ geometry. 2D simulations suggest that shape effects do not have a significant impact on plasma evolution and observed yield trends are primarily the result of 1D kinetic mix mechanisms. Simulations are also compared against available experimental data recorded at the OMEGA laser facility. In particular, synthetic x-ray emission images show good qualitative agreement with experimental results, albeit with an apparent timing discrepancy for the two-sided vacuum target. More generally, we demonstrate the potential of hybrid-PIC simulations for full-system modeling and experimental design, including collisional absorption of laser energy, plasma evolution, mix, and fusion burn.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Kinetic Plasma Simulation Capabilities in the MOOSE Framework: Verification of Particle-Particle Collisions

High-fidelity simulations of complex plasma systems allow researchers to gain key insights into and understanding of these systems. To facilitate massively parallel high-fidelity plasma simulations, finite-element-based particle-in-cell capabilities are being developed within the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) based framework called Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER). While SALAMANDER’s primary objective is modeling edge plasmas and plasma-facing components in fusion devices, the particle-in-cell capabilities being developed are general and will support modeling low-temperature plasmas as well. Previously, collisionless magnetostatic simulation capabilities have been verified with the two-stream and Dorey-Guest-Harris instabilities, and single particle motion. Collisions were implemented using the direct simulation Monte Carlo method, and verification of this capability will be presented here several verification problems: relaxation of a randomly initialized gas to a Maxwellian distribution, Fourier heat flow, and comparison of reaction rates to both analytic calculations and those calculated using a multi-term Boltzmann solver.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Orbit-averaging and deposition accuracy for runaway electron beams in hybrid kinetic-MHD simulations of the runaway plateau

We develop a new procedure that combines the kinetic orbit runaway electron code (KORC) and the NIMROD extended-magnetohydrodynamic code to simulate runaway electrons (REs) in the post-disruption plateau. KORC integrates guiding-center orbits, with a barycentric-based binary search strategy providing initial guesses for the Newton–Raphson logical-to-physical coordinate inversion, ensuring reliable particle-to-mesh mapping in NIMROD, whose fields remain static for the present study. Samples are drawn in accord with experimental parallel current profiles of RE beams during the plateau phase. Deposition in NIMROD is verified through comparison with a Python-based finite-element code that ensures periodicity in the poloidal direction and continuity at the magnetic axis. Accurate representation of near-axis fields requires finer mesh resolution to prevent under- and overshoots in current density from orbit inaccuracies. Yet, at a fixed particle count, increasing mesh resolution amplifies statistical noise in the deposited fields. An orbit-averaging method accumulates partial current deposits over multiple kinetic steps and reduces the statistical noise with little added computational cost. By coupling kinetic routines from KORC directly into the NIMROD codebase, these developments lay essential groundwork for future self-consistent KORC–NIMROD coupling.

Algorithms and data structure↗
DOE OSTI2025

Nonlinear evolution, propagation, electron-trapping, and damping effects of ion-acoustic solitons using fully kinetic PIC simulations

We investigate ion acoustic solitary waves (solitons) of varying amplitudes in a one-dimensional plasma using fully kinetic particle-in-cell simulations. The initial soliton conditions are based on the Korteweg–de Vries (KdV) equation, treating ions as a cold species and electrons with finite temperature. Our findings reveal that KdV solitons evolve nonlinearly to a saturated state at higher amplitude, deviating from KdV predictions for ion density and electric potential, and from the Boltzmann relation for electron density. At this saturated state, the KdV model cannot accurately describe the soliton behavior. For small amplitudes, Sagdeev's model describes the saturated state, but not the soliton width; for larger amplitudes, it models the width accurately, but not the amplitude. These discrepancies arise from assuming a Boltzmann relation for electron density, while electron trapping creates non-Boltzmann densities—a deviation that increases with soliton amplitude. Additionally, we observe that the soliton amplitude oscillates roughly at the electron bounce frequency. The soliton is better described by Schamel's electron density formulation and a modified KdV equation incorporating electron trapping. The soliton velocity matches best with predictions from Sagdeev's and Schamel's models. Moreover, the soliton speed–amplitude relationship differs from existing theoretical predictions. Finally, we find minimal ion and electron Landau damping effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Influences of δB contribution and parallel inertial term of energetic particles on MHD-kinetic hybrid simulations: a case study of the 1/1 internal kink mode

The magnetohydrodynamic-kinetic (MHD-kinetic) hybrid model (Park et al 1992 Phys. Fluids B 4 2033–7) has been widely applied in studying energetic particles (EPs) problems in fusion plasmas for past decades. The pressure-coupling scheme or the current-coupling scheme is adopted in this model. However, two noteworthy issues arise in the model application: firstly, the coupled term introduced in the pressure-coupling scheme, (∇•P h ) ⟂ , is often simplified by ∇•P h , which is equivalent to neglecting the parallel inertial term of EPs; secondly, besides the $δf$ contribution caused by changing in the EP distribution function, the magnetic field perturbation (the $δB$ contribution) generated during development of the instabilities should also be considered, but it is often ignored in existing hybrid simulations. In this paper, we derive the analytical formulations under these two coupling schemes and then numerically study the representative case of the linear stability of the $m/n$ = $1/1$ internal kink mode (IKM) (Fu et al 2006 Phys. Plasmas 13 052517) by using the CLT-K code. Further, it is found that the approximated models can still yield reasonable results when EPs are isotopically distributed. But it fails completely in cases with anisotropic EP distributions. In addition, we further investigate the influence of EP's orbit width on the stability of IKM and verify the equivalence between pressure-coupling scheme and the current-coupling scheme.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

Self-Consistent Simulation of the Magnetotail

Kinetic effects may be critically important at various stages of magnetotail evolution prior to and after the break up of substorms. The quasi-adiabaticity of ion dynamics results in the appearance of pronounced off diagonal terms in the pressure tensor (an effect that cannot be reproduced using the standard MHD (magnetohydrodynamic) approach). Only the appearance of these terms could support the formation of thin current sheets that recent experimental data have repeatedly indicated. We present in this paper the results of a self-consistent large-scale kinetic simulation of the Earth's magnetotail. We launch ions from the plasma mantle and self-consistently adjust the background magnetic field by computing the perturbations magnetic field from the particle current. We find that a thin current sheet with a half-thickness of approximately 0.3 R(sub E) quickly develops in the simulation. Also, the magnetotail evolves into a periodic state where the location of the X-line in the tail oscillates between x approximately 40 R(sub E) and x approximately 70 R(sub E) on a 5-minute time scale. Evidence of this periodic motion is also seen in the auroral precipitation profile of ions.

Peroomian, Vahe↗
DOE OSTI2026

Capturing Secondary Kinetic Instabilities in Three‐Dimensional Dayside Reconnection Using an Improved Gradient‐Based Closure

Magnetic reconnection is a highly dynamic process that excites a wide variety of kinetic waves and instabilities. Transverse current sheet instabilities such as the lower-hybrid drift and secondary drift-kink instabilities in particular have been shown by kinetic simulations to modify the reconnection and introduce significant turbulence and mixing to the reconnection layer. Past studies using the ten-moment fluid model to capture important kinetic physics such as the electron inertia and full representation of the pressure tensor proved advantageous to a two-fluid representation of reconnection, but the model struggled when using a local relaxation closure for the heat flux to replicate the current sheet instabilities and subsequent mixing seen in kinetic simulations. This work uses the Gkeyll software framework to perform simulations of asymmetric reconnection based on the 16 October 2015 MMS crossing of a diffusion region, the Burch event. An improved gradient-based heat flux closure is implemented, showing significant improvement in secondary kinetic instabilities that grow in the current sheet. These instabilities generate turbulence which leads to growth of secondary magnetic islands and flux ropes.

Bradshaw, K. [Princeton University, NJ (United Sta↗
NASA NTRS2004

A Massively Parallel Hybrid Dusty-Gasdynamics and Kinetic Direct Simulation Monte Carlo Model for Planetary Applications

In order to understand the global structure, dynamics, and physical and chemical processes occurring in the upper atmospheres, exospheres, and ionospheres of the Earth, the other planets, comets and planetary satellites and their interactions with their outer particles and fields environs, it is often necessary to address the fundamentally non-equilibrium aspects of the physical environment. These are regions where complex chemistry, energetics, and electromagnetic field influences are important. Traditional approaches are based largely on hydrodynamic or magnetohydrodynamic (MHD) formulations and are very important and highly useful. However, these methods often have limitations in rarefied physical regimes where the molecular collision rates and ion gyrofrequencies are small and where interactions with ionospheres and upper neutral atmospheres are important. At the University of Michigan we have an established base of experience and expertise in numerical simulations based on particle codes which address these physical regimes. The Principal Investigator, Dr. Michael Combi, has over 20 years of experience in the development of particle-kinetic and hybrid kinetichydrodynamics models and their direct use in data analysis. He has also worked in ground-based and space-based remote observational work and on spacecraft instrument teams. His research has involved studies of cometary atmospheres and ionospheres and their interaction with the solar wind, the neutral gas clouds escaping from Jupiter s moon Io, the interaction of the atmospheres/ionospheres of Io and Europa with Jupiter s corotating magnetosphere, as well as Earth s ionosphere. This report describes our progress during the year. The contained in section 2 of this report will serve as the basis of a paper describing the method and its application to the cometary coma that will be continued under a research and analysis grant that supports various applications of theoretical comet models to understanding the inner comae of comets (grant NAGS- 13239 from the Planetary Atmospheres program).

Combi, Michael R.↗
NASA NTRS2022

Probing the High-energy Gamma-Ray Emission Mechanism in the Vela Pulsar via Phase-resolved Spectral and Energy-dependent Light-curve Modeling

Recent kinetic simulations sparked a debate regarding the emission mechanism responsible for pulsed GeV γ-ray emission from pulsars. Some models invoke curvature radiation, while other models assume synchrotron radiation in the current sheet. We interpret the curved spectrum of the Vela pulsar as seen by H.E.S.S. II (up to ∼100 GeV) and the Fermi Large Area Telescope to be the result of curvature radiation due to primary particles in the pulsar magnetosphere and current sheet. We present phase-resolved spectra and energy-dependent light curves using an extended slot gap and current-sheet model, invoking a step function for the accelerating electric field as motivated by kinetic simulations. We include a refined calculation of the curvature radius of particle trajectories in the lab frame, impacting the particle transport, predicted light curves, and spectra. Our model reproduces the decrease of the flux of the first peak relative to the second one, evolution of the bridge emission, near-constant phase positions of peaks, and narrowing of pulses with increasing energy. We can explain the first of these trends because we find that the curvature radii of the particle trajectories in regions where the second γ-ray light-curve peak originates are systematically larger than those associated with the first peak, implying that the spectral cutoff of the second peak is correspondingly larger. However, an unknown azimuthal dependence of the E field, as well as uncertainty in the precise spatial origin of the GeV emission, precludes a simplistic discrimination of emission mechanisms.

Monica Barnard↗
NASA NTRS2008

Plasma Sheet Circulation Pathways

Global simulations of Earth's magnetosphere in the solar wind compute the pathways of plasma circulation through the plasma sheet. We address the pathways that supply and drain the plasma sheet, by coupling single fluid simulations with Global Ion Kinetic simulations of the outer magnetosphere and the Comprehensive Ring Current Model of the inner magnetosphere, including plasmaspheric plasmas. We find that the plasma sheet is supplied with solar wind plasmas via the magnetospheric flanks, and that this supply is most effective for northward IMF. For southward IMF, the innermost plasma sheet and ring current region are directly supplied from the flanks, with an asymmetry of single particle entry favoring the dawn flank. The central plasma sheet (near midnight) is supplied, as expected, from the lobes and polar cusps, but the near-Earth supply consists mainly of slowly moving ionospheric outflows for typical conditions. Work with the recently developed multi-fluid LFM simulation shows transport via plasma "fingers" extending Earthward from the flanks, suggestive of an interchange instability. We investigate this with solar wind ion trajectories, seeking to understand the fingering mechanisms and effects on transport rates.

Moore, Thomas E.↗
NASA NTRS2015

The Role of the Hall Effect in Global Structure and Dynamics of Planetary Magnetospheres: Ganymede as a Case Study

We present high-resolution Hall MHD simulations of Ganymede's magnetosphere demonstrating that Hall electric fields in ion-scale magnetic reconnection layers have significant global effects not captured in resistive MHD simulations. Consistent with local kinetic simulations of magnetic reconnection, our global simulations show the development of intense field-aligned currents along the magnetic separatrices. These currents extend all the way down to the moon's surface, where they may contribute to Ganymede's aurora. Within the magnetopause and magnetotail current sheets, Hall J x B forces accelerate ions to the local Alfven speed in the out-of-plane direction, producing a global system of ion drift belts that circulates Jovian magnetospheric plasma throughout Ganymede's magnetosphere. We discuss some observable consequences of these Hall-induced currents and ion drifts: the appearance of a sub-Jovian 'double magnetopause' structure, an Alfvenic ion jet extending across the upstream magnetopause, and an asymmetric pattern of magnetopause Kelvin-Helmholtz waves.

Magnetic Reconnection↗
NASA NTRS2015

Mission Oriented Support and Theory (MOST) for MMS -- The Goddard Space Flight Center/University of California Los Angeles Interdisciplinary Science Program

The MOST IDS team was tasked with focusing on two general areas: The first was to participate with the Fast Plasma Investigation (FPI) team in the development of virtual detectors that model the instrument responses of the MMS FPI sensors. The virtual instruments can be 'flown through' both simulation data (from magnetohydrodynamic, hybrid, and kinetic simulations) and Cluster and THEMIS spacecraft data. The goal is to determine signatures of magnetic reconnection expected during the MMS mission. Such signatures can serve as triggers for selection of burst mode downloads. The chapter contributed by the FPI team covers that effort in detail and, therefore, most of that work has not been included here. The second area of emphasis, and the one detailed in this chapter, was to build on past and present knowledge of magnetic reconnection and its physical signatures. Below we describe intensive analyses of Cluster and THEMIS data together with theoretical models and simulations that delineate the plasma signatures that surround sites of reconnection, including the effects of turbulence as well as the detailed kinetic signatures that indicate proximity to reconnection sites. In particular, we point out that particles are energized in several regions, not only at the actual site of reconnection.

Goldstein, Melvyn L.↗
DOE OSTI2026

A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes

Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Constraining the Low-Temperature Oxidation Mechanism of n -Hexanol through the Detection and Identification of C 6 Elusive Intermediates

Alcohol-based fuels are currently considered to be viable energy carriers for the transportation sector. Consequently, a comprehensive mechanistic understanding of the low-temperature oxidation of alcohols is essential for application in advanced low-temperature compression engines. Here, in this work, a multidimensional approach involving experimental investigations, kinetic modeling, and theoretical calculations was used to provide new insights into the low-temperature oxidation mechanism of a C 6 alcohol, n -hexanol (CH 3 (CH 2 ) 5 OH), through the detection and identification of elusive C 6 intermediates. The oxidation of n-hexanol was investigated in a jet-stirred reactor under stoichiometric conditions (ϕ = 1.0), an initial fuel concentration of 2%, a residence time of 2 s, a temperature range between 500 and 660 K, and a pressure of 700 Torr. The reactants, intermediates, and final products were detected and identified by means of molecular-beam mass spectrometry coupled with single-photon ionization employing tunable synchrotron-generated vacuum ultraviolet radiation. Chemical kinetic simulations were performed using a previously published kinetic model (Togbé et al., Energy Fuels 2010, 11, 5859−5875) to predict the reactivity of n-hexanol and elucidate the predominant formation pathways of the observed low-temperature species. Experimental photoionization efficiency curves in conjunction with ab initio calculations, enabled the identification of important low-temperature species, such as C 6 unsaturated alcohols, C 6 olefinic hydroperoxides, C 6 cyclic ethers, C 6 diones, and C 6 ketohydroperoxides. The results of this study provide valuable insight into the mechanism of the low-temperature oxidation chemistry of n -hexanol, contributing to the development of kinetic models for the low-temperature oxidation of n-hexanol and other long-chain linear alcohols.

alcohols↗
DOE OSTI2025

Planar laser-induced fluorescence system for spatiotemporal ion velocity distribution function measurements

In this work, we present a planar laser-induced fluorescence (PLIF) system for measuring two-dimensional (2D), spatiotemporally resolved ion velocity distribution functions (IVDFs). A continuous-wave tunable diode laser produces a laser sheet that irradiates the plasma, and the resulting fluorescence is captured by an intensified CCD (ICCD) camera. Fluorescence images recorded at varying laser wavelengths are converted into 2D IVDFs using the Doppler shift principle. The developed diagnostic is implemented in an electron beam generated E × B plasma with a bulk plasma density of $\sim\!\!{10^{10}}{\text{c}}{{\text{m}}^{ - 3}}$ . The developed diagnostic is validated against a conventional single-point LIF method using photomultiplier tube-based detection, while significantly reducing the total measurement time by the number of spatial positions measured. The time-resolving capability of this diagnostic is tested by oscillating the plasma between two nominal operating modes with different density profiles and triggering the ICCD camera by the externally driven plasma oscillation. The measured 2D IVDF maps reveal several signatures of ion dynamics in this plasma source, including radially outflowing ions across the electric field and anomalous ion heating at the periphery, consistent with recent kinetic simulations and theoretical studies. A possible correlation between these ion kinetic features and rotating spoke structures is discussed.

E×B plasma↗
DOE OSTI2025

Planar Collisionless Shock Simulations with the Semi-implicit Particle-in-cell Model FLEKS

This study investigates the applicability of the semi-implicit particle-in-cell code FLexible Exascale Kinetic Simulator (FLEKS) to heliospheric shock simulations. We examine one- and two-dimensional local planar shock simulations, initialized using MHD states with upstream conditions representative of plasmas in the hypersonic, β ∼ 1 regime, for both quasi-perpendicular and quasi-parallel configurations. The refined algorithm in FLEKS proves robust, enabling accurate shock simulations with a grid resolution on the order of the electron inertial length d e . Our simulations successfully capture key shock features, including shock structures (foot, ramp, overshoot, and undershoot), upstream and downstream waves (fast magnetosonic, whistler, Alfvén ion-cyclotron, and mirror modes), and non-Maxwellian particle distributions. Crucially, we find that at least two spatial dimensions are critical for accurately reproducing downstream-wave physics in quasi-perpendicular shocks and capturing the complex dynamics of quasi-parallel shocks, including surface rippling, shocklets, short, large-amplitude magnetic structures, magnetic reconnection, and jets. Furthermore, our parameter studies demonstrate the impact of mass ratio and grid resolution on shock physics. This work provides valuable guidance for selecting appropriate physical and numerical parameters for shock simulations using a semi-implicit PIC method, paving the way for incorporating kinetic shock processes into large-scale collisionless plasma simulations with the MHD-AEPIC model.

plasma astrophysics↗