Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic 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.

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

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

3D simulation of lower-hybrid drift waves in strong guide field asymmetric reconnection in laboratory experiments

The magnetic reconnection experiment has recently seen short wavelength ([Formula: see text]) lower-hybrid waves near the electron diffusion region in strong guide field reconnection. Based on plasma parameters from the experiment, we perform a three-dimensional fully kinetic simulation in order to investigate the generation of the lower-hybrid waves and their effects on the reconnection process. We find that the low-beta regions around the reconnection site are unstable to the lower-hybrid drift instability propagating in the outflow direction, driven by the difference between the electron and ion outflows. The waves modify the electron distributions, leading to periodic opening and closing of gaps in electron velocity space, and provide a small contribution to the anomalous resistivity. Finally, the simulation results are discussed in the context of space observations and laboratory experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Interpenetrating Plasma Simulations

This project investigated multifluid models for simulating interpenetrating plasmas. These models solve the Euler equations for each plasma stream along with an appropriate model for electron dynamics. The project developed a novel code, Euclid, that simulates these models using high-order conservative spatial discretization methods and explicit and semi-implicit time integration methods. Euclid was used to solve several multifluid test cases representative of plasma dynamics in high-energy density physics experiments. Kinetic simulations of interpenetrating flows using the LOKI code revealed novel physical effects of the ion-ion streaming instability. A new multi-flow collision operator was developed for the LOKI code. A simple model of the ion-ion streaming instability was also developed. Several presentations were given on the project results; three journal articles have been published, one has been submitted, and two are planned.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effect of Initial Microstructure on the Toughness of Coarse-Grained Heat-Affected Zone in a Microalloyed Steel

Toughness of the coarse-grained-heat-affected-zone (CGHAZ) strongly depends on the prior austenite grain size. The prior austenite grain size is affected not only by chemical composition, thermal cycle, and dissolution of second-phase particles, but also by the initial microstructure. The effect of base metal microstructure (ferrite/pearlite obtained by air cooling and martensite obtained by water-quenching) on Charpy impact toughness of the CGHAZ has been investigated for different heat inputs for high-heat input welding of a microalloyed steel. A welding thermal cycle with a heat input of 100 kJ/cm and 400 kJ/cm were simulated on the MMS-300 system. Despite a similar microstructure in the CGHAZ of both the base metals, the average Charpy impact energy for the air-cooled base metal was found to be higher than the water-quenched base metal. Through thermo-kinetic simulations, it was found that a higher enrichment of Mn/C at the ferrite/austenite transformation interface of the CGHAZ of water-quenched base metal resulted in stabilizing austenite at a lower A1 temperature, which resulted in a coarser austenite grain size and eventually lowering the toughness of the CGHAZ.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Collisionless relativistic magnetic reconnection driven by electron vortices in laser-plasma interaction

Magnetic reconnection (MR) is a fundamental process in space and laboratory plasmas. The appearance of high power lasers opens a new way to investigate MR under the relativistic condition. In this paper, relativistic collisionless MR driven by two ultra-intense lasers and a pair of asymmetric targets is studied numerically via the kinetic simulations. The static magnetic fields produced by the electron vortex structures with opposite magnetic polarities approach each other driven by the magnetic pressure and the density gradient. The antiparallel magnetic fields annihilate accompanied with the topological variation and the corresponding magnetic field energy is being dissipated to the kinetic energy of the nonthermal charged particles. Besides the outflows along the current sheet, a fast particle bunch is accelerated perpendicularly contributed by the displacement current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A moment-conserving discontinuous Galerkin representation of the relativistic Maxwellian distribution

Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell–Jüttner distribution, on a discrete simulation grid is challenging. Difficulties arise from the finite velocity bounds of the domain, which may not capture the entire distribution function, as well as errors introduced by projecting the function onto a discrete grid. Here, we present a novel scheme for iteratively correcting the moments of the projected distribution applicable to all grid-based discretizations of the relativistic kinetic equation. In addition, we describe how to compute the needed nonlinear quantities, such as Lorentz boost factors, in a discontinuous Galerkin scheme through a combination of numerical quadrature and weak operations. The resulting method accurately captures the distribution function and ensures that the moments match the desired values to machine precision.

astrophysical plasmas↗
DOE OSTI2023

Kinetic characterization of 3D magnetic reconnection: A transformative step (Final Report)

We investigated magnetic reconnection in the Earth’s magnetopause, magnetotail, and the bow shock, by means of fully kinetic simulations and analysis of space data obtained by NASA’s Magnetospheric Multiscale mission. We elucidated signatures of particle energization in particle distribution functions in the vicinity of reconnection X-line, properties of reconnection in shock-driven turbulence, and wave excitations and nonlinear structures. We advanced the understanding of magnetic reconnection from the viewpoint of particle kinetics in the diffusion region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Kinetic corrections to heat-flow and Nernst advection for laser heated plasmas

