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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Ionization wave propagation in a He plasma jet in a controlled gas environment

Characterizing ionization wave propagation in low temperature plasma jets is critical to predicting production of reactive species and plasma–surface interactions for biomedical applications and surface functionalization. In this work, results from optical emission and laser induced fluorescence measurements of the ionization wave in a He plasma jet operating in a controlled gas environment are discussed and used for comparison with numerical modeling. Additionally, the ionization wave was observed using ICCD (Intensified Charge Coupled Device) imaging and characterized by time and spatially resolved electron density measurements using laser-collision-induced fluorescence. The plasma jet was initially characterized using pure He (nominally at 200 Torr), while varying pressure and voltage. When operating in pure He, the ionization wave broadly expands exiting the plasma tube. Increasing the operating pressure reduces the speed and isotropic expansion of the ionization wave. Furthermore, the jet operated with a humid He shroud was also studied. The humid He shroud results in the electron density increasing and having an annular profile due to the lower ionization potential of H2O compared to He and localized photoionization in the mixing region. Numerical modeling highlighted the importance of resonance radiation emitted by excited states of He, photoelectron emission from the quartz tube, and the kinetic behavior of the electrons produced by photoionization ahead of the ionization front.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

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

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

ITG↗
DOE OSTI2021

Statistical description of coalescing magnetic islands via magnetic reconnection

The physical picture of interacting magnetic islands provides a useful paradigm for certain plasma dynamics in a variety of physical environments, such as the solar corona, the heliosheath and the Earth's magnetosphere. In this work, we derive an island kinetic equation to describe the evolution of the island distribution function (in area and in flux of islands) subject to a collisional integral designed to account for the role of magnetic reconnection during island mergers. This equation is used to study the inverse transfer of magnetic energy through the coalescence of magnetic islands in two dimensions. We solve our island kinetic equation numerically for three different types of initial distribution: Dirac delta, Gaussian and power-law distributions. The time evolution of several key quantities is found to agree well with our analytical predictions: magnetic energy decays as $\tilde {t}^{-1}$ , the number of islands decreases as $\tilde {t}^{-1}$ and the averaged area of islands grows as $\tilde {t}$ , where $\tilde {t}$ is the time normalised to the characteristic reconnection time scale of islands. General properties of the distribution function and the magnetic energy spectrum are also studied. Finally, we discuss the underlying connection of our island-merger models to the (self-similar) decay of magnetohydrodynamic turbulence.

Physics↗
DOE OSTI2022

The influence of surface impurities on photoelectric currents driven by intense soft x rays

In an x-ray driven cavity experiment, an intense flux of soft x rays on the emitting surface produces significant emission of photoelectrons having several kiloelectronvolts of kinetic energy. At the same time, rapid heating of the emitting surface occurs, resulting in the release of adsorbed surface impurities and subsequent formation of an impurity plasma. This numerical study explores a simple model for the photoelectric currents and the impurity plasma. In this work, attention is given to the effect of varying the composition of the impurity plasma. The presence of protons or hydrogen molecular ions leads to a substantially enhanced cavity current, while heavier plasma ions are seen to have a limited effect on the cavity current due to their lower mobility. Additionally, it is demonstrated that an additional peak in the current waveform can appear due to the impurity plasma. A correlation between the impurity plasma composition and the timing of this peak is elucidated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Alfvén eigenmode classification based on ECE diagnostics at DIII-D using deep recurrent neural networks

Modern tokamaks have achieved significant fusion production, but further progress towards steady-state operation has been stymied by a host of kinetic and MHD instabilities. Control and identification of these instabilities is often complicated, warranting the application of data-driven methods to complement and improve physical understanding. In particular, Alfvén eigenmodes are a class of ubiquitous mixed kinetic and MHD instabilities that are important to identify and control because they can lead to loss of confinement and potential damage to the walls of a plasma device. In the present work, we use reservoir computing networks to classify Alfvén eigenmodes in a large labeled database of DIII-D discharges, covering a broad range of operational parameter space. However, despite the large parameter space, we show excellent classification and prediction performance, with an average hit rate of 91% and false alarm ratio of 7%, indicating promise for future implementation with additional diagnostic data and consolidation into a real-time control strategy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

