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

Impact of gasoline composition on the effects of nitric oxide on autoignition and knock in a DISI engine

Modern spark-ignition engines use exhaust gas recirculation (EGR) to dilute the charge and suppress knock, enabling the use of higher compression ratios and/or more optimum combustion phasing for higher efficiency. The effectiveness of EGR is affected by the composition of the fuel and its chemical-kinetic interactions with combustion products. Among those, nitric oxide (NO) has been shown to strongly affect autoignition reactivity. However, the impact of fuel composition of the effect of NO on reactivity is not well-understood. Here, in this study, engine experiments were conducted to assess the impact of NO seeded to the intake on knock-limited operation of two gasoline fuels (high cycloalkane content, or HCA, and high olefin content, or HO). Results showed that compositionally-different fuels responded differently to NO. HCA, which was less knock-limited than HO for NO < 200 ppm, became more knock-limited for NO > 200 ppm. Moreover, it was found that differences in knock between fuels were caused by differences in autoignition chemistry and not in the sequential autoignition process of the end gas that occurs due to thermal stratification. Chemical kinetic simulations were performed to better understand the experimental results. For HCA, intermediate-temperature heat release had a greater impact on autoignition reactivity than low-temperature heat release, while the opposite was observed for HO. For both fuels, NO enhances the magnitude of low-temperature heat release via NO + HO 2 → NO 2 + OH. The effect of NO on reactivity was stronger for HCA because OH produced from NO helped to overcome the OH quenching effect of cyclopentane, a main species in HCA. In contrast, HO had relatively strong inherent low-temperature chemistry arising from iso-octane, which reduced the impact of NO on reactivity. For the range of NO mole fractions tested in this study, in-cylinder NO increased fuel’s knock propensity, especially for fuels with mild low-temperature chemistry.

Autoignition↗
DOE OSTI2024

Investigating the potential of a higher reactivity fuel to achieve faster heat-up of aftertreatment systems

Shortening the aftertreatment system heat-up period will be crucial to meet upcoming emissions regulations, which mandate significant reductions in tail-pipe emissions for on-road heavy-duty vehicles. This work studies the impact of increasing fuel reactivity on the combustion of post injections during the expansion stroke. The work demonstrates an operating strategy using a blend of di-butyl ether (a potential low-carbon biofuel component) and #2 diesel that achieves higher exhaust enthalpy, while maintaining similar engine-out emissions (NOx and combustion efficiency) in comparison to the baseline #2 diesel fuel (cetane number 41). On average, a 12% increase in exhaust enthalpy was observed for the blend (cetane number 60) relative to the baseline fuel. The ability to operate the higher reactivity fuel at a more retarded post-injection timing, where combustion with #2 diesel was found to be more unstable, and less complete was key to achieving higher exhaust enthalpy. However, a small fuel penalty (3%–5%) was also associated with the higher reactivity fuel with this operating strategy. The formation of fuel lean regions, due to the long ignition delays at late injection conditions, was predicted to be the leading contributor to the increased emissions of products of partial combustion. Here, chemical kinetic simulations using surrogates for the fuels used in the experiments, demonstrated the ability of a higher cetane fuel to outperform the baseline fuel under fuel lean conditions with equivalence ratios less than 0.7.

09 BIOMASS FUELS↗
DOE OSTI2025

Directed energy deposition of functionally graded V-4Cr-4Ti to Fe-9Cr transition for fusion power systems

This study proposes a graded structure via additive manufacturing for divertor and first wall blanket applications in fusion reactors. Materials were selected based on thermodynamic calculations to operate from 1100 °C at the plasma-facing level to 550 °C at the structural steel level. Conventional joining methods often lead to failures due to discrete reaction layers with significant mechanical property differences. Using laser beam-directed energy deposition (LB-DED), this study demonstrates the fabrication of a VCrTi-Gr91 steel functionally graded component through a novel process parameter optimization framework. A systematic approach included powder characterization, single-track depositions, and construction of printability maps. Near full-density specimens of each interlayer were additively manufactured, and a transition from V-based alloys to reduced activation ferritic martensitic steels was achieved. Computational material selection of interlayer alloys and thermodynamic/diffusion kinetics simulations prevented most interface incompatibilities. A brittle intermetallic formed at one interface, causing cracking, which was not predicted by current thermodynamic models. Transition alloy design approach was updated with a more recent database and a mitigation strategy has been proposed to eliminate the formation of deleterious intermetallic phases. Ultimately, LB-DED has proven effective for producing multi-material graded systems for fusion applications, with the demonstrated process parameter optimization framework applicable to various materials.

Additive manufacturing↗
DOE OSTI2024

