Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma models”

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

Determination of confinement regime boundaries via separatrix parameters on Alcator C-Mod based on a model for interchange-drift-Alfvén turbulence

The separatrix operational space (SepOS) model (Eich and Manz 2021 Nucl. Fusion 61 086017) is shown to predict the L–H transition, the L-mode density limit, and the ideal magnetohydrodynamic ballooning limit in terms of separatrix parameters for a wide range of Alcator C-Mod plasmas. The model is tested using Thomson scattering measurements across a wide range of operating conditions on C-Mod, spanning $\overline{n}_{\mathrm{e}}$ = 0.3–5.5 $\times 10^{20}\,$m −3 , $B_{\mathrm{t}} = 2.5$–8.0 T, and $B_{\mathrm{p}}$= 0.1–1.2 T. An empirical regression for the electron pressure gradient scale length, $\lambda_{{p}_{\mathrm{e}}}$, against a turbulence control parameter, $\alpha_{\mathrm{t}}$, and the poloidal fluid gyroradius, $\rho_{\mathrm{s,p}}$, is constructed for H-modes and found to require positive exponents for both regression parameters, indicating turbulence widening of near-scrape-off layer widths at high $\alpha_{\mathrm{t}}$ and an inverse scaling with $B_{\mathrm{p}}$, consistent with results on ASDEX Upgrade. The SepOS model is also tested in the unfavorable drift direction and found to apply well to all three boundaries, including the L–H transition as long as a correction to the Reynolds energy transfer term, $\alpha_\mathrm{RS} \lt 1$ is applied. I-modes typically exist in the unfavorable drift direction for values of $\alpha_{\mathrm{t}} \lesssim 0.3$. Finally, an experiment studying the transition between the Type-I ELMy and EDA H-mode is analyzed using the same framework. It is found that a recently identified boundary $\alpha_{\mathrm{t}} = 0.55$ at excludes most EDA H-modes but that the balance of wavenumbers responsible for the L-mode density limit, namely $k_\mathrm{EM} = k_\mathrm{RBM}$, may better describe the transition on C-Mod. The ensemble of boundaries validated and explored is then applied to project regime access and limit avoidance for the SPARC primary reference discharge parameters.

ELM suppression↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ∼1/Np with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly non-equilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ~1/$\sqrt{N_p}$ with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly nonequilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Lethality is Local, but Survival is Systemic: Temporal and Multi-Organ Responses to Chlorine Gas Exposure in a Murine Model

Chlorine gas (Cl2) is a highly toxic chemical associated with both localized lung injury and systemic health effects. While pulmonary damage has been well characterized, the systemic inflammatory and metabolic responses remain poorly understood. We aimed to define the temporal and multi-organ responses to Cl2 exposure in a murine model, with a focus on identifying spatiotemporal inflammation and its impact on survival and lethality. SKH1 mice were exposed for 10 min to varying concentrations of Cl2 (94.4–810 ppm, representative of non-lethal, LD10, and LD50 doses) and monitored for respiratory function, perfusion, and acidosis using organ-specific imaging. At multiple time points (40 min, 6 h, 24 h, and 7 d), we measured phosphoproteins, cytokines, chemokines, growth factors, and metabolic hormones in the lungs, heart, cortex, and plasma. Statistical modeling and logistic regression were used to identify biomarkers associated with lethality and survival. We found that lung injury was the primary cause of potential lethality, particularly via early phosphoprotein signaling disruptions. However, survival correlated with early systemic coordination of inflammatory and metabolic signals across organs. Perfusion and acidosis imaging were strongly associated with chemokine and hormone responses. Key survival-associated plasma biomarkers included decreased insulin, increased ghrelin, and decreased eotaxin. While potential lethality from Cl2 exposure is locally driven by pulmonary injury, survival depends on systemic, multi-organ responses that occur rapidly post-exposure. Within this model, our findings identify a potential therapeutic window to enhance survival and suggest candidate biomarkers that may be explored translationally for both triage and treatment of chlorine-related incidents.

chlorine gas↗
DOE OSTI2021

Novel 2D velocity estimation method for large transient events in plasmas

