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

Role of electrostatic perturbation on kinetic resistive wall mode with application to spherical tokamak

Abstract A more complete non-perturbative magnetohydrodynamic (MHD)-kinetic hybrid formulation is developed by including the perturbed electrostatic potential δφ in the particle Lagrangian. The fluid-like counter-parts of the hybrid equations, in the Chew-Goldberger-Low high-frequency limit, are also derived and utilized to test the new toroidal implementation in the MARS-K code. Application of the updated non-perturbative hybrid model for a high- β spherical tokamak plasma in MAST finds that the perturbed electrostatic potential generally plays a minor role in the n = 1 ( n is the toroidal mode number) resistive wall mode instability. The effect of δφ is largely destabilizing, with the growth rate of the instability increased by several (up to 20) percent as compared to the case without including δφ . A similar relative change is also obtained for the kinetic-induced resonant field amplification effect at high- β in the MAST plasma considered. The updated capability of the MARS-K code allows quantitative exploration of drift kinetic effects on various MHD instabilities and the antenna-driven plasma response where the electrostatic perturbation, coupled to magnetic perturbations, may play important roles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Physics-guided machine learning approaches to predict the ideal stability properties of fusion plasmas

One of the biggest challenges to achieve the goal of producing fusion energy in tokamak devices is the necessity of avoiding disruptions of the plasma current due to instabilities. The disruption event characterization and forecasting (DECAF) framework has been developed in this purpose, integrating physics models of many causal events that can lead to a disruption. Two different machine learning approaches are proposed to improve the ideal magnetohydrodynamic (MHD) no-wall limit component of the kinetic stability model included in DECAF. First, a random forest regressor (RFR), was adopted to reproduce the DCON computed change in plasma potential energy without wall effects for a large database of equilibria from the national spherical torus experiment (NSTX). This tree-based method provides an analysis of the importance of each input feature, giving an insight into the underlying physics phenomena. Secondly, a fully-connected neural network has been trained on sets of calculations with the DCON code, to get an improved closed form equation of the no-wall β limit as a function of the relevant plasma parameters indicated by the RFR. The neural network has been guided by physics theory of ideal MHD in its extension outside the domain of the NSTX experimental data. The estimated value has been incorporated into the DECAF kinetic stability model and tested against a set of experimentally stable and unstable discharges. Moreover, the neural network results were used to simulate a real-time stability assessment using only quantities available in real-time. Finally, the portability of the model was investigated, showing encouraging results by testing the NSTX-trained algorithm on the mega ampere spherical tokamak (MAST).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Nonlinear susceptibilities for weakly turbulent magnetized plasma: Electrostatic approximation

The plasma weak turbulence theory is a perturbative nonlinear theory, which has been proven to be quite valid in a number of applications. However, the standard weak turbulence theory found in the literature is fully developed for highly idealized unmagnetized plasmas. As many plasmas found in nature and laboratory are immersed in a background static magnetic field, it is necessary to extend the existing discussions to include the effects of ambient magnetic field. Such a task is quite formidable, however, which has prevented fundamental and significant progresses in the subject matter. The central difficulty lies in the formulation of the complete nonlinear response functions for magnetized plasmas. The present paper derives the nonlinear susceptibilities for weakly turbulent magnetized plasmas up to the third order nonlinearity, but in doing so, a substantial reduction in mathematical complexity is achieved by the use of Bessel function addition theorem (or sum rule). The present paper also constructs the weak turbulence wave kinetic equation in a formal sense. For the sake of simplicity, however, the present paper assumes the electrostatic interaction among plasma particles. Fully electromagnetic generalization is a subject of a subsequent paper.

Physics↗
DOE OSTI2021

Bolt: A Fast Solver for Kinetic Theories Using a High-Resolution Constrained Transport Scheme [Slides]

Our understanding of collisionless and semi-collisional plasmas in the nonlinear regime is limited by the expense of computing solutions numerically. Bolt is a fast, GPU-accelerated code for rapidly computing such solutions with accurate transport and an approximate collision operator. Such calculations are relevant to both problems in astrophysics, such as heat conduction and magnetic reconnection in accretion disks around black holes, and also to programmatic interests at LANL. The bolt code paper, demonstrating accuracy via a suite of test problems calculated on kodiak, is currently in preparation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Transport coefficients of quasiparticle models within a new relaxation time approximation of the Boltzmann equation

