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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Bohm Criterion of Plasma Sheaths away from Asymptotic Limits

The plasma exit flow speed at the sheath entrance is constrained by the Bohm criterion. The so-called Bohm speed regulates the plasma particle and power exhaust fluxes to the wall, and it is commonly deployed as a boundary condition to exclude the sheath region in quasineutral plasma modeling. In this work, the Bohm criterion analysis is performed in the intermediate plasma regime away from the previously known limiting cases of adiabatic laws and the asymptotic limit of infinitesimal Debye length in a finite-size system, using the transport equations of an anisotropic plasma. The resulting Bohm speed has explicit dependence on local plasma heat flux, temperature isotropization, and thermal force. Comparison with kinetic simulations demonstrates its accuracy over the plasma-sheath transition region in which quasineutrality is weakly perturbed and the Bohm criterion applies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Neutralization of ion beam by electron injection: Accumulation of cold electrons

Ion beam charge neutralization by electron injection is a complex kinetic process. Recent experiments show that the resulting self-potential of the ion beam after neutralization by plasma is much lower than the temperature of plasma electrons [Stepanov et al., Phys. Plasmas 23, 043113 (2016)], indicating that kinetic effects are important and may affect the neutralization of the ion beam. We performed a numerical study of the charge neutralization process of an ion beam making use of a two-dimensional electrostatic particle-in-cell code. The results show that the process of charge neutralization by electron injection is composed of two stages. During the first stage, the self-potential of the beam is higher than the temperature of injected electrons (T e /e) and all injected electrons are captured by the ion beam. During the second stage, hot electrons escape from the ion beam and the beam self-potential (φ) decreases because cold electrons slowly accumulate resulting in the beam self-potential φ to become much lower than Te/e in agreement with previous experimental observations at Princeton Advanced Teststand. We also determined that the resulting φ scales as φ~√T e , in agreement with previous experimental observations from Gabovich's group. In addition, the results show that the transverse position of the electron source has a great impact on ion beam neutralization. A slight shift of the electron source as relevant to the ion thrusters leads to a large increase in the beam self-potential because of an increase in potential energy of injected electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

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

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

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Production and Persistence of Extreme Two-temperature Plasmas in Radiative Relativistic Turbulence

Turbulence is a predominant process for energizing electrons and ions in collisionless astrophysical plasmas, and thus is responsible for shaping their radiative signatures (luminosity, spectra, and variability). To better understand the kinetic properties of a collisionless radiative plasma subject to externally driven turbulence, we investigate particle-in-cell simulations of relativistic plasma turbulence with external inverse Compton cooling acting on the electrons. We find that ions continuously heat up while electrons gradually cool down (due to the net effect of radiation), and hence the ion-to-electron temperature ratio T i /T e grows in time. Furthermore, we show that Ti/Te is limited only by the size and duration of the simulations (reaching ${T}_{i}/{T}_{e}\sim {10}^{3}$), indicating that there are no efficient collisionless mechanisms of electron–ion thermal coupling. This result has implications for models of radiatively inefficient accretion flows, such as observed in the Galactic center and in M87, for which so-called two-temperature plasmas with ${T}_{i}/{T}_{e}\gg 1$ have been invoked to explain their low luminosity. Additionally, we find that electrons acquire a quasi-thermal distribution (dictated by the competition of turbulent particle energization and radiative cooling), while ions undergo efficient nonthermal acceleration (acquiring a harder distribution than in equivalent nonradiative simulations). There is a modest nonthermal population of high-energy electrons that are beamed intermittently in space, time, and direction; these beamed electrons may explain rapid flares in certain high-energy astrophysical systems (e.g., in the Galactic center). These numerical results demonstrate that extreme two-temperature plasmas can be produced and maintained by relativistic radiative turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An investigation of shock formation vs shock mitigation of colliding plasma jets

