Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic plasmas”

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 577 records · Page 32

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↗

bolt: A GPU-Accelerated Solver for Marginally Collisional Kinetic Physics Using a High- Resolution Constrained Transport Scheme

Collisionless and marginally collisional plasmas are very expensive to simulate due to high dimensionality and large scale separations, especially in global problems with additional length scales. bolt addresses this challenge by achieving high performance on GPUs using the arrayfire performance portability library. We have demonstrated accuracy on a range of test problems, and the bolt code paper is currently in preparation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enabling Low-Temperature (LTP) Ignition Technologies for Multi-Mode Engines through the Development of a Validated High-Fidelity LTP Model for Predicative Simulations Tools

The goal of multi-mode engine architectures is to extend current lean-burn dilution limits with renewable fuels, which requires spark plugs to deposit high energies (hundreds of mJ) in order to initiate ignition and complete combustion. At elevated energy deposition rates, spark plugs experience increased electrode erosion and thermal losses, which ultimately shortens the spark-plug lifetime and lowers ignition efficiency. As such, in order to safeguard the efficiency gains of multi-mode concepts, new and improved ignition technologies are required. Recently, non-equilibrium low-temperature plasmas (LTP) have been shown to promote energy-efficient ignition via quenching and transport of electronically excited atoms and molecules, selective radical production and fast heating of hydrocarbon/air mixtures [1-2]. Thus, LTP is seen as a technology that can potentially improve the energy extraction efficiency of fuels, while enabling kinetically controlled combustion modes towards fuel leaner conditions to realize current DOE VTO goals of improving the sustainability of future mobility [3]. Although many previous studies have demonstrated the efficacy of plasma-assisted ignition to enhance combustion, the detailed enhancement mechanisms remain largely unknown, especially for oxygenated fuels and at elevated pressures that are most relevant to practical engine conditions. These barriers hinder the development of accurate and comprehensive numerical models that seek to describe LTP-based ignition in existing engine design software tools and methods. Current state-of-the-art simulation capabilities for LTP ignition systems are in need of improvements since they deliver qualitative results only due to important limitations of existing approaches. Firstly, validated kinetic models with elementary steps for plasma discharges in oxygenated fuel/air mixtures of relevance to the transportation sector are required. Such kinetic models do not exist at present and will be developed and validated within this project. Secondly, plasma discharges and reactive mixture ignition are multi-scale, unsteady processes requiring high-performance numerical methods and software that execute efficiently on DOE supercomputers. Such software does not exist at present and will be developed and applied to practical LTP ignition scenarios as part of this project. Thirdly, experimental databases that are tailored to serve as benchmark in support of the development of predictive computational models of LTP ignition do not exist and will be part of this project.

33 ADVANCED PROPULSION SYSTEMS↗

Measurements and kinetic modeling of O 2 vibrational kinetics in O 2 –Ar mixtures partially dissociated by a Ns pulse discharge

Vibrational kinetics of O 2 is studied during the O atom recombination in an O 2 –Ar mixture, partially dissociated by a burst of ns discharge pulses in a heated plasma flow reactor. The time-resolved temperature in the discharge afterglow is determined by Rayleigh scattering. Time-resolved O atom number density is measured by ps Two-Photon absorption Laser Induced Fluorescence, calibrated in xenon. Time-resolved vibrational level populations of molecular oxygen, O 2 (v= 8–20), are measured by ps Laser Induced Fluorescence (LIF), with the absolute calibration by NO LIF. Time-resolved ozone number density is monitored by broadband UV absorption. The results are compared with the predictions of a state-specific kinetic model. The experimental data indicate a rapid initial decay of O 2 (v) populations generated by electron impact in the discharge, due to the vibration-translation (V–T) relaxation by O atoms. This is followed by a slower population reduction, on the time scale much longer compared to that for V–T relaxation or vibration-vibration (V–V) exchange. Both O atoms and the O 2 (v) populations decay on the same time scale, indicating that chemical reactions initiated by the O atom recombination result in the generation of vibrationally excited O 2 molecules. These trends are reproduced by the kinetic model, which shows that the reaction of O atoms with ozone is the dominant pathway of O 2 (v) generation at the present conditions. The predicted relative O 2 (v) populations are close to the experimental results, but absolute number densities differ from the experimental data. This is likely due to uncertainties in the absolute calibration of LIF measurements and in the spectroscopic model used in the data reduction. The present work demonstrates the capability for the absolute, time-resolved measurements of vibrationally excited O 2 in recombining gas flows, to quantify the energy partition in the recombination reactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Critical magnetic Reynolds number of the turbulent dynamo in collisionless plasmas

