Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA THEORY”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Degradation of performance in ICF implosions due to Rayleigh–Taylor instabilities: A Hamiltonian perspective

The Rayleigh–Taylor instability (RTI) is an ubiquitous phenomenon that occurs in inertial-confinement-fusion (ICF) implosions and is recognized as an important limiting factor of ICF performance. To analytically understand the RTI dynamics and its impact on ICF capsule implosions, we develop a first-principle variational theory that describes an imploding spherical shell undergoing RTI. The model is based on a thin-shell approximation and includes the dynamical coupling between the imploding spherical shell and an adiabatically compressed fluid within its interior. Using a quasilinear analysis, we study the degradation trends of key ICF performance metrics (e.g., stagnation pressure, residual kinetic energy, and areal density) as functions of initial RTI parameters (e.g., the initial amplitude and Legendre mode), as well as the 1D implosion characteristics (e.g., the convergence ratio). We compare analytical results from the theory against nonlinear results obtained by numerically integrating the governing equations of this reduced model. Our findings emphasize the need to incorporate polar flows in the calculation of residual kinetic energy and demonstrate that higher convergence ratios in ICF implosions lead to significantly greater degradation of key performance metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Quasilinear theory: the lost ponderomotive effects and why they matter

Quasilinear theory (QLT) has been used for modeling wave–plasma interactions for decades but remains largely heuristic. Plasma inhomogeneity, ponderomotive effects, microscopic fluctuations, and collisions are not easily accommodated from first principles in QLT, and typically are ignored entirely, due to the limitations of the standard Fourier–Laplace global-mode approach. This results in inconsistencies, for example, violation of the action conservation for nonresonant waves. However, these issues can be avoided, and the theory can be substantially generalized and corrected, if QLT is formulated using more suitable analytical tools, particularly, the Weyl symbol calculus. Here, an attempt is made to deliver an accessible review of this modern formulation, provide intuitive calculations for special cases, and elaborate on the connection with the ‘oscillation-center QLT’ originally proposed by Dewar (Phys Fluids 16:1102, 1973). A Fokker–Planck equation for a ‘dressed’ distribution is derived from the Klimontovich equation and captures quasilinear diffusion, ponderomotive forces, and interactions with background fields for a generic Hamiltonian, so many known formulations of QLT for specific plasma models become corollaries of a single unifying theory. Also, waves are allowed to be off-shell (not constrained by a dispersion relation), which allows them to accommodate microscopic fluctuations. This leads to a collision integral of the Balescu–Lenard type that has all the usual properties but is not restricted to any specific plasma model. For on-shell waves, a generalized version of the classic oscillation-center QLT is obtained. Finally, combined with the wave-kinetic equation, this formulation not only conserves particles, momentum, and energy, like the classic QLT but also reinstates the action conservation for nonresonant waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Alfvén eigenmode-driven zonal modes saturate and heat thermal ions by cross-scale interactions

In scenarios where a sustained energetic particle source strongly drives toroidal Alfvén eigenmodes (TAE), and phase-space transport is insufficient to saturate TAE, this novel theory of TAE-zonal mode (ZM)-turbulence—self-regulated by cross-scale interactions (including collisionless ZF damping) – merits consideration. Zonal modes are driven by Reynolds and Maxwell stresses, without the onset of modulational instability. TAE evolution in the presence of ZMs conserves energy and closes the system feedback loop. The saturated zonal shears can be sufficient to suppress ambient drift-ion temperature gradient (ITG) turbulence, achieving an enhanced core confinement regime. The saturated state is regulated by linear and turbulent zonal flow drag. This regulation leads to bursty TAE spectral oscillations, which overshoot while approaching saturation. Heating by both collisional and collisionless ZM damping deposits alpha particle energy into the thermal plasma, achieving effective alpha channeling. This theory offers a mechanism for EP-induced transport barrier formation, and predicts a novel thermal ion heating mechanism.

ITB↗
DOE OSTI2025

Measurement and Modeling of Electron Temperature in Laboratory Photoionized Plasmas Relevant to Astrophysics