R-Value Measurements Performed on Uranium Targets Irradiated with Fission Spectrum Neutrons FY 2021 for F2019 Project

The separation and characterization of two irradiated uranium targets, a depleted uranium (DU) and a highly enriched uranium (HEU) target, was conducted in April of 2021. The two targets were assembled at Los Alamos National Laboratory (LANL) and irradiated using the critical assembly at the National Criticality Experiments Research Center (NCERC). Splits of the dissolved targets were received by Pacific Northwest National Laboratory (PNNL) after which the PNNL and LANL teams chemically separated the solutions using independent separation schemes and analyzed the separated fractions for short lived actinides and fission products. Chemical separations at PNNL were traced with stable or radioactive tracers to allow for the determination of chemical yields, analyzing using either inductively coupled plasma optical emission spectroscopy (ICP-OES), inductively coupled plasma mass spectrometry (ICP-MS) or gamma emission analysis (GEA) depending on the nature of the tracer. Several other analytical techniques were used by PNNL including kinetic phosphorescence analysis (KPA) and thermal ionization mass spectrometry (TIMS) depending on the analyte’s need. Comparisons were made between current and historical PNNL and LANL, as well as literature values. Overall, there was agreement between the two laboratories for the bulk of analytes, with some notable exceptions such as 111 Ag, and 141,143,144 Ce. Included in these comparisons were the short-lived actinides 237 U, 239 Np, the fission products 89 Sr, 91 Y, 95/97 Zr, 99 Mo, 111 Ag, 115/115 mCd, 136/137 Cs, 140 Ba, 141/143/144 Ce, 147 Nd, 153 Sm, 156 Eu, and 161 Tb, providing both total atoms as well as the R-values. The data presented in this report represents the sixth NCERC irradiation of HEU and DU (FY13, FY14, FY15, FY17, FY18, FY21), and their subsequent separation and analysis.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Intermittent Electron-Only Reconnection at Lunar Mini-Magnetospheres

Lunar crustal magnetic anomalies (LCMA) are sub-ion-gyroradius structures that have been shown to stand off the solar wind (SW) plasma from the Moon’s surface, forming shock-like discontinuities and reflecting incident SW protons. In this Letter, the results of high-resolution, 2D fully kinetic simulations show a bursty electron-only magnetic reconnection in the SW-LCMA interaction region, characterized by the quasiperiodic formation and ejection of magnetic islands and strong parallel electron flows along the X-point separator lines. The islands are observed to modify the magnetic pressure pileup and Hall electric field above the LCMA, leading to sharp increases in reflected protons that drive electromagnetic fluctuations downstream and short distances upstream in the SW.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Lagrangian particle model for 3D simulation of pellets and SPI fragments in tokamaks

A 3D numerical model for the ablation of pellets and shattered pellet injection (SPI) fragments in tokamaks in the plasma disruption mitigation and fueling parameter space has been developed based on the Lagrangian particle code [R. Samulyak, X. Wang, H.-S. Chen, Lagrangian Particle Method for Compressible Fluid Dynamics, J. Comput. Phys., 362 (2018), 1-19]. The pellet code implements the low magnetic Reynolds number MHD equations, kinetic models for the electronic heating, a pellet surface ablation model, an equation of state that supports multiple ionization states, radiation, and a model for grad-B drift of the ablated material across the magnetic field. The Lagrangian particle algorithm is highly adaptive, capable of simulating a large number of fragments in 3D while eliminating numerical difficulties of dealing with the tokamak background plasma. The code has achieved good agreement with theory for spherically symmetric ablation flows. Axisymmetric simulations of neon and deuterium pellets in magnetic fields ranging from 1 to 6 Tesla have been compared with previous simulations using the FronTier code, and very good agreement has also been obtained. Furthermore, the main physics contribution of the paper is a detailed study of the influence of 3D effects, in particular grad-B drift, on pellet ablation rates and properties of ablation clouds. Smaller reductions of ablation rates in magnetic fields compared to axially symmetric simulations have been demonstrated because the ablated material is not confined to narrowing channels in the presence of grad-B drift. Contribution of various factors in the grad-B drift model has also been quantified.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Dynamics of dust particles confined in imposed potential structures in strongly magnetized, low-temperature plasmas