Here, this work studies the interaction between colliding plasma jets to understand regimes in which jet merging results in shock formation vs regimes in which the shock formation is mitigated due to the collisionless interpenetration of the jets. A kinetic model is required for this study because fluid models will always produce a shock upon the collision of plasma jets. The continuum-kinetic, Vlasov–Maxwell–Dougherty model with one velocity dimension is used to accurately capture shock heating, along with a novel coupling with a moment equation to evolve perpendicular temperature for computational efficiency. As a result, this relatively inexpensive simulation can be used for detailed scans of the parameter space toward predictions of shocked vs shock-mitigated regimes, which is of interest for several fusion concepts such as plasma-jet-driven magneto-inertial fusion, high-energy-density plasmas, astrophysical phenomena, and other laboratory plasmas. The initial results obtained using this approach are in agreement with the preliminary outcomes of the plasma liner experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Higher Harmonics in Multipactor Induced Plasma Ionization Breakdown near a Dielectric Surface

In this Letter, the novel physics of higher harmonic (HH) generation of the normal electric field near a dielectric surface is reported in multipactor induced plasma ionization breakdown, as determined by kinetic particle-in-cell simulations. The observed HH frequency is around ten times the fundamental rf driving frequency, but lower than the electron plasma frequency. A theory is constructed which indicates that stream plasma interaction-induced instability is the mechanism of HH generation in the collisional regime. The HH frequency and its corresponding growth rate of the HH oscillation amplitude from the theory are in good agreement with kinetic particle-in-cell simulations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

The scaling of kinetic and transport behaviors in the solution-phase chemistry of a plasma–liquid interface

The reactions at a plasma–liquid interface often involve species such as the solvated electron or the hydroxyl radical, which initiate the reduction or oxidation of solution-phase reactants (so-called scavengers) or are consumed by their own second-order recombination. Here, the mathematical scaling of the reaction–diffusion equations at the interface is used to obtain a characteristic time that can be used to determine the transition from highly efficient scavenger reduction or oxidation to lower efficiencies due to transport limitations. The characteristic time (tc) is validated using numerical solutions of the reaction–diffusion equations. When the scavenger kinetics are faster than second-order recombination, this characteristic transition time scales proportionally with the scavenger diffusivity (Ds) and the square of the scavenger bulk concentration (SB) and inversely proportional to the electron flux (J) squared; that is, t c = D s S B 2 F 2 /J 2 , where F is Faraday’s constant. However, when the scavenger kinetics are comparable or slower than second-order recombination, this scaling does not hold. Extending this analysis to three-dimensional systems shows that the profile of the electron flux on the surface affects the spatial location where reactions are most effective. Finally, the assessment of the implications of these behaviors for the reactor design highlights how effectively controlling the electron flux and solution transport may be necessary to improve the efficiency of scavenger reactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Low-pressure inductively coupled plasmas in hydrogen: impact of gas heating on the spatial distribution of atomic hydrogen and vibrationally excited states

Non-equilibrium inductively coupled plasmas (ICPs) operating in hydrogen are of significant interest for applications including large-area materials processing. Increasing control of spatial gas heating, which drives the formation of neutral species density gradients and the rate of gas-temperature-dependent reactions, is critical. In this study, we use 2D fluid-kinetic simulations with the Hybrid Plasma Equipment Model to investigate the spatially resolved production of atomic hydrogen in a low-pressure planar ICP operating in pure hydrogen (10–20 Pa or 0.075–0.15 Torr, 300 W). The reaction set incorporates self-consistent calculation of the spatially resolved gas temperature and 14 vibrationally excited states. We find that the formation of neutral-gas density gradients, which result from spatially non-uniform electrical power deposition at constant pressure, can drive significant variations in the vibrational distribution function and density of atomic hydrogen when gas heating is spatially resolved. This highlights the significance of spatial gas heating on the production of reactive species in relatively high-power-density plasma processing sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Relativistic Alfvén Turbulence at Kinetic Scales