The electron temperature of photoionized plasmas characterizes the thermalization of photoelectrons, impacts the charge-state distribution, emissivity, and opacity through atomic recombination processes, and is needed to perform detailed comparisons with theory predictions. We discuss temperature measurements in laboratory photoionized plasmas and a comparison with model calculations done with several theory approximations and codes. These include a radiation-hydrodynamics simulation and two nonequilibrium heating and ionization models that tracked the evolution of the internal energy of the electrons. Furthermore, for the same physics model and X-ray flux time history, calculations were performed assuming steady-state or time-dependent conditions. The time history of steady-state results correlates with that of the X-ray flux, while that of the time-dependent cases does not, and it is qualitatively and quantitatively different from the steady-state case. Steady-state results significantly overestimated temperature measurements, while time-dependent results produced better approximations, which suggests the importance of transient effects in the experiment and also the need for time-resolved measurements.

Laboratory astrophysics↗
DOE OSTI2025

Dynamics of heavy quarks in strongly coupled $\mathcal{N}$ = 4 SYM plasma

We calculate the probability distribution P(k) for a heavy quark with velocity v propagating through strongly coupled N = 4 SYM plasma in the ’t Hooft limit (N c → ∞, λ = g 2 N c → ∞) at a temperature T to acquire a momentum k due to interactions with the plasma. This distribution encodes the well-known drag coefficient η D and the transverse and longitudinal momentum diffusion coefficients κ T and κ L . The jet quenching parameter $\hat{q}$ can be extracted from P(k) for v = 1. Going beyond these known Gaussian characteristics of P(k), our calculation determines all of the higher order and mixed moments to leading order in 1/$\sqrt{λ}$ for the first time. These non-Gaussian features of P(k) include qualitatively novel correlations between longitudinal energy loss and transverse momentum broadening at nonzero v. We show that all higher moments scale characteristically with an effective temperature of the boosted plasma in the heavy quark rest frame, and we demonstrate that these non-Gaussian characteristics can be sizable in magnitude and even dominant in physically relevant situations. We use these results to derive a Kolmogorov equation for the evolution of the probability distribution for the total momentum of a heavy quark that propagates through strongly coupled plasma. This evolution equation accounts for all higher order correlations between transverse momentum broadening and longitudinal energy loss, which we have calculated from first principles. It reduces to a Fokker-Planck equation when truncated to only include the effects of η D , κ T and κ L . Remarkably, while heavy quarks do not reach kinetic equilibrium with the plasma if evolved with this Fokker-Planck equation, by showing that the Boltzmann distribution is a static solution of the all-order Kolmogorov equation that we have derived we demonstrate that heavy quarks do reach kinetic equilibrium if evolved with this equation. Our results thus provide a dynamically complete framework for understanding the thermalization of a heavy quark that may be initially far from equilibrium in the strongly coupled N = 4 SYM plasma — as well as new insight into heavy quark transport and equilibration in quark-gluon plasma.

Holography and Hydrodynamics↗
DOE OSTI2024

Automated Bayesian high-throughput estimation of plasma temperature and density from emission spectroscopy

Here, this paper introduces a novel approach for automated high-throughput estimation of plasma temperature and density using atomic emission spectroscopy, integrating Bayesian inference with sophisticated physical models. We provide an in-depth examination of Bayesian methods applied to the complexities of plasma diagnostics, supported by a robust framework of physical and measurement models. Our methodology is demonstrated using experimental observations in the field of magneto-inertial fusion, focusing on individual and sequential shot analyses of the Plasma Liner Experiment at LANL. The results demonstrate the effectiveness of our approach in enhancing the accuracy and reliability of plasma parameter estimation and in using the analysis to reveal the deep hidden structure in the data. This study not only offers a new perspective of plasma analysis but also paves the way for further research and applications in nuclear instrumentation and related domains.

Bayesian inference↗
DOE OSTI2025

Quantifying fission yields at the National Ignition Facility using depleted uranium foil experiments

There are programs for high-Z shell experiments at the National Ignition Facility (NIF). For shells made of actinide material, a quantitative fission diagnostic is needed in order to determine how much fission took place and whether the fission was sufficient to produce a non-negligible heat source in the burning capsule. Here, we present a viable coupled experimental and theoretical technique for making quantitative fission measurements possible. The proposed scheme involves using a well-characterized set of depleted uranium foils outside an NIF capsule to verify the conversion of xenon and krypton fission fragments collected at the Radiochemical Analysis of Gaseous Samples (RAGS) facility into total fission yield. We present the calculations needed for this conversion, including the decays in and out of fission fragment chains during the RAGS pump-down of the NIF chamber.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