A new phenomenon called imposed ordered structures has been of particular interest to the dusty plasma community within the past several years. These structures are a new type of pattern formation within a dusty plasma in which the dust particles become fixed to a background confining potential whose spatial structure is determined by some conducting element present in the plasma. In previous works, this element has typically been a conducting wire mesh. One of the unanswered questions is whether the dust particles become trapped beneath the conducting surface of the mesh or in the gaps (holes) between the wires. This work makes use of a new electrode whose shape is designed to mimic that of wire mesh, but with much larger dimensions to facilitate in situ diagnostic measurements. Observations of the dust show that particles become confined to the regions beneath the holes of this new electrode. Measurements of the dust particle velocities allow for a determination of the kinetic energy within the dust cloud which show that the particles are in an energetic steady state. Further, comparisons of spatial profiles of the particle velocities also show that the dust particles typically reside in areas of increased plasma glow, possibly being trapped by plasma filaments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Plasma Dynamics and Nonthermal Particle Acceleration in 3D Nonrelativistic Magnetic Reconnection

Understanding plasma dynamics and nonthermal particle acceleration in 3D magnetic reconnection has been a long-standing challenge. In this paper, we explore these problems by performing large-scale fully kinetic simulations of multi-X-line plasmoid reconnection with various parameters in both the weak- and strong-guide-field regimes. In each regime, we have identified its unique 3D dynamics that lead to field-line chaos and efficient acceleration, and we have achieved nonthermal acceleration of both electrons and protons into power-law spectra. The spectral indices agree well with a simple Fermi acceleration theory that includes guide-field dependence. In the low-guide-field regime, the flux rope kink instability governs the 3D dynamics for efficient acceleration. The weak dependence of the spectra on the ion-to-electron mass ratio and β ($\ll$1) implies that the particles are sufficiently magnetized for Fermi acceleration in our simulations. While both electrons and protons are injected at reconnection exhausts, protons are primarily injected by perpendicular electric fields through Fermi reflections and electrons are injected by a combination of perpendicular and parallel electric fields. The magnetic power spectra agree with in situ magnetotail observations, and the spectral index may reflect a reconnection-driven size distribution of plasmoids instead of the Goldreich–Sridhar vortex cascade. As the guide field becomes stronger, the oblique flux ropes of large sizes capture the main 3D dynamics for efficient acceleration. Intriguingly, the oblique flux ropes can also experience flux rope kink instability, to drive extra 3D dynamics. This work has broad implications for 3D reconnection dynamics and particle acceleration in heliophysics and astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Hybrid-VPIC: An open-source kinetic/fluid hybrid particle-in-cell code

Hybrid-VPIC is an extension of the open-source high-performance particle-in-cell (PIC) code VPIC incorporating hybrid kinetic ion/fluid electron solvers. This paper describes the models that are available in the code and gives an overview of applications of the code to space and laboratory plasma physics problems. Particular choices in how the hybrid solvers were implemented are documented for reference by users. A few solutions for handling numerical complications particular to hybrid codes are also described. Finally, special emphasis is given to the computationally taxing problem of modeling mix in collisional high-energy-density regimes, for which more accurate electron fluid transport coefficients have been implemented for the first time in a hybrid PIC code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Precise stellarator quasi-symmetry can be achieved with electromagnetic coils

