Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma confinement”

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 307 records · Page 17

An ideal MHD δW stability analysis that bypasses the Newcomb equation

In the following work, we demonstrate the efficacy of a Riccati δW ideal MHD stability analysis that bypasses the numerically intractable integration of the Newcomb equation. By transforming the linear Newcomb equation into a quadratic Riccati equation, an accurate and equivalent δW analysis is performed that is shown to enjoy some numerical advantages. Additionally, we demonstrate that the Riccati approach is better conditioned than its Newcomb counterpart at the magnetic axis, and we apply dynamical systems insights to examine its behavior at singular points. We further discuss the constraints involved in solving for the extremal admissible perturbations of a δW analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Toroidal modeling of thermal particle drift kinetic effects and sub-sonic plasma flow on internal kink mode

The stability of the n = 1 ( n is the toroidal mode number) internal kink mode in a tokamak plasma is numerically investigated, utilizing the full toroidal, ideal magnetohydrodynamic (MHD) code MARS-F [Y. Q. Liu et al. , Phys. Plasmas 7, 3681 (2000)] and the nonperturbative MHD-kinetic hybrid code MARS-K [Y. Liu et al. , Phys. Plasmas 15, 112503 (2008)]. This study focuses on two physics effects: the sheared toroidal flow of plasma and the drift kinetic effects from thermal particles (ions and electrons) on the internal kink instability. Within the subsonic flow assumption, the flow and flow shear effects on the internal kink are generally weak, but the kinetic effects can be strong. A significant reduction in the mode growth rate is predicted by the nonperturbative MHD-kinetic hybrid computations, when the precessional drift resonance contributions from both particle species are included and when the mode is not too strongly unstable according to the fluid theory. The stabilization/destabilization of sheared toroidal flow depends on the radial location of the local flow shear. Strong destabilization (stabilization) occurs when a negative (positive) flow shear is located near the q = 1 rational surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Current channel evolution in ideal Z pinch for general velocity profiles

Recent diagnostic advances in gas-puff Z pinches at the Weizmann Institute for the first time allow the reconstruction of the current flow as a function of time and radius. These experiments show an unexpected radially-outward motion of the current channel, as the plasma moves radially-inward. In this paper, a mechanism that could explain this current evolution is described. Here, we examine the impact of advection on the distribution of current in a cylindrically symmetric plasma. In the case of metric compression, |v r | ∝ r, the current enclosed between each plasma fluid element and the axis is conserved, and so the current profile maintains its shape. We show that for more general velocity profiles, this simple behavior quickly breaks down, allowing for non-conservation of current in a compressing conductor, rapid redistribution of the current density, and even for the formation of reverse currents. In particular, a specific inward radial velocity profile is shown to result in radially-outward motion of the current channel, recovering the surprising current evolution discovered at the Weizmann Institute.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation of supersonic molecular beam injection fueling into H-mode plasmas on EAST using BOUT++

Transport dynamic 2D simulations for supersonic molecular beam injection (SMBI) fueling into H-mode deuterium plasmas on the Experimental Advanced Superconducting Tokamak (EAST) is first simulated using a seven-field two-fluid model in the BOUT++ framework. The SMB is assumed to be injected into plasma from the midplane at the low field side with a fixed width and constant molecular flux. The different densities and injection velocities of SMBI are investigated within the upper single-null geometry of EAST. The simulations indicate that the SMBI has a self-shielding effect on molecules' inward transport into the plasma, and the deposition of SMBI leads to a large increase in plasma density and decrease in plasma temperature. There is a velocity threshold for SMB penetrating the pedestal and depositing at the top of density pedestal. However, the deposition point would be back toward the plasma boundary after SMB arriving at the deepest penetration position. Comparing the different molecular injection velocities and densities, the outcomes demonstrate that the depth of deposition is closely related to the injection velocity rather than the molecular density. The simulated results show good agreement with the EAST experiments by comparing the electron density profiles obtained by simulation and experiment, respectively. These results will be helpful for guiding future experiments, as well as the design of the SMBI system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coherence imaging spectroscopy at Wendelstein 7-X for impurity flow measurements