A compact x-ray spectrometer for measurements of electron temperature distributions in inertial confinement fusion implosions at OMEGA

The Wedge Range Filter (WRF), commonly used for proton spectroscopy at the OMEGA Laser Facility and National Ignition Facility, is adapted to measure the x-ray continuum spectrum through transmission measurement using a continuous-gradient filter. Continuum x rays emitted from the hotspot of an implosion contain information about the plasma composition and electron temperature. The WRF data are leveraged to probe this distribution, specifically the electron temperature distribution. In this work, the data recorded with the WRF are forward modeled using a temperature distribution model folded with the WRF response function. An uncertainty analysis is conducted through a Bayesian regression algorithm using a Hamiltonian Monte Carlo sampler. This analysis enables the uncertainties in the instrument response to be folded into the uncertainty estimation of the electron temperature and absolute x-ray emission. Data analysis for a series of OMEGA implosions is presented and compared with radiation hydrodynamic simulations.

Lasers↗
DOE OSTI2025

Orbit-averaging and deposition accuracy for runaway electron beams in hybrid kinetic-MHD simulations of the runaway plateau

We develop a new procedure that combines the kinetic orbit runaway electron code (KORC) and the NIMROD extended-magnetohydrodynamic code to simulate runaway electrons (REs) in the post-disruption plateau. KORC integrates guiding-center orbits, with a barycentric-based binary search strategy providing initial guesses for the Newton–Raphson logical-to-physical coordinate inversion, ensuring reliable particle-to-mesh mapping in NIMROD, whose fields remain static for the present study. Samples are drawn in accord with experimental parallel current profiles of RE beams during the plateau phase. Deposition in NIMROD is verified through comparison with a Python-based finite-element code that ensures periodicity in the poloidal direction and continuity at the magnetic axis. Accurate representation of near-axis fields requires finer mesh resolution to prevent under- and overshoots in current density from orbit inaccuracies. Yet, at a fixed particle count, increasing mesh resolution amplifies statistical noise in the deposited fields. An orbit-averaging method accumulates partial current deposits over multiple kinetic steps and reduces the statistical noise with little added computational cost. By coupling kinetic routines from KORC directly into the NIMROD codebase, these developments lay essential groundwork for future self-consistent KORC–NIMROD coupling.

Algorithms and data structure↗
DOE OSTI2026

An ultra-fast method for generating synthetic down-scattered neutron data for inertial confinement fusion implosions

In inertial confinement fusion experiments at the National Ignition Facility, asymmetries are probed by a variety of neutron diagnostics, including neutron imaging systems, real-time neutron activation diagnostics (RTNADs), and neutron spectrometers. It is often useful to generate synthetic data based on these diagnostics to validate and tune models. However, current methods of doing so using Monte Carlo particle tracing are time-consuming. In this paper, an ultra-fast method is presented for generating synthetic neutron images, RTNAD data, and spectrometry data using line integrals and 3D convolutions. While it does not contain as much physics as particle tracing codes, it is thousands of times faster and produces nearly identical data. This enables analysis techniques that depend on generating large amounts of synthetic data, which will prove very useful for the study of asymmetries going forward.

Deuterium↗
DOE OSTI2026

Chaos as the cause of randomness in tearing mode onset times in DIII-D ITER baseline scenario plasmas

The first evidence of widespread chaotic plasma dynamics relevant to the onset of disruptive tearing modes is found in DIII-D tokamak discharges. This evidence is obtained by calculating positive maximal Lyapunov exponents from Mirnov signals using the Rosenstein algorithm, implying the chaotic growth and/or rotation of magnetohydrodynamic modes in DIII-D [Luxon, Nucl. Fusion 42, 614 (2002)] plasmas. Such chaotic dynamics offer an explanation for the random tearing mode onset times previously observed in ITER baseline scenario discharges. Further, the Lyapunov timescales associated with this chaos could inhibit reliable prediction of tearing mode onset beyond the angular momentum confinement timescale; this constraint should inform the design of plasma control systems on ITER.