We investigate the transport properties of a kinetic theory model that is tuned to describe the thermodynamic properties of QCD at zero chemical potential using a new formulation of the relaxation time approximation. In contrast to previous approaches, the latter is constructed to preserve the fundamental properties of the collision term of the Boltzmann equation for any energy dependence of the relaxation time. A novel choice of matching conditions is implemented to ensure that the background mean-field depends only on the temperature even when the system is out of equilibrium. We provide a consistent analysis of how the transport coefficients of relativistic Navier-Stokes theory vary with the energy dependence of the relaxation time. We also show that the entropy production of this theory is consistent with the second law of thermodynamics and verify that it is independent of the matching conditions employed. We used this fact to calculate the matching independent combination of transport coefficients.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Role of edge-localized neoclassical tearing modes in quiescent H-mode plasmas in the DIII-D tokamak

The edge-harmonic oscillations (EHOs) in standard quiescent H-mode (QH-mode) plasmas in DIII-D are consistent with edge-localized neoclassical tearing modes (NTMs) based on nonlinear two-fluid MHD simulations. Using kinetic equilibria constrained by edge profile measurements, the MHD simulations show that the n = 1 NTM and its harmonics can be destabilized at the pedestal top of QH-mode plasma by the edge bootstrap current. The simulations further show that the unstable NTMs can saturate either at small (<2% ψ N ) or large (>4% ψ N ) island width depending on the magnitude of the edge bootstrap current, where ψ N is the normalized radius in poloidal flux. The onset of the EHO also results in a prompt decrease in the pedestal width and height, consistent with simulation results for the onset of the NTM at the top of the QH-mode pedestal. This suggests that the avoidance of edge-localized modes (ELMs) in QH-mode can be attributed to the enhanced local transport induced by the NTM that is sufficient to prevent the expansion of the pedestal to an unstable width, analogous to the mechanism explored for ELM suppression by resonant magnetic perturbations. Nonlinear MHD simulations scanning the E × B frequency and the ratio of parallel and perpendicular thermal diffusivity (χ | /χ ⊥ ) at the pedestal top show that edge-localized NTMs are destabilized for conditions of high E × B frequency, high pedestal temperature, and low pedestal density, qualitatively consistent with experimental conditions required for observing the EHO.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Linear gyrokinetic stability of a high β non-inductive spherical tokamak

Spherical tokamaks (STs) have been shown to possess properties desirable for a fusion power plant such as achieving high plasma β and having increased vertical stability. To understand the confinement properties that might be expected in the conceptual design for a high β ST fusion reactor, a 1 GW ST plasma equilibrium was analyzed using local linear gyrokinetics to determine the type of micro-instabilities that arise. Kinetic ballooning modes and micro-tearing modes are found to be the dominant instabilities. The parametric dependence of these linear modes was determined and, from the insights gained, the equilibrium was tuned to find a regime marginally stable to all micro-instabilities at θ 0 = 0.0. This work identifies the most important micro-instabilities expected to generate turbulent transport in high β STs. The impact of such modes must be faithfully captured in first-principles-based reduced models of anomalous transport that are needed for predictive simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Planar Collisionless Shock Simulations with the Semi-implicit Particle-in-cell Model FLEKS

This study investigates the applicability of the semi-implicit particle-in-cell code FLexible Exascale Kinetic Simulator (FLEKS) to heliospheric shock simulations. We examine one- and two-dimensional local planar shock simulations, initialized using MHD states with upstream conditions representative of plasmas in the hypersonic, β ∼ 1 regime, for both quasi-perpendicular and quasi-parallel configurations. The refined algorithm in FLEKS proves robust, enabling accurate shock simulations with a grid resolution on the order of the electron inertial length d e . Our simulations successfully capture key shock features, including shock structures (foot, ramp, overshoot, and undershoot), upstream and downstream waves (fast magnetosonic, whistler, Alfvén ion-cyclotron, and mirror modes), and non-Maxwellian particle distributions. Crucially, we find that at least two spatial dimensions are critical for accurately reproducing downstream-wave physics in quasi-perpendicular shocks and capturing the complex dynamics of quasi-parallel shocks, including surface rippling, shocklets, short, large-amplitude magnetic structures, magnetic reconnection, and jets. Furthermore, our parameter studies demonstrate the impact of mass ratio and grid resolution on shock physics. This work provides valuable guidance for selecting appropriate physical and numerical parameters for shock simulations using a semi-implicit PIC method, paving the way for incorporating kinetic shock processes into large-scale collisionless plasma simulations with the MHD-AEPIC model.