Dynamics of fast transient events are challenging to be analyzed with high time resolution. Such events can occur in fusion plasmas such as the filaments during edge-localized modes (ELMs). Here, we present a robust method—the spatial displacement estimation—for estimating the displacements of structures with fast dynamics from high spatial and time resolution imaging diagnostics [e.g., gas-puff imaging (GPI)] with sampling time temporal resolution. First, a background suppression method is shown, which suppresses the slowly time-evolving and spatially non-uniform background in the signal. In the second step, a two-dimensional polynomial trend subtraction method is presented to tackle the remaining polynomial order trend in the signal. After performing these pre-processing steps, the spatial displacement of the propagating structure is estimated from the two-dimensional spatial cross-correlation coefficient function calculated between consecutive frames. The method is tested for its robustness and accuracy by simulated Gaussian events and spatially displaced random noise. An example application of the method is presented on propagating ELM filaments measured by the GPI system on the National Spherical Torus Experiment spherical tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Deep Electric Field Predictions by Drift-Reduced Braginskii Theory with Plasma-Neutral Interactions Based on Experimental Images of Boundary Turbulence

We present 2-dimensional turbulent electric field calculations via physics-informed deep learning consistent with (i) drift-reduced Braginskii theory under the framework of an axisymmetric fusion plasma with purely toroidal field and (ii) experimental estimates of the fluctuating electron density and temperature on open field lines obtained from analysis of gas puff imaging of a discharge on the Alcator C-Mod tokamak. The inclusion of effects from the locally puffed atomic helium on particle and energy sources within the reduced plasma turbulence model are found to strengthen correlations between the electric field and electron pressure. The neutrals are also directly associated with broadening the distribution of turbulent field amplitudes and increasing ${\bf E \times B}$ shearing rates. As a result, this demonstrates a novel approach in plasma experiments by solving for nonlinear dynamics consistent with partial differential equations and data without encoding explicit boundary nor initial conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Weibel instability beyond bi-Maxwellian anisotropy

The shape of the anisotropic velocity distribution function, beyond the realm of strict Maxwellians can play a significant role in determining the evolution of the Weibel instability dictating the dynamics of self-generated magnetic fields. For non-Maxwellian distribution functions, we show that the direction of the maximum growth rate wave vector changes with shape. In this work, we investigate different laser-plasma interaction model distributions which show that their Weibel generated magnetic fields may require closer scrutiny beyond the second moment (temperature) anisotropy ratio characterization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Impurity transport in the pedestal of H-mode plasmas with resonant magnetic perturbations

Experiments on DIII-D show that resonant magnetic perturbations (RMPs) reduce the impurity confinement time in the H-mode pedestal below that of naturally ELMy plasmas. STRAHL modeling of discharges with aluminum laser blow off (LBO) injection show increased diffusion in the pedestal of H-mode discharges with RMPs. The increased diffusion from the RMPs provides improved impurity control compared to the naturally ELMy discharges and prevents the buildup of impurities near the edge of the plasma. Tungsten LBO injection into discharges with and without RMPs led to striking differences in the evolution of the discharge. The naturally ELMy discharge became ELM free and showed an increase in edge radiation, leading to the eventual radiative collapse of the plasma. Injecting a similar number ofWparticles into discharges with RMPs, either with or without ELMs, led to a measurable accumulation of core W but did not significantly impact the evolution of the discharge. Furthermore, these results indicate that RMPs are beneficial for controlling the buildup of impurities in the pedestal region of the plasma.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Plasma Dynamics in Low-Electron-Beta Environments

Recent in situ measurements by the MMS and Parker Solar Probe missions bring interest to small-scale plasma dynamics (waves, turbulence, magnetic reconnection) in regions where the electron thermal energy is smaller than the magnetic one. Examples of such regions are the Earth’s magnetosheath and the vicinity of the solar corona, and they are also encountered in other astrophysical systems. In this brief review, we consider simple physical models describing plasma dynamics in such low-electron-beta regimes, discuss their conservation laws and their limits of applicability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

An extended Vlasov–Fokker–Planck approach for kinetic simulations of laser plasmas

Vlasov–Fokker–Planck simulation codes occupy an important niche in modeling laser-produced plasmas, since they are well suited to studying the effect of collisions on electron kinetic phenomena, especially energy transport. One of the most important elements of energy transport is the absorption of laser light by the plasma; however, simulating this in detail requires resolving oscillations of the laser light, whose characteristic timescale is orders of magnitude shorter than the simulation time needed to study transport physics. For this reason, most Vlasov–Fokker–Planck codes used to study electron transport in laser plasmas rely on simplified models of the laser–plasma coupling. Their underlying assumptions nominally preclude their use for modeling laser light having short-scale structure in space or time, such as broadband lasers. In this work, we derive a more general computational framework suitable for arbitrarily structured laser fields. Furthermore, our approach is based on an extended set of Vlasov–Fokker–Planck equations that separately solve for the low- and high-frequency plasma response. We implement these extended Vlasov–Fokker–Planck equations in the spherical harmonic code K2 and demonstrate the performance of the method on several laser absorption test problems, with particular attention to the judicious selection of time steps, time integrators, and spherical harmonic truncation, according to the intensity and spectrum of the laser light under consideration. Comparison with the widely used Langdon absorption operator shows the Langdon operator performs remarkably well for predicting laser heating in the simple cases considered here, even in situations that would seem to violate its underlying assumptions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Nonlinear gyrokinetic predictions of SPARC burning plasma profiles enabled by surrogate modeling