Magnetic fields with quasi-symmetry are known to provide good confinement of charged particles and plasmas, but the extent to which quasi-symmetry can be achieved in practice has remained an open question. Recent work [M. Landreman and E. Paul, Phys. Rev. Lett. 128, 035001, 2022] reports the discovery of toroidal magnetic fields that are quasi-symmetric to orders-of-magnitude higher precision than previously known fields. We show that these fields can be accurately produced using electromagnetic coils of only moderate engineering complexity, that is, coils that have low curvature and that are sufficiently separated from each other. Our results demonstrate that these new quasi-symmetric fields are relevant for applications requiring the confinement of energetic charged particles for long time scales, such as nuclear fusion. The coils’ length plays an important role for how well the quasi-symmetric fields can be approximated. For the longest coil set considered and a mean field strength of 1 T, the departure from quasi-symmetry is of the order of Earth’s magnetic field. Additionally, we find that magnetic surfaces extend far outside the plasma boundary used by Landreman and Paul, providing confinement far from the core. Simulations confirm that the magnetic fields generated by the new coils confine particles with high kinetic energy substantially longer than previously known coil configurations. In particular, when scaled to a reactor, the best found configuration loses only 0.04% of energetic particles born at midradius when following guiding center trajectories for 200 ms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Frequency multiplication with toroidal mode number of kink/fishbone modes on a static HL-2A-like tokamak

In the presence of energetic particles (EPs), the long-lived mode (LLM) frequency multiplication with n = 1, 2, 3, or higher is often observed on HL-2A, where n is the toroidal mode number. Hybrid kinetic-MHD model simulations of the energetic particle (EP) driven kink/fishbone modes on a static HL-2A-like tokamak using NIMROD code find that when the background plasma pressure is relatively high, and the EP pressure and the beam energy are relatively low, the mode frequency increases almost linearly with EP pressure, and the frequency is proportional to n (‘frequency multiplication’), even in the absence of any equilibrium plasma rotation. In addition, the frequency multiplication persists as the safety factor at the magnetic axis q 0 varies. In the absence of EPs, the growth rate of the 1/1 mode is the largest; however, as the EP pressure increases, the growth rate of 2/2 modes or 3/3 modes becomes dominant, suggesting that higher-n modes are more vulnerable to EPs. These results may shed light on the understanding of the toroidal mode number dependence of kink/fishbone modes in the advanced scenarios of tokamaks with weak or reversed central magnetic shear.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Exploring extreme magnetization phenomena in directly driven imploding cylindrical targets

Abstract This paper uses extended-magnetohydrodynamics (MHD) simulations to explore an extreme magnetized plasma regime realizable by cylindrical implosions on the OMEGA laser facility. This regime is characterized by highly compressed magnetic fields (greater than 10 kT across the fuel), which contain a significant proportion of the implosion energy and induce large electrical currents in the plasma. Parameters governing the different magnetization processes such as Ohmic dissipation and suppression of instabilities by magnetic tension are presented, allowing for optimization of experiments to study specific phenomena. For instance, a dopant added to the target gas-fill can enhance magnetic flux compression while enabling spectroscopic diagnosis of the imploding core. In particular, the use of Ar K-shell spectroscopy is investigated by performing detailed non-LTE atomic kinetics and radiative transfer calculations on the MHD data. Direct measurement of the core electron density and temperature would be possible, allowing for both the impact of magnetization on the final temperature and thermal pressure to be obtained. By assuming the magnetic field is frozen into the plasma motion, which is shown to be a good approximation for highly magnetized implosions, spectroscopic diagnosis could be used to estimate which magnetization processes are ruling the implosion dynamics; for example, a relation is given for inferring whether thermally driven or current-driven transport is dominating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Neutron generation dynamics inside a MA-class dense plasma focus Z-pinch

Dense plasma focii (DPFs) are appealing as energy efficient sources of short pulses of ions, neutrons, and x rays. The output of these sources is expected to scale with input current (I4), but has been shown to drop at the MA level [S. K. H. Auluck, “On the failure of neutron yield scaling in the dense plasma focus,” Phys. Plasmas 30, 080701 (2023)]. New results on the MegaJOuLe Neutron Imaging Radiography DPF showed neutron yield production in agreement with the input current scaling beyond the previously observed drop. This work provides insight into the pinch formation on a DPF and reports on the two different mechanisms leading to neutron generation inside a DPF using a combination of kinetic simulations and experimental data. A combination of particle-in-cell (PIC) and 1D shock theory results are used to describe the pinch formation and disassembly and the corresponding thermonuclear and beam-target mechanisms. The temporal evolution of the pinch column predicted by the PIC simulations shows qualitative agreement with the experimental data from plasma photon emission as well as temporal neutron pulse shapes. In MJ-class DPFs, both thermonuclear and beam-target mechanisms can occur over the course of the implosion and contribute to the total neutron production. Hence the neutron source size of a DPF will change throughout the implosion. Experimental neutron radiographs show the increase in source size as the pinch breaks apart, in agreement with simulation's prediction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Statistical data analysis of x-ray spectroscopy data enabled by neural network accelerated Bayesian inference

