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 451 records · Page 25

Computational general relativistic force-free electrodynamics: II. Characterization of numerical diffusivity

Scientific codes are an indispensable link between theory and experiment; in (astro-)plasma physics, such numerical tools are one window into the universe’s most extreme flows of energy. The discretization of Maxwell’s equations – needed to make highly magnetized (astro)physical plasma amenable to its numerical modeling – introduces numerical diffusion. It acts as a source of dissipation independent of the system’s physical constituents. Understanding the numerical diffusion of scientific codes is the key to classifying their reliability. It gives specific limits in which the results of numerical experiments are physical. We aim at quantifying and characterizing the numerical diffusion properties of our recently developed numerical tool for the simulation of general relativistic force-free electrodynamics by calibrating and comparing it with other strategies found in the literature. Our code correctly models smooth waves of highly magnetized plasma. We evaluate the limits of general relativistic force-free electrodynamics in the context of current sheets and tearing mode instabilities. We identify that the current parallel to the magnetic field (j ∥ ), in combination with the breakdown of general relativistic force-free electrodynamics across current sheets, impairs the physical modeling of resistive instabilities. We find that at least eight numerical cells per characteristic size of interest (e.g., the wavelength in plasma waves or the transverse width of a current sheet) are needed to find consistency between resistivity of numerical and of physical origins. High-order discretization of the force-free current allows us to provide almost ideal orders of convergence for (smooth) plasma wave dynamics. The physical modeling of resistive layers requires suitable current prescriptions or a sub-grid modeling for the evolution of j ∥ .

79 ASTRONOMY AND ASTROPHYSICS↗

Dissipation in the magnetic turbulence of reversed field pinch plasmas

Reversed field pinch plasmas are subject to tearing instability that creates a broad spectrum of magnetic fluctuations. The dominant fluctuations have poloidal and toroidal mode numbers (m,n) = (1,6 - 10) and can grow to 2 - 3% of the mean magnetic field. Through non-linear coupling this growth culminates in a strong reconnection event and broadening of the magnetic spectrum extending to the ion gyroradius scale. Multiple developments occur during the reconnection stage: ions and electrons are energized, magnetic fluctuation amplitudes increase, plasma flow is halted, and the toroidal magnetic flux increases in a sawtooth-like fashion as the RFP dynamo becomes stronger. Magnetic fluctuations are measured in the plasma edge at multiple radial locations from r/a = 0.75 to 0.96 to assess and characterize the magnetic turbulence. The measured spectrum perpendicular to the mean field, S(k ⊥ ), can be fit to a model spectrum consisting of power law and exponential component with one free parameter that characterizes dissipation. The measured dissipation is much larger than estimated from classical viscous or resistive dissipation, but it is consistent with a flow damping measurement of anomalous viscosity. The measurements show an evolution of the spectrum during which fluctuation power builds up in the smallest wavenumbers and cascades to the larger wavenumber due to the non-linear coupling between the linear (m,n) = (1, > 6) and the non-linear (m,n) = (0, 1) tearing mo

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

L-H transition trigger physics in ITER-similar plasmas with applied n = 3 magnetic perturbations

