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

A new approach to the evaluation and solution of the relativistic kinetic dispersion relation and verification with continuum kinetic simulation

Here, the present work describes a new approach to evaluation and root finding for the kinetic dispersion relation of Langmuir waves, which is central to the analytical understanding of collisionless damping in plasmas. The plasma dispersion function is solved to machine precision using direct integration in the complex plane in combination with an analytic evaluation of the residue to account for the deformation along the Landau contour. To efficiently attain machine precision, the contour is displaced in the complex plane prior to integration, and numerical subtleties related to the placement of the contour are discussed. The approach is generic in that it applies to arbitrary distribution functions, with the present manuscript focused on relativistic cases. Detailed verification of results via direct kinetic simulation in a variety of configuration space dimensions is also presented. Finally, the technique is applied to the challenging case of highly relativistic (i.e. extremely hot) plasmas. Here we show both qualitative agreement with prior work, as well as the disappearance of the Landau root which would have significant implication for real-life observation or experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Development of a synthetic phase contrast imaging diagnostic for turbulence studies at Wendelstein 7-X

Abstract We present a synthetic phase contrast imaging (PCI) diagnostic for studying turbulence at the Wendelstein 7-X (W7-X) stellarator. We first describe the implemented instrument response model, which captures diffraction effects, detector noise, and the long-wavelength cutoff due to the phase plate of the PCI system. To verify the instrument response model, we show that it is capable of reproducing the PCI signal generated by the sound wave speaker used for calibration at W7-X. Next, we discuss the calculation of synthetic PCI signals based on the global, nonlinear gyrokinetic codes GENE-3D and EUTERPE, including results from some of the first stellarator simulations of this type with kinetic electrons (KEs) in GENE-3D. While the simulations used in this work lack a neoclassical radial electric field, which is crucial for reproducing experimental PCI signals, they do indicate that the dominant rotation direction and velocities of the turbulent fluctuations can be inferred from the wave number-frequency spectra of the PCI signals, as expected. The synthetic PCI wave number spectra are further shown to be similar to those of the line-integrated fluctuating electron density, with distinct differences between adiabatic and KE simulations, explainable by previously published turbulence models. For example, the wave number spectra of all adiabatic electron simulations analyzed here follow a power law with an exponent close to −5 for sufficiently large wave numbers. This indicates that universal features of electron density turbulence at W7-X may be studied using the PCI system.

EUTERPE↗
DOE OSTI2020

X-ray heating and electron temperature of laboratory photoionized plasmas

We discuss experimental and modeling results for the x-ray heating and temperature of laboratory photoionized plasmas. A method is used to extract electron temperature based on the analysis of transmission spectroscopy data that is independent of atomic kinetics modeling. Boltzmann kinetic and radiation-hydrodynamics simulations emphasized the critical interplay between atomic physics and heating in photoionized plasmas. They also demonstrated the dramatic impact of photoexcitation on excited state populations, line emissivity, and radiation cooling. Temperature calculations performed with astrophysical codes significantly overestimated the measurement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Towards the first plasma-electron screening experiment

The enhancement of fusion reaction rates in a thermonuclear plasma by electron screening of the Coulomb barrier is an important plasma-nuclear effect that is present in stellar models but has not been experimentally observed. Experiments using inertial confinement fusion (ICF) implosions may provide a unique opportunity to observe this important plasma-nuclear effect. Herein, we show that experiments at the National Ignition Facility (NIF) have reached the relevant physical regime, with respect to the density and temperature conditions, but the estimated impacts of plasma screening on nuclear reaction rates are currently too small and need to be increased to lower the expected measurement uncertainty. Detailed radiation hydrodynamics simulations show that practical target changes, like adding readily available high-Z gases, and significantly slowing the inflight implosion velocity, while maintaining inflight kinetic energy, might be able to push these conditions to those where plasma screening effects may be measurable. We also perform synthetic data exercises to help understand where the anticipated experimental uncertainties will become important. But challenges remain, such as the detectability of the reaction products, non-thermal plasma effects, species separation, and impacts of spatial and temporal gradients. This work lays the foundation for future efforts to develop an important platform capable of the first plasma electron screening observation.

inertial confinement fusion↗
DOE OSTI2022

