Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “collisional plasma”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The kinetic analog of the pressure–strain interaction

Energy transport in weakly collisional plasma systems is often studied with fluid models and diagnostics. However, the applicability of fluid models is limited when collisions are weak or absent, and using a fluid approach can obscure kinetic processes that provide key insights into the physics of energy transport. Kinetic diagnostics retain all of the information in 3D-3V phase space and thereby reach beyond the insights of fluid models to elucidate the mechanisms responsible for collisionless energy transport. In this work, we derive the Kinetic Pressure–Strain (KPS): a kinetic analog of the pressure–strain interaction, which is the channel between flow energy density and internal energy density in fluid models. Through two case studies of electron Landau damping, we demonstrate that the KPS diagnostic can elucidate kinetic mechanisms that are responsible for energy transport in this channel, just as the related field–particle correlation is known to identify kinetic mechanisms of transport between electromagnetic field energy density and kinetic energy density in particle flows. In addition, we show that resonant electrons play a major role in transferring energy between fluid flows and internal energy during the process of Landau damping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A self-consistent hybrid model of kinetic striations in low-current argon discharges

A self-consistent hybrid model of standing and moving striations was developed for low-current DC discharges in noble gases. We introduced the concept of surface diffusion in phase space ( r , u ) (where u denotes the electron kinetic energy) described by a tensor diffusion in the nonlocal Fokker–Planck kinetic equation for electrons in the collisional plasma. Electrons diffuse along surfaces of constant total energy ε = u - eφ ( r ) between energy jumps in inelastic collisions with atoms. Numerical solutions of the 1d1u kinetic equation for electrons were obtained by two methods and coupled to ion transport and Poisson solver. We studied the dynamics of striation formation in Townsend and glow discharges in argon gas at low discharge currents using a two-level excitation-ionization model and a ‘full-chemistry’ model, which includes stepwise and Penning ionization. Standing striations appeared in Townsend and glow discharges at low currents, and moving striations were obtained for the discharge currents exceeding a critical value. These waves originate at the anode and propagate towards the cathode. We have seen two types of moving striations with the two-level and full-chemistry models, which resemble the s and p striations previously observed in the experiments. Simulations indicate that processes in the anode region could control moving striations in the positive column plasma. The developed model helps clarify the nature of standing and moving striations in DC discharges of noble gases at low discharge currents and low gas pressures.

Physics↗
DOE OSTI2022

Quantum-inspired method for solving the Vlasov-Poisson equations

Kinetic simulations of collisionless (or weakly collisional) plasmas using the Vlasov equation are often infeasible due to high-resolution requirements and the exponential scaling of computational cost with respect to dimension. Recently, it has been proposed that matrix product state (MPS) methods, a quantum-inspired but classical algorithm, can be used to solve partial differential equations with exponential speed-up, provided that the solution can be compressed and efficiently represented as a MPS within some tolerable error threshold. Here, in this work, we explore the practicality of MPS methods for solving the Vlasov-Poisson equations for systems with one coordinate in space and one coordinate in velocity, and find that important features of linear and nonlinear dynamics, such as damping or growth rates and saturation amplitudes, can be captured while compressing the solution significantly. Furthermore, by comparing the performance of different mappings of the distribution functions onto the MPS, we develop an intuition of the MPS representation and its behavior in the context of solving the Vlasov-Poisson equations, which will be useful for extending these methods to higher-dimensional problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Dynamic Phase Alignment in Inertial Alfvén Turbulence

In weakly collisional plasma environments with sufficiently low electron beta, Alfvénic turbulence transforms into inertial Alfvénic turbulence at scales below the electron skin depth, k ⊥ d e ≳ 1. We argue that, in inertial Alfvénic turbulence, both energy and generalized kinetic helicity exhibit direct cascades. We demonstrate that the two cascades are compatible due to the existence of a strong scale dependence of the phase alignment angle between velocity and magnetic field fluctuations, with the phase alignment angle scaling as cos α k ∝ $k$$^{–1}_{⊥}$. The kinetic and magnetic energy spectra scale as ∝ $k$$^{–5/3}_{⊥}$ and ∝ $k$$^{–11/3}_{⊥}$, respectively. As a result of the dual direct cascade, the generalized helicity spectrum scales as ∝ $k$$^{–5/3}_{⊥}$, implying progressive balancing of the turbulence as the cascade proceeds to smaller scales in the k ⊥ d e >> 1 range. Turbulent eddies exhibit a phase-space anisotropy k ∥ ∝ $k$$^{5/3}_{⊥}$, consistent with critically balanced inertial Alfvén fluctuations. Furthermore, our results may be applicable to a variety of geophysical, space, and astrophysical environments, including the Earth’s magnetosheath and ionosphere, solar corona, and nonrelativistic pair plasmas, as well as to strongly rotating nonionized fluids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Bolt: A Fast Solver for Kinetic Theories Using a High-Resolution Constrained Transport Scheme [Slides]