plasma astrophysics↗
DOE OSTI2020

Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Low temperature hydrogen plasma permeation in palladium and its alloys for fuel recycling in fusion systems

Superpermeation of hydrogen isotopes through metal foils is a critical component for efficient fuel recycling in fusion power systems. In that context hydrogen permeation through foils of palladium and its alloys with silver and copper was studied at low temperature (60 - 200°C) under plasma exposure. These alloys differ significantly in both bulk and surface properties, and comparisons can provide mechanistic insights. Here, permeation was observed only during plasma operation, confirming the negligible contribution of molecular hydrogen to the observed flux. As-received foils required surface treatment to achieve top performance. For Pd and Pd75Ag25 an oxidation treatment increased permeation an order of magnitude, but proved unstable as this desirable surface was reduced under hydrogen plasma exposure. In contrast, an Ar plasma cleaning step provided both high and stable flux. As-received FCC phase Pd60Cu40 foils required annealing to transform it into the high permeability BCC phase that delivered top performance. All foils displayed similar temperature dependence with flux declining with temperature, suggesting that the primary rate-limiting step is absorption of superthermal hydrogen. Among these foils the hydrogen flux through PdCu was 3 - 5X greater than that of Pd or PdAg, which were similar. The superiority of BCC PdCu is attributed to its superior hydrogen desorption kinetics. Using PdCu 100% permeation of supplied hydrogen was achieved, and the flux saturated with increasing plasma power at values >10 –2 mol H•m –2 •s– 1 , and under these conditions permeation rates are equivalent with or without the membrane present. The fluxes achieved are the highest reported to date at these conditions, and the results highlight the important roles of both surface and bulk properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Cross-Beam Energy Transfer Saturation by Ion Heating

Cross-beam energy transfer (CBET) saturation by ion heating was measured in a gas-jet plasma characterized using Thomson scattering. Here, a wavelength-tunable ultraviolet (UV) probe laser beam was interacted with four intense UV pump beams to drive large-amplitude ion-acoustic waves. For the highest-intensity interactions, the power transferred to the probe laser dropped, demonstrating ion-acoustic wave saturation. Over this time, the ion temperature was measured to increase by a factor of 7 during the 500-ps interaction. Particle-in-cell simulations show ion trapping and a subsequent ion heating consistent with measurements. Linear kinetic CBET models were found to agree well with the observed energy transfer when the measured plasma conditions were used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On the minimum transport required to passively suppress runaway electrons in SPARC disruptions

In Izzo et al (2022 Nucl. Fusion 62 096029), state-of-the-art modeling of thermal and current quench (CQ) magnetohydrodynamics (MHD) coupled with a self-consistent evolution of runaway electron (RE) generation and transport showed that a non-axisymmetric (n = 1) in-vessel coil could passively prevent RE beam formation during disruptions in SPARC, a compact high-field tokamak projected to achieve a fusion gain Q > 2 in DT plasmas. However, such suppression requires finite transport of REs within magnetic islands and re-healed flux surfaces; conservatively assuming zero transport in these regions leads to an upper bound of RE current ~1MA compared to ~8.7MA of pre-disruption plasma current. Further investigation finds that core-localized electrons, within r/a < 0.4 and with kinetic energies ~0.2–15 MeV, contribute most to the RE plateau formation. Yet only a relatively small amount of transport, i.e. a diffusion coefficient ~18 m 2 s -1 , is needed in the core to fully mitigate these REs. Properly accounting for (a) the CQ electric field's effect on RE transport in islands and (b) the contribution of significant RE currents to disruption MHD may help achieve this.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Time-resolved CO 2 , CO, and N 2 vibrational population measurements in Ns pulse discharge plasmas