Abstract In a strongly magnetized, magnetically dominated relativistic plasma, Alfvénic turbulence can extend to scales much smaller than the particle inertial scales. It leads to an energy cascade somewhat analogous to inertial- or kinetic-Alfvén turbulent cascades existing in nonrelativistic space and astrophysical plasmas. Based on phenomenological modeling and particle-in-cell numerical simulations, we propose that the energy spectrum of such relativistic kinetic-scale Alfvénic turbulence is close to k −3 or slightly steeper than that due to intermittency corrections or Landau damping. We note the analogy of this spectrum with the Kraichnan spectrum corresponding to the enstrophy cascade in 2D incompressible fluid turbulence. Such turbulence strongly energizes particles in the direction parallel to the background magnetic field, leading to nearly one-dimensional particle momentum distributions. We find that these distributions have universal log-normal statistics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

How transient effects can change the inferred plasma temperatures when using K-shell sulfur line ratios in a time dependent iron sulfide plasma

This paper describes how the steady state atomic kinetics approximation can underestimate the electron temperature determined from K-shell lines in ps-time-scale transient plasmas. In particular, we model the temperature determination of solid FeS targets used in opacity experiments at the Orion laser facility from the ratio of sulfur He-α to Ly-α lines. Such experiments use short-pulse lasers to heat a thin microdot of FeS buried in a plastic target to temperatures of more than 1 keV and densities of approximately 1–2 g/cm 3 . Using atomic kinetics calculations based on a temperature history from a radiation hydrodynamic simulation of the target evolution, the peak temperature inferred from the sulfur line ratios is 1.29 keV at 3.1 ps as compared with the input peak temperature of 1.41 keV at 2.0 ps. There is a time lag of 1.2 ps at the peak, and an overall 0.5 ps time lag in the temporal history of the temperature as the plasma cools over the next 10 ps.

Rilloraza, Francisco↗
DOE OSTI2021

Nonlinear kinetic simulation study of the ion–ion streaming instability in single- and multi-ion species plasmas

The nonlinear evolution of the ion–ion streaming instability (IISI) is studied using numerical techniques novel to this problem that afford direct insight into the evolution of the particle distributions of each species. Here, during the linear phase of the instability, we demonstrate quantitative agreement with linear kinetic theory. Subsequently, the electrostatic field generated by the IISI causes ring-like velocity distributions of ions to form that are both heated and slowed to varying degrees relative to their initial flows. Due to variation in the trapping conditions for ion species of differing charge-to-mass ratio, when flows of multiple species interact, the nonlinear evolution of each species can be starkly different: we show a case where a lighter ion species is completely stopped by a heavier ion species via the IISI alone (i.e., without collisions) and, for the first time, demonstrate how the IISI can introduce a relative flow between ion species that initially have the same flow velocities, thereby separating them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Numerical modeling of plasma assisted deflagration to detonation transition in a microscale channel

Here, this work numerically studies the plasma assisted deflagration to detonation transition (DDT) of H 2 /O 2 mixtures in a microscale channel with detailed chemistry and transport. The results show that the DDT onset time is non-monotonically dependent on the discharge pulse number. The DDT is accelerated with small pulse numbers, whereas retarded with large ones. Two different DDT regimes, respectively at a small and large plasma discharge number, via acoustic choking of the burned gas and plasma-enhanced reactivity gradient without acoustic choking, are observed. Without plasma discharge, pronounced pressure and temperature gradients in front of the flame are generated by acoustic compression after the choking of the burned gas, triggering DDT via autoignition. With small plasma pulse numbers, the plasma-generated species enhance the ignition kinetics and lead to an increased reactivity in the boundary layer. After the choking of the burned gas, the plasma-enhanced reactivity advances the sequence of autoignition near the wall, strengthens ignition-shock wave coupling, and accelerates DDT. However, with a large discharge pulse number, a direct autoignition initiating DDT can occur without the acoustic choking of the burned gas due to the strongly accelerated reactivity and elevated temperature. In this case, DDT onset is retarded because the elevated temperature increases sonic velocity and the increased reactivity accelerates fuel oxidation in front of the flame, decelerating the formation of a leading shock and subsequent pressure buildup ahead of the flame. The present modeling reveals that no matter with or without plasma discharge, DDT is initiated by autoignition in thermal, pressure, and reactivity gradient fields via the Zel'dovich gradient mechanism. The acoustic choking of the burned gas may not be the necessary condition of DDT with strong plasma-enhanced reactivity gradient. This work provides an answer to the experimentally observed non-monotonic DDT onset time by plasma, which provides guidance to control DDT in advanced detonation engines and fire safety of hydrogen-fueled catalytic reactors in microchannels by non-equilibrium plasma discharge.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Higher-order nonequilibrium term: Effective power density quantifying evolution towards or away from local thermodynamic equilibrium

