Search NASA⌕ Search

SEARCH · Search NASA

Results for “Plasma stability”

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

Toroidal modeling of anisotropic thermal transport and energetic particle effects on stability of resistive plasma resistive wall mode

Effects of anisotropic thermal transport on the linear stability of the resistive plasma resistive wall mode (RPRWM) are investigated by the magnetohydrodynamic-kinetic hybrid code MARS-K [Liu et al., Phys. Plasmas 15, 112503 (2008)], including the kinetic contribution from energetic particles (EPs). It is found that thermal transport can further stabilize the RPRWM in the presence of drift kinetic contributions from EPs. This is different from the fluid model, which always predicts destabilization of the mode by thermal transport. Furthermore, the thermal transport effect is found to amplify the role played by an adiabatic term, associated with the radial distribution of EPs’ birth energy, in modifying the mode stability as well as the mode eigenfunction. The shape of the equilibrium profiles of EPs, in particular that of the temperature, is also found to strongly affect the mode stabilization. This profile effect is more pronounced in the absence of thermal transport. MARS-K computations show that the stabilizing effect by thermal transport is more likely to occur at slower plasma rotation and lower EP energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Study of the Alfven eigenmodes stability in CFQS plasma using a Landau closure model

The aim of this study is to analyze the stability of the Alfven eigenmodes (AE) in the Chinese First Quasi-axisymmetric Stellarator (CFQS). The AE stability is calculated using the code FAR3d that solves the reduced MHD equations to describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moment for the energetic particles (EP) species including the effect of the helical couplings and acoustic modes. The Landau damping and resonant destabilization effects are added in the model by a given closure relation. The simulation results indicate the destabilization of n = 1 to 4 AEs by EP during the slowing down process, particularly n = 1 and n = 2 toroidal AEs (TAE), n = 3 elliptical AE (EAE) and n = 4 non circular AE (NAE). If the resonance is caused by EPs with an energy above 17 keV (weakly thermalized EP), n = 2 EAEs and n = 3 NAEs are unstable. On the other hand, EPs with an energy below 17 keV (late thermalization stage) lead to the destabilization of n = 3 and n = 4 TAEs. The simulations for an off-axis NBI injection indicate the further destabilization of n = 2 to 4 AEs although the growth rate of the n = 1 AEs slightly decreases, so no clear optimization trend with respect to the NBI deposition region is identified. In addition, n = 2, 4 helical AE (HAE) are unstable above an EP β threshold. Also, if the thermal β of the simulation increases (higher thermal plasma density) the AE stability of the plasma improves. Overall, the simulations including the effect of the finite Larmor radius and electron-ion Landau damping show the stabilization of the n = 1 to 4 EAE/NAEs as well as a decrease of the growth rate and frequency of the n = 1 to 4 BAE/TAEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Parallelized real-time physics codes for plasma control on DIII-D

A real-time safe multi-threading library was developed on the DIII-D plasma control system to optimize the real-time TORBEAM and real-time STRIDE physics codes. These physics codes are crucial for future fusion power plant operation as they provide information about electron cyclotron wave propagation and heating as well as inform about ideal plasma stability limits. The real-time TORBEAM code executed consistently in under 20 ms while the real-time STRIDE code computes in 100 ms. The multi-threading library developed in this work can be applied to other real-time physics-based codes that will be crucial for the next generation of fusion devices.

DIII-D↗
DOE OSTI2024

Implementation of a real-time MSE system

Motional Stark effect polarimetry is a key diagnostic for plasma fusion research since its usage on PBX-M. The MSE diagnostic measures the radial magnetic pitch angle profile in a plasma from a neutral beam by observation of Stark split D-alpha emission from atoms excited by collision with ions and electrons in the plasma. The pitch angle measurement is used with equilibrium reconstruction codes to determine the q-profile for studies of plasma stability, confinement, and transport. Historically, the algorithm was used in a post-processing fashion. The goal of our work was to apply this method in real time and pass the results to the plasma control system computer for real-time equilibrium reconstruction and control.

Instruments & Instrumentation↗
DOE OSTI2025