Our understanding of collisionless and semi-collisional plasmas in the nonlinear regime is limited by the expense of computing solutions numerically. Bolt is a fast, GPU-accelerated code for rapidly computing such solutions with accurate transport and an approximate collision operator. Such calculations are relevant to both problems in astrophysics, such as heat conduction and magnetic reconnection in accretion disks around black holes, and also to programmatic interests at LANL. The bolt code paper, demonstrating accuracy via a suite of test problems calculated on kodiak, is currently in preparation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

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↗
DOE OSTI2021

Dynamic Phase Alignment in Inertial Alfven Turbulence

In weakly-collisional plasma environments with sufficiently low electron beta, Alfvenic turbulence transforms into inertial Alfvenic turbulence at scales below the electron skin-depth. We argue that, in inertial Alfvenic turbulence, both energy and generalized kinetic helicity exhibit direct cascades. We demonstrate that the two cascades are compatible due to the existence of a strong scale-dependence of the phase alignment angle between velocity and magnetic field fluctuations, with the phase alignment angle scaling as proportional to the perpendicular wavenumber to the power of -1. The kinetic and magnetic energy spectra scale as the perpendicular wavenumber to the power of -5/3 and -11/3, respectively. As a result of the dual direct cascade, the generalized-helicity spectrum scales as the perpendicular wavenumber to the power of -5/3, implying progressive balancing of the turbulence as the cascade proceeds to smaller scales in the sub-electron-skin-depth range. Turbulent eddies exhibit a phase-space anisotropy, with the parallel wavenumber scaling as the perpendicular wavenumber to the power of 5/3, consistent with critically-balanced inertial Alfven fluctuations. Our results may be applicable to a variety of geophysical, space, and astrophysical environments, including the Earth's magnetosheath and ionosphere, solar corona, non-relativistic pair plasmas, as well as to strongly rotating non-ionized fluids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Closure theory for high-collisionality multi-ion plasmas

A general formalism is developed to construct and solve a system of linearized moment equations for parallel and perpendicular closures in high-collisionality plasmas. It is applicable for multiple ion species with arbitrary masses, temperatures, charges, and densities. The convergence of closure coefficients is evaluated by increasing the number of moments from 2 to 32 for scalar, vector, and rank-2 tensor moments. As an example, the complete set of closure coefficients for a deuterium-carbon plasma over the entire Hall parameter range is presented. Furthermore, the closure coefficients at various temperature ratios show that the one-temperature closure coefficients can differ significantly from the two-temperature coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Resolving the mystery of electron perpendicular temperature spike in the plasma sheath

A large family of plasmas has collisional mean-free-path much longer than the non-neutral sheath width, which scales with the plasma Debye length. The plasmas, particularly the electrons, assume strong temperature anisotropy in the sheath. The temperature in the sheath flow direction (Te∥) is lower and drops toward the wall as a result of the decompressional cooling by the accelerating sheath flow. The electron temperature in the transverse direction of the flow field (Te⊥) not only is higher but also spikes up in the sheath. This abnormal behavior of Te⊥ spike is found to be the result of a negative gradient of the parallel heat flux of transverse degrees of freedom (qes) in the sheath. The non-zero heat flux qes is induced by pitch-angle scattering of electrons via either their interaction with self-excited electromagnetic waves in a nearly collisionless plasma or Coulomb collision in a collisional plasma, or both in the intermediate regime of plasma collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On the collisional damping of plasma velocity space instabilities

For plasma velocity space instabilities driven by particle distributions significantly deviated from a Maxwellian, weak collisions can damp the instabilities by an amount that is significantly beyond the collisional rate itself. This is attributed to the dual role of collisions that tend to relax the plasma distribution toward a Maxwellian and to suppress the linearly perturbed distribution function. The former effect can dominate in cases where the unstable non-Maxwellian distribution is driven by collisionless transport on a timescale much shorter than that of collisions, and the growth rate of the ideal instability has a sensitive dependence on the distribution function. The whistler instability driven by electrostatically trapped electrons is used as an example to elucidate such a strong collisional damping effect of plasma velocity space instabilities, which is confirmed by first-principles kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Kinematics of femtosecond laser-generated plasma expansion: Determination of sub-micron density gradient and collisionality evolution of over-critical laser plasmas