A common approach to assess the nature of energy conversion in a classical fluid or plasma is to compare power densities of the various possible energy conversion mechanisms. A leading research area is quantifying energy conversion for systems that are not in local thermodynamic equilibrium (LTE), as is common in a number of fluid and plasma systems. Here we introduce the “higher-order nonequilibrium term” (HORNET) effective power density, which quantifies the rate of change of departure of a phase space density from LTE. It has dimensions of power density, which allows for quantitative comparisons with standard power densities. We employ particle-in-cell simulations to calculate HORNET during two processes, magnetic reconnection and decaying kinetic turbulence in collisionless magnetized plasmas, that inherently produce non-LTE effects. We investigate the spatial variation of HORNET and the time evolution of its spatial average. By comparing HORNET with power densities describing changes to the internal energy (pressure dilatation, Pi − D , and divergence of the vector heat flux density), we find that HORNET can be a significant fraction of these other measures (8% and 35% for electrons and ions, respectively, for reconnection; up to 67% for both electrons and ions for turbulence), meaning evolution of the system towards or away from LTE can be dynamically important. Applications to numerous plasma phenomena are discussed. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Electrical and Spectroscopic Diagnostics as Real‐Time Metallization Indicators During Hydrogen Plasma Smelting Reduction

This study investigates the hydrogen (H 2 ) plasma reduction process of direct-reduced-iron-grade hematite ore at different arc currents (100–200 A) in an Ar–5% H 2 atmosphere at 0.9 bar. Iron ore samples (10 g) were exposed to a plasma arc, and the reduction/metallization kinetics were analyzed over fixed time intervals. Electrical diagnostics revealed that the arc voltage exhibited takeover-mode oscillations which were suppressed at higher currents due to stronger electromagnetic coupling. The voltage dropped significantly as metallization approached ∼95%, linked to increased electrical conductivity of the metallic iron (Fe) in the ore as well as Fe evaporation into the arc, lowering the arc resistance. A simplified Elenbaas–Heller model supported this explanation and confirmed that Fe vapor concentration enhances plasma conductivity. Optical emission spectroscopy focused on the plasma–metal interface revealed the plasma's optically thick nature, as the primary Fe I 526.95 nm line experienced self-absorption. However, weaker Fe I lines (404.58, 438.35 nm) normalized to Ar I 696.5 nm provided a reliable proxy for metallization. These diagnostics, electrical and spectroscopic, effectively track metallization in real-time during H 2 plasma smelting reduction.

08 HYDROGEN↗
DOE OSTI2025

Structured light approaches in laser-based plasma diagnostics

There is a growing demand for plasma diagnostics suitable for industrial plasma reactors employed in semiconductor nanofabrication, especially relevant to microelectronics and quantum information systems. Such reactors typically have limited optical access and pose considerable diagnostic challenges, including intense background emission, significant thermal loads, and contamination of optical viewports. In this study, we outline research into structured light techniques (laser beams with tailored spatial, temporal, or phase characteristics) that effectively overcome these issues using laser-induced fluorescence (LIF) as an example. The focus of presented diagnostics is on ion kinetics analysis within an industrial plasma source, although this approach is broadly applicable to other plasma systems and diagnostic contexts. We present a confocal LIF implementation using an axicon-generated Bessel annular beam, achieving spatial resolutions of approximately 5 mm at a focal distance of 300 mm, with potential improvements to about 1 mm. This approach matches conventional orthogonal LIF performance but requires only one optical port. Wavelength-modulation LIF employs nonlinear laser wavelength tuning to measure spectral line derivatives, suppressing background emission and enhancing details of spectral line shape. Additionally, we present new results on applying vortex beams (laser beams carrying orbital angular momentum, OAM) for LIF measurements in an industrial plasma device. These measurements enable simultaneous axial and tangential velocity determination using a single laser beam and have been tested with xenon ion transition. Initial quantification of results was performed. Together, these structured-light approaches provide robust, background-resilient, multi-dimensional diagnostics for complex plasma environments.