Impact of the Out‐Of‐Plane Flow Shear on Magnetic Reconnection at the Flanks of Earth's Magnetopause

Abstract Magnetic reconnection changes the magnetic field topology and facilitates the energy and particle exchange at magnetospheric boundaries such as the Earth's magnetopause. The flow shear perpendicular to the reconnecting plane prevails at the flank magnetopause under southward interplanetary magnetic field conditions. However, the effect of the out‐of‐plane flow shear on asymmetric reconnection is an open question. In this study, we utilize kinetic simulations to investigate the impact of the out‐of‐plane flow shear on asymmetric reconnection. By systematically varying the flow shear strength, we analyze the flow shear effects on the reconnection rate, the diffusion region structure, and the energy conversion rate. We find that the reconnection rate increases with the upstream out‐of‐plane flow shear, and for the same upstream conditions, it is higher at the dusk side than at the dawn side. The diffusion region is squeezed in the outflow direction due to magnetic pressure which is proportional to the square of the Alfvén Mach number of the shear flow. The out‐of‐plane flow shear increases the energy conversion rate , and for the same upstream conditions, the magnitude of is larger at the dusk side than at the dawn side. This study reveals that out‐of‐plane flow shear not only enhances the reconnection rate but also significantly boosts energy conversion, with more pronounced effects on the dusk‐side flank than on the dawn‐side flank. These insights pave the way for better understanding the solar wind‐magnetosphere interactions.

Liang, Haoming↗
DOE OSTI2026

Time-resolved X-ray imaging of the current filamentation instability in solid-density plasmas

The streaming of energetic charged particles can magnetize astrophysical and laboratory plasmas via the current filamentation instability. Despite its importance, the experimental characterization of this instability has remained a challenge. Here, we report an experiment combining a high-intensity optical laser with a high-brightness X-ray free electron laser that successfully images the instability in solid-density plasmas with 200 nm spatial and 50 fs temporal resolution. We characterize the development of μm-scale filamentary structures and their evolution over tens of picoseconds through a non-linear merging process. The measured plasma density modulations and long merging time reveal the critical importance of space-charge effects and ion motion on this electron-driven instability. Supporting theoretical analysis and kinetic simulations help distinguish the relative role of space-charge and resistive effects. Our findings indicate that magnetic fields on the order of 10 megagauss are produced, with important implications for transport and radiation emission of energetic particles in plasmas.

