Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Magnetized plasma”

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

Isotope effects and Alfvén eigenmode stability in JET H, D, T, DT, and He plasmas

While much about Alfvén eigenmode (AE) stability has been explored in previous and current tokamaks, open questions remain for future burning plasma experiments, especially regarding exact stability threshold conditions and related isotope effects; the latter, of course, requiring good knowledge of the plasma ion composition. In the JET tokamak, eight in-vessel antennas actively excite stable AEs, from which their frequencies, toroidal mode numbers, and net damping rates are assessed. The effective ion mass can also be inferred using measurements of the plasma density and magnetic geometry. Thousands of AE stability measurements have been collected by the Alfvén Eigenmode Active Diagnostic in hundreds of JET plasmas during the recent Hydrogen, Deuterium, Tritium, DT, and Helium-4 campaigns. In this novel AE stability database, spanning all four main ion species, damping is observed to decrease with increasing Hydrogenic mass, but increase for Helium, a trend consistent with radiative damping as the dominant damping mechanism. These data are important for confident predictions of AE stability in both non-nuclear (H/He) and nuclear (D/T) operations in future devices. In particular, if radiative damping plays a significant role in overall stability, some AEs could be more easily destabilized in D/T plasmas than their H/He reference pulses, even before considering fast ion and alpha particle drive. Active MHD spectroscopy is also employed on select HD, HT, and DT plasmas to infer the effective ion mass, thereby closing the loop on isotope analysis and demonstrating a complementary method to typical diagnosis of the isotope ratio.

Alfvén eigenmodes↗
DOE OSTI2025

Nyström type exponential integrators for strongly magnetized charged particle dynamics

Solving for charged particle motion in electromagnetic fields (i.e. the particle pushing problem) is a computationally intensive component of particle-in-cell (PIC) methods for plasma physics simulations. This task is especially challenging when the plasma is strongly magnetized due numerical stiffness arising from the wide range of time scales between highly oscillatory gyromotion and long term macroscopic behavior. A promising approach to solve these problems is by a class of methods known as exponential integrators that can solve linear problems exactly and are A-stable. This work extends the standard exponential integration framework to derive Nyström-type exponential integrators that integrates the Newtonian equations of motion as a second-order differential equation directly. In particular, we derive second-order and third-order Nyström-type exponential integrators for strongly magnetized particle pushing problems. Numerical experiments show that the Nyström-type exponential integrators exhibit significant improvement in computation speed over the standard exponential integrators.

general physics↗
DOE OSTI2021

Modeling Snowflake Divertors in MAST-U Tokamak

We report that in a snowflake (SF) divertor, two magnetic field nulls are placed close to each other, creating four strike points (SPs) compared to two in a standard X-point divertor. In preparation for MAST-U experiments, magnetic configurations with the standard and SF divertors with various locations and separation distances of the nulls were modeled using the two-dimensional multi-fluid code UEDGE with a full plasma transport model featuring charge-state-resolved sputtered carbon impurities. The complex interplay of the plasma transport and magnetic configurations was comprehensively studied using a simple model for the theoretically predicted fast plasma mixing driven by the 'churning' mode instability in the two-null SF region. The modeling results show that (1) all SF-plus configurations and SF-minus configuration with closely located nulls produce the same plasma parameters and heat fluxes at the same SPs; (2) SF divertors approach the outer and inner SP detachment conditions at lower upstream density w.r.t. the standard divertor; (3) heat flux profiles at primary SPs are substantially broadened and peak values are reduced in SF configurations w.r.t. SN divertors; this broadening becomes more pronounced with the fast plasma mixing increase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Instability of an electron-plasma shear layer in an externally imposed strain flow

The E x B shear instability of a two-dimensional (2D) filament (i.e., a thin, rectangular strip perpendicular to the magnetic field) of magnetized pure electron plasma is investigated experimentally in the presence of an externally imposed strain flow. Data are acquired using a specialized Penning–Malmberg trap in which strain flows can be applied in 2D by biasing segmented electrodes surrounding the plasma. The E x B drift dynamics are well-described by the Drift-Poisson equations, which are isomorphic to the 2D Euler equations describing ideal fluids. Thus, the experimental results correspond to the Rayleigh instability of a shear layer in a 2D ideal fluid, where the electron density is analogous to the fluid vorticity. Shear layers are prepared by stretching initially axisymmetric electron vortices using a strong, applied strain flow. The data at early times are in quantitative agreement with a linear model which extends Rayleigh’s work to account for the influence of an external strain flow. In the presence of weak strain, the system approximately maintains a phase relationship that corresponds to an instantaneous Rayleigh eigenmode. The instability develops into the nonlinear regime later in time and at smaller spatial scales as the strain rate is increased. A secondary vortex pairing instability is observed, but it is suppressed when the strain-to-vorticity ratio exceeds roughly 0.025. In this way, vorticity transport perpendicular to the filament is diminished due to the applied strain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