In the last decade, Coherence Imaging Spectroscopy (CIS) has shown distinctive results in measuring ion flow velocities in the edge of magnetically confined plasma devices. Its 2D spatially resolved measurement capabilities and its high optical throughput are ideal for investigating the impurity behavior in the complex 3D magnetic island topology edge of Wendelstein 7-X (W7-X). However, a highly precise and stable calibration method is required for a reliable diagnostic operation. A new level of precision and stability has been achieved for the two CIS systems installed at W7-X with the use of a new calibration source, a continuous tunable laser commercially available only since 2015. A specific prototype model was successfully adapted to the challenging requirements of W7-X, granting high accuracy (±0.01 pm) and flexibility (spectral range: 450–650 nm) in the wavelength calibration required for measuring low-Z impurity ion flow velocities. These features opened up new investigation possibilities on temperature stability and wavelength response of the CIS components, allowing to fully characterize and validate the W7-X systems. The CIS diagnostic was operational throughout the last W7-X experimental campaign. Measured velocities on the order of ~20–30 km/s were observed, corroborated by comparisons with measurements with Mach probes.

47 OTHER INSTRUMENTATION↗

Gyrokinetic particle simulations of interactions between energetic particles and magnetic islands induced by neoclassical tearing modes

Interactions between energetic particles (EPs) and neoclassical tearing mode (NTM) islands in the DIII-D tokamak are studied using the global gyrokinetic toroidal code (GTC). GTC simulations find that the EP radial profile is partially flattened within the magnetic island regions and that there are stochastic regions in the particle phase space. Radial particle flux is induced mainly around the magnetic island regions and decreases with time to almost zero when the initial EP distribution achieves a new steady-state in the absence of EP sources. Stochastic regions of magnetic field lines induced by the superposition of multiple islands have weak effects on the particle flux when the width of stochastic regions is smaller than the EP drift orbit width. In conclusion, the perturbed parallel EP current induced by the magnetic islands has weak stabilizing effects on the linear growth rate of the NTM instability in this DIII-D experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental characterization of a section of a spherically imploding plasma liner formed by merging hypersonic plasma jets

In this work, we report experimental results on merging of hypersonic plasma jets, which is the fundamental building block for forming spherically imploding plasma liners as a potential standoff compression driver for magneto-inertial fusion. Jets are formed and launched by contoured-gap coaxial plasma guns mounted at the six vertices and the center of a hexagon covering approximately one-tenth of the surface area of a 9-ft-diameter spherical chamber. First, from experiments with two and three merging jets of four different species (N, Ar, Kr, and Xe), we show that (1) density spatial non-uniformities can be large (with electron-density jumps ranging from 2.9 for N to 6.6 for Xe) when shocks form upon jet merging, but smaller (density jumps <2) when shocks do not form; (2) jet impurities (20% Ti in these experiments) can increase the level of density spatial non-uniformity by increasing the collisionality of jet merging, leading to shock formation rather than potentially more desirable shockless jet merging; and (3) the liner Mach number can remain high (≳10), as required for plasma liners to be an effective compression driver. Secondly, from experiments with six and seven merging jets using Ar, we present results with improved jet-to-jet mass balance of <2% across jets, including (1) evidence of substantially increased balance in the jet merging and symmetry of the liner structure and (2) potentially favorable changes in the jet-merging morphology with the addition of the seventh jet. Ultimately, for both experiments, we present comparisons between experimental and synthetic data from three-dimensional hydrodynamic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Axisymmetric simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

A benchmark exercise for the modeling of vertical displacement events (VDEs) is presented and applied to the 3D nonlinear magnetohydrodynamic codes M3D-C 1 , JOREK, and NIMROD. The simulations are based on a vertically unstable NSTX equilibrium enclosed by an axisymmetric resistive wall with a rectangular cross section. A linear dependence of the linear VDE growth rates on the resistivity of the wall is recovered for sufficiently large wall conductivity and small temperatures in the open field line region. We report the benchmark results show good agreement between the VDE growth rates obtained from linear NIMROD and M3D-C 1 simulations and from the linear phase of axisymmetric nonlinear JOREK, NIMROD, and M3D-C 1 simulations. Axisymmetric nonlinear simulations of a full VDE performed with the three codes are compared, and an excellent agreement is found regarding the plasma location and plasma currents, as well as eddy and halo currents in the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computational analysis of ion orbital loss in diverted positive- and negative-triangularity tokamaks