Error field measurements with rotating RMP fields for DIII-D H-mode

3D magnetic sensors are employed to identify the amplitude and toroidal phase of error fields (EF) by analyzing the torque balance for magnetic islands entrained by rotating resonant magnetic perturbations (RMPs) in DIII-D H-mode plasmas. This technique of torque balance allows for efficient error field identification, offering a valuable tool for scenario-specific and optimized error field compensation (EFC) and requires only magnetic diagnostics. The torque balance used in this work includes the contribution from electromagnetic torque due to error fields, wall response, RMP fields, and viscous. Results show that viscous torque plays a crucial role, particularly during locked modes and H-mode plasmas, ensuring accurate data fits with lower residuals. The torque balance technique reveals that the L- and H-mode plasmas have distinct EF configurations, and consistent with a model-based EF analysis including MHD response in IPEC and the SURFMN EF simulation. This technique shows great robustness in measuring the intrinsic EF amplitude regardless of its amplitude or toroidal phase. Repeated discharges with EFC disparities exhibit consistent results of intrinsic error field within a reasonable range near the “standard” error field compensation. Additionally, the use of a rotating n = 1 resonant magnetic perturbation offers the advantage of reducing disruption risks by entraining saturated magnetic islands. These findings are instrumental for optimizing EF correction in fusion devices, thereby enhancing tearing mode suppression and overall plasma stability.

3D magnetic sensors↗
DOE OSTI2023

Machine learning-based real-time kinetic profile reconstruction in DIII-D

Abstract Kinetic equilibrium reconstruction plays a vital role in the physical analysis of plasma stability and control in fusion tokamaks. However, the traditional approach is subjective and prone to human biases. To address this, the consistent automatic kinetic equilibrium reconstruction (CAKE) method was introduced, providing objective results. Nonetheless, its offline nature limits its application in real-time plasma control systems (PCSs). To address this limitation, we present RTCAKENN, a machine learning model that approximates 7 CAKE-level output profiles, namely pressure, inverse q , toroidal current density, electron temperature and density, carbon ion impurity temperature and rotation profiles, using real-time available inputs. The deep neural network consists of an encoder layer, where the scalars and interdependent inputs such as plasma boundary coordinates and motional Stark effect data are encoded using multi-layer perceptrons (MLPs), while profile inputs are encoded by 1D convolutional layers. The encoded data is passed through a MLP for latent feature extraction, before being decoded in the decoding layers, which consist of upsampling and convolutional layers. RTCAKENN has been implemented in the DIII-D PCS and our model achieves accuracy comparable to CAKE and surpasses existing real-time alternatives. Through clever dropout training, RTCAKENN exhibits robustness and can operate even in the absence of Thomson scattering data or charge exchange recombination data. It executes in under 8 ms in the real-time environment, enabling future application in real-time control and analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Routes to high-performance operation in Wendelstein 7-X: turbulence suppression with shaping of the density profile

Steep density gradients generally lead to improved plasma performance in the neoclassically optimized stellarator W7-X. This is evident in the global energy confinement time as well as in the ion temperature and can be explained by a strong reduction of the ion temperature gradient turbulence. Such conditions can experimentally be realized by several methods: injection of cryogenic hydrogen pellets, appropriate combination of neutral beam and electron cyclotron resonance heating (ECRH) and, in some cases, with low power of ECRH after preconditioning of the first wall. The duration of the improved phases is determined by the ability to sustain the steep density gradient, by technical limitations of the involved systems and, eventually, by the plasma stability. This paper gives an overview of relevant experimental results and presents example discharges where the improved confinement conditions could be extended to multiple seconds: up to 4 s using neutral beam injection and from 14 to 40 s with steady state pellet injection. In these plasmas the turbulent thermal diffusivity is reduced by a factor of 3 to 4 in a broad radial range, which allows high ion temperatures of up to 3 keV at the densities of about 1.5 • 10 20 m −3 .

high performance plasma↗
DOE OSTI2023

Conceptual design of the LHCD system on CFETR