Chaotic dynamics↗
DOE OSTI2026

Effects of screening and pressure ionization on the electron broadening of spectral lines in dense plasmas

Collisions between electrons and radiating atoms broaden spectral absorption and emission lines in dense plasmas. High densities also introduce screening and pressure ionization effects that distort the wave functions of both bound and free electrons. In order to study how dense plasma effects influence the electron broadening of spectral lines, this paper incorporates electron wave functions from an average-atom (AA) model to calculate the linewidth of the B III 2⁢𝑝−2⁢𝑠 transition at 𝑇 = 10 eV for mass densities ranging from 𝜌 = 10 −4 to 0.4⁢ g⁡/cc. The calculation method uses the impact approximation, allowing the linewidth to be written in terms of electron-collision cross sections and an interference term. Compared to an otherwise identical calculation that uses Coulomb free wave functions, the AA method is found to modify both the cross sections and the resulting linewidth at sufficiently high density by introducing screening and pressure-ionized bound states. Screening lowers the cross sections at low energies and near electron excitation thresholds, while pressure-ionized bound states introduce resonances into the continuum. Thus, as the density increases, the relative linewidth between the AA and Coulomb calculations follows a general decrease because of screening, with sharp increases at various intervals due to pressure ionization. Finally, the AA results are also compared with a common approach to introduce screening through the interaction potential and reduced models that use the Bethe formula for the inelastic electron-collision cross sections.

electronic excitation & ionization↗
DOE OSTI2025

Quantum computing for energy correlators

In recent years, energy correlators have emerged as powerful observables for probing the fragmentation dynamics of high-energy collisions. We introduce the first numerical strategy for calculating energy correlators using the Hamiltonian lattice approach, providing access to the intriguing nonperturbative dynamics of these observables. Furthermore, motivated by rapid advances in quantum computing hardware and algorithms, we propose a quantum algorithm for calculating energy correlators in quantum field theories. This algorithm includes ground state preparation, the application of source, sink, energy flux and real-time evolution operators, and the Hadamard test. We validate our approach by applying it to the SU(2) pure gauge theory in 2 + 1 dimensions on 3 × 3 and 5 × 5 honeycomb lattices with 𝑗 max = $\frac{1}{2}$ at various couplings, utilizing both classical methods and the quantum algorithm, the latter tested using the IBM emulator for specific configurations. The results are consistent with the expected behavior of the strong coupling regime and motivate a more comprehensive study to probe the confinement dynamics across the weak and strong coupling regimes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

LLNL FESP Theory Highlights: November 2024

This talk presented a refined understanding of plasma-facing thermionic emitter cathodes by examining the limitations of the conventional space charge-limited model, which overlooks ion trapping dynamics within the virtual cathode. While conventional theory distinguishes between temperature-limited and space charge-limited regimes based on electron current constraints, recent findings reveal that trapped ions can form a quasi-neutral "second plasma," expanding upstream and interacting with the primary plasma. These effects are amplified in non-planar geometries, such as cylindrical or spherical sheaths, where trapped ions create a positive feedback loop with the upstream plasma. To address these phenomena, a new "Aid-and-Compete" model was developed, offering insights into current enhancement and sheath mode transitions relevant for plasma devices, including emissive cathodes in electric propulsion, fusion applications, thin film deposition systems, and plasma diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Taylor limit studies for local helicity injection plasma startup