The L–H transition power threshold PLH is observed to increase with applied n = 3 resonant magnetic perturbations (RMP) in ITER-similar-shape plasmas with balanced neutral beam torque injection in DIII-D. The increase is most pronounced with added electron–cyclotron heating: PLH increases with decreasing edge plasma collisionality as PLH/PLH-08 ~ (ν*)−0.5, where PLH-08 is the 2008 ITPA multi-machine power threshold scaling. This result raises concerns for H-mode access at low edge collisionality in ITER, where RMP may have to be applied before the L–H transition to safely suppress the first edge-localized mode. Non-axisymmetric modifications with RMP include a simultaneous reduction of the radial electric field (Er) well depth and E × B shear. This can be attributed to increasing edge toroidal co-current rotation, and is consistent with substantially increased local long-wavelength turbulence (measured via beam emission spectroscopy). At high RMP perturbation strength the edge electric field Er reverses sign locally (becomes positive), with changes in dominant turbulence modes. Edge magnetic stochasticity provides an attractive explanation of the observed modifications, and the observed changes in toroidal rotation and Er are consistent with a simple fluid model describing radial electron current flow along stochastic fieldlines. The observed collisionality dependence of the L-mode edge electric field with applied RMP is also qualitatively consistent with this model. Reflectometry data indicate a significant reduction of the normalized L-mode radial density gradient a/Ln at high RMP field with simultaneous increase in radial particle flux and electron thermal flux from power balance analysis. We conjecture that the increase of PLH with RMP results from the combined effects of reduced E × B flow shear (increasing turbulent transport levels) and toroidal/poloidal flow modulation due to edge stochasticity. Initial experiments indicate that non-resonant n = 3 magnetic perturbations lead only to relatively small changes in Er, E × B shear and fluctuation characteristics, and have less impact on the L–H transition power threshold. This motivates further exploration of the RMP spectrum dependence of PLH for possible mitigation of the observed threshold increase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Comptonization by reconnection plasmoids in black hole coronae – III. Dependence on the guide field in pair plasma

ABSTRACT We perform non-radiative two-dimensional particle-in-cell simulations of magnetic reconnection for various strengths of the guide field (perpendicular to the reversing field), in magnetically dominated electron–positron plasmas. Magnetic reconnection under such conditions could operate in accretion disc coronae around black holes. There, it has been suggested that the transrelativistic bulk motions of reconnection plasmoids containing inverse-Compton-cooled electrons could Compton-upscatter soft photons to produce the observed non-thermal hard X-rays. Our simulations are performed for magnetizations 3 ≤ σ ≤ 40 (defined as the ratio of enthalpy density of the reversing field to plasma enthalpy density) and guide field strengths 0 ≤ Bg/B0 ≤ 1 (normalized to the reversing field strength B0). We find that the mean bulk energy of the reconnected plasma depends only weakly on the flow magnetization but strongly on the guide field strength – with Bg/B0 = 1 yielding a mean bulk energy twice smaller than Bg/B0 = 0. Similarly, the dispersion of bulk motions around the mean – a signature of stochasticity in the plasmoid chain’s motions – is weakly dependent on magnetization (for σ ≳ 10) but strongly dependent on the guide field strength – dropping by more than a factor of two from Bg/B0 = 0 to Bg/B0 = 1. In short, reconnection in strong guide fields (Bg/B0 ∼ 1) leads to slower and more ordered plasmoid bulk motions than its weak guide field (Bg/B0 ∼ 0) counterpart.

Astronomy & Astrophysics↗

Fusion gain from cylindrical liner-driven implosions of field reversed configurations

MagLIF experiments [M.R. Gomez et al., Phys. Plasmas 22, 056306 (2015)] on Z have demonstrated the basic principles of Magneto-Inertial Fusion (MIF) for wall confined plasmas. Other MIF schemes have been proposed based on the liner implosion of closed field magnetically confined plasmas such as Field Reversed Configurations (FRCs) [T. P. Intrator et al., Phys. Plasmas 15, 042505 (2008)]. We present a semi-analytical model of liner driven FRC implosions that predicts the fusion gain of such systems. The model predicts a fusion gain near unity for an FRC imploded by a liner driven by the Z Machine. We show that FRCs could be formed and imploded at the Z facility using the AutoMag liner concept [S. A. Slutz et al., Phys. Plasmas 24, 012704 (2017)]. Furthermore, an initial bias magnetic field can be supplied by the external magnets used in MagLIF experiments. The reverse field is then supplied by an AutoMag liner, which has helical conducting paths imbedded in an insulating substance. Experiments [Shipley et al., Phys. Plasmas 26, 052705 (2019)] have demonstrated that AutoMag can generate magnetic fields greater than 30 Tesla inside of the liner. We have performed 2D Radiation MHD simulations of the formation and implosion of an FRC, which are in good agreement with the analytical model. The FRC formation process could be studied on small pulsed power machines delivering about 1 MA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transport from electron-scale turbulence in toroidal magnetic confinement devices