With the capability of saving flux consumption in the current ramp up phase, controlling the safety factor (q) profile, providing the required off-axis current drive (CD), a lower hybrid current drive (LHCD) system was eventually determined to be used on the China Fusion Engineering Test Reactor (CFETR). Here, in this paper, a 20 MW/4.6 GHz system composed of 40 units of 500 kW continuous wave (CW) klystron amplifier is preliminarily designed. The feasibility of the frequency choice of 4.6 GHz in the physical aspect is discussed. In order to minimize the transmission loss, a TE01 over-mode of circular waveguide will be used in the transmission system. Although previous calculations indicates that high field side (HFS) launcher location results in waves damping at inner region than low field side (LFS) launcher, which is favorable for plasma stability, HFS launcher will be extremely difficult in engineering due to the limited space. Calculations show that the tritium breeding ratio (TBR) will be decreased by ~ 1.5% with HFS launcher, while it will be decreased by ~ 0.22% only with LFS launcher. As a result, the LH power will be coupled to plasma from the top port at LFS. The designed passive active multi-junction (PAM) launcher is arranged in an array of 5 rows and 8 columns with the height of 1181 mm and the width of 794 mm. A shielding block with a thickness of 20 cm will be equipped around the feeding waveguides to protect the feeding waveguides and to prevent neutron leakage. A preliminary scheme of the remote maintenance for the LFS launcher is given. The modeling results shows that the density for optimum coupling is ~ 1.8 × 10 17 /m 3 , lower than the 4.6 GHz cut-off density n e_co = 2.6 × 10 17 /m 3 . Around this density, the power directivity (D p ) can be high as 71%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The dual-electron cyclotron emission based measurement of 3D structures on DIII-D tokamak

The measurements of the locked mode and other 3D magnetohydrodynamic modes are essential for real-time control of plasma stability in future fusion reactors. Current magnetic measurements present non-negligible challenges for the application on fusion pilot plants (FPPs) due to neutron radiation and pose challenging design issues. In this paper, we discuss an alternative solution with microwave-based measurements. The dual-electron cyclotron emission (ECE) method has been developed and applied on DIII-D. The dual-ECE measurement method demonstrates clear identification of magnetic islands compared to single ECE, which is highly beneficial for locked mode feedback control on the FPPs. Toroidally separated ECE measurements provide accurate three-dimensional structure of magnetic islands and tearing modes. This technique also offers deeper physical insights into the penetration of resonant magnetic perturbations. Numerical calculations indicate that small locked modes (magnetic island width > 2 cm) can be detected using the dual-ECE measurement method in ITER plasma. This will enhance measurement capabilities for the safe operation of FPP plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Transients using low-high pulsed power in inductively coupled plasmas

Pulsed inductively coupled plasmas (ICPs) are widely deployed in the fabrication of semi-conductor devices. Pulse repetition frequencies (PRFs) of up to tens of kHz are commonly used during plasma etching for the high power densities they generate during the pulse-on period, and for their unique chemistries during the pulse-off period. The use of highly attaching halogen gases produces low electron densities during the pulse-off period, and these low density can result in instabilities, E-H transitions and ignition delays when applying power on the next pulse. To mitigate these possibilities, a low-level power environment could be maintained during “pulse-off” to moderate the minimum plasma density, therefore reducing ignition delays and enhancing plasma stability. In this work, ICPs sustained by 5 kHz pulsed power using Ar/Cl 2 mixtures at 20 mTorr were computationally investigated. For these conditions, the computed electron temperature (T e ) reaches a quasi-steady state during both the high- and low-power excitation. The model predicts that within the electromagnetic skin-depth, T e spikes to a high value during a low-to-high power modulation, and to a low value during a high-to-low power modulation. At the same time, a few cm above the substrate, there is little modulation in T e , as electron power convected from the skin depth disperses in traversing the reactor. Finally, the positive and negative spikes, and convection of transients across the reactors, are functions of power ramping time and gas mixtures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

MCTrans++: a 0-D model for centrifugal mirrors