Abstract Time-resolved CO 2 and N 2 vibrational populations and translational-rotational temperature are measured in a CO 2 –N 2 plasma sustained by a ns pulse discharge burst in plane-to-plane geometry. Time-resolved, absolute number density of CO generated in the plasma is also inferred from the experimental data. CO 2 and CO vibrational populations are measured by mid-IR, tunable quantum cascade laser absorption spectroscopy, and N 2 vibrational populations are measured by the ns broadband vibrational CARS. Transient excitation of N 2 and CO 2 asymmetric stretch vibrational energy modes is detected during the discharge burst. The time-resolved rate of CO generation does not correlate with N 2 or CO 2 ( ν 3 ) vibrational temperatures, indicating that CO 2 dissociation via the vibrational excitation is insignificant at the present conditions. The rate of CO generation decreases gradually during the discharge burst. The estimated specific energy cost of the CO product is close to that of N atoms in pure nitrogen, measured previously at similar operating conditions. Comparison of the experimental data with the kinetic modeling analysis indicates that CO 2 dissociation in collisions with electronically excited N 2 molecules is the dominant channel of CO generation at the present conditions, although the inferred CO yield in these processes is significantly lower than 1. The effect of vibrational energy transfer between N 2 and CO 2 on the plasma chemical processes is insignificant. The kinetic model underpredicts a rapid reduction of the N 2 and CO 2 ( ν 3 ) vibrational temperatures during the later half of the discharge burst and in the afterglow. V–T relaxation of N 2 by N and O atoms generated in the ns pulse discharge plasma does not affect the vibrational relaxation rate in a significant way. However, rapid V–T relaxation of CO 2 by O atoms has a significant effect on the relaxation rate. The difference between the experimental data and the modeling predictions may be due to the unknown scaling of the CO 2 –O V–T rates with the vibrational quantum number.

Physics↗
DOE OSTI2024

A physics-informed deep learning description of Knudsen layer reactivity reduction

A physics-informed neural network (PINN) is used to evaluate the fast ion distribution in the hot spot of an inertial confinement fusion target. The use of tailored input and output layers to the neural network is shown to enable a PINN to learn the parametric solution to the Vlasov–Fokker–Planck equation in the absence of any synthetic or experimental data. As an explicit demonstration of the approach, the specific problem of Knudsen layer fusion yield reduction is treated. Here, the predictions from the Vlasov–Fokker–Planck PINN are used to provide a non-perturbative solution of the fast ion tail in the vicinity of the hot spot, thus allowing the spatial profile of the fusion reactivity to be evaluated for a range of collisionalities and hot spot conditions. Excellent agreement is found between the predictions of the Vlasov–Fokker–Planck PINN and the results from traditional numerical solvers with respect to both the energy and spatial distribution of fast ions and the fusion reactivity profile, demonstrating that the Vlasov–Fokker–Planck PINN provides an accurate and efficient means of determining the impact of Knudsen layer yield reduction across a broad range of plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The effect of rigid electron rotation on the Grad–Shafranov equilibria of a class of FRC devices

Rigid electron rotation of a Rotamak-FRC produces a pressure flux function that is more peaked than the Solov'ev flux function. This paper explores the implications of this peaked pressure flux function, including the isothermal case, which appear when the temperature profile is broader than the density profile, creating both benefits and challenges to a Rotamak-FRC based fusion reactor. In this regime, the density distribution becomes very peaked, enhancing the fusion power. The separatrix has a tendency to become oblate, which can be mitigated by flux conserving current loops. Plasma extends outside the separatrix, notably in the open field line region. Furthermore, this model does not apply to very kinetic FRCs or FRCs in which there are significant ion flows, but it may have some applicability to their outer layers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Isotope impact on Alfvén eigenmodes and fast ion transport in DIII-D