The intracluster medium of galaxy clusters is an extremely hot and diffuse, nearly collisionless plasma, which hosts dynamically important magnetic fields of ~μG strength. Seed magnetic fields of much weaker strength of astrophysical or primordial origin can be present in the intracluster medium. In collisional plasmas, which can be approximated in the magnetohydrodynamical (MHD) limit, the turbulent dynamo mechanism can amplify weak seed fields to strong dynamical levels efficiently by converting turbulent kinetic energy into magnetic energy. However, the viability of this mechanism in weakly collisional or completely collisionless plasma is much less understood. In this study, we explore the properties of the collisionless turbulent dynamo using three-dimensional hybrid-kinetic particle-in-cell simulations. We explore the properties of the collisionless turbulent dynamo in the kinematic regime for different values of the magnetic Reynolds number, Rm, initial magnetic-to-kinetic energy ratio, (E mag /E kin ) i , and initial Larmor ratio, (r Larmor /L box ) i , i.e. the ratio of the Larmor radius to the size of the turbulent system. We find that in the ‘un-magnetized’ regime, (r Larmor /L box ) i > 1, the critical magnetic Reynolds number for the dynamo action Rm crit ≈ 107 ± 3. In the ‘magnetized’ regime, (r Larmor /L box ) i ≲ 1, we find a marginally higher Rm crit = 124 ± 8. We find that the growth rate of the magnetic energy does not depend on the strength of the seed magnetic field when the initial magnetization is fixed. We also study the distribution and evolution of the pressure anisotropy in the collisionless plasma and compare our results with the MHD turbulent dynamo.

79 ASTRONOMY AND ASTROPHYSICS↗

Multi-messenger dynamic imaging of laser-driven shocks in water using a plasma wakefield accelerator

Understanding dense matter hydrodynamics is critical for predicting plasma behavior in environments relevant to laser-driven inertial confinement fusion. Traditional diagnostic sources face limitations in brightness, spatiotemporal resolution, and in their ability to detect relevant electromagnetic fields. In this work, we present a dual-probe, multi-messenger laser wakefield accelerator platform combining ultrafast X-rays and relativistic electron beams at 1 Hz, to interrogate a free-flowing water target in vacuum, heated by an intense 200 ps laser pulse. This scheme enables high-repetition-rate tracking the evolution of the interaction using both particle types. Betatron X-rays reveal a cylindrically symmetric shock compression morphology assisted by low-density vapor, resembling foam-layer-assisted fusion targets. The synchronized electron beam detects time-evolving electromagnetic fields, uncovering charge separation and ion species differentiation during plasma expansion – phenomena not captured by photons or hydrodynamic simulations. We show that combining both probes provides complementary insights spanning kinetic to hydrodynamic regimes, highlighting the need for hybrid physics models to accurately predict fusion-relevant plasma behavior.

imaging techniques↗

Effect of energetic ions on edge-localized modes in tokamak plasmas

The most efficient and promising operational regime for the International Thermonuclear Experimental Reactor tokamak is the high-confinement mode. In this regime, however, periodic relaxations of the plasma edge can occur. These edge-localized modes pose a threat to the integrity of the fusion device. Here we reveal the strong impact of energetic ions on the spatio-temporal structure of edge-localized modes in tokamaks using nonlinear hybrid kinetic–magnetohydrodynamic simulations. A resonant interaction between the fast ions at the plasma edge and the electromagnetic perturbations from the edge-localized mode leads to an energy and momentum exchange. Energetic ions modify, for example, the amplitude, frequency spectrum and crash timing of edge-localized modes. The simulations reproduce some observations that feature abrupt and large edge-localized mode crashes. The results indicate that, in the International Thermonuclear Experimental Reactor, a strong interaction between the fusion-born alpha particles and ions from neutral beam injection, a main heating and fast particle source, is expected with predicted edge-localized mode perturbations. This work advances the understanding of the physics underlying edge-localized mode crashes in the presence of energetic particles and highlights the importance of including energetic ion kinetic effects in the optimization of edge-localized mode control techniques and regimes that are free of such modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The effects of pre-plasma scale length on the relativistic electron beam directionality