On the study of hydrodynamic instabilities in the presence of background magnetic fields in high-energy-density plasmas

Blast-wave-driven hydrodynamic instabilities are studied in the presence of a background B-field through experiments and simulations in the high-energy-density (HED) physics regime. In experiments conducted at the Laboratoire pour l’utilisation des lasers intenses (LULI), a laser-driven shock-tube platform was used to generate a hydrodynamically unstable interface with a prescribed sinusoidal surface perturbation, and short-pulse x-ray radiography was used to characterize the instability growth with and without a 10-T B-field. The LULI experiments were modeled in FLASH using resistive and ideal magnetohydrodynamics (MHD), and comparing the experiments and simulations suggests that the Spitzer model implemented in FLASH is necessary and sufficient for modeling these planar systems. These results suggest insufficient amplification of the seed B-field, due to resistive diffusion, to alter the hydrodynamic behavior. Although the ideal-MHD simulations did not represent the experiments accurately, they suggest that similar HED systems with dynamic plasma-β (=2μ0ρv2/B2) values of less than ∼100 can reduce the growth of blast-wave-driven Rayleigh–Taylor instabilities. These findings validate the resistive-MHD FLASH modeling that is being used to design future experiments for studying B-field effects in HED plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Machine learning methods for probabilistic locked-mode predictors in tokamak plasmas

A rotating tokamak plasma can interact resonantly with the external helical magnetic perturbations, also known as error fields. This can lead to locking and then to disruptions. We leverage machine learning (ML) methods to predict the locking events. We use a coupled third-order nonlinear ordinary differential equation model to represent the interaction of the magnetic perturbation and the plasma rotation with the error field. This model is sufficient to describe qualitatively the locking and unlocking bifurcations. Here, we explore using ML algorithms with the simulation data and experimental data, focusing on the methods that can be used with sparse datasets. These methods lead to the possibility of the avoidance of locking in real-time operations. We describe the operational space in terms of two control parameters: the magnitude of the error field and the rotation frequency associated with the momentum source that maintains the plasma rotation. The outcomes are quan- tified by order parameters that completely characterize the state, whether locked or unlocked. We use unsupervised ML methods to classify locked/unlocked states and note the usefulness of a certain normalization of the order parameters. Three supervised ML classifiers are used in suite to estimate the probability of locking in the region of control parameter space with hysteresis, i.e., the set of control parameters for which both locked and unlocked states can exist. The results show that a neural network gives the best estimate of the locking probability. An analogy of the present locking model with the van der Waals equation of state is also provided.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Energetic particle physics: Chapter 7 of the special issue: on the path to tokamak burning plasma operation

We review the physics of energetic particles (EPs) in magnetically confined burning fusion plasmas with focus on advances since the last update of the ITER Physics Basis (Fasoli et al 2007 Nucl. Fusion 47 S264). Topics include basic EP physics, EP generation, diagnostics of EPs and instabilities, the interaction of EPs and thermal plasma instabilities, EP-driven instabilities, energetic particle modes (EPMs), and turbulence, linear and nonlinear stability and simulation of EP-driven instabilities and EPMs, 3D effects, scenario optimization strategies based on EP phase-space control, EPs in reduced field scenarios in ITER before DT, and the physics of runaway electrons. We describe the simulation and modeling of EPs in fusion plasmas, including instability drive and damping as well as EP transport, with a range of approaches from first-principles to reduced models, including gyrokinetic simulations, kinetic-MHD models, gyrofluid models, reduced models, and semi-analytical approaches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Electron-scale current sheets and energy dissipation in 3D kinetic-scale plasma turbulence with low electron beta