An optical diagnostic based on resonant absorption of laser light in a plasma is introduced and is used for the determination of density scale lengths in the range of 10 nm to >1 μm at the critical surface of an overdense plasma. This diagnostic is also used to extract the plasma collisional frequency, allowing inference of the temporally evolving plasma composition on the tens of femtosecond timescale. This is found to be characterized by two eras: the early time and short scale length expansion (L < 0.1λ), where the interaction is highly collisional and target material dependent, followed by a period of material independent plasma expansion for longer scale lengths (L > 0.1λ); this is consistent with a hydrogen plasma decoupling from the bulk target material. Density gradients and plasma parameters on this scale are of importance to plasma mirror optical performance and comment is made on this theme.

Scott, G. G. (ORCID:000000032772372X)↗
DOE OSTI2025

Microtearing stability and turbulence in the pedestal: Linear gyrokinetics, reduced models, and nonlinear turbulent transport

Microtearing modes can play a crucial role in electron heat transport in tokamak plasmas, affecting both energy confinement and overall performance. This study investigates microtearing modes (MTM) stability and turbulence in a JET pedestal through gyrokinetic simulations using the Gene code, complemented by a reduced eigenvalue model. The focus is on how MTM properties depend on key plasma parameters, including collisionality and plasma beta β—the ratio of plasma pressure to magnetic pressure—the normalized toroidal wavenumber k y ρ s ⁠, where ρ s denotes the ion sound gyroradius (typically a few millimeters in edge plasmas) and isotope mass. Collisionality enhances MT growth rates, while increasing β leads to a shift from MTMs to kinetic-ballooning modes, typically for k y ρ s ⁠, where ρ s ≲ 0.2⁠. A purely collisionless branch of MTMs persists at low k y ρ s ⁠, where ρ s with distinctive properties including non-negligible particle flux and ion thermal transport. Isotope mass scans reveal modest reduction of MTM growth rates as ion mass decreases. Nonlinear simulations produce experimentally relevant transport levels. Numerical experiments turning off zonal flows and fields identify the critical role of zonal flows and zonal fields in regulating MTM turbulence. Their removal leads to a significant increase in electron heat flux. These findings provide new insight into MTM-driven transport and its impact on tokamak confinement and lay a foundation for reduced modeling and predictive capabilities.

Electrostatics↗
DOE OSTI2023

Observation of long-radial-range-correlation in turbulence in high-collisionality high-confinement plasmas on DIII-D

Abstract We report on the observation of spatially asymmetric turbulent structures with a long radial correlation length in the core of high-collisionality H -mode plasmas on DIII-D tokamak. These turbulent structures develop from shorter wavelength turbulence and have a radially elongated structure. The envelope of turbulence spans a broad radial range in the mid-radius region, leading to streamer-like transport events. The underlying turbulence is featured by intermittency, long-term memory effect, and the characteristic spectrum of self-organized criticality. The amplitude and the radial scale increase substantially when the shearing rate of the mean flow is reduced below the turbulent scattering rate. The enhanced long-radial-range-correlated (LRRC) transport events are accompanied by apparent degradation of normalized energy confinement time. The emergence of such LRRC transport events may serve as a candidate explanation for the degrading nature of H-mode core plasma confinement at high collisionality on DIII-D tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Global gyro-kinetic ion temperature gradient and trapped electron mode turbulence modelling in X -point geometry in negative and positive triangularity