Equilibrium reconstruction of DIII-D plasmas using predictive modeling of the pressure profile

New workflows have been developed for predictive modeling of magnetohydrodynamic (MHD) equilibrium in tokamak plasmas. The goal of this work is to predict the MHD equilibrium in tokamak discharges without having measurements of the kinetic profiles. The workflows include a cold start tool, which constructs all the profiles and power flows needed by transport codes; a Grad–Shafranov equilibrium solver; and various codes for the sources and sinks. For validation purposes, a database of DIII-D tokamak discharges has been constructed that is comprised of scans in the plasma current, toroidal magnetic field, and triangularity. Initial efforts focused on developing a workflow utilizing an empirically derived pressure model tuned to DIII-D discharges with monotonic safety factor profiles. This workflow shows good agreement with experimental kinetic equilibrium calculations, but is limited in that it is a single fluid (equal ion and electron temperatures) model and lacks H-mode pedestal predictions. The best agreement with the H-mode database is obtained using a theory-based workflow utilizing pressure profile predictions from a coupled TGLF turbulent transport and EPED pedestal models together with external magnetics and Motional Stark Effect (MSE) data to construct the equilibrium. Here, we obtain an average root mean square error of 5.1% in the safety factor profile when comparing the predicted and experimental kinetic equilibrium. We also find good agreement with the plasma stored energy, internal inductance, and pressure profiles. Including MSE data in the theory-based workflow results in noticeably improved agreement with the q-profiles in high triangularity discharges in comparison with the results obtained with magnetic data only. The predictive equilibrium workflow is expected to have wide applications in experimental planning, between-shot analysis, and reactor studies.

Kinsey, J. E. (ORCID:0000000193347473)↗
DOE OSTI2024

Studies of High Energy Density Discharge and Laser-Driven Deflagrating Plasma Stagnations

The goal of this proposal is to develop an understanding of the thermalization and mixing of high energy density (HED) interpenetrating and stagnating plasma deflagrations on two different platforms. The first is opposing gas-fed Z-pinch plasma expansions in a facility at Stanford University. The second is the convergence of multiple deflagration fronts generated by inverted corona laser ablations through experiments carried out by the P.I. as part of a larger team at the University of Rochester Laboratory for Laser Energetics. The link between the two platforms is the role played by kinetic effects that impact mixing and the conversion of directed to thermal energy. Such kinetic effects are important because of the relatively high temperatures despite high plasma densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Absorption of relativistic multi-picosecond laser pulses in wire arrays

We study the interaction of intense multi-picosecond laser pulses with arrays of carbon wires attached to solid substrates. Here, we find that laser absorption in wire arrays resembles that in flat targets with very large uniform plasma density gradients. Performing two-dimensional particle-in-cell simulations, we optimize target parameters like wire thickness and -distance for energy absorption of a 2picosecond laser pulse with a large focal spot; this has implications for x-ray- and charged particle source development.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Fully Kinetic Simulations of Proton-beam-driven Instabilities from Parker Solar Probe Observations

The expanding solar wind plasma ubiquitously exhibits anisotropic nonthermal particle velocity distributions. Typically, proton velocity distribution functions (VDFs) show the presence of a core and a field-aligned beam. Novel observations made by the Parker Solar Probe (PSP) in the innermost heliosphere have revealed new complex features in the proton VDFs, namely anisotropic beams that sometimes experience perpendicular diffusion. In this study, we use a 2.5D fully kinetic simulation to investigate the stability of proton VDFs with anisotropic beams observed by PSP. Our setup consists of a core and an anisotropic beam population that drift with respect to each other. This configuration triggers a proton beam instability from which nearly parallel fast magnetosonic modes develop. Our results demonstrate that before this instability reaches saturation, the waves resonantly interact with the beam protons, causing perpendicular heating at the expense of the parallel temperature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A deep learning approach to fast analysis of collective Thomson scattering spectra