ABSTRACT Three-dimensional kinetic-scale turbulence is studied numerically in the regime where electrons are strongly magnetized (the ratio of plasma species pressure to magnetic pressure is βe = 0.1 for electrons and βi = 1 for ions). Such a regime is relevant in the vicinity of the solar corona, the Earth’s magnetosheath, and other astrophysical systems. The simulations, performed using the fluid-kinetic spectral plasma solver (sps) code, demonstrate that the turbulent cascade in such regimes can reach scales smaller than the electron inertial scale, and results in the formation of electron-scale current sheets (ESCS). Statistical analysis of the geometrical properties of the detected ESCS is performed using an algorithm based on the medial axis transform. A typical half-thickness of the current sheets is found to be on the order of electron inertial length or below, while their half-length falls between the electron and ion inertial length. The pressure–strain interaction, used as a measure of energy dissipation, exhibits high intermittency, with the majority of the total energy exchange occurring in current structures occupying approximately 20 per cent of the total volume. Some of the current sheets corresponding to the largest pressure–strain interaction are found to be associated with Alfvénic electron jets and magnetic configurations typical of reconnection. These reconnection candidates represent about 1 per cent of all the current sheets identified.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Effective Resistivity in Relativistic Collisionless Reconnection

Magnetic reconnection can power spectacular high-energy astrophysical phenomena by producing nonthermal energy distributions in highly magnetized regions around compact objects. By means of two-dimensional fully kinetic particle-in-cell (PIC) simulations, we investigate relativistic collisionless plasmoid-mediated reconnection in magnetically dominated pair plasmas with and without a guide field. In X-points, where diverging flows result in a nondiagonal thermal pressure tensor, a finite residence time for particles gives rise to a localized collisionless effective resistivity. Here, for the first time for relativistic reconnection in a fully developed plasmoid chain, we identify the mechanisms driving the nonideal electric field using a full Ohm law by means of a statistical analysis based on our PIC simulations. We show that the nonideal electric field is predominantly driven by gradients of nongyrotropic thermal pressures. We propose a kinetic physics motivated nonuniform effective resistivity model that is negligible on global scales and becomes significant only locally in X-points. It captures the properties of collisionless reconnection with the aim of mimicking its essentials in nonideal magnetohydrodynamic descriptions. This effective resistivity model provides a viable opportunity to design physically grounded global models for reconnection-powered high-energy emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Overview of the SPARC tokamak

The SPARC tokamak is a critical next step towards commercial fusion energy. SPARC is designed as a high-field ( $B_0 = 12.2$ T), compact ( $R_0 = 1.85$ m, $a = 0.57$ m), superconducting, D-T tokamak with the goal of producing fusion gain $Q>2$ from a magnetically confined fusion plasma for the first time. Currently under design, SPARC will continue the high-field path of the Alcator series of tokamaks, utilizing new magnets based on rare earth barium copper oxide high-temperature superconductors to achieve high performance in a compact device. The goal of $Q>2$ is achievable with conservative physics assumptions ( $H_{98,y2} = 0.7$ ) and, with the nominal assumption of $H_{98,y2} = 1$ , SPARC is projected to attain $Q \approx 11$ and $P_{\textrm {fusion}} \approx 140$ MW. SPARC will therefore constitute a unique platform for burning plasma physics research with high density ( $\langle n_{e} \rangle \approx 3 \times 10^{20}\ \textrm {m}^{-3}$ ), high temperature ( $\langle T_e \rangle \approx 7$ keV) and high power density ( $P_{\textrm {fusion}}/V_{\textrm {plasma}} \approx 7\ \textrm {MW}\,\textrm {m}^{-3}$ ) relevant to fusion power plants. SPARC's place in the path to commercial fusion energy, its parameters and the current status of SPARC design work are presented. This work also describes the basis for global performance projections and summarizes some of the physics analysis that is presented in greater detail in the companion articles of this collection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electromagnetic total- f algorithm for gyrokinetic particle-in-cell simulations of boundary plasma in XGC