Comparative modelling of ion temperature gradient (ITG)/trapped electron mode (TEM) turbulence in negative (NT) and positive (PT) triangularity plasma shapes was done using the nonlinear global full-f gyrokinetic particle code JOREK-GK in the realistic X-point tokamak geometry including the Scrape Off Layer (SOL) for TCV and DIII-D parameters. A comparison of JOREK-GK code with the gyrokinetic codes GS2 and GENE was done using NT/PT triangularity TCV L-modes parameters showing good agreement between codes in linear growth rates and clear beneficial effect of NT as compared to PT. Global non-linear modelling of the ITG/TEM saturated turbulence for realistic DIII-D NT pulses was done and compared with numerically constructed PT equilibrium with the same plasma profiles. Existence of longer correlation length of density fluctuations in PT compared to NT was demonstrated. Stronger and more sheared zonal flows are generated via Reynolds stress in NT compared to PT. These factors are stabilizing for TEM/ITG turbulence in NT and lead to smaller heat fluxes and heat conductivities in NT compared to PT configuration. Weak dependence of plasma confinement on collisionality and plasma rotation was found in modelling of DIII-D NT shots similar to the experiment. The confinement scaling with normalized ion gyro-radius ρ* was estimated both for NT and PT. Bohm-like scaling was obtained in both configurations, however with better confinement for NT compared to PT which could be favourable factor for reactor size machines.

ITG↗
DOE OSTI2021

Generalized impurity pinch in partially magnetized multi-ion plasma

In a two-ion-species plasma with disparate ion masses, heavy ions tend to concentrate in the low-temperature region of collisionally magnetized plasma and in the high-temperature region of collisionally unmagnetized plasma, respectively. Moreover, collisional magnetization can be determined as the ratio of the light ion gyrofrequency to the collision frequency of light and heavy ion species, and the behavior of this effect in the intermediate regime of partially magnetized plasma is predominantly dependent on this Hall parameter. Multi-ion cross-field transport has been described before in the collisionally magnetized plasma regime, and generalized pinch relations, which describe densities of ion species in equilibrium in that plasma, are found in the literature. Here, the role of collisional magnetization and Larmor magnetization in multi-ion collisional transport is clarified, and generalized pinch relations are extended to the partially magnetized regime in which the ion Hall parameter may be small, as long as electrons remain collisionally magnetized. Equilibrium ion density profiles have the same dependence on external forces and on each other regardless of collisional magnetization of ions. The expansion of the range of validity of multi-ion collisional transport models makes them applicable to a wider range of laboratory plasma conditions. In particular, ion density profiles evolve sufficiently fast for radial impurity transport to be observable around stagnation on MagLIF, leading to expulsion of heavy ion impurities from the hotspot as long as plasma becomes sufficiently collisionally magnetized during the implosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Evaluation of a collisional radiative model for electron temperature determination in hydrogen plasma

A collisional-radiative (CR) model that extracts the electron temperature, T e , of hydrogen plasmas from Balmer-line-ratio measurements is examined for the plasma electron density, n e , and T e ranges of 10 10 –10 15 cm –3 and 5–500 eV, respectively. The CR code, developed and implemented in Python, has a forward component that computes the densities of excited states up to n = 15 as functions of T e , n e , and the molecular-to-atomic neutral ratio r(H 2 /H). The backward component provides ne and r(H 2 /H) as functions of the Balmer ratios to predict the T e . The model assumes Maxwellian electrons. Furthermore, the density profiles of the electrons and of the molecular and atomic hydrogen neutrals are shown to be of great importance, as is the accuracy of the line-ratio measurement method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Feasibility Study: High-k Temperature (HiT) Fluctuation Diagnostic

The objective of the proposed research was to explore the feasibility of making measurements of high-k (electron-scale, ETG-relevant) electron temperature fluctuations in the core plasma of tokamaks and stellarators. We found that soft X-ray imaging combined with our novel pseudolocal tomographic algorithm enables local measurements of low-k (ion-scale, ITG/TEM and KBM/MTM relevant) electron temperature fluctuations in spherical tokamaks. The design parameters needed to consider deploying this new diagnostic are presented, and predictions for what the diagnostic would measure in L-mode and H-mode plasmas at NSTX-U, a high performance spherical tokamak in the USA, suggest that the new diagnostic can be used to study turbulence and transport over a wide range of plasma parameters, such as collisionality and plasma beta.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An implicit particle code with exact energy and charge conservation for electromagnetic studies of dense plasmas

A collisional particle code based on implicit energy- and charge-conserving methods is presented. A modified version of the particle-suppressed Jacobian-Free Newton-Krylov method that can enhance the solver efficiency is introduced. Mathematically, it is shown that this new approach can be viewed as a fixed-point iteration method for the particle positions. The model can exactly conserve global energy and local charge and can efficiently use time steps larger than the plasma period. In conclusion, the algorithm's ability to simulate dense plasmas accurately and efficiently is quantified by simulating the dynamic compression of a plasma slab via a magnetic piston in 1D planar geometry.

97 MATHEMATICS AND COMPUTING↗