Fast analysis of collective Thomson scattering ion acoustic wave features using a deep convolutional neural network model is presented. The network was trained from spectra to predict the plasma parameters, including ion velocities, population fractions, and ion and electron temperatures. A fully kinetic particle-in-cell simulation was used to model a laboratory astrophysics experiment and simulate a diagnostic image of the ion acoustic wave feature. Network predictions were compared with Bayesian inference of the plasma model parameters for both the simulated and experimentally measured images. Both approaches were fairly accurate predicting the simulated image and the network predictions matched a good portion of the Bayesian results for the experimentally measured image. The Bayesian approach is more robust to noise and motivates future work to train deep learning models with realistic noise. The advantage of the deep learning model is making thousands of predictions in a few hundred milliseconds, compared to a few seconds to minutes per prediction for the optimization and Bayesian approaches presented here. The results demonstrate promising capabilities of deep learning models to analyze Thomson data orders of magnitude faster than conventional methods when using the neural network for standalone analysis. If more rigorous analysis is needed, neural network predictions can be used to quickly initialize other optimization methods and increase chances of success. This is especially useful when the dataset becomes very large or highly dimensional and manually refining initial conditions for the entire dataset are no longer tractable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Fusion burn-propagation simulations using the collisional and radiative particle-in-cell code TRIFORCE

The ability to accurately model burn propagation in inertial confinement fusion plasmas is crucial for advancing fusion energy research. This work presents enhancements to the triforce hybrid fluid-kinetic multiphysics code, focusing on its kinetic half, which employs the particle-in-cell (PIC) method with Monte Carlo collisions (MCC). We use a moment-preserving collision model that mitigates numerical noise, particularly in spherical geometries where particle weights vary significantly. Additionally, we refine the treatment of inverse bremsstrahlung to account for electron–ion collision frequency reductions in degenerate plasmas and incorporate a blackbody radiation source to enable realistic photon injection. These improvements enable the simulation of 1-dimensional (1D) spherical fusion burn propagation in deuterium–tritium plasmas. Benchmark comparisons with the hydra radiation-hydrodynamics code confirm that triforce accurately captures the dynamics of hot-spot expansion and burn propagation, demonstrating sensitivity to ignition thresholds consistent with theoretical models. Findings show the ignition cliff to be less steep in our work compared to radiation-hydrodynamic modeling. These results highlight the role of kinetic effects in fusion ignition physics and underscore the necessity of hybrid fluid-kinetic models for advancing predictive capabilities in high-energy-density plasma systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Kinetic Landau-fluid closures of non-Maxwellian distributions

New kinetic Landau-fluid closures, based on the cutoff Maxwellian distribution, are derived. A special static case is considered (the frequency ω=0). In the strongly collisional regime, our model reduces to Braginskii's heat flux model, and the transport is local. In the weak collisional regime, our model indicates that the heat flux is non-local and recovers the Hammett–Perkins model while the value of the cutoff velocity approaches to infinity. We compare the thermal transport coefficient χ of Maxwellian, cutoff Maxwellian and super-Gaussian distribution. The results show that the reduction of the high-speed tail particles leads to the corresponding reduction of the thermal transport coefficient χ across the entire range of collisionality, more reduction of the free streaming transport toward the weak collisional regime. In the collisionless limit, χ approaches to zero for the cutoff Maxwellian and the super-Gaussian distribution but remains finite for Maxwellian distribution. χ is complex if the cutoff Maxwellian distribution is asymmetric. The Im(χ) approaches to different convergent values in both collisionless and strongly collisional limit, respectively. It yields an additional streaming heat flux in comparison with the symmetric cutoff Maxwellian distribution. Furthermore, due to the asymmetric distribution, there is a background heat flux q0 though there is no perturbation. Furthermore, the derived Landau-fluid closures are general for fluid moment models, and applicable for the cutoff Maxwellian distribution in an open magnetic field line region, such as the scape-off-layer of Tokamak plasmas, in the thermal quench plasmas during a tokamak disruption, and the super-Gaussian electron distribution function due to inverse bremsstrahlung heating in laser-plasma studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Implicit highly-coupled single-ion Hall-MHD formulation for hybrid particle-in-cell codes