We report the simplified δf mixed-variable/pullback electromagnetic simulation algorithm implemented in XGC for core plasma simulations by Cole et al. [Phys. Plasmas 28, 034501 (2021)] has been generalized to a total- f electromagnetic algorithm that can include, for the first time, the boundary plasma in diverted magnetic geometry with neutral particle recycling, turbulence, and neoclassical physics. The δf mixed-variable/pullback electromagnetic implementation is based on the pioneering work by Kleiber and Mischenko et al. [Kleiber et al., Phys. Plasmas 23, 032501 (2016); Mishchenko et al., Comput. Phys. Commun. 238, 194 (2019)]. An electromagnetic demonstration simulation is performed in a DIII-D-like, H-mode boundary plasma, including a corresponding comparative electrostatic simulation, which confirms that the electromagnetic simulation is necessary for a higher fidelity understanding of the electron particle and heat transport even at the low- β pedestal foot in the vicinity of the magnetic separatrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Spherical compression of an applied magnetic field in inertial confinement fusion

Applying an external magnetic field to laser-driven inertial confinement fusion implosions is a promising approach for enhancing fusion yield. The field is compressed with the plasma, producing a magnetized hotspot that anisotropically suppresses thermal losses and traps alpha particles, making performance sensitive to the compressed field orientation. Here, we derive a simple, readily applicable analytic model that enables rapid evaluation of the compressed field topology and show that ablation into the hotspot amplifies the central field, while the ablated ice near the hotspot edge develops a decaying, radially bent field, with a discontinuity in the field direction. The radially bent field renders thermal insulation at the hotspot edge negligible and largely independent of the applied field strength, whereas insulation in the hotspot core still depends strongly on the applied field. Applying the model to non-axial initial field configurations, we find that an initially applied mirror field provides the greatest suppression, followed by the standard axial field.

Physics - Plasma physics↗
DOE OSTI2023

New linear stability parameter to describe low-β electromagnetic microinstabilities driven by passing electrons in axisymmetric toroidal geometry

In magnetic confinement fusion devices, the ratio of the plasma pressure to the magnetic field energy, β, can become sufficiently large that electromagnetic microinstabilities become unstable, driving turbulence that distorts or reconnects the equilibrium magnetic field. In this paper, a theory is proposed for electromagnetic, electron-driven linear instabilities that have current layers localised to mode-rational surfaces and binormal wavelengths comparable to the ion gyroradius. The model retains axisymmetric toroidal geometry with arbitrary shaping, and consists of orbit-averaged equations for the mode-rational surface layer, with a ballooning space kinetic matching condition for passing electrons. The matching condition connects the current layer to the large scale electromagnetic fluctuations, and is derived in the limit that β is comparable to the square root of the electron-to-ion-mass ratio. Electromagnetic fluctuations only enter through the matching condition, allowing for the identification of an effective β that includes the effects of equilibrium flux surface shaping. The scaling predictions made by the asymptotic theory are tested with comparisons to results from linear simulations of micro-tearing and electrostatic microinstabilities in MAST discharge #6252, showing excellent agreement. In particular, it is demonstrated that the effective β can explain the dependence of the local micro-tearing mode (MTM) growth rate on the ballooning parameter θ 0 – possibly providing a route to optimise local flux surfaces for reduced MTM-driven transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On the role of mode resonances in regulating zonal-flow-moderated plasma microturbulence

Abstract The onset of turbulent heat transport at a higher temperature gradient than the critical gradient of linear instability, known as the Dimits shift, is a recurring feature of nonlinear simulations for magnetically confined fusion plasmas. Resonance in the nonlinear coupling between the modes that dominate energy transfer can lead to suppression of turbulence and transport above the linear critical gradient. As an expression of this resonance, gyrokinetic simulations show a quasi-coherent interaction between streamers and sidebands coupled through the zonal flow within the Dimits regime. This mechanism is further confirmed by use of artificial complex frequencies which break the resonance. By incorporating corresponding saturation physics, the standard quasilinear model for rapid head flux prediction is improved, which can now predict reduced heat flux in the Dimits regime. In particular, the triplet correlation time, the lifetime of the nonlinear interaction, is shown to be well-approximated by combinations of linear eigenvalues, and yields good representations of the heat flux variation both in and above the Dimits regime. Thus, a reduced but predictive model for transport near the critical gradient of zonal-flow saturated turbulence now exists.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Sheath constraints on turbulent magnetised plasmas