Bayesian inference applied to x-ray spectroscopy data analysis enables uncertainty quantification necessary to rigorously test theoretical models. However, when comparing to data, detailed atomic physics and radiation transfer calculations of x-ray emission from non-uniform plasma conditions are typically too slow to be performed in line with statistical sampling methods, such as Markov Chain Monte Carlo sampling. Furthermore, differences in transition energies and x-ray opacities often make direct comparisons between simulated and measured spectra unreliable. Here, we present a spectral decomposition method that allows for corrections to line positions and bound–bound opacities to best fit experimental data, with the goal of providing quantitative feedback to improve the underlying theoretical models and guide future experiments. In this work, we use a neural network (NN) surrogate model to replace spectral calculations of isobaric hot-spots created in Kr-doped implosions at the National Ignition Facility. The NN was trained on calculations of x-ray spectra using an isobaric hot-spot model post-processed with Cretin, a multi-species atomic kinetics and radiation code. The speedup provided by the NN model to generate x-ray emission spectra enables statistical analysis of parameterized models with sufficient detail to accurately represent the physical system and extract the plasma parameters of interest.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Electron Energization by Inertial Alfvén Waves in Density Depleted Flux Tubes at Jupiter

Juno satellite observations have illustrated substantial examples of broadband electron energization up to 10 5 –10 6 eV levels. In order to explain these observations, we use a hybrid gyrofluid kinetic-electron model in an untilted dipolar topology to illustrate energization to high levels in weak current conditions by inertial Alfvén waves close to the Jupiter ionosphere for ambient plasma densities and magnetic field perturbations inferred from Juno satellite observations. The key to the high energization is the extremely low densities evident in the observations which necessitates the acceleration of electrons to very high velocities in order to carry the field-aligned current.

58 GEOSCIENCES↗
DOE OSTI2021

Three dimensional low-mode areal-density non-uniformities in indirect-drive implosions at the National Ignition Facility

To achieve hotspot ignition, an inertial confinement fusion implosion must achieve high hotspot pressure that is inertially confined by a dense shell of DT fuel. This requires a symmetric implosion having high in-flight shell velocity and high areal density at stagnation. The size of the driver and scale of the capsule required can be minimized by maintaining a high efficiency of energy coupling from the imploding shell to the hotspot. Significant 3D low mode asymmetries, however, are commonly observed in indirect-drive implosions and reduce the coupling of shell kinetic energy to the hotspot. To better quantify the magnitudes and impacts of shell density asymmetries, we have developed new analysis techniques and analytic models [Hurricane et al., Phys. Plasmas 27(6), 062704 (2020)]. To build confidence in the underlying data, we have also developed an analytic neutron transport model to cross-compare two independent measurements of asymmetry, which shows excellent agreement across shots for mode-1 (ℓ = 1). This work also demonstrates that asymmetry can introduce potential sampling bias into down-scattered ratio measurements causing the solid-angle-average and uncertainty-weighted-average down-scattered ratios to differ significantly. Diagnosing asymmetries beyond mode-1 (ℓ > 1) presents significant challenges. Using new diagnostic instruments and analysis techniques, however, evidence of significant Legendre mode P2 (ℓ = 2, m = 0) and additional 3D asymmetries (ℓ > 1, m ≠ 0) are beginning to emerge from the high precision activation diagnostic data (real-time nuclear activation detectors) and down-scattered neutron imaging data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