Plasma transport driven by turbulence ultimately determines the energy confinement performance of controlled fusion devices regardless of their confinement schemes and configurations. A large variety of plasma instabilities have been proposed for driving turbulence responsible for anomalous plasma transport beyond classical/neoclassical transport due to collisions. Although ion-scale turbulence usually dominates due to its large eddy size and saturation level, electron-scale turbulence has been recognized to be important in regions where ion-scale turbulence is suppressed (e.g., in internal transport barrier and in spherical tokamak H-mode plasmas) or is close to marginality. Electron-scale turbulence has been shown to nonlinearly interact with ion-scale turbulence, which modifies the dynamics of both and affects the resulting plasma transport, particularly when ion-scale instability is weakly driven. In this review paper, we focus on electron-scale turbulence that is believed to operate in magnetically confinement fusion devices and aim to provide a review of theoretical, numerical, and experimental developments in understanding electron-scale turbulence and its role in driving anomalous plasma turbulence. In particular, we focus on the electrostatic electron temperature gradient (ETG) mode which is the most widely recognized plasma instability underlying electron-scale turbulence observed in magnetically confined plasmas. Finally, we note that there are other less studied instabilities that might be responsible for observed electron-scale turbulence, most notably ubiquitous mode, and short-wavelength ion temperature gradient (SWITG) mode, which will be briefly touched on in this review.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Uncovering turbulent plasma dynamics via deep learning from partial observations

One of the most intensely studied aspects of magnetic confinement fusion is edge plasma turbulence which is critical to reactor performance and operation. Drift-reduced Braginskii two-fluid theory has for decades been widely applied to model boundary plasmas with varying success. Towards better understanding edge turbulence in both theory and experiment, we demonstrate that a novel multi-network physics-informed deep learning framework constrained by partial differential equations can accurately learn turbulent fields consistent with the two-fluid theory from partial observations of electron pressure which is not otherwise possible using conventional equilibrium models. Furthermore, this technique presents a novel paradigm for the advanced design of plasma diagnostics and validation of magnetized plasma turbulence theories in challenging thermonuclear environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electromagnetic torque measurements in DIII-D using internal/external magnetic field decomposition

Given spatially resolved measurements of normal and tangential components of the magnetic field just outside the surface of a magnetically confined plasma, the field at the measurement location can be uniquely decomposed into contributions from the plasma and from external sources. This principle allows direct measurement of the electromagnetic torque on the plasma without knowledge of the distribution of the internal and external currents, similar to the more well-known formalism using the Maxwell stress tensor. The internal/external field decomposition also enables a mixed approach that incorporates any explicitly known current distributions (e.g., from non-axisymmetric coils). We discuss the requirements and limitations of such an approach to torque measurements. Experimental measurements of the torque evolution as a rotating tearing mode locks to the wall in the DIII-D tokamak are consistent with a simple model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measurement of Zero-Frequency Fluctuations Generated by Coupling between Alfvén Modes in the JET Tokamak

We report the first experimental detection of a zero-frequency fluctuation that is pumped by an Alfvén mode in a magnetically confined plasma. Core-localized Alfvén modes of frequency inside the toroidicity-induced gap (and its harmonics) exhibit three-wave coupling interactions with a zero-frequency fluctuation. The observation of the zero-frequency fluctuation is consistent with theoretical and numerical predictions of zonal modes pumped by Alfvén modes, and is correlated with an increase in the deep core ion temperature, temperature gradient, confinement factor H 89,P , and a reduction in the main ion heat diffusivity. Despite the energetic particle transport induced by the Alfvén eigenmodes, the generation of a zero-frequency fluctuation that can suppress the turbulence leads to an overall improvement of confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Efficient analysis of magnetic field line behavior in toroidal plasmas

