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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Validation of the GFS model for gyrokinetic stability of NSTX pedestal data

This study presents a large database validation of the gyro fluid system (GFS) model for linear gyrokinetic stability for high-mode (H-mode) edge transport barrier conditions in the national spherical torus experiment (NSTX) tokamak. The database of linear stability calculations with the CGYRO gyrokinetic code was produced using plasma profile measurements from NSTX discharges to identify kinetic ballooning modes (KBM), trapped electron modes (TEM), and micro-tearing modes (MTM) that limit the pressure profile gradient in the H-mode barrier. A novel Bayesian optimization approach determines optimal resolution parameters for GFS specifically for spherical tokamak pedestal conditions. Our results demonstrate that GFS, with optimized resolution, can achieve accurate linear stability analysis in NSTX pedestal conditions for reduced resolution compared to CGYRO. GFS can accurately find the KBM, TEM, and MTM instability branches. Parametric analysis reveals that GFS accuracy in this extreme pedestal parameter range is degraded for low magnetic shear and near the separatrix conditions. These findings establish GFS as a fast linear eigenmode solver for spherical tokamak pedestal gyrokinetic stability and demonstrate a systematic methodology for determining the optimum resolution settings.

Yang, Minglei [Oak Ridge National Laboratory (ORNL↗
DOE OSTI2024

Synchrotron-driven Instabilities in Relativistic Plasmas of Arbitrary Opacity

Recent work has shown that synchrotron emission from relativistic plasmas leads the electron distribution to form an anisotropic ring in momentum space, which can be unstable to both kinetic and hydrodynamic instabilities. Fundamental to these works was the assumption that the plasma was optically thin, allowing all emitted radiation to escape. Here, we examine the behavior of these instabilities as the plasma becomes more optically thick. To do this, we extend a recently developed Fokker–Planck operator for synchrotron emission and absorption in mildly relativistic plasmas to ultrarelativistic plasmas. For a given set of plasma parameters, photons emitted by higher-energy electrons tend to be higher frequency, and thus more easily escape the plasma. As a result, the ratio of the photon emission rate (radiative drag) to absorption rate (radiative diffusion) for a given electron is extremely energy dependent. Given this behavior, we determine the critical parameters that control the opacity, and show how the plasma gradually transitions to become more isotropic and stable at higher opacity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Plasma decay of nanosecond pulsed laser-produced Ar and Ar–H 2 O sparks at atmospheric pressure

Time-resolved diagnostics were applied to investigate free-electron properties in nanosecond laser-produced discharges generated in atmospheric pressure Ar and in Ar–3%H 2 O. The discharges were generated using 23 ns, 1064 nm laser pulses. Broadband plasma imaging and laser Thomson scattering were combined with optical emission spectroscopy, with particular emphasis on the Stark broadening of the H α and H β lines. The plasma exhibited a bright emission that persists for up to 30–40 µs after breakdown. Plasma emission was then followed by a very weak glow emission that persisted for up to 19 ms after breakdown. Peak electron number density of ∼2 × 10 17 cm −3 and electron temperature of ∼7 eV were measured. An excellent agreement between both techniques was obtained regarding absolute electron number densities. The inferred free-electron temporal decay dynamics are consistent with processes dominated by hydrodynamic expansion and two- and three-body electron–ion recombination. These results provide benchmark data for modeling nanosecond laser discharges and demonstrate the reliability of combining Thomson scattering with Stark broadening in atmospheric laser sparks.

Thomson scattering↗
DOE OSTI2024

Final Report for Center for Tokamak Transient Simulations at USU

Providing plasma fluid codes like NIMROD with continuum drift kinetic (CDK) physics that is quantitatively valid and computationally feasible throughout the spatial domain is difficult. Work at Utah State University (USU), in collaboration with the Center for Tokamak Transient Simulations (CTTS), focused on applying CDK closures in disruption-related calculations. Three examples where kinetic physics is paramount are (1) the electron stress tensor closure in Ohms law for accurately describing neoclassical tearing mode (NTM) evolution, (2) runaway electron (RE) density (nRE) and current (jRE) moments in NIMROD’s extended MHD model for self-consistent evolution of RE populations during disruptions and, (3) energetic ion effects on a myriad of MHD instabilities. While NTM simulations and continuum and PIC approaches to energetic ions in NIMROD have been a major goals of USU’s closure work for several years, the development of self-consistent CDK RE capability in NIMROD was started and extended considerably during the CTTS effort. Some goals of CDK RE in NIMROD are to explore the effects of the 2D relativistic phase space in 4D simulations and compare with NIMROD’s fluid RE model. Four publications and two PhD theses came out of the USU CTTS effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Cross-beam energy transfer in conditions relevant to direct-drive implosions on OMEGA

In cross-beam energy transfer (CBET), the interference of two laser beams ponderomotively drives an ion-acoustic wave that coherently scatters light from one beam into the other. This redirection of laser beam energy can severely inhibit the performance of direct-drive inertial confinement fusion (ICF) implosions. To assess the role of nonlinear and kinetic processes in direct-drive-relevant CBET, the energy transfer between two laser beams in the plasma conditions of an ICF implosion at the OMEGA laser facility was modeled using particle-in-cell simulations. For typical laser beam intensities, the simulations are in excellent agreement with linear kinetic theory, indicating that nonlinear processes do not play a role in direct-drive implosions. At higher intensities, CBET can be modified by pump depletion, backward stimulated Raman scattering, or ion trapping, depending on the plasma density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Collisionless conduction in a high-beta plasma: a collision operator for whistler turbulence

The regulation of electron heat transport in high-β, weakly collisional, magnetized plasma is investigated. A temperature gradient oriented along a mean magnetic field can induce a kinetic heat-flux-driven whistler instability (HWI), which back-reacts on the transport by scattering electrons and impeding their flow. Previous analytical and numerical studies have shown that the heat flux for the saturated HWI scales as β$^{-1}_{e}$. These numerical studies, however, had limited scale separation and consequently large fluctuation amplitudes, which calls into question their relevance at astrophysical scales. To this end, we perform a series of particle-in-cell simulations of the HWI across a range of β$_e$ and temperature-gradient length scales under two different physical set-ups. The saturated heat flux in all of our simulations follows the expected β$^{-1}_{e}$ scaling, supporting the robustness of the result. We also use our simulation results to develop and implement several methods to construct an effective collision operator for whistler turbulence. The results point to an issue with the standard quasi-linear explanation of HWI saturation, which is analogous to the well-known 90° scattering problem in the cosmic-ray community. Despite this limitation, the methods developed here can serve as a blueprint for future work seeking to characterize the effective collisionality caused by kinetic instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Ionization disequilibrium in K- and L-shell ions

Time-gated Sc K-shell and Ge L-shell spectra are presented from a range of characterized thermodynamic states spanning ion densities of 1019–1020cm−3 and plasma temperatures around 2000 eV. For the higher densities studied and temperatures from 1000 to 3000 eV, the Sc and Ge x-ray emission spectra are consistent with steady-state calculations from the modern atomic kinetics model SCRAM. At the lower ion densities achieved through plasma expansion, however, the model calculations require a higher plasma temperature to reproduce the observed Ge spectrum. We attribute this to ionization disequilibrium of the Sc because the ionization time scales exceed the hydrodynamic timescale when the inferred temperatures diverge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Magnetic Field Amplification by a Nonlinear Electron Streaming Instability

Magnetic field amplification by relativistic streaming plasma instabilities is central to a wide variety of high-energy astrophysical environments as well as to laboratory scenarios associated with intense lasers and electron beams. We report on a new secondary nonlinear instability which arises for relativistic dilute electron beams after the saturation of the linear Weibel instability. This instability grows due to the transverse magnetic pressure associated with the beam current filaments, which cannot be quickly neutralized due to the inertia of background ions. We show that it can amplify the magnetic field strength and spatial scale by orders of magnitude, leading to large-scale plasma cavities with strong magnetic field and to very efficient conversion of the beam kinetic energy into magnetic energy. The instability growth rate, saturation level, and scale length are derived analytically and shown to be in good agreement with fully-kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Stability of perpendicular magnetohydrodynamic shocks in materials with ideal and nonideal equations of state

Magnetized target fusion approach to inertial confinement fusion involves the formation of strong shocks that travel along a magnetized plasma. Shocks, which play a dominant role in thermalizing the upstream kinetic energy generated in the implosion stage, are seldom free from perturbations, and they wrinkle in response to upstream or downstream disturbances. In Z-pinch experiments, significant plasma instability mitigation was observed with pre-embedded axial magnetic fields. To isolate effects, in this work we theoretically study the impact of perpendicular magnetic fields on the planar shock dynamics for different equations of state. For fast magnetosonic shocks in ideal gases, it was found that the magnetic field amplifies the intensity of the perturbations when $\gamma$ > 2 or it weakens them when $\gamma$ < 2. Weak shocks have been found to be stable regardless of the magnetic plasma intensity and gas compressibility; however, for sufficiently strong shocks the magnetic fields can promote a neutral stability/SAE at the shock if the adiabatic index is higher than 1 + $\sqrt{2}$. Finally, results have been validated with numerical simulations performed with the FLASH code.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Plasma instabilities dominate radioactive transients magnetic fields: the self-confinement of leptons in Type Ia and core-collapse supernovae, and kilonovae

The light curves of radioactive transients, such as supernovae and kilonovae, are powered by the decay of radioisotopes, which release high-energy leptons through $\beta ^+$ and $\beta ^-$ decays. These leptons deposit energy into the expanding ejecta. As the ejecta density decreases during expansion, the plasma becomes collisionless, with particle motion governed by electromagnetic forces. In such environments, strong or turbulent magnetic fields are thought to confine particles, though the origin of these fields and the confinement mechanism have remained unclear. Using fully kinetic particle-in-cell (PIC) simulations, we demonstrate that plasma instabilities can naturally confine high-energy leptons. These leptons generate magnetic fields through plasma streaming instabilities, even in the absence of pre-existing fields. The self-generated magnetic fields slow lepton diffusion, enabling confinement, and transferring energy to thermal electrons and ions. Our results naturally explain the positron trapping inferred from late-time observations of thermonuclear and core-collapse supernovae. Furthermore, they suggest potential implications for electron dynamics in the ejecta of kilonovae. We also estimate synchrotron radio luminosities from positrons for Type Ia supernovae and find that such emission could only be detectable with next-generation radio observatories from a Galactic or local-group supernova in an environment without any circumstellar material.

instabilities↗
DOE OSTI2021

Accurate numerical, integral methods for computing drift-kinetic Trubnikov-Rosenbluth potentials

A novel numerical method is employed to compute the integral form of the axi-symmetric Trubnikov-Rosenbluth potentials. Two methods for quadrature in pitch-angle are described and their convergence properties are studied. Careful attention is given to quadrature over a singular Green's function. Here it is shown that an infinite series representation of the Green's function can be used more efficiently than its closed form involving complete elliptic integrals. Then a collocation method in speed, with its associated quadrature scheme, is laid out and its convergence properties are studied. Using the proposed scheme, accurate low-order moments of the field collision operator are obtained using relatively few velocity space degrees of freedom. The scheme is showcased by solving for the equilibrium, axi-symmetric bootstrap current in tokamaks. A C 0 Gauss-Lobatto-Legendre finite element pitch-angle basis with vertex nodes at the trapped/passing boundary is shown, in the context of the integral methods used, to be much more efficient than the more common Legendre polynomial expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Conservative DG method for the micro-macro decomposition of the Vlasov–Poisson–Lenard–Bernstein model

The micro-macro (mM) decomposition approach is considered for the numerical solution of the Vlasov– Poisson–Lenard–Bernstein (VPLB) system, which is relevant for plasma physics applications. In the mM approach, the kinetic distribution function is decomposed as $f$ = $\mathscr{E}$[$ρ_f$] + g , where $\mathscr{E}$ is a local equilibrium distribution, depending on the macroscopic moments $ρ_f$ = $∫_{\mathbb{R}}$ $efdv$ = $\langle$$ef$$\rangle$$_{\mathbb{R}}$, where e = (1, $v$, $\frac{1}{2}$ $v^2$) T , and $g$, the microscopic distribution, is defined such that $\langle$$eg$$\rangle$$_{\mathbb{R}}$ = 0. We aim to design numerical methods for the mM decomposition of the VPLB system, which consists of coupled equations for $ρ_f$ and $g$. To this end, we use the discontinuous Galerkin (DG) method for phase-space discretization, and implicit-explicit (IMEX) time integration, where the phase-space advection terms are integrated explicitly and the collision operator is integrated implicitly. We give special consideration to ensure that the resulting mM method maintains the $\langle$$eg$$\rangle$$_{\mathbb{R}}$ = 0 constraint, which may be necessary for obtaining (i) satisfactory results in the collision dominated regime with coarse velocity resolution, and (ii) unambiguous conservation properties. The constraint-preserving property is achieved through a consistent discretization of the equations governing the micro and macro components. Here, we present numerical results that demonstrate the performance of the mM method. The mM method is also compared against a corresponding DG-IMEX method solving directly for $f$.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Simulation of a strong steady-state plasma shock in a warm dense matter regime

The structure of collisional plasma shocks has been subject to an extensive, multi-decadal investigation—in the hydrodynamic, hybrid kinetic ion/electron fluid, and fully kinetic ion/electron limits. Despite this thoroughness, all of these studies apply exclusively to classical, weakly coupled plasmas. Here, we report the first results for a planar hydrodynamic simulation of a strong, steady-state shock in a subspace of the warm dense matter (WDM) regime. Specifically, we consider a plasma of fully degenerate electrons with moderate-to-strongly coupled ions. Since the WDM ion and electron transport coefficients and equation of state differ markedly from their non-degenerate, weak-coupling equivalents, we find that the structure of a WDM plasma shock notably deviates from the ideal plasma picture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Runaway electron seed formation at reactor-relevant temperature

Systematic variation of the pre-disruption core electron temperature (T e ) from 1 to 12 keV using an internal transport barrier scenario reveals a dramatic increase in the production of ‘seed’ runaway electrons (REs), ultimately accessing near-complete conversion of the pre-disruption current into sub-MeV RE current. Injected Ar pellets are observed to ablate more intensely and promptly as T e rises. At high T e , the observed ablation exceeds predictions from published thermal ablation models. Simultaneously, the thermal quench (TQ) is observed to significantly shorten with increasing T e —a surprising result. While the reason for the shorter TQ is not yet understood, candidate mechanisms include: insufficiently accurate thermal ablation models, enhanced ablation driven by the seed RE population, or significant parallel heat transport along stochastic fields. Kinetic modeling that self-consistently treats the plasma cooling via radiation, the induced electric field, and the formation of the seed RE is performed. Including the combined effect of the inherent dependence of hot-tail RE seeding on T e together with the shortened TQ, modeling recovers the progression towards near-complete conversion of the pre-disruption current to RE current as T e rises. Measurement of the HXR spectrum during the early current quench (CQ) reveals a trend of decreasing energy with pre-disruption T e . At the very highest T e (≈ 12 keV), ≈ 100% conversion of the thermal current to runaway current is found. The energy of this peculiar RE beam is inferred to be sub-MeV as it emits vanishingly few MeV hard x-rays (HXRs). These measurements demonstrate novel TQ dynamics as Te is varied and illustrate the limitations of treating the RE seed formation problem without considering the inter-related dependencies of the pellet ablation, radiative energy loss, and resultant variations of the TQ duration. Finally, if the observed shortening of the TQ with increasing T e extends to fusion-grade plasmas, than their propensity to form large quantities of RE seeds at high T e may be far worse than previously thought. Positively, the high T e scenario in DIII-D produces REs so prodigiously that it can serve as a meaningful new platform for demonstrating RE avoidance techniques.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Exact and locally implicit source term solvers for multifluid-Maxwell systems

Recently, a family of models that couple multifluid systems to the full Maxwell equations have been used in laboratory, space, and astrophysical plasma modeling. These models are more complete descriptions of the plasma than reduced models like magnetohydrodynamic (MHD) since they are derived more closely from the full kinetic Vlasov-Maxwell system, without assumptions like quasi-neutrality, negligible electron mass, etc. Thus these models naturally retain non-ideal MHD effects like electron inertia, Hall term, pressure anisotropy/nongyrotropy, displacement current, among others. One obstacle to broader application of these model is that an explicit treatment of their source terms leads to the need to resolve rapid processes like plasma oscillation and electron cyclotron motion, even when these are not important. In this paper, we suggest two ways to address this issue. First, we derive the analytic solutions to the source update equations, which can be implemented as a practical, but less generic solver. We then develop a time-centered, locally implicit algorithm to update the source terms, allowing stepping over the fast kinetic time-scales. For a plasma with S species, the locally implicit algorithm involves inverting a local (3 S + 3) × (3 S + 3) matrix only, thus is very efficient. The performance can be further increased by using the direct update formulas to skip null calculations. In this paper, we present benchmarks illustrating the exact energy-conservation of the locally implicit solver, as well as its efficiency and robustness for both small-scale, idealized problems and largescale, complex systems. The locally implicit algorithm can be also easily extended to include other local sources, like collisions and ionization, which are difficult to solve analytically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Effects of equilibrium pressure on plasma response to RMPs in a spherical tokamak

This study presents a comprehensive analysis of the equilibrium pressure on the plasma response to resonant magnetic perturbations (RMPs) in the spherical tokamak (ST) MAST-U, employing both single-fluid and MHD-kinetic hybrid models (implemented via the MARS-F/K codes). As a key finding, the study identifies two different pressure-driven eigenmodes, exhibiting Sturmian property, that affect the Troyon no-wall limits for the onset of the n = 1 and n = 2 ( n is the toroidal mode number) ideal external kink instabilities as well as the corresponding plasma response to the applied RMP. With increasing equilibrium pressure, the plasma response to RMPs is significantly enhanced in the ST plasma, particularly in the high-pressure regime where kinetic effects strongly stabilize the external kink instability. The Troyon no-wall limit divides the plasma response into two regions: well below the limit, the response amplitudes and trends (versus pressure) are similar between the fluid and kinetic models; as the equilibrium pressure approaches the Troyon limit, the kinetic model predicts significant amplification of the RMP field, up to 30 times for cases considered. A relatively weak dependence of the optimal coil phasing on the equilibrium pressure is computed in this ST plasma, similar to the trend obtained for the conventional aspect ratio devices. These findings underscore the importance of incorporating kinetic effects in accurate prediction of the plasma response to RMPs in high-pressure ST tokamak plasmas and provide a theoretical basis for optimizing RMP-based control of the edge-localized modes in future ST devices.

RMP↗
DOE OSTI2021

Geometric electrostatic particle-in-cell algorithm on unstructured meshes

We present a geometric particle-in-cell (PIC) algorithm on unstructured meshes for studying electrostatic perturbations with frequency lower than electron gyrofrequency in magnetized plasmas. In this method, ions are treated as fully kinetic particles and electrons are described by the adiabatic response. The PIC method is derived from a discrete variational principle on unstructured meshes. To preserve the geometric structure of the system, the discrete variational principle requires that the electric field is interpolated using Whitney 1-forms, the charge is deposited using Whitney 0-forms and the electric field is computed by discrete exterior calculus. The algorithm has been applied to study the ion Bernstein wave (IBW) in two-dimensional magnetized plasmas. The simulated dispersion relations of the IBW in a rectangular region agree well with theoretical results. In a two-dimensional circular region with fixed boundary condition, the spectrum and eigenmode structures of the IBW are obtained from simulations. We compare the energy conservation property of the geometric PIC algorithm derived from the discrete variational principle with that of previous PIC methods on unstructured meshes. The comparison shows that the new PIC algorithm significantly improves the energy conservation property.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