The rudiments of a particle-based single-fluid two-temperature magnetohydrodynamic (MHD) algorithm have been outlined in Thoma et al. (2013). The extension of this algorithm to include the effect of Hall physics is described in this paper. An implicit leapfrog version of the algorithm, which allows timesteps large compared to the resistive decay time and other relevant timescales, has recently been added to a hybrid particle-in-cell code. In standard MHD the Hall term in the generalized Ohm’s law can often be neglected when the Hall parameter is small. This term must, however, be retained in regimes where it is non-negligible. The retention of displacement current in Maxwell’s equations avoids the numerical difficulties associated with the whistler mode, which are encountered in standard explicit Hall-MHD codes, and allows the algorithm to be incorporated into hybrid particle-in-cell codes, for which particles may migrate from a kinetic to fluid to MHD description based upon local ambient plasma conditions. A highly-coupled implicit Hall-MHD formalism is presented, in which displacement current can either be retained or neglected. Even when displacement current is neglected, the highly-coupled implicit formalism avoids the restrictive timesteps for the whistler mode in explicit Hall-MHD codes. A comparison of numerical and analytic dispersion analysis demonstrates the feasibility of this approach and establishes relevant constraints to assure numerical stability. The implementation of the algorithm is described, and test simulation results in 1D and 2D in both linear and nonlinear regimes are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Data-driven, multi-moment fluid modeling of Landau damping

Deriving governing equations of complex physical systems based on first principles can be quite challenging when there are certain unknown terms and hidden physical mechanisms in the systems. In this work, we apply a deep learning architecture to learn fluid partial differential equations (PDEs) of a plasma system based on the data acquired from a fully kinetic model. Here, the learned multi-moment fluid PDEs are demonstrated to incorporate kinetic effect such as Landau damping. Based on the learned fluid closure, the data-driven, multi-moment fluid modeling can well reproduce all the physical quantities derived from the fully kinetic model. The calculated damping rate of Landau damping is consistent with both the fully kinetic simulation and the linear theory. The data-driven fluid modeling of PDEs for complex physical systems may be applied to improve the fluid closure and reduce the computational cost of multi-scale modeling of global systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Mitigation of mode-one asymmetry in laser-direct-drive inertial confinement fusion implosions

Nonuniformities present in the laser illumination and target in laser-driven inertial confinement fusion experiments lead to an asymmetric compression of the target, resulting in an inefficient conversion of shell kinetic energy to thermal energy of the hot-spot plasma. In this paper, the effects of asymmetric compression of cryogenic deuterium tritium laser-direct-drive implosions are examined using a suite of nuclear and x-ray diagnostics on the OMEGA laser. The neutron-averaged hot-spot velocity (u → hs ) and apparent ion temperature (T i ) asymmetry are determined from neutron time-of-flight measurements of the primary deuterium tritium fusion neutron energy spectrum, while the areal density (ρR) of the compressed fuel surrounding the hot spot is inferred from measurements of the scattered neutron energy spectrum. Here, the low-mode perturbations of the hot-spot shape are characterized from x-ray self-emission images recorded along three quasi-orthogonal lines of sight. Implosions with significant mode-1 laser drive asymmetries show large hot-spot velocities (>100 km/s) in a direction consistent with the hot-spot elongation observed in x-ray images, measured T i asymmetry, and ρR asymmetry. Laser drive corrections have been applied through shifting the initial target location in order to mitigate the observed asymmetry. With the asymmetry corrected, a more-symmetric hot spot is observed with reduced u → hs , T i asymmetry, ρR asymmetry, and a 30% increase in the fusion yield.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Weak magnetohydrodynamic turbulence theory revisited

Two recent papers, P. H. Yoon and G. Choe, Phys. Plasmas 28, 082306 (2021) and Yoon et al., Phys. Plasmas 29, 112303 (2022), utilized in the derivation of the kinetic equation for the intensity of turbulent fluctuations the assumption that the wave spectra are isotropic, that is, the ensemble-averaged magnetic field tensorial fluctuation intensity is given by the isotropic diagonal form, ⟨δBiδBj⟩k=⟨δB2⟩kδij. However, it is more appropriate to describe the incompressible magnetohydrodynamic turbulence involving shear Alfvénic waves by modeling the turbulence spectrum as being anisotropic. That is, the tensorial fluctuation intensity should be different in diagonal elements across and along the direction of the wave vector, ⟨δBiδBj⟩k=12 ⟨δB⊥2⟩k(δij−kikj/k2)+⟨δB∥2⟩k(kikj/k2). In the present paper, we thus reformulate the weak magnetohydrodynamic turbulence theory under the assumption of anisotropy and work out the form of nonlinear wave kinetic equation.