Abstract A solid target in contact with a plasma charges (negatively) to reflect the more mobile species (electrons) and thus keep the bulk plasma quasineutral. To shield the bulk plasma from the charged target, there is an oppositely (positively) charged sheath with a sharp electrostatic potential variation on the Debye length scale λ D . In magnetised plasmas where the magnetic field is inclined at an oblique angle α with the target, some of the sheath potential variation occurs also on the ion sound gyroradius length scale ρ S cos ⁡ α , caused by finite ion gyro-orbit distortion and losses. We consider a collisionless and steady-state magnetised plasma sheath whose thickness l ms ∼ max ( λ D , ρ S cos ⁡ α ) is smaller than the characteristic length scale L of spatial fluctuations in the bulk plasma, such that the limit l ms / L → 0 is appropriate. Spatial structures are assumed to be magnetic field-aligned. In the case of small magnetic field angle α ∼ δ ≡ ρ S / L ≪ 1 , electric fields tangential to the target transport ions towards the target via E × B drifts at a rate comparable to the one from parallel streaming. A generalised form of the kinetic Bohm–Chodura criterion at the sheath entrance is derived by requiring that the sheath electric field have a monotonic spatial decay far from the target. The criterion depends on tangential gradients of potential and ion distribution function, with additional nontrivial conditions.

kinetic Bohm–Chodura condition↗
DOE OSTI2023

A reduced-model (nSOLT) simulation of neutral recycling effects on plasma turbulence in the divertor region of MAST-U

The 2D scrape-off-layer turbulence code (nSOLT), which includes 1D kinetic neutral–plasma interactions, is applied to study effects of neutral recycling on plasma turbulence for parameters illustrative of the MAST-U divertor region. Neutral recycling is modeled by injecting a fraction of the parallel plasma flux to the divertor back into the simulation domain as a source of Franck–Condon-distributed neutrals. Stationary sources, concentrated at the magnetic separatrix, model plasma streaming into the divertor region from the upstream scrape-off-layer and sustain plasma turbulence absent neutral recycling. Starting from one such no-neutrals equilibrium, we initiate recycling in a numerical experiment designed to diagnose and identify the effects of various neutral–plasma interactions on the divertor plasma, divertor turbulence, and plasma exhaust. The onset of recycling triggers an initial burst of enhanced cross field plasma transport that is quelled by ionization cooling and charge–exchange (CX) friction, with growing neutral pressure, leading to a quiescent, turbulence-free state. Diagnosis of this transient burst reveals that (1) the sudden increase in plasma density due to ionization dominates the onset of the burst, (2) electron cooling due to ionization increases collisionality and disconnects blob filaments from the sheath, and (3) CX friction drives tripole polarization of a blob that can dominate the curvature-driven dipole polarization, leading to the stagnation of blob propagation and reduced radial turbulent transport. It is shown that CX friction is negligible compared to sheath physics in determining equilibrium mean flow shearing rates, for parameters considered herein (specifically a short connection length to the divertor target), while it can significantly reduce interchange-instability growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Microtearing stability and turbulence in the pedestal: Linear gyrokinetics, reduced models, and nonlinear turbulent transport

Microtearing modes can play a crucial role in electron heat transport in tokamak plasmas, affecting both energy confinement and overall performance. This study investigates microtearing modes (MTM) stability and turbulence in a JET pedestal through gyrokinetic simulations using the Gene code, complemented by a reduced eigenvalue model. The focus is on how MTM properties depend on key plasma parameters, including collisionality and plasma beta β—the ratio of plasma pressure to magnetic pressure—the normalized toroidal wavenumber k y ρ s ⁠, where ρ s denotes the ion sound gyroradius (typically a few millimeters in edge plasmas) and isotope mass. Collisionality enhances MT growth rates, while increasing β leads to a shift from MTMs to kinetic-ballooning modes, typically for k y ρ s ⁠, where ρ s ≲ 0.2⁠. A purely collisionless branch of MTMs persists at low k y ρ s ⁠, where ρ s with distinctive properties including non-negligible particle flux and ion thermal transport. Isotope mass scans reveal modest reduction of MTM growth rates as ion mass decreases. Nonlinear simulations produce experimentally relevant transport levels. Numerical experiments turning off zonal flows and fields identify the critical role of zonal flows and zonal fields in regulating MTM turbulence. Their removal leads to a significant increase in electron heat flux. These findings provide new insight into MTM-driven transport and its impact on tokamak confinement and lay a foundation for reduced modeling and predictive capabilities.

Electrostatics↗