Multi-channel, nonlinear predictions of core temperature and density profiles are performed for the SPARC tokamak accounting for both kinetic neoclassical and fully nonlinear gyro-kinetic turbulent fluxes. A series of flux-tube, nonlinear, electromagnetic simulations using the CGYRO code with six gyrokinetic species are coupled to a nonlinear optimizer using Gaussian process regression techniques. The simultaneous evolution of energy sources, including alpha heat, radiation, and energy exchange, coupled with these high fidelity models and techniques, leads to a converged solution in electron temperature, ion temperature and electron density channels with a minimal number of expensive gyrokinetic simulations without compromising accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Light Dark Matter from Entropy Dilution

We show that a thermal relic which decouples from the standard model (SM) plasma while relativistic can be a viable dark matter (DM) candidate, if the decoupling is followed by a period of entropy dilution that heats up the SM, but not the dark sector. Such diluted hot relics can be as light as few keV, while accounting for the entirety of the DM, and not conflicting with cosmological and astrophysical measurements. The requisite dilution can be achieved via decays of a heavy state that dominates the energy budget of the universe in the early matter dominated era. The heavy state decays into the SM particles, heats up the SM plasma, and dilutes the hidden sector. The interaction required to equilibrate the two sectors in the early universe places a bound on the maximum possible dilution as a function of the decoupling temperature. As an example of diluted hot relic DM we consider a light Dirac fermion with a heavy dark photon mediator. We present constraints on the model from terrestrial experiments (current and future), astrophysics, and cosmology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Mean-field dynamo as a quantum-like modulational instability

Presented here is a novel formulation of the mean-field dynamo as a modulational instability of magnetohydrodynamic (MHD) turbulence. This formulation, termed mean-field wave kinetics (MFWK), is based on the Weyl symbol calculus and allows describing the interaction between the mean fields (magnetic field and fluid velocity) and turbulence without requiring scale separation that is commonly assumed in the literature. The turbulence is described by the Wigner–Moyal equation for the spectrum of the two-point correlation matrix (Wigner matrix) of magnetic-field and velocity fluctuations and depicts the turbulence as an effective plasma of quantum-like particles that interact via the mean fields. Eddy–eddy interactions, which serve as ‘collisions’ in this effective plasma, are modelled within the standard minimal tau approximation to aid comparison with existing theories. Using MFWK, the non-local electromotive force is calculated for generic turbulence from first principles, modulo the limitations of MFWK. This result is then used to study, both analytically and numerically, the modulational modes of MHD turbulence, which appear as linear instabilities of the said effective quantum-like plasma of fluctuations. The standard α 2 -dynamo and other known results are reproduced as special cases. A new dynamo effect is predicted that is driven by correlations between the turbulent flow velocity and the turbulent current.

astrophysical plasmas↗
DOE OSTI2023

Physics and applications of dusty plasmas: The Perspectives 2023

Dusty plasmas are electrically quasi-neutral media that, along with electrons, ions, neutral gas, radiation, and electric and/or magnetic fields, also contain solid or liquid particles with sizes ranging from a few nanometers to a few micrometers. These media can be found in many natural environments as well as in various laboratory setups and industrial applications. As a separate branch of plasma physics, the field of dusty plasma physics was born in the beginning of 1990s at the intersection of the interests of the communities investigating astrophysical and technological plasmas. An additional boost to the development of the field was given by the discovery of plasma crystals leading to a series of microgravity experiments of which the purpose was to investigate generic phenomena in condensed matter physics using strongly coupled complex (dusty) plasmas as model systems. Finally, the field has gained an increasing amount of attention due to its inevitable connection to the development of novel applications ranging from the synthesis of functional nanoparticles to nuclear fusion and from particle sensing and diagnostics to nano-contamination control. The purpose of the present perspectives paper is to identify promising new developments and research directions for the field. As such, dusty plasmas are considered in their entire variety: from classical low-pressure noble-gas dusty discharges to atmospheric pressure plasmas with aerosols and from rarefied astrophysical plasmas to dense plasmas in nuclear fusion devices. Both fundamental and application aspects are covered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Modeling of injected current stream-induced 3D perturbations in local helicity injection plasmas