Reduced models for approximating the impact of kinetic electron behavior on the transport of thermal energy and magnetic field are investigated. The thermal flux limiter has improved agreement with Vlasov–Fokker–Planck data when a harmonic form is used that adjusts the electron mean free path to account for electron–electron collisions; these results apply to both unmagnetized and magnetized plasmas. Once a magnetic field is incorporated, the mean free path should also be modified using the electron gyroradius. A flux limiter on Nernst advection of magnetic fields is also required; a form that limits Nernst by the same fraction as the thermal heat-flow best reproduces kinetic simulations. A flux limiter form for the cross terms (Righi–Leduc and cross-gradient-Nernst) is also suggested. Hohlraum simulations relevant to fusion experiments on the National Ignition Facility are found to be sensitive to all of these details.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Kinetic Monte Carlo simulations of structural evolution during anneal of additively manufactured materials

Our experiments indicated that upon a post-processing anneal, an additively manufactured 316L stainless steel exhibits cubic grains rather than the conventional equiaxed grains. In this work, we have used kinetic Monte Carlo simulations to explore the origin of these cubic grains. First, we implemented a new kinetic Monte Carlo model in parallel code SPPARKS to simulate grain growth and recrystallization under a residual energy distribution. Our model incorporates physical properties and real-time, as opposed to generic properties and relative time. We further validated that our SPPARKS simulations reproduced the expected kinetic behavior of single-grain evolution. We then used the validated approach to simulate the anneal of an additively manufactured material under the same conditions used in our experiments. We found that the cubic grains can origin from a periodically varying residual energy that may be present in additively manufactured materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Transport physics dependence of Bohm speed in presheath–sheath transition

We report the ion exit flow speed at the sheath entrance is constrained by the Bohm criterion, which is used as a boundary condition for simulations that do not resolve the sheath region. Traditional Bohm criterion analysis invokes the equation of state and, thus, ignores transport physics in the sheath transition problem. An expression for the Bohm speed away from the asymptotic limit is derived from a set of anisotropic plasma transport equations. The thermal force, collisional temperature isotropization, and heat flux enter into the evaluation of the Bohm speed. By comparison with kinetic simulation results, this expression is shown to be accurate in the presheath–sheath transition region rather than a single point at the sheath entrance over a broad range of collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Constraining the Reaction Rate of Criegee Intermediates with Carboxylic Acids during the Multiphase Ozonolysis of Aerosolized Alkenes

Criegee intermediates (CIs) are central to the ozonolysis of unsaturated organic compounds, controlling the formation of small fragmentation reaction products and higher molecular weight oligomeric compounds through competing unimolecular and bimolecular pathways. In particular, the reaction of thermalized CI with carboxylic acids (RCOOH) to produce oligomeric α-acyloxyalkyl hydroperoxides (AAHPs) is of interest to the chemical transformation of organic compounds in the atmosphere. In the gas phase, the CI + RCOOH reaction proceeds at supercollisional rates, consistent with a barrierless reaction mechanism. However, the rate of this reaction in condensed phases remains uncertain. Here, we report the results of aerosol flow tube studies of the ozonolysis of submicron organic aerosol containing mixtures of an internal n-alkene (Z-9-tricosene, tri) and a saturated acid (2-hexyldecanoic acid, HDA). The decay of the reactants and formation of reaction products are monitored using atmospheric pressure chemical ionization mass spectrometry (APCI–MS). A 6-fold decrease in the reactive uptake of ozone is observed with increasing acid concentration, consistent with a previous study using alcohol additives. The decay of HDA with ozone exposure shows that HDA scavenges CIs at rates comparable to competing unimolecular pathways involving CIs. The decay kinetics of HDA are analyzed using a simple kinetic model in order to constrain the ratio of the CI + RCOOH and unimolecular CI rate constants. Empirical fitting to this model, combined with theoretical estimates of the unimolecular loss rate, yields a CI + RCOOH rate constant of 1.85 ± 0.27 × 10 –19 cm 3 molec –1 s –1 , which is six orders of magnitude slower than the expected diffusion limit in a liquid organic matrix. Stochastic kinetic simulations are conducted with different values of the CI + RCOOH rate constant to validate this result and show that the diffusion-limited rate is indeed much too rapid to reproduce the experimental results.

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

Flow Channels and the Generation of Alfvénic Turbulence Along Storm‐Time Inner Magnetospheric Field‐Lines

Abstract A feature of Earth's storm‐time magnetosphere outside the plasmapause is the occurrence of broad‐spectrum Alfvénic fluctuations. In this letter observations from the Van Allen Probes are compared with 3‐D fluid‐kinetic simulations of an evolving convective flow channel to investigate the mechanisms generating the observed spectrum. It is shown how narrow channels of fast convection are unstable to the Kelvin‐Helmholtz instability which on closed field‐lines initiates a cascade to small scales. Sustained driving of the flow combined with reflection from the topside ionosphere leads to the generation of an intensified spectrum of electromagnetic structures having similar spectral and morphological characteristics to those observed. This process couples enhanced magnetospheric convection to kinetic scale electromagnetic fluctuations that drive particle transport, scattering and energization through the outer radiation belt and ring current during geomagnetic storms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Thermodynamic Phase Transition in Magnetic Reconnection

By examining the entropy production in fully kinetic simulations of collisional plasmas, it is shown that the transition from collisional Sweet-Parker reconnection to collisionless Hall reconnection may be viewed as a thermodynamic phase transition. The phase transition occurs when the reconnection electric field satisfies E = $E_{D}\sqrt{m{e}/{m}_i}$ , where ${m{e}/{m}_i}$ is the electron-to-ion mass ratio and E D is the Dreicer electric field. This condition applies for all ${m{e}/{m}_i}$, including ${m{e}/{m}_i}$ = 1 , where the Hall regime vanishes and a direct phase transition from the collisional to the kinetic regime occurs. In the limit ${m{e}/{m}_i}$ → 0, this condition is equivalent to there being a critical electron temperature $Te ≈ m_iΩ^{2}_{i}δ^{2}$, where $Ω_i$ is the ion cyclotron frequency and δ is the current sheet half-thickness. The heat capacity of the current sheet changes discontinuously across the phase transition, and a critical power law is identified in an effective heat capacity. A model for the time-dependent evolution of an isolated current sheet in the collisional regime is derived.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Data-driven, multi-moment fluid modeling of Landau damping

Deriving governing equations of complex physical systems based on first principles can be quite challenging when there are certain unknown terms and hidden physical mechanisms in the systems. In this work, we apply a deep learning architecture to learn fluid partial differential equations (PDEs) of a plasma system based on the data acquired from a fully kinetic model. Here, the learned multi-moment fluid PDEs are demonstrated to incorporate kinetic effect such as Landau damping. Based on the learned fluid closure, the data-driven, multi-moment fluid modeling can well reproduce all the physical quantities derived from the fully kinetic model. The calculated damping rate of Landau damping is consistent with both the fully kinetic simulation and the linear theory. The data-driven fluid modeling of PDEs for complex physical systems may be applied to improve the fluid closure and reduce the computational cost of multi-scale modeling of global systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Off-Equilibrium Reactivity of Boron-Enriched Metal Diboride Surfaces in Electroreduction Conditions

Boron-based materials, featuring B-dependent reactivity and diverse phases, are emerging as promising catalyst systems. However, the catalytic mechanism on many borides remains poorly understood due to complex surface reconstructions under reaction conditions. Here, we investigate the MoB 2 surface in conditions of hydrogen evolution reaction in acidic media, using grand canonical global optimization, grand canonical density functional theory, ab initio molecular dynamics, free energy surface sampling, and an analytical model for electrochemical barrier evaluation. We propose a boron-enrichment strategy to tune the surface reactivity of the hexagonal face of MoB 2 . We reveal the dynamic nature of the B-enriched surface under H coverage and kinetic trapping of the system in the metastable regime with an extensive examination of the deactivation pathways. The metastable center B site on B-enriched surfaces, featuring buckled-up configuration and a usual relaxation effect, is found to be highly active toward HER via the Volmer–Heyrovsky mechanism. In conclusion, this work demonstrates how off-equilibrium behaviors can arise from the interplay between adsorbate coverage and surface reconstruction on a seemingly simple surface, and we present a theoretical framework and computational workflows to address these behaviors, along with other realistic complexities, in kinetics simulations.

Adsorption↗
DOE OSTI2025

Ohm’s Law, the Reconnection Rate, and Energy Conversion in Collisionless Magnetic Reconnection

Magnetic reconnection is a ubiquitous plasma process that transforms magnetic energy into particle energy during eruptive events throughout the universe. Reconnection not only converts energy during solar flares and geomagnetic substorms that drive space weather near Earth, but it may also play critical roles in the high energy emissions from the magnetospheres of neutron stars and black holes. In this review article, we focus on collisionless plasmas that are most relevant to reconnection in many space and astrophysical plasmas. Guided by first-principles kinetic simulations and spaceborne in-situ observations, we highlight the most recent progress in understanding this fundamental plasma process. We start by discussing the non-ideal electric field in the generalized Ohm’s law that breaks the frozen-in flux condition in ideal magnetohydrodynamics and allows magnetic reconnection to occur. We point out that this same reconnection electric field also plays an important role in sustaining the current and pressure in the current sheet and then discuss the determination of its magnitude (i.e., the reconnection rate), based on force balance and energy conservation. This approach to determining the reconnection rate is applied to kinetic current sheets with a wide variety of magnetic geometries, parameters, and background conditions. We also briefly review the key diagnostics and modeling of energy conversion around the reconnection diffusion region, seeking insights from recently developed theories. Finally, future prospects and open questions are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Artificial correlation heating in PIC simulations

The Particle-in-Cell (PIC) method, a cornerstone in plasma modeling, is widely employed for its ability to simulate kinetic phenomena in device-scale domains. Part of what makes this possible is that computational macroparticles represent many physical particles. It converges under certain constraints, including a grid spacing that resolves the Debye length and a time step small enough to respect the Courant–Friedrichs–Lewy condition and plasma frequency stability limit. Here, we introduce a new constraint necessary to avoid Artificial Correlation Heating (ACH). This requires that the macroparticle coupling strength be smaller than one, Γ ω < 1, where Γ ω ≡ Γω 2/3 ⁠, Γ = Z 2 e 2 /(4πε 0 ak B T) is the physical coupling strength, and w is the macroparticle weight. This is particularly relevant to 3D simulations of dense plasmas, which are becoming common with modern computing power. If this condition is violated, the finite macroparticle weight artificially enhances the coupling strength and causes the plasma to heat until the macroparticle coupling strength is near unity, depending on the grid resolution. A comprehensive model of ACH is developed that incorporates electron density, temperature, macroparticle weight, and grid resolution. It is then tested using PIC simulations, delineating the boundaries of the method's applicability and offering a predictive framework for ACH. Moreover, the research explores a runaway heating process induced by ACH in the presence of ionization, which can lead to numerical instability. A conclusion of this study is that the onset of ACH can impose a more stringent constraint on the macroparticle weight and average number of macroparticles per cell than what is typically expected, particularly in 3D simulations of dense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Energy Partition in Collisionless Counterstreaming Plasmas

Fast, counterstreaming plasma outflows drive magnetic field amplification, plasma heating, and particle acceleration in numerous astrophysical environments, from supernova remnant shocks to active galactic nuclei jets. Understanding how, in the absence of Coulomb collisions, energy is redistributed between the different plasma species remains a fundamental open question. We use three-dimensional, fully kinetic simulations to investigate energy partition in weakly magnetized counterpropagating plasmas. Our results reveal a complex interplay between different processes, where at early times the Weibel instability drives a first stage of magnetic field amplification and at late times the kinking of current filaments drives a second amplification stage via a dynamo-type mechanism. Electrons are heated primarily during the latter phase through magnetic pumping. By the time the flows thermalize, we observe that the final temperature ratio T e /T i and energy partition depend on the ion-to-electron mass ratio. For electron–proton flows, the electron thermal energy only reaches up to a few percent of the initial ion kinetic energy.

Marret, Alexis [SLAC National Accelerator Laborato↗
DOE OSTI2025

Electron-scale Magnetic Holes Generation Driven by Whistler-to-Bernstein Mode Conversion in Fully Kinetic Plasma Turbulence

Magnetic holes (MHs) are coherent structures characterized by a strong and localized magnetic field amplitude dip, commonly observed in the heliosphere. These structures come in different sizes, from magnetohydrodynamic to kinetic scales. Subion-scale MHs are usually sustained by an electron current vortex and exhibit a strong electron temperature anisotropy, with higher temperatures perpendicular to the background magnetic field. Magnetospheric multiscale observations (MMSs) have revealed electron-scale MHs to be ubiquitous in the turbulent Earth’s magnetosheath and the solar wind, potentially playing an important role in the energy cascade and dissipation. Despite abundant observations, the origin of electron-scale MHs is still unclear and debated. In this work, we use fully kinetic simulations to investigate the role of plasma turbulence in generating electron-scale MHs. We find that the turbulence spontaneously produces electron-scale MHs via the following mechanism: first, large-scale turbulent velocity shears produce regions with high electron temperature anisotropy; these localized regions become unstable, generating oblique electron-scale whistler waves; as they propagate over the inhomogeneous turbulent background, whistler fluctuations develop an electrostatic component, turning into Bernstein-like modes; the strong electrostatic fluctuations produce current filaments that merge into an electron-scale current vortex; the resulting electron vortex locally reduces the magnetic field amplitude, finally evolving into an electron-scale MH. We show that MHs generated by this mechanism have properties consistent with MMSs and nontrivial kinetic features with a “mushroom”-shaped electron velocity distribution function. Our results have potential implications for understanding the formation and occurrence of electron-scale MHs in astrophysical turbulent and space environments, such as the Earth’s magnetosheath and the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