Abstract Measurements of beam driven Alfvén Eigenmode (AE) activity in matched deuterium (D) and hydrogen (H) DIII-D plasmas show a dramatic difference in unstable mode activity and fast ion transport for a given injected beam power. The dependence of the unstable AE spectrum in reversed magnetic shear plasmas on beam and thermal species is investigated in the current ramp by varying beam power in a sequence of discharges for fixed thermal and beam species at fixed density. In general, a spectrum of Reversed Shear Alfvén Eigenmodes (RSAEs) and Toroidal Alfvén Eigenmodes (TAEs) are driven unstable with sub-Alfvénic D beam injection while primarily only RSAEs are driven unstable for the H beam cases investigated. Further, for a given beam power, the driven AE amplitude is always reduced with H beams relative to D and for H thermal plasma relative to pure D or mixed D/H plasmas. Estimates of the fast ion stored energy combined with modeling using the hybrid kinetic-MHD code MEGA indicate that the dominant mechanism contributing to the difference between H and D beam drive is the faster classical slowing down of H beam ions relative to D and the resultant lower beam ion pressure. Calculations of the AE induced stored energy deficits using the reduced critical gradient model TGLFEP show quantitative agreement with the observed dependencies on injected power, isotope and minimum safety factor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Encapsulation of Monolayer 2D Materials Using Kinetic Energy-Controlled Pulsed Laser Deposition

The integration of monolayer (ML) two-dimensional (2D) materials into next-generation microelectronics, optoelectronics, and sensors is hindered by their sensitivity to environmental exposure. Deposition of additional layers for encapsulation or growth on ML 2D materials by versatile but energetic plasma techniques such as pulsed laser deposition (PLD) has not been considered at the monolayer level because of potential damage caused by hyperthermal species with kinetic energies (KEs) exceeding the threshold displacement energy (TDE) of the ML. Here, we describe a general strategy to understand and mitigate damage during PLD by reducing the incident KE of ablated species below the TDE of the 2D monolayer using background gas collisions. Ion flux diagnostics, combined with in situ Raman spectroscopy of monolayer graphene during PLD of amorphous boron nitride (a-BN) as a dielectric encapsulation layer, show that damage is primarily correlated with fast ions that penetrate the background gas in accordance with Beer’s Law and are often overlooked in ICCD imaging due to the dominance of the bright, delayed plasma luminescence. Significantly, if fast ions are eliminated and a ∼2 nm-thick a-BN layer is “soft landed”, the monolayer graphene is effectively protected from damage by high KE species in the boron nitride plasma plume. Deposited a-BN films display a characteristic dielectric constant of 3.6 at 100 kHz and tunable charge injection properties. Our results enable PLD as a viable option for encapsulation and thin film growth onto ML 2D materials, with implications for both fundamental research and device integration.

2D materials↗
DOE OSTI2021

Kinetic studies of excited singlet oxygen atom O( 1 D) reactions with ethanol

In this work, themultichannel reaction of excited singlet oxygen atom with ethanol, O( 1 D) + C 2 H 5 OH (1), was studied in a photolysis flow reactor coupled with mid-infrared Faraday rotation spectroscopy (FRS) and UV-IR direct absorption spectroscopy (DAS) at 297 K with reactor pressures of 60, 120, and 150 Torr (bath He). The excited singlet oxygen atom was generated through the photolysis of O 3 at 266 nm. The photon flux and O( 1 D) concentrations were determined by in situ actinometry based on O 3 depletion. Temporal profiles of OH and H 2 O were monitored via DAS signals at ca. 3568.62 and 3568.29 cm –1 , while temporal profiles of HO 2 were measured via FRS signals at ca. 1396.90 cm –1 . The branching ratios of the target reaction (1) were determined by fitting temporal profiles to simulations from an in-house reaction mechanism. Two major reaction channels were identified as CH 3 CHOH + OH and CH 3 O + CH 2 OH, and their branching ratios were determined as 0.46 ± 0.12 and 0.42 ± 0.11, respectively. A specific HO 2 + RO 2 reaction between HO 2 and O 2 CH 2 CH 2 OH (β-RO 2 ) at the low-temperature range is estimated in this work as HO 2 + O 2 CH 2 CH 2 OH → products with a rate constant of 7 × 10 –12 cm 3 molecule –1 s –1 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