The centrifugal mirror confinement scheme incorporates supersonic rotation of a plasma into a magnetic mirror device. This concept has been shown experimentally to drastically decrease parallel losses and increase plasma stability as compared with prior axisymmetric mirrors. MCTrans++ is a dimensionless (0-D) scoping tool which rapidly models experimental operating points in the Centrifugal Mirror Fusion Experiment (CMFX) at the University of Maryland. In the low-collisionality regime, parallel losses can be modelled analytically. A confining potential is set up that is partially ambipolar and partially centrifugal. Due to the stabilizing effects of flow shear, the perpendicular losses can be modelled as classical. Radiation losses such as bremsstrahlung and cyclotron emission are taken into account. A neutrals model is included, and, in some circumstances, charge-exchange losses are found to exceed all other loss mechanisms. We use the SUNDIALS ARKODE library to solve the underlying equations of this model; the resulting software is suitable for scanning large parameter spaces, and can also be used to model time-dependent phenomena such as a capacitive discharge. MCTrans++ has been used to verify results from prior centrifugal mirrors, create an experimental plan for CMFX and find configurations for future reactor-scale fusion devices.

Schwartz, Nick R. (ORCID:0000000296846435)↗
DOE OSTI2020

Hamiltonian formulations for perturbed dissipationless plasma equations

The Hamiltonian formulations for the perturbed Vlasov–Maxwell equations and the perturbed ideal magnetohydrodynamics (MHD) equations are expressed in terms of the perturbation derivative ∂ F / ∂ ϵ ≡ [ F , S ] of an arbitrary functional F [ ψ ] of the Vlasov–Maxwell fields ψ = ( f , E , B ) or the ideal MHD fields ψ = ( ρ , u , s , B ) , which are assumed to depend continuously on the (dimensionless) perturbation parameter ϵ . In this study, [ , ] denotes the functional Poisson bracket for each set of plasma equations and the perturbation action functional S is said to generate dynamically accessible perturbations of the plasma fields. The new Hamiltonian perturbation formulation introduces a framework for functional perturbation methods in plasma physics and highlights the crucial roles played by polarization and magnetization in Vlasov–Maxwell and ideal MHD perturbation theories. One application considered in this paper is a formulation of plasma stability that guarantees dynamical accessibility and leads to a natural generalization to higher-order perturbation theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Nonlinear coupling of tearing modes in reversed field pinch plasmas with stepped pressure profiles

A theory of the three-wave coupling of triplets of tearing modes in toroidal pinches [i.e., either reversed field pinches (RFPs) or tokamaks] was proposed by R. Fitzpatrick [Phys. Plasmas 6, 1168 (1999)]. However, this theory only applies to toroidal pinches with negligible equilibrium plasma pressure gradients. Such a limitation is particularly inappropriate to RFPs. This paper generalizes the analysis of R. Fitzpatrick [Phys. Plasmas 6, 1168 (1999)] in order to take the equilibrium pressure gradient into account. However, for the sake of simplicity, a stepped pressure profile, rather than a continuous profile, is employed. In the limit in which the number of steps becomes very large, the results obtained from the generalized theory are presumably equivalent to those that would have been achieved using a continuous pressure profile. The generalized theory is used to investigate the formation of the characteristic toroidally localized pattern of phase-locked m = 1 and m = 0 tearing modes in RFP plasmas that is known as the “slinky” pattern. The incorporation of the equilibrium plasma pressure into the analysis is found to be of crucial importance when determining the properties of the pattern. This is the case because the plasma pressure controls the number of unstable m = 1 and m = 0 tearing modes, and also significantly affects the strength of three-wave coupling, as well as the phase relation between the phase-locked m = 1 and m = 0 modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Nonlinear ablative Rayleigh-Taylor instability: Increased growth due to self-generated magnetic fields

The growth rate of the nonlinear ablative Rayleigh-Taylor (RT) instability is enhanced by magnetic fields self-generated by the Biermann battery mechanism; a scaling for this effect with perturbation height and wavelength is proposed and validated with extended-magnetohydrodynamic simulations. The magnetic flux generation rate around a single RT spike is found to scale with the spike height. The Hall parameter, which quantifies electron magnetization, is found to be strongly enhanced for short-wavelength spikes due to Nernst compression of the magnetic field at the spike tip. Furthermore, the impact of the magnetic field on spike growth is through both the suppressed thermal conduction into the unstable spike and the Righi-Leduc heat flow deflecting heat from the spike tip to the base. Righi-Leduc is found to be the dominant effect for small Hall parameters, while suppressed thermal conduction dominates for large Hall parameters. These results demonstrate the importance of considering magnetic fields in all perturbed inertial confinement fusion hot spots.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electron Modulational Instability in the Strong Turbulent Regime for an Electron Beam Propagating in a Background Plasma