The mechanism of ion orbital loss in diverted negative-triangularity (NT) tokamaks is significantly different from that of positive-triangularity (PT) tokamaks. The orbital loss from the edge of NT tokamaks can be quite large, because the separatrix X-point is located on the outboard side, meaning that trapped ions crossing the separatrix are lost to the inboard divertor plate after a mirror reflection on the side of the strong magnetic field. NT tokamaks have difficulties in building up steep pressure gradients in the edge regions. However, the majority of trapped ions crossing the separatrix in PT tokamaks survive (through complete banana motions) without reaching the X-point region. For NT tokamaks, topological changes in the trapped ion trajectories can be achieved in the presence of the radial electric field, thus reducing the orbital loss, whereas for PT tokamaks, the orbital loss is relatively insensitive to the radial electric field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear global gyrokinetic delta- f turbulence simulations in a quasi-axisymmetric stellarator

We use the global gyrokinetic stellarator code XGC-S to look at turbulence physics in geometry of the proposed quasi-axisymmetric stellarator, QUASAR. In this article, the XGC-S code is first verified for nonlinear simulations of ion temperature gradient-driven turbulence by a successful benchmark with the established nonlinear tokamak code XGC1 in a CYCLONE-like tokamak case. We then apply the XGC-S code to QUASAR for nonlinear turbulence physics. We find that the highly poloidally localized linear mode gives way to a more evenly distributed nonlinear turbulent structure. Furthermore, the radial extent of the mode increases and occupies a large fraction of the volume, despite the initially localized temperature profile, as the profiles relax.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhanced particle flux due to localized divertor MHD instability in DIII-D tokamak

An MHD mode with a frequency of <10 kHz has been identified near the inner strike point from various diagnostics, i.e., divertor Langmuir probes, magnetics sensors, and interferometers, but does not appear in the upstream and core diagnostics. This MHD mode is associated with magnetic oscillations of ≳5 G, has a long wavelength in the toroidal direction with toroidal mode number n=1, but is localized in a narrow radial region of a few cm. The mode appears when the outer strike point remains attached and the inner strike point nearly detaches, grows with increasing density, and eventually weakens and vanishes as the outer strike point detaches. This mode results in particle flux with an order of magnitude higher than the background plasmas near the inner strike point. The mode characteristics are consistent with the Current-Convective-Instability theory prediction. Initial simulations based on experimental input have found oscillations with similar frequencies but weaker amplitude.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Resolving ECRH deposition broadening due to edge turbulence in DIII-D

Microwave heat pulse propagation experiments have demonstrated a correlation between millimeter-scale turbulence and deposition profile broadening of electron cyclotron (EC) waves on the DIII-D tokamak. In a set of discharges in DIII-D, a variation in edge density fluctuations on the mm-scale is associated with 40%–150% broader deposition profiles, expressed in terms of normalized minor radius, as compared with equilibrium ray tracing. The 1D power profile is determined from transport analysis of the electron temperature response to EC power modulation using perturbative analysis with a square wave power modulation at 20–70 Hz, producing a series of Fourier harmonics that are fit collectively to resolve transport. Fitting an integrated heat flux expressed in the Fourier basis of the modulation to diffusive, convective, and coupled transport terms in a linear model can resolve the broadened EC deposition width from the power perturbation to resolve a broadening in each case. The best fit degree of beam broadening observed scales approximately linearly with the Doppler backscattering measured fluctuation level in the steep gradient region. Quantifying the effect of edge fluctuation broadening on EC current drive power needs of future devices will require 3D full-wave codes that can be validated on the current generation of machines. These DIII-D experiments provide a quantitative measure of fluctuation effects and a dataset to benchmark full-wave simulations that can model and eventually predict nonlinear effects neglected by 1D equilibrium beam and ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear δf particle in cell gyrokinetic simulations of the microtearing mode