Romadanov, Ivan [Princeton Plasma Physics Laborato↗
DOE OSTI2020

Precursor selection in hybrid molecular beam epitaxy of alkaline-earth stannates

One of the challenges of oxide molecular beam epitaxy (MBE) is the synthesis of oxides containing metals with high electronegativity (metals that are hard to oxidize). The use of reactive organometallic precursors can potentially address this issue. To investigate the formation of radicals in MBE, we explored three carefully chosen metal-organic precursors of tin for SnO 2 and BaSnO 3 growth: tetramethyltin (TMT), tetraethyltin (TET), and hexamethylditin (HMDT). All three precursors produced single-crystalline, atomically smooth, and epitaxial SnO 2 (101) films on r-Al 2 O 3 (101¯2) in the presence of oxygen plasma. The study of growth kinetics revealed reaction-limited and flux-limited regimes except for TET, which also exhibited a decrease in the deposition rate with increasing temperature above ~800 °C. Contrary to these similarities, the performance of these precursors was dramatically different for BaSnO 3 growth. TMT and TET were ineffective in supplying adequate tin, whereas HMDT yielded phase-pure, stoichiometric BaSnO 3 films. Significantly, HMDT resulted in phase-pure and stoichiometric BaSnO 3 films even without the use of an oxygen plasma (i.e., with molecular oxygen alone). Furthermore, these results are discussed using the ability of HMDT to form tin radicals and therefore assisting with Sn → Sn 4+ oxidation reaction. Structural and electronic transport properties of films grown using HMDT with and without oxygen plasma are compared. This study provides guideline for the choice of precursors that will enable the synthesis of metal oxides containing hard-to-oxidize metals using reactive radicals in MBE.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Energy-angle distribution of the ions in the RF sheath of ICRH antennas

Radio-frequency sheaths forming at the surface of ion cyclotron resonance heating (ICRH) actuators in fusion experiments are associated with enhanced impurity sputtering from ICRH plasma-facing components (PFCs). The minimization of impurity fluxes from the ICRH PFC is a vital task for the usage of ICRH systems. Capturing the ion kinetics is important for an accurate description of plasma–material interactions, because the ion dynamics plays a crucial role in RF sheaths. Here, we present a hybrid particle-in-cell (hPIC) model able to capture the kinetic behavior of the ions. We analyze the kinetic ion energy-angle distributions (IEADs) impacting the RF antenna and its dependence on different plasma and RF sheath parameters. In particular, the IEAD dependence on RF frequency and magnetic field alignment is investigated. Using hPIC, we simulated a case emulating the latest experimental campaign from JET. The simulation showed that under specific plasma and RF parameters, the kinetic motion of ions results in a cusp formation in the IEAD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Experimental Validation of a Kinetic Ballooning Mode in High-Performance High-Bootstrap Current Fraction Fusion Plasmas

We report the observation of a set of coherent high frequency electromagnetic fluctuations that leads to a turbulence induced self-regulating phenomenon in the DIII-D high bootstrap current fraction plasma. The fluctuations have frequency of 130~220kHz, the poloidal wave length and phase velocity are 16~30 m -1 and ~30 km/s, respectively in the outboard midplane with the estimated toroidal mode number n~5- 9. The fluctuations are located in the internal transport barrier (ITB) region at large radius and are experimentally validated to be kinetic ballooning modes (KBM). Furthermore, quasilinear estimation predicts the KBM to be able to drive experimental particle flux and non-negligible thermal flux, suggesting its significant role in regulating the ITB saturation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