Solenoid-free tokamak startup techniques are essential for spherical tokamaks and offer a pathway to cost reduction and design simplification in fusion energy systems. Local helicity injection (LHI) is one such approach, employing compact edge current sources to drive open field line current that initiates and sustains tokamak plasmas. The recently commissioned Pegasus-III Experiment provides a platform for advancing this and other solenoid-free startup methods. This study investigates the effect of LHI on magnetic topology in Pegasus-III plasmas. A helical filament model represents the injected current, and the linear plasma response to its three-dimensional field is calculated with M3D-C1. Poincaré mapping reveals substantial flux surface degradation in all modeled cases. The onset of overlapping magnetic structures and large-scale surface deformation begins near Ψ N ≈ 0.37, indicating a broad region of perturbed topology extending toward the edge. In rotating plasmas, both single-fluid and two-fluid models exhibit partial screening of the n = 1 perturbation, with two-fluid calculations showing stronger suppression near the edge. In contrast, the absence of rotation leads to strong resonant field amplification in the single-fluid case, while the two-fluid case with zero electron rotation mitigates this amplification and preserves edge screening. Magnetic probe measurements indicate that modeling the stream with spatial spreading—representing distributed current and/or oscillatory motion—better reproduces measured magnetic power profiles than a rigid filament model. The results underscore the role of rotation and two-fluid physics in screening stream perturbations and point to plasma flow measurements and refined stream models as key steps toward improving predictive fidelity.

Schaefer, Carolyn E. [Univ. of Wisconsin, Madison,↗
DOE OSTI2022

Deep electric field predictions by drift-reduced Braginskii theory with plasma-neutral interactions based upon experimental images of boundary turbulence

We present 2-dimensional turbulent electric field calculations via physics-informed deep learning consistent with (i) drift-reduced Braginskii theory under the framework of an axisymmetric fusion plasma with purely toroidal field and (ii) experimental estimates of the fluctuating electron density and temperature obtained from analysis of gas puff imaging of a discharge on the Alcator C-Mod tokamak. The inclusion of effects from the locally puffed atomic helium on particle and energy sources within the reduced plasma turbulence model are found to strengthen correlations between the electric field and electron pressure. The neutrals are also directly associated with an observed broadening in the distribution of turbulent field amplitudes and increased E×B shearing rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Variations in lithium vapor cave performance predictions due to radial transport and recycling assumptions

The lithium vapor cave is a detached divertor design that uses a single private flux region baffle to contain a dense cloud of lithium vapor to dissipate heat flux. Plasma flows are created via fuel gas puffing in order to minimize lithium contamination of the main plasma. Significant modeling using the 2D edge code SOLPS-ITER has already been performed, predicting that sufficient target heat flux reductions (q$^{max}_{Target}$ MW m −2 ) with ( $n$ Li /$n$ e ) LCFS < 0.05 is possible in a case with 90 MW m −2 unmitigated heat flux in NSTX-U. Low heat flux and low upstream concentration was found with a variety of combinations of divertor geometries, target recycling coefficients, upstream plasma parameters, lithium evaporation locations and deuterium fueling locations, with variations in performance found for each design choice. However, the most universal uncertainty of SOLPS-ITER simulations has until now remained unaddressed systematically, namely the cross-field anomalous particle and heat diffusivities. This article aims to bound the uncertainty in lithium concentration prediction as a result of the assumed deuterium cross-field transport and recycling. For the simulations presented here, a factor of 2.1 increase in the upstream lithium density prediction across a factor of four decrease to the assumed deuterium radial particle diffusivity is found. This result is compared across different assumed deuterium recycling coefficients, known to reduce with lithium injection. Upstream lithium density is found to vary by a factor of 2.4 across a feasible range of recycling coefficients at similar $n$$^{OMP,sep}_{e}$.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Quantum Information for Fusion Energy Sciences (Final Technical Report)

The simulation of plasma dynamics is a critical area of Fusion Energy Sciences (FES) due to it’s usefulness in predicting, controlling, and confining plasmas in the context of potential fusion reactors. The simulation of plasmas is a computationally difficult problem in both classical and quantum physics, motivating investigation into the potential of quantum computers to simulate these systems. This project took several concrete steps towards this goal by developing tools for improving the control, characterization, and calibration of quantum gates on a superconducting quantum computer, developing error suppression and mitigation tools to reduce errors on the quantum computer, and utilizing these advancements to simulate reduced models of plasma dynamics on the quantum computer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