We study collective processes for an electron beam propagating through a background plasma using simulations and analytical theory. A new regime where the instability of a Langmuir wave packet can grow locally much faster than ion frequency is clearly identified. The key feature of this new regime is an electron modulational instability that rapidly creates a local Langmuir wave packet, which in its turn produces local charge separation and strong ion density perturbations because of the action of the ponderomotive force, such that the beam-plasma wave interaction stops being resonant. Three evolution stages of the process and observed periodic burst features are discussed. In conclusion, different physical regimes in the plasma and beam parameter space are demonstrated for the first time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effect of negative triangularity on tearing mode stability in tokamak plasmas

Abstract The influence of negative triangularity (NT) of the plasma shape on the n = 1 (n is the toroidal mode number) tearing mode (TM) stability has been numerically investigated, with results compared to that of the positive triangularity (PT) counterpart. By matching the safety factor profile for a series of toroidal equilibria, several important plasma parameters, including the triangularity, the plasma equilibrium pressure, the plasma resistivity as well as the toroidal rotation, have been varied. The TM localized near the plasma edge is found to be more unstable in the NT plasmas as compared to the PT counterpart. The fundamental reason for this difference is the lack of favorable average curvature stabilization in NT configurations. Direct comparison of the Mercier index corroborates this conclusion. For the core-localized mode, where the difference in the local triangularity between NT and PT becomes small and the curvature stabilization is significantly reduced, larger Shafranov shift in the plasma core associated with the NT configuration results in more stable TM. The plasma toroidal flow generally stabilizes the TM in plasmas with both NTs and PTs. The flow stabilization is however weaker in the case of negative triangularity with finite plasma pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A generalized effective potential for differentially rotating plasmas

Global stability of differentially rotating plasma is investigated using a generalized effective potential. We first, for a current-free system, obtain a general form of an effective potential in terms of the free energies of global curvature and gradients of rotation for non-axisymmetric disturbances. We then examine the stability of differentially rotating disks for several rotation profiles and present the associated effective potential for the onset of these instabilities in the MHD regime. In particular, results for global axisymmetric magnetorotational instability as well as local and global non-axisymmetric modes are presented. The latter constitute two distinct non-axisymmetric modes, a high frequency local MRI and a global low-frequency non-axisymmetric mode (the magneto-curvature mode, introduced in Ebrahimi and Pharr [Astrophys. J. 936, 145 (2022)]), confined either between two Alfvénic resonances or an Alfvénic resonance and a boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Kinetic equilibrium reconstruction and the impact on stability analysis of KSTAR plasmas

High fidelity kinetic plasma equilibrium reconstructions are an essential requirement for accurate stability and disruption prediction analyses to support continuous operation of high beta tokamak plasmas. The present kinetic equilibrium reconstructions of plasmas in the KSTAR device include plasma density and temperature profiles from Thomson scattering and ion temperature from charge exchange spectroscopy diagnostics, and allowance for fast particle pressure. In addition, up to 25 channels of motional Stark effect diagnostic data are used to constrain the magnetic field pitch angle profile in the plasma to produce a reliable computation of the safety factor, q, profile. In this work, H-mode plasmas exhibit clear pedestal characteristics in the reconstructed pressure profile compared to internal transport barrier or L-mode plasmas. The plasma configuration and vertical position of inner strike points are validated by CCD and infrared camera images. Ideal and resistive magnetohydrodynamic (MHD) stability analyses using the DCON and resistive DCON codes utilize these kinetic equilibrium reconstructions to compare to experimental plasma stability. Equilibria with sufficiently low convergence error can provide reliable computation of ideal and resistive MHD stability analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