Solenoid-free tokamak startup techniques can simplify the design and reduce the cost of tokamak-based fusion energy systems. The newly commissioned Pegasus-III spherical tokamak provides a dedicated platform for developing a scalable solenoid-free startup approach. Local helicity injection (LHI) is one promising technique being developed on Pegasus-III, which uses compact edge current sources to drive open field line current I inj that is redistributed via helicity-conserving instabilities to initiate tokamak plasmas. These instabilities relax the system toward a minimum energy state described by Taylor relaxation theory. This relaxation process imposes a global current limit for LHI plasmas, referred to as the Taylor limit I TL ⁠, which can dictate the plasma current I p evolution. Experiments on Pegasus and Pegasus-III show that increasing this limit leads to more effective utilization of the available helicity input and higher I p plasmas. Initial Pegasus-III experiments have verified the expected I TL $\varpropto$ ( I inj B T ) 1/2 scaling holds for B T and I inj up to 0.3 T on axis and 12 kA ⁠, respectively. Additionally, data from different injector arrays show that the width of the injector aperture directly alters I TL and can be leveraged for future injector design optimization. Throughout these scaling studies, discharges were overdriven at I p = I TL by providing excess helicity input. Intermittent bursts of n = 1 magnetic activity consistent with large-scale reconnection events occurred more often as the system was increasingly overdriven leading to a flattening of the λ = μ 0 J ∥ / |B| profile, where J ∥ is the parallel current density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Overview of recent experimental results on the EAST Tokamak

Since the last IAEA-FEC in 2021, significant progress on the development of long pulse steady state scenario and its related key physics and technologies have been achieved, including the reproducible 403 s long-pulse steady-state H-mode plasma with pure radio frequency (RF) power heating. A thousand-second time scale (~1056 s) fully non-inductive plasma with high injected energy up to 1.73 GJ has also been achieved. The EAST operational regime of high β P has been significantly extended (H 98y2 > 1.3, β P ~ 4.0, β N ~ 2.4 and n e /n GW ~ 1.0) using RF and neutral beam injection (NBI). The full edge localized mode suppression using the n = 4 resonant magnetic perturbations has been achieved in ITER-like standard type-I ELMy H-mode plasmas with q 95 ≈ 3.1 on EAST, extrapolating favorably to the ITER baseline scenario. The sustained large ELM control and stable partial detachment have been achieved with Ne seeding. The underlying physics of plasma-beta effect for error field penetration, where toroidal effect dominates, is disclosed by comparing the results in cylindrical theory and MARS-Q simulation in EAST. Breakdown and plasma initiation at low toroidal electric fields (<0.3 V m -1 ) with EC pre-ionization is developed. A beneficial role on the lower hybrid wave injection to control the tungsten concentration in the NBI discharge is observed for the first time in EAST suggesting a potential way toward steady-state H-mode NBI operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Experimental and computational study of phase space dynamics in strongly coupled plasmas with steep density gradients

Understanding how plasmas thermalize when density gradients are steep remains a fundamental challenge in plasma physics, with direct implications for fusion experiments and astrophysical phenomena. Standard hydrodynamic models break down in these regimes, and kinetic theories make predictions that have never been directly tested. Here, we present the first detailed phase-space measurements of a strongly coupled plasma as it evolves from sharp density gradients to thermal equilibrium. Using laser-induced fluorescence imaging of an ultracold calcium plasma, we track the complete ion distribution function f(x,v,t). We discover that commonly used kinetic models (Bhatnagar–Gross–Krook and Lenard–Bernstein) overpredict thermalization rates, even while correctly capturing the initial counterstreaming plasma formation. Our measurements reveal that the initial ion acceleration response scales linearly with electron temperature, and that the simulations underpredict the initial ion response. In our geometry we demonstrate the formation of well-controlled counterpropagating plasma beams. This experimental platform enables precision tests of kinetic theories and opens new possibilities for studying plasma stopping power and flow-induced instabilities in strongly coupled systems.

Bergeson, Scott (ORCID:0000000231249226)↗
DOE OSTI2024

Interplay between the non-resonant streaming instability and self-generated pressure anisotropies

The non-thermal particles escaping from collision-less shocks into the surrounding medium can trigger a non-resonant streaming instability that converts parts of their drift kinetic energy into large amplitude magnetic field perturbations, and promote the confinement and acceleration of high energy cosmic rays. We present simulations of the instability using an hybrid-Particle-in-Cell approach including Monte Carlo collisions, and demonstrate that the development of the non-resonant mode is associated with important ion pressure anisotropies in the background plasma. Depending on the initial conditions, the anisotropies may act on the instability by lowering its growth and trigger secondary micro-instabilities. Introducing collisions with neutrals yield a strong reduction of the magnetic field amplification as predicted by linear fluid theory. In contrast, Coulomb collisions in fully ionized plasmas are found to mitigate the self-generated pressure anisotropies and promote the growth of the magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