The effects of pre-plasmas on the electron beam directionality was experimentally and numerically investigated. Single material and layered targets made of Ti and/or CH were used to simultaneously measure high-energy (≥3 MeV) electrons along two directions, pre-pulse energy and pre-plasma density. The electron directionality is quantified by using a new parameter, the electron energy ratio of the total kinetic energies along the two directions. Measurements and radiation–hydrodynamic (RH) simulations show that a large (≥3.5 μm) plasma scale length at the critical surface enhances electrons along the laser axis, and such pre-plasma conditions could only be achieved with the CH targets. Particle-in-cell simulations were performed on the RH generated pre-plasmas from Ti and CH targets, and the results show that the CH target provided conditions for higher forward momentum gains by electrons. First, the CH target allowed longer distances for electrons to interact with laser. Second, the intense laser pulse modified the critical surface, but the resulting surface differed. The CH target resulted in a smooth surface where a retro-reflection was observed while the Ti target resulted in a rippled surface that scattered the reflected light. As results, the CH electrons gained higher forward momentum via a direct-laser-acceleration in the counter propagating laser fields. So the results presented in this article show a way of controlling the high-energy electron directionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Direct evidence of the effect of a moderate external magnetic field on stimulated Raman scattering in the kinetic regime

We present results from an experiment carried out at the OMEGA-EP laser facility that investigated the effect of a perpendicular magnetic field on stimulated Raman scattering (SRS) and report the first direct measurement of magnetic mitigation of SRS. A 13-T magnetic field generated by pulsed-power coils was imposed on a gas jet plasma, and a novel three-picket interaction beam was used to explore SRS reflectivity for several plasma conditions within single shots. The time-resolved backscattered light shows that SRS was mitigated by the external 13-T magnetic field in the kinetic regime (kλ D ∼ 0.3, where k is the electron plasma wave's wavenumber and λ D is electron Debye length) and at a density of n e /n cr ∼ 0.10. On the other hand, we also measured an enhancement of SRS reflectivity at lower density (n e /n cr < 0.08). We discuss experimental results in the context of magnetohydrodynamic and particle-in-cell simulations that scan SRS dynamics for a variety of plasma conditions. While the experimental evidence of SRS mitigation validates prior work on the kinetic simulation of SRS in an external magnetic field, we also find that other mechanisms such as SRS rescatter can lead to an enhancement of measured SRS reflectivity in simulations of parameters relevant to our experiment.

Gas jet↗

Radial electric field and density fluctuations measured by Doppler reflectometry during the post-pellet enhanced confinement phase in W7-X

Radial profiles of density fluctuations and the radial electric field, Er, have been measured using Doppler reflectometry during the post-pellet enhanced confinement phase achieved, under different heating power levels and magnetic configurations, during the 2018 W7-X experimental campaign. A pronounced Er-well is measured with local values as high as -40 kV m -1 in the radial range ρ ~ 0.7–0.8 during the post-pellet enhanced confinement phase. The maximum Er intensity scales with both the plasma density and electron cyclotron heating power level, following a similar trend to the plasma energy content. A good agreement is found when the experimental Er profiles are compared to simulations carried out using the neoclassical codes, the drift kinetic equation solver (DKES) and kinetic orbit-averaging solver for stellarators (KNOSOS). The density fluctuation level decreases from the plasma edge toward the plasma core and the drop is more pronounced in the post-pellet enhanced confinement phase than in reference gas-fuelled plasmas. Besides, in the post-pellet phase, the density fluctuation level is lower in the high iota magnetic configuration than in the standard one. Therefore, to determine whether this difference is related to the differences in the plasma profiles or to the stability properties of the two configurations, gyrokinetic simulations have been carried out using the codes stella and EUTERPE. The simulation results point to the plasma profile evolution after the pellet injection and the stabilization effect of the radial electric field profile as the dominant players in the stabilization of the plasma turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Beta-dependent Properties of Solar Wind Current Sheets at 1 au

We utilize plasma turbulence in interplanetary coronal mass ejections as a natural laboratory to analyze more than 4600 current sheets observed over a broad range of electron and proton beta values, 10 −2 ≲ β e ≲ 10 and 10 −3 ≲ β p ≲ 10. The analysis shows that the shear angle Δθ and the thickness λ CS of the current sheets depend on both electron and proton beta. Specifically, the dependencies on electron beta are Δθ ≈ 12.6° $β^{0.41}_{e}$ and, λ CS /λ p ≈ 1.9 $β^{0.26}_{e}$ where λ p is the local proton inertial length. We argue that the beta dependence of the shear angle is an intrinsic feature of solar wind turbulence arising from the natural correlation between turbulence intensity and plasma beta. In contrast, the beta dependence of current sheet scales is likely a general characteristic of plasma turbulence. According to recent theory, current sheets formed in turbulence are disrupted by the electron tearing instability once their thickness is smaller than a critical scale. We demonstrate that normalizing the thickness of current sheets by this critical scale eliminates the dependence on beta across the entire range of plasma beta. This observation not only supports the theory suggesting that current sheets mediate the transition from the inertial to the kinetic cascade through the electron tearing instability, but also demonstrates the plausibility of similar processes in astrophysical plasma. The presented results also indicate that the theory may be applicable beyond its formal limits, m e /m i ≪ β ≪ 1.

79 ASTRONOMY AND ASTROPHYSICS↗

Benchmark for two-dimensional large scale coherent structures in partially magnetized E × B plasmas—community collaboration & lessons learned

Low-temperature plasmas (LTPs) are essential to both fundamental scientific research and critical industrial applications. As in many areas of science, numerical simulations have become a vital tool for uncovering new physical phenomena and guiding technological development. Code benchmarking remains crucial for verifying implementations and evaluating performance. This work continues the Landmark benchmark initiative, a series specifically designed to support the verification of LTP codes. In this study, seventeen simulation codes from a collaborative community of nineteen international institutions modeled a partially magnetized E × B Penning discharge. The emergence of large scale coherent structures, or rotating plasma spokes, endows this configuration with an enormous range of time scales, making it particularly challenging to simulate. The codes showed excellent agreement on the rotation frequency of the spoke as well as key plasma properties, including time-averaged ion density, plasma potential, and electron temperature profiles. Achieving this level of agreement came with challenges, and we share lessons learned on how to conduct future benchmarking campaigns. Comparing code implementations, computational hardware, and simulation runtimes also revealed interesting trends, which are summarized with the aim of guiding future plasma simulation software development.

benchmarking↗

Thermal quench of open field plasma intercepting with recycling walls

When a fusion plasma suddenly intercepts a solid surface (wall or pellet), thermal collapse is distinctly kinetic & has novel physics. VPIC simulations & theory revealed the fascinating dynamics of four propagating fronts controlling parallel electron temperature cooling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fluid-kinetic modeling of a high power density radio frequency inductively coupled positive hydrogen ion source

High power density radio-frequency (RF) inductively coupled positive ion sources are attractive candidates for next-generation neutral beam injection (NBI) systems, where higher injected power and longer pulse lengths are desired without sacrificing source reliability. Operating at absorbed power densities of order $\gt 1~\mathrm{W\,cm}^{-3}$ places these sources in a regime with stronger gas heating, higher dissociation, and non-Maxwellian electron energy distributions. The Large Uniform Plasma for Ionizing Neutrals (LUPIN) is an RF inductively coupled plasma source designed to explore this high power density regime and to provide guidance for a positive ion source upgrade for the DIII-D NBI system. LUPIN is designed to operate at up to 20 kW of RF power at 2 MHz, coupling energy through a cylindrical quartz vessel to achieve target ion current densities of $2100\,\mathrm{A\,m}^{-2}$ . This paper presents fluid-kinetic modeling of LUPIN using the hybrid plasma equipment model where electrons are treated kinetically, and the simulations reveal that electron energy distribution function transitions from nearly Maxwellian in the core to bi-Maxwellian towards the edge. Parametric simulations investigate the effects of RF power, gas pressure, and frequency on plasma density, ion flux, and uniformity. Parametric sweeps reveal that increasing power shifts the primary ionization channel from molecular to atomic with diminishing flux gains due to skin-depth contraction and gas rarefaction. Higher frequency localizes heating and increases $\mathrm{H}_2^+$ and $\mathrm{H}_3^+$ delivery to the grid, while elevated pressure boosts ionization yet hinders ion transport due to increase in collisionality.

inductively coupled plasma↗

On the local measurement of electric currents and magnetic fields using Thomson scattering in Weibel-unstable plasmas

Here, we demonstrate the capability of the Thomson Scattering (TS) diagnostic to measure locally the microscopic electron and ion currents in counter-streaming plasmas unstable to the Weibel or current-filamentation instability. Synthetic TS spectra are calculated with particle distribution functions obtained from particle-in-cell simulations and used to accurately reproduce the simulated currents. We show that this technique allows accurate local measurements of the magnetic field, thus opening the way for the complete experimental characterization of the growth rate, saturation, and nonlinear dynamics of electromagnetic instabilities in plasmas. We illustrate the application of this diagnostic to experimental TS data, which yields local measurements of the magnetic field in Weibel-unstable plasmas and indicates that the magnetic energy density reaches ~1% of the kinetic energy density of the flows, in agreement with previous numerical studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Extreme Heating of Minor Ions in Imbalanced Solar-wind Turbulence

Minor ions in the solar corona are heated to extreme temperatures, far in excess of those of the electrons and protons that comprise the bulk of the plasma. These highly nonthermal distributions make minor ions sensitive probes of the collisionless processes that heat the corona and power the solar wind. The recent discovery of the "helicity barrier" offers a mechanism in which imbalanced Alfvénic turbulence in low-β plasmas preferentially heats protons over electrons, generating high-frequency, proton-cyclotron-resonant fluctuations. We use the hybrid-kinetic particle-in-cell code Pegasus++ to drive imbalanced Alfvénic turbulence in a 3D low-β plasma with additional passive ion species, He 2+ and O 5+ . A helicity barrier naturally develops, followed by clear phase-space signatures of oblique proton-cyclotron-wave heating and Landau-resonant heating from the imbalanced Alfvénic fluctuations. The former results in characteristically arced ion velocity distribution functions, whose non-bi-Maxwellian features are shown by linear ALPS calculations to be critical to the heating process. Additional features include a steep transition-range electromagnetic spectrum, proton-cyclotron waves propagating in the direction of the imbalance, significantly enhanced proton-to-electron heating ratios, ion temperatures that are considerably more perpendicular with respect to magnetic field, and extreme heating of heavier species in a manner consistent with mass scalings inferred from spacecraft measurements. None of these features are realized in an otherwise equivalent simulation of balanced turbulence. If seen simultaneously in the fast solar wind, these signatures of the helicity barrier would testify to the necessity of incorporating turbulence imbalance in a complete theory for the evolution of the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗

Electron Dynamics Within a MITL Containing a Load

In this paper, we derive the vacuum electric fields within specific cylindrically symmetric magnetically insulated transmission lines (MITLs) in the limit of an infinite speed of light for an arbitrary time-dependent current. We focus our attention on two types of MITLs: the radial MITL and a spherically curved MITL. We then simulate the motion of charged particles, such as electrons, present in these MITLs due to the vacuum fields. In general, the motion of charged particles due to the vacuum fields is highly nonlinear since the fields are nonlinear functions of spatial coordinates and depend on an arbitrary time-dependent current drive. Using guiding center theory, however, one can describe the gross particle kinetics using a combination of $\textbf {E} \times \textbf {B}$ and $\nabla B$ drifts. In addition, we compare our approximate inner MITL field models and particle kinetics with those from a fully electromagnetic simulation code. We find that the agreement between the approximate model and the electromagnetic simulations is excellent.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Decoupling plasma, catalyst, and gaseous mechanisms for non-oxidative methane conversion

Direct non-oxidative methane (CH 4 ) conversion to value-added hydrogen (H 2 ) and C 2 products remains hindered by fundamental catalytic scaling constraints and rapid surface deactivation at elevated temperatures. Plasma-enabled catalysis offers a promising route to overcome the thermodynamic and kinetic barriers of direct non-oxidative methane upgrading at mild conditions, yet control over C–C product selectivity and catalyst stability remains elusive. Here, we establish a unified mechanistic framework including Langmuir–Hinshelwood (L–H) and Langmuir–Rideal (L–R) mechanisms that disentangles the roles of plasma excitation (including vibrationally activated species and radicals), surface temperature (T sur ), and catalyst binding energy in steering CH 4 conversion to H 2 and C 2 hydrocarbons. Through a combination of density functional theory (DFT) informed microkinetic modeling, in situ and ex situ surface characterization, and product quantification under dielectric barrier discharge conditions, we show that vibrationally excited CH 4 lowers activation barriers selectively for dissociative chemisorption, enabling surface activation across a wide range of transition metal catalysts at low thermal energy input. We find that once CH 4 is dissociatively chemisorbed, the branching between C 2 H 2 , C 2 H 4 , and C 2 H 6 is governed by surface properties (carbon binding energy, T sur , etc), regardless of plasma excitation. The DFT informed microkinetic model decouples the effects of molecular activation from surface properties and indentifies operating windows that maximize target yields while suppressing carbon accumulation and subsequent catalytic inactivation. Experiments on polycrystalline Cu/Al 2 O 3 , Ni/Al 2 O 3 , and Pt/Al 2 O 3 validate these predictions, revealing catalyst-dependent branching toward ethane or ethylene and distinct deactivation profiles. We unify these trends into a generalized three-dimensional plasma-thermal-catalytic design space, from which reduced descriptors such as T vib /T sur in the limit of vibrationally excited L–H pathways emerge as predictive metrics. These results enable rational tuning of methane conversion pathways and unlock selective C 2 formation using earth-abundant metals under mild plasma conditions.

catalyst inactivation↗