Physics↗
DOE OSTI2023

Parametric study of Alfvénic instabilities driven by runaway electrons during the current quench in DIII-D

Abstract To avoid or mitigate runaway electron (RE) beams in ITER, RE-driven instabilities are actively studied as a complimentary technique to massive material injection. In this work we report experimental dependencies of Alfvénic instabilities driven by REs during the current quench in DIII-D on plasma and material injection parameters. These instabilities, observed in the frequency range of 0.1–3 MHz, correlate with increased RE loss and thus may play a role in non-sustained RE beams. It was found that as the toroidal magnetic field ( B T ) decreases, the RE population becomes more energetic, the energy of instabilities increases, and no RE beam is observed when the maximum energy of REs exceeds 15 MeV (or when B T is below 1.8 T). Analysis of disruptions at plasma core temperature ( T e ) of 1 keV and 8 keV shows that the RE population is much less energetic (with the maximum energy of only about 3 MeV) when T e is high, and no instabilities are observed in this case. Besides disruptions above caused by Ar injection, cases with Ne and D 2 injections were also studied. Both Ne and D 2 injections cause no sustained RE beams, however, for different reasons. Measurements of the instability polarization indicate that it is of predominantly compressional nature at the edge, which is consistent with modeling suggesting excitation of compressional Alfvén eigenmodes. However, drive of global Alfvén eigenmodes is also possible at low frequencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Mitigation of mode-one asymmetry in laser-direct-drive inertial confinement fusion implosions

Nonuniformities present in the laser illumination and target in laser-driven inertial confinement fusion experiments lead to an asymmetric compression of the target, resulting in an inefficient conversion of shell kinetic energy to thermal energy of the hot-spot plasma. In this paper, the effects of asymmetric compression of cryogenic deuterium tritium laser-direct-drive implosions are examined using a suite of nuclear and x-ray diagnostics on the OMEGA laser. The neutron-averaged hot-spot velocity (~uhs) and apparent ion temperature (Ti) asymmetry are determined from neutron time-of-flight measurements of the primary deuterium tritium fusion neutron energy spectrum, while the areal density (rhoR) of the compressed fuel surrounding the hot spot is inferred from measurements of the scattered neutron energy spectrum. The low-mode perturbations of the hot-spot shape are characterized from x-ray self-emission images recorded along three quasi-orthogonal lines of sight. Implosions with significant mode-1 laser drive asymmetries show large hot-spot velocities (>100 km/s) in a direction consistent with the hot-spot elongation observed in x-ray images, measured Ti asymmetry, and rhoR asymmetry. Laser drive corrections have been applied through shifting the initial target location in order to mitigate the observed asymmetry. With the asymmetry corrected, a more-symmetric hot spot is observed with reduced ~uhs, Ti asymmetry, rhoR asymmetry, and a 30% increase in the fusion yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

1D Simulation of Avalanche to Streamer to Spark Transition of Plasma Discharge in Ammonia-Air Combustion

The backward problem of plasma assisted combustion emphasizes evaluating the effect of the evolving thermochemical state on the plasma discharge. This paper investigates the dependence of avalanche to streamer to spark formation dynamics and kinetics on the gas composition and temperature at different points in an ammonia-air premixed laminar flame using a self-consistent multigrid-based 1D plasma solver. Different values of ..alpha.., the coefficient for effective ionization events per unit length, have been reported for electron avalanches in air and stoichiometric NH3-air mixtures. The streamer inception has been shown to obey the Meek's criterion. An exponential reduction in streamer and spark formation time has been observed from plasma simulations at different points in the unburnt, pre-heat zone, reaction zone and the fully burnt regions of the premixed flame. While the enhancement of the reduced electric field with increasing temperature affects effective ionization, there exists a minimum breakdown field for streamer formation, which does not vary proportionally with the changing number density of the gas. The change in the mixture from reactants (NH3, O2, N2) to products of complete combustion of ammonia in air (N2, H2O) has also been shown to affect the streamer and spark formation. Finally, the major pathways during the streamer and spark phases which are responsible for producing important radicals used in combustion of NH3 are also discussed.

ADVANCED PROPULSION SYSTEMS,INORGANIC, ORGANIC, PH↗