Nonlinear simulations are carried out for the microtearing mode using particle-based δf gyrokinetic simulations for parameters relevant to spherical tokamaks. The present study finds that the microtearing mode can generate significant electron heat flux, which is predominantly carried out by the electromagnetic component of the heat flux with a negligible contribution from the electrostatic component. The mode sustains without the electrostatic component. We observe that the electron heat flux increases with the electron temperature gradient. The heat flux exhibits a very weak dependence on the collisions. It increases with electron β initially; however, at very high β, the electron heat flux is reduced.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear error field response in the presence of plasma rotation and real frequencies due to favorable curvature

This works describes the interaction of an phenomenon known to arise in tokamak fusion experiments with the magnetic perturbations that can arise at the edge of such experiments. We present nonlinear NIMROD resistive MHD simulations of the response of a rotating plasma to an error field when the plasma has weakly damped linear tearing modes (TMs), stabilized by a pressure gradient and favorable curvature. The favorable curvature leads to the Glasser effect: the occurrence of real frequencies and stabilization with positive stability index Δ'.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Direct imaging of plasma waves using ultrafast electron microscopy

A femtosecond plasma imaging modality based on a new development of ultrafast electron microscope is introduced. We investigated the laser-induced formation of high-temperature electron microplasmas and their subsequent non-equilibrium evolution. Based on a straightforward field imaging principle, we directly retrieve detailed information about the plasma dynamics, including plasma wave structures, particle densities, and temperatures. We discover that directly subjected to a strong magnetic field, the photo-generated microplasmas manifest in novel transient cyclotron echoes and form new wave states across a broad range of field strengths and different laser fluences. Intriguingly, the transient cyclotron waves morph into a higher frequency upper-hybrid wave mode with the dephasing of local cyclotron dynamics. The quantitative real-space characterizations of the non-equilibrium plasma systems demonstrate the feasibilities of a new microscope system in studying the plasma dynamics or transient electric fields with high spatiotemporal resolutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

RF current condensation in the presence of turbulent enhanced transport

Sharp temperature gradients in a magnetically confined plasma can lead to turbulent motion of the plasma. This turbulence in turn enhances the transport of heat across magnetic field lines. The enhanced transport impacts the temperature differential that can be sustained in magnetic islands between the island center and its periphery. It is shown in this work that, by limiting this temperature differential, this enhanced transport can have a profound influence on the extent to which the RF current condensation effect stabilizes the island growth. Interestingly, because the heat transport is no longer simply linear in the temperature gradient, the RF current condensation effect also exhibits entirely new hysteresis phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theory of edge localized mode suppression by static resonant magnetic perturbations in the DIII-D tokamak

According to a recent paper [Hu et al ., Phys. Plasmas 26, 120702 (2019)], mode penetration at the top of the pedestal is a necessary and sufficient condition for the suppression of edge localized modes (ELMs) by a resonant magnetic perturbation (RMPs) in an H-mode tokamak discharge. This paper employs asymptotic matching theory to model a particular DIII-D discharge in which ELMs were suppressed by an externally generated, static, n ¼ 2, RMP whose amplitude was modulated at a frequency of 1 Hz. It is demonstrated that the response of the plasma to the applied RMP, in the immediate vicinities of the rational (i.e., resonant) surfaces, is governed by nonlinear, rather than by linear, physics. This is the case because the magnetic island widths associated with driven reconnection exceed the linear layer widths, even in cases where driven reconnection is strongly suppressed by plasma rotation. The natural frequency at a given rational surface (i.e., the helical frequency at which the locally resonant component of the RMP would need to propagate in order to maximize driven reconnection) is found to be offset from the local E B frame in the ion diamagnetic direction. The size of the offset is mostly determined by neoclassical poloidal rotation. Finally, the predictions of a fully nonlinear plasma response model are found to be broadly consistent with the DIII-D experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