The confinement of plasmas in tokamaks and stellarators depends on magnetic field lines lying in nested toroidal surfaces. The transition near the plasma edge away from the lines lying in magnetic surfaces defines properties of divertors. The transition in time defines properties of disruptions. Divertor design and disruption analyses require a detailed understanding of these transitions. The use of a Fourier transform coupled with a Gaussian window function allows far more information to be extracted about these transitions using far shorter field line integrations than can be obtained using traditional methods based on Poincaré plots. The physics of divertors and disruptions is reviewed to clarify why the type of information that can be gained from more efficient methods of analysis is of central importance to the fusion program based on magnetic confinement.

Boozer, Allen H. [Columbia Univ., New York, NY (Un↗

A Brief Overview of VPIC [Slides]

VPIC 2.0 is an open source, fully kinetic, fully relativistic, general purpose particle-in-cell (PIC) code for modeling a broad range of plasmas—from micron-scale short-pulse laser-plasma interactions to Earth-scale magnetospheric plasmas, including laser-plasma instability, magnetic flux ropes, and magnetic reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Parallel Algebraic Multigrid for Fusion and Higher-Order PDEs

Multigrid methods play a key role in large-scale scientific simulation because they are among the fastest and most scalable approaches for solving the underlying sparse linear systems of equations that arise from a wide array of Partial Differential Equation (PDE) discretizations. Algebraic multigrid (AMG) is a special type of multigrid method that depends only on the description of the linear system, giving it better portability and broader applicability than geometric multigrid, as it requires no explicit knowledge of the problem geometry. Even though these methods are widely used today, there are still applications where further development is needed. In this report, we focus on PDEs with higher-order terms (e.g., fourth order), concentrating on a PDE that arises in tokamak edge plasma simulations (a tokamak is a machine that confines a plasma using magnetic fields and is believed to be the leading plasma confinement concept for future fusion power plants). General multigrid relaxes a linear system on coarser grids and reverses this process with interpolation, but standard AMG methods struggle with the aforementioned higher-order PDEs. We investigate cyclic coarsening and interpolation heuristics, as well as new iterative approximation methods of refining the solution at each grid to improve the existing multigrid approach. To this end, we ensure that these techniques are transferable to a parallelized setting with LLNL’s supercomputers.

97 MATHEMATICS AND COMPUTING↗

Non-linear Plasma Wake Growth of Electron Holes

An object’s wake in a plasma with small Debye length that drifts across the magnetic field is subject to electrostatic electron instabilities. Such situations include, for example, the moon in the solar wind wake and probes in magnetized laboratory plasmas. The instability drive mechanism can equivalently be considered drift down the potential-energy gradient or drift up the density gradient. The gradients arise because the plasma wake has a region of depressed density and electrostatic potential into which ions are attracted along the field. The non-linear consequences of the instability are analysed in this paper. At

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radiative damping of toroidal Alfvén eigenmode in low-shear plasmas

Instabilities of Alfvén eigenmodes (AEs) are of significant concern because they can enhance the cross-field transport of fusion-born alpha particles beyond the neoclassical level in magnetic fusion plasmas. The threshold value of alpha-particle pressure for exciting AEs depends critically on the damping rate of AEs. The damping mechanisms include kinetic damping due to interactions with thermal particles, continuum damping due to AE frequency crossing Alfvén continuum, and radiative damping due to emitting kinetic Alfvén waves (KAWs). The radiative damping is substantial and can even prevail in high-temperature burning plasmas [1]. We revisit the radiative damping analytic theory for TAE in plasmas with low positive magnetic shear, considering TAE with an eigenfrequency near the bottom of TAE-gap and with poloidal harmonics of the same sign (even TAE). In contrast to earlier papers, we provide the damping calculations in real space rather than Fourier space. This approach is straightforward technically and more enlightening from a physics standpoint for benchmarking numerical calculations of radiative damping. The parametric dependence of the resulting damping rate agrees with that of Refs. [2-5], but it has a smaller numerical factor in front of it.

Alpha-particle driven instability↗