Schoenwaelder, Christopher [SLAC National Accelera↗
DOE OSTI2025

An analytical model of “Electron-Only” magnetic reconnection rates

Abstract “Electron-only” reconnection, which is both uncoupled from the surrounding ions and much faster than standard reconnection, is arguably ubiquitous in turbulence. One critical step to understanding the rate in this novel regime is to model the outflow speed that limits the transport of the magnetic flux, which is super ion Alfvénic but significantly lower than the electron Alfvén speed based on the asymptotic reconnecting field. Here we develop a simple model to determine this limiting speed by taking into account the multiscale nature of reconnection, the Hall-mediated electron outflow speed, and the pressure buildup within the small system. The predicted scalings of rates and various key quantities compare well with fully kinetic simulations and can be useful for interpreting the observations of NASA’s Magnetospheric-Multiscale (MMS) mission and other ongoing missions.

Physics↗
DOE OSTI2025

The role of the droplet interface in controlling the multiphase oxidation of thiosulfate by ozone

Predicting reaction kinetics in aqueous microdroplets, including aerosols and cloud droplets, is challenging due to the probability that the underlying reaction mechanism can occur both at the surface and in the interior of the droplet. Additionally, few studies directly measure the surface activities of doubly charged anions, despite their prevalence in the atmosphere. Here, deep-UV second harmonic generation spectroscopy is used to probe surface affinities of the doubly charged anions thiosulfate, sulfate, and sulfite, key species in the thiosulfate ozonation reaction mechanism. Thiosulfate has an appreciable surface affinity with a measured Gibbs free energy of adsorption of -7.3 ± 2.5 kJ mol -1 in neutral solution, while sulfate and sulfite exhibit negligible surface propensity. The Gibbs free energy is combined with data from liquid flat jet ambient pressure X-ray photoelectron spectroscopy to constrain the concentration of thiosulfate at the surface in our model. Stochastic kinetic simulations leveraging these novel measurements show that the primary reaction between thiosulfate and ozone occurs at the interface and in the bulk, with the contribution of the interface decreasing from ~65% at pH 5 to ~45% at pH 13. Additionally, sulfate, the major product of thiosulfate ozonation and an important species in atmospheric processes, can be produced by two different pathways at pH 5, one with a contribution from the interface of >70% and the other occurring predominantly in the bulk (>98%). The observations in this work have implications for mining wastewater remediation, atmospheric chemistry, and understanding other complex reaction mechanisms in multiphase environments. Future interfacial or microdroplet/aerosol chemistry studies should carefully consider the role of both surface and bulk chemistry.

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

Effects of wave damping and finite perpendicular scale on three-dimensional Alfvén wave parametric decay in low-beta plasmas

Shear Alfvén wave parametric decay instability (PDI) provides a potential path toward significant wave dissipation and plasma heating. However, fundamental questions regarding how PDI is excited in a realistic three-dimensional (3D) open system and how the finite perpendicular wave scale—as found in both laboratory and space plasmas—affects the excitation remain poorly understood. Here, we present the first 3D, open-boundary, hybrid kinetic-fluid simulations of kinetic Alfvén wave PDI in low-beta plasmas. Key findings are that the PDI excitation is strongly limited by the wave damping present, including electron–ion collisional damping (represented by a constant resistivity) and geometrical attenuation associated with the finite-scale Alfvén wave, and ion Landau damping of the child acoustic wave. The perpendicular wave scale alone, however, plays no discernible role: waves of different perpendicular scales exhibit similar instability excitation as long as the magnitude of the parallel ponderomotive force remains unchanged. These findings are corroborated by theoretical analysis and estimates. This new understanding of 3D kinetic Alfvén wave PDI physics is essential for laboratory study of the basic plasma process and may also aid future evaluation of the relevance/role of PDI in low-beta space plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Advances in laser-based bremsstrahlung x-ray sources. I. Optimizing laser-accelerated electrons

In this work, we have performed a suite of kinetic simulations of relativistic laser–plasma interaction under settings relevant to recent and planned experiments on a variety of laser systems. The goal of the study is to illuminate the physics of laser–target coupling and to provide guidance for how to optimize these sources for applications. It is shown that the production of relativistic electrons is maximized when conditions of relativistic induced transparency (RIT) in dense plasmas can be achieved over a large interaction volume at the time of arrival of most intense part of the laser pulse. RIT is shown to enhance both the numbers of relativistic electrons and the energies of the electrons, leading to an increased x-ray dose. A variety of approaches to enhancing laser–target coupling are considered. These include optimizing the effects of low-density pre-plasma (arising either from finite laser pedestal or from the use of foam coatings) and of modifying the laser focusing geometry to reduce effects of filamentation and self-focusing. Evidence of a novel approach to achieving stable laser propagation over distances of tens of micrometers in a plasma gradient is also presented. These conditions coincide with plasma and laser conditions explored in recent experiments on the Omega EP laser system and compare favorably with an analytic criterion for stable laser propagation in relativistically underdense plasma obtained from a nonlinear Wentzel–Kramers–Brillouin analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Neutron generation dynamics inside a MA-class dense plasma focus Z-pinch

Dense plasma focii (DPFs) are appealing as energy efficient sources of short pulses of ions, neutrons, and x rays. The output of these sources is expected to scale with input current (I4), but has been shown to drop at the MA level [S. K. H. Auluck, “On the failure of neutron yield scaling in the dense plasma focus,” Phys. Plasmas 30, 080701 (2023)]. New results on the MegaJOuLe Neutron Imaging Radiography DPF showed neutron yield production in agreement with the input current scaling beyond the previously observed drop. This work provides insight into the pinch formation on a DPF and reports on the two different mechanisms leading to neutron generation inside a DPF using a combination of kinetic simulations and experimental data. A combination of particle-in-cell (PIC) and 1D shock theory results are used to describe the pinch formation and disassembly and the corresponding thermonuclear and beam-target mechanisms. The temporal evolution of the pinch column predicted by the PIC simulations shows qualitative agreement with the experimental data from plasma photon emission as well as temporal neutron pulse shapes. In MJ-class DPFs, both thermonuclear and beam-target mechanisms can occur over the course of the implosion and contribute to the total neutron production. Hence the neutron source size of a DPF will change throughout the implosion. Experimental neutron radiographs show the increase in source size as the pinch breaks apart, in agreement with simulation's prediction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Fast ion stabilization of tilt in large radius FRCs

The field reversed configuration (FRC) has been a curious case in plasma physics research in that early MHD analysis suggested FRCs should be grossly unstable, while experimental results contradicted that prediction. Later, this theory was able to resolve this contradiction by understanding that finite Larmor radius effects largely negated the MHD predictions. Similarly, previous theoretical studies of beam driven FRCs predicted that such system would be unstable to beam driven modes while, again, experimental results indicated the contradiction. In this paper, we reconcile the theoretical understanding of beam driven modes with experimental observations of stability in these systems. By self-consistently capturing fast ion generation from neutral beam injection and its impact on the plasma equilibrium, we show that low amplitude perturbations in the magnetic field, driven by betatron particles, modify the precession frequencies of the betatron particles such that the drive for compressional Alfvén waves in the thermal plasma is reduced. Finally, we are able to demonstrate, for the first time, stable beam driven FRC evolution at high S*/E in 3D kinetic simulations.

Betatrons↗
DOE OSTI2026

Understanding cold electron impact on parallel-propagating whistler chorus waves via moment-based quasilinear theory

Earth's magnetosphere hosts a wide range of collisionless particle populations that interact through various wave-particle processes. Among these, cold electrons, with energies below 100 eV, often dominate the plasma density but remain poorly characterized due to measurement challenges such as spacecraft charging and photoelectron contamination. Understanding the contribution of these cold populations to wave–particle interaction is of significant interest. Recent kinetic simulations identified a secondary drift-driven instability, in which parallel-propagating whistler-mode chorus waves excite oblique electrostatic whistler waves near the resonance cone and Bernstein-mode turbulence. These secondary modes enable a new channel of energy transfer from the parallel-propagating whistler wave to the cold electrons. In this work, we develop a moment-based quasilinear theory of the secondary instabilities to quantify such energy exchange. Our results show that these secondary instabilities persist for a wide range of parameters and, in many cases, lead to nearly complete damping of the primary wave. Such secondary instability might limit the amplitude of parallel-propagating whistler waves in Earth's magnetosphere and might explain why high-amplitude oblique whistler or electron Bernstein waves are rarely observed simultaneously with high-amplitude field-aligned whistler waves in the inner magnetosphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

The role of focusing geometry in MeV x-ray production from petawatt laser–solid interaction

Relativistic laser–plasma interactions provide a compact and flexible route to generating bright, ultrashort pulses of MeV x rays, with applications in high-energy density science, nuclear physics, and radiography. Despite extensive study of intensity scaling in laser-driven electron acceleration, the role of focusing geometry and focal-volume effects in MeV radiation production remains insufficiently understood. Here, we present experimental and kinetic simulation results from the Texas Petawatt Laser (120 J, 140 fs) in which the focusing geometry (f/3 or f/1.5) is varied, while the laser energy and pulse duration are fixed. Experimentally, the f/3 geometry produces approximately three times more MeV radiation than the f/1.5 geometry, despite its twofold lower nominal vacuum intensity. This result is based on the time-integrated radiation per unit solid angle, as measured along the diagnostic line of sight. Three-dimensional particle-in-cell simulations reproduce this trend when a modest (10s of μm) effective focal plane shift is introduced, demonstrating that relativistic laser–plasma coupling is highly sensitive to focal geometry in the presence of a preplasma. This behavior is consistent with differences in interaction length and effective intensity at the critical surface between the two configurations. The sensitivity of the tightly focused f/1.5 configuration ($z_R ≈$⁠ 5 μ m) reflects the combined influence of thermal lensing and an extended preplasma, while the f/3 geometry (⁠$z_R ≈$ 80 μ m) remains comparatively robust. These results demonstrate a breakdown of conventional intensity scaling and identify focusing geometry as a critical control parameter for MeV electron and x-ray generation at petawatt powers. Implications for laser–plasma-based radiographic facilities are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Tearing Stable Stationary ITER Baseline Operation in DIII-D

We present a tearing stable, stationary and reproducible operating solution for DIII-D plasmas characterized by the ITER Baseline Scenario normalized parameter set and shape. In these plasmas low differential rotation (Δf1,2) between the core and edge is identified as the single root cause of the onset of the primary limiting instability, the 2,1 Neoclassical Tearing Mode (NTM), while a group of other examined parameters have a much weaker impact on the stability. Explanation is offered by drift kinetic simulations, which show a six-fold reduction in the NTM onset threshold due to diminished stabilizing ion polarization currents when the seed island drift frequency ceases in the local plasma frame. Based on the experimental observation and the explanation provided by the theory, a new method is developed to sustain Δf1,2 by modifying the edge neoclassical potential through neutral gas fueling. This method enables stationary ITER baseline operation free of disruptive 2,1 tearing modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Predicted thresholds for RMP ELM suppression access in double-null configurations

Modeling of DIII-D plasmas spanning shapes from single to double-null (DN) reveals new insights into the nature of resonant magnetic perturbation (RMP) conditions necessary for edge-localized mode (ELM) suppression. The suppression of ELMs with RMPs has proven difficult in DN configurations, where no device has thus far reported any hints of suppression. Modeling using the GPEC code finds a reduced high-field side response closer to DN shaping. The resulting synthetic diagnostic measurements are consistent with what has been observed in experiments on DIII-D, validating the plasma response model in this regime. While common metrics for suppression do not illustrate a clear distinction between single-null (SN) and DN cases, the pedestal top resonant field does show a ∼20% decrease at DN shaping in modeling. Field penetration is assessed using linear tearing theory, which demonstrates a lack of sufficient pedestal top resonant flux in the DN shape, requiring at least 1.5 × greater RMP coil currents than what was used in experiment. Analysis from drift kinetic simulations further indicate up to 2 × larger critical island widths are required at the pedestal top for tearing mode growth compared to SN cases. Effective island widths inferred from 3D ideal MHD are also analyzed, where maximum widths in lab coordinates indicate a threshold of ∼18–24× the ion gyroradius for sufficient profile flattening for ELM suppression. These results suggest that ELM suppression may be possible in DN with sufficiently large RMP coil amplitude. Future prospects of achieving RMP ELM suppression in the DN configuration may involve going to lower triangularity, which is also highlighted in this work.

RMP ELM suppression↗
DOE OSTI2025

Towards modelling AR Sco: calibration – reproducing high-energy pulsar emission and testing convergence to Aristotelian electrodynamics

In recent years, kinetic simulations have been crucial to further our understanding of pulsar electrodynamics. Yet, due to the large-scale separation between the gyro-period and the stellar rotation period, resolving the particle gyration has been computationally unfeasible for realistic pulsar parameters. The main aim of this work is comparing our gyro-phase-resolved model with a gyro-centric pulsar model, where our model solves the general equations of motion with included radiation reaction using a higher order numerical solver with adaptive time-steps. Specifically, we aim to (i) reproduce a pulsar’s high-energy emission maps, namely one with 10 per cent of the surface B-field strength of Vela, and the spectra produced by an independent gyro-centric pulsar emission model; and (ii) test convergence of these results to the radiation-reaction limit of Aristotelian electrodynamics. (iii) Additionally, we identify the effect that a large $E_{\parallel }$-field has on the trajectories and radiation calculations. We find that we can reproduce the curvature radiation emission maps and spectra well, using 10 per cent field strengths of the Vela pulsar and injecting our particles at a higher altitude in the magnetosphere. Using sufficiently large $E_{\parallel }$-fields, our numeric results converge to the analytic radiation-reaction limit trajectories. Additionally, we illustrate the importance of accounting for the $\mathbf {E}\times \mathbf {B}$-drift in the particle trajectories and radiation calculations, validating the Harding and collaborators’ model approach. Lastly, we found that our model deals very well with the high-radiation-reaction and high-field regimes present in pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Kinetic study of strong shock waves in fully ionized plasmas

Fully kinetic simulations are used in this work to study strong shock waves in fully ionized plasmas for charge state 𝑍 ranging from 1 to 18. The ion stress tensor amplitude agrees well with classical predictions within the compression layer; however, its spatial structure is noticeably shaped by nonlocal transport effects. The electron stress tensor remains significantly smaller than that of the ions for all values of 𝑍 considered. Ion overheating in the compression layer is observed to scale approximately linearly with 𝑍. Both findings diverge from hydrodynamic models based on classical electron stress formulations. In conclusion, the electrons exhibit strong nonlocal behavior in the compression layer, with Knudsen numbers approaching unity for 𝑍 ≥ 6.

Physics - Plasma physics↗
DOE OSTI2025

Combined Influence of Rotation and Scrape-Off Layer Drifts on Recycling Asymmetries in Tokamak Plasmas

Coupled 2D fluid-kinetic simulations of a DIII-D high confinement tokamak plasma show that plasma rotation coupled with drift effects near the plasma edge play a significant role in the creation of the observed poloidal distribution of neutrals. It is observed that including either drift or rotation effects enhances particle flux at the inner target in the case of ion 𝐵×∇𝐵 drift toward the 𝑋-point. However, the particle flux asymmetry is significantly higher with the combination of drifts and rotation than either effect alone. The heightened particle flux asymmetry allows for improved simulation of the strong in-out asymmetry of the Lyman-𝛼 brightness profiles measured in the experiment. Enhancement of radial transport of parallel momentum changes the upstream scrape-off layer flow pattern, increasing the fraction of deuterium flux that reaches the inboard divertor entrance while lowering that which arrives at the outboard. In conclusion, this Letter indicates that by combining drifts, rotation, and viscous coupling, existing boundary plasma models can achieve a satisfactory agreement with experimentally measured neutral asymmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