SEARCH · Search NASA

Results for “peeling ballooning modes”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A study of resistive peeling–ballooning modes across spherical tokamaks

We investigate how non-ideal-magnetohydrodynamics (MHD) effects, in particular plasma resistivity, impact the peeling–ballooning stability thresholds in spherical tokamaks. This analysis follows the discovery of resistive kink-peeling modes in ELMing National Spherical Torus Experiment (NSTX) discharges. In the present study we extend this modeling to ELMing pulses in the Mega Ampere Spherical Tokamak (MAST) and MAST—Upgrade (MAST-U), where we find a clear resistive scaling for peeling–ballooning modes. While in NSTX ideal-MHD predicts stability for ELMing discharges, in MAST-U we find that the plasma is slightly unstable to peeling–ballooning modes, but is fully stabilized once diamagnetic effects are considered in terms of a growth rate normalization. A resistive power law scaling is calculated for these modes on MAST-U, which lies in between that of tearing modes and resistive interchange modes. A comparison between M3D-C1 and NIMROD shows reasonable agreement for this scaling. Resistivity destabilizes the modes and the peeling–ballooning unstable domain is considerably expanded in both, MAST and MAST-U. In addition to the MAST/-U pulses we also analyze resistive PB stability in a NSTX-similarity discharge on DIII-D. While having a different aspect ratio from NSTX, this discharge uses NSTX-like shaping parameters, toroidal field and plasma current. By considering these discharges alongside NSTX cases, we identify conditions influencing the onset of resistive peeling–ballooning modes. Furthermore, our findings indicate that magnetic shear in the pedestal region is closely linked to the emergence of resistive edge modes.

MAST-U↗

DOE OSTI2025

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗

DOE OSTI2025

Data for "A study of resistive peeling-ballooning modes across spherical tokamaks"

The data in this repository correspond to the figures of the paper "A study of resistive peeling-ballooning modes across spherical tokamaks" published in the journal Plasma Physics and Controlled Fusion. The figures are created based on simulations with the M3D-C1 and NIMROD codes using Python post-processing tools.

edge localized modes↗

DOE OSTI2026

Impact of impurities on peeling–ballooning modes and turbulence in tokamak plasmas

This study investigates the impact of various impurity species on peeling–ballooning (PB) modes and microturbulence in tokamak plasmas through the extension of traditional two-fluid and gyro-landau-fluid (GLF) models. By incorporating finite Larmor radius (FLR) effects, the analysis provides a comprehensive understanding of impurity-driven impact and its interaction with plasma turbulence. Depending on charge state and local plasma conditions, heavy impurities may exhibit gyro-radii larger than those of main ions, which are captured in the extended GLF model presented. Following the presentation of modified two-fluid equations incorporating impurity effects, we systematically analyze the distinctions between impurity and main ion dynamics and their resultant feedback mechanisms on plasma behavior. Derivation of the linear dispersion relation enables quantification of impurity-mediated modifications to: plasma vorticity, diamagnetic drift and gyroviscous effects, electron Hall physics, and FLR effects. BOUT++ – based linear simulations corroborate this formalism, demonstrating systematic stabilization of PB modes upon impurity seeding. And then operational implications for practical impurity control strategies in tokamak devices are proposed. The results underscore the necessity of impurity management to maintain stability and optimize plasma confinement, with specific focus on how FLR effects contribute to transport dynamics. This work paves the way for enhanced modeling and simulation efforts, supporting the development of strategies to control impurity-induced turbulence and improve overall reactor performance.

BOUT++ simulation↗

DOE OSTI2022

Impact of pedestal density gradient and collisionality on ELM dynamics

BOUT++ turbulence simulations are conducted to capture the underlying physics of small ELM characteristics achieved by increasing separatrix density via controlling strike points from vertical to horizontal divertor plates for three EAST discharges. BOUT++ linear simulations show that the most unstable modes change from high-n ideal ballooning modes to intermediate-n peeling–ballooning modes and eventually to peeling–ballooning stable plasmas in the pedestal. Nonlinear simulations show that the fluctuation is saturated at a high level for the lowest separatrix density. The ELM size decreases with increasing separatrix density, until the fraction of this energy lost during the ELM crash becomes less than 1% of the pedestal stored energy, leading to small ELMs. Simulations indicate that small ELMs can be triggered either by the marginally peeling–ballooning instability near the peak pressure gradient position inside the pedestal or by a local instability in the pedestal foot with a larger separatrix density gradient. The pedestal collisionality scan for type-I ELMs with steep pedestal density gradient shows that both linear growth rate and ELM size decrease with increasing collisionality. On the contrary, the pedestal collisionality and pedestal density width scan with a weak pedestal density gradient indicate small ELMs can either be triggered by a high-n ballooning mode or by a low-n peeling mode in a low collisionality region 0.04–0.1. We report the simulations indicate the weaker the linear unstable modes near marginal stability with small linear growth rate, the lower nonlinearly saturated fluctuation intensity and the smaller turbulence spreading from the linear unstable zone to stable zone in the nonlinear saturation phase, leading to small ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Critical role of current-driven instabilities for ELMs in NSTX

The impact of different extended-magnetohydrodynamic (MHD) contributions on the stability thresholds of peeling-ballooning modes in ELMing and ELM-free plasmas in the spherical tokamak NSTX is investigated with the initial value code M3D-C1. We show that ELMing discharges in NSTX are limited by resistive current-driven peeling modes, whereas non-ELMing wide-pedestal H-mode discharges are located near the ideal pressure-driven ballooning threshold. It is demonstrated that extended-MHD can lead to more reliable edge stability predictions than existing ideal-MHD models. Resistive peeling-ballooning modes are found to exist well before the ideal stability threshold is met, and kink-peeling modes exhibit considerable sensitivity to plasma resistivity. Other effects not considered in ideal-MHD models affect PB modes in NSTX in a weaker way. Gyroviscous stress appears stabilizing such that the stability boundary lies closer to the experimental point. Equilibrium rotation can suppress ideal core modes and thus isolate edge modes. These results are important for the development of a predictive pedestal model for low-aspect ratio tokamaks.

pedestal↗

DOE OSTI2024

MHD stability analysis against pressure and current-driven modes in the SMall Aspect Ratio Tokamak

Linear magnetohydrodynamic (MHD) simulations for the SMall Aspect Ratio Tokamak (SMART) have been carried out for the first time, for both positive (PT) and negative triangularity (NT) shaped plasmas using the MARS-F code. The MHD stability of projected SMART plasmas against internal kinks, infernal modes and edge peeling-ballooning modes have been analyzed for a wide range of realistic equilibria. A stabilization of internal kinks and infernal modes is observed when increasing the safety factor profile and reducing plasma beta. PT shaped plasmas are more stable against both internal kinks and infernal modes than their counterpart NT shaped plasmas. Toroidal flows have little impact on the MHD stability of the internal kinks, but they have a strong stabilizing effect on infernal modes, which can be further mitigated in NT shaped plasmas. The MHD stability of peeling-ballooning modes is reduced in NT shaped plasmas, as observed in conventional tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Pedestal stability analysis on MAST in preparation for MAST-U

In preparation for the upcoming MAST-U campaign, pedestal stability of spherical tokamaks is revisited by investigating standard H-mode discharges on MAST. As a step beyond previous studies, both ion and electron profiles are used for obtaining equilibria and a diverse set of pedestals is evaluated. Stability analysis with the ELITE and CGYRO codes shows that MAST pedestals are constrained by kinetic ballooning modes and medium toroidal mode number peeling-ballooning modes, with most unstable modes ranging from n = 25 to n = 45. In discharges with a steep q profile at the edge a larger number of poloidal harmonics is excited for each toroidal mode. In conclusion, a comparison with discharges on DIII-D with matched shape and similar non-dimensional parameters indicates that the increased shear at lower aspect ratio stabilizes low n peeling modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Nonlinear MHD modeling of n = 1 RMP -induced pedestal transport and mode coupling effects on ELM suppression in KSTAR

We report fully suppressing edge-localized modes (ELMs), e.g., with resonant magnetic perturbations (RMPs), is essential to reach and sustain high-performance steady-state H-mode plasmas because large ELMs can significantly reduce the lifetime of divertor components in future tokamak reactors. RMP-driven ELM suppression in KSTAR has been modeled by coupling the neoclassical transport code PENTRC to the nonlinear 3D MHD code JOREK. We have found that the radial transport from the combined effects of the kink-peeling, tearing response, and neoclassical toroidal viscosity can explain the pedestal degradation observed in experiments. In addition, it has been found that the RMP response can increase the inter-ELM heat flux on the lower outer divertor by redistributing the heat transport between the divertor plates. In addition to the degraded pedestal, ELM suppression is also attributable to the RMP-induced mode interactions. While the linear stability of peeling-ballooning mode (PBMs) improves owing to the degraded pedestal, the PBM and RMP interaction increases the spectral transfer between edge harmonics, preventing catastrophic growth and the crash of unstable modes. Here, it turns out that the magnetic islands near the pedestal top can play a vital role in mediating the mode interactions.

mode coupling↗

DOE OSTI2023

Nonlinear second order electromagnetic gyrokinetic theory for a tokamak plasma

The steep plasma pressure gradient that forms at the edge of the high confinement, H-mode regime of tokamak operation provides free energy to drive electromagnetic micro-instabilities that are widely believed to influence the transport processes in this so-called pedestal region. This high pressure gradient also provides a high current density (bootstrap current), known to influence ballooning mode stability and to be important for driving kink modes in the ideal magneto-hydrodynamic plasma model (so-called peeling-ballooning modes). Furthermore, efficient, steady state future tokamak power plants must operate with a large bootstrap current in the core and especially concerning spherical tokamaks, confinement will be influenced by electromagnetic turbulence. To accommodate these important situations, conventional electromagnetic gyrokinetic theory is extended to incorporate neoclassical effects in the equilibrium drives, allowing $B_{\vartheta} ~ B_0$ (B0 is the confining magnetic field, and $B_{\vartheta}$ is its poloidal component). This provides a global gyrokinetic model that self-consistently captures the consequences of large bootstrap current fractions on the equilibrium distribution functions.

kink modes↗

DOE OSTI2022

Different roles of elongation and triangularity in the H-mode pedestal height enhancement through strong plasma shaping

While the peeling–ballooning mode is known to be stabilized by strong plasma shaping in a high-confinement mode (H-mode), a significant difference is found in the roles of elongation and triangularity. The elongation mainly stabilizes the high-n ballooning branch, while the triangularity of the low-n peeling one with some destabilization of the high-n ballooning one, where n is the toroidal mode number. This difference is conjectured mainly due to the magnetic shear effect which is substantially enhanced in a high triangularity regime, giving a destabilization on the ballooning branch, while a stabilization on the peeling one. As shown here, an important consequence of these different roles of the two shape parameters is that they should work together to stabilize simultaneously both branches of the PBM, implying the necessity of synergy between them for the large pedestal height enhancement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Characteristics of grassy ELMs and their impact on the divertor heat flux width

We report BOUT++ turbulence simulations are conducted for a 60 s steady-state long pulse high β p EAST grassy ELM discharge. BOUT++ linear simulations show that the unstable mode spectrum covers a range of toroidal mode numbers from low-n (n = 10–15) peeling–ballooning modes (P–B) to high-n (n = 40–80) drift-Alfvén instabilities. Nonlinear simulations show that the ELM crash is triggered by low-n peeling modes and fluctuation is generated at the peak pressure gradient position and radially spread outward into the scrape-off-layer, even though the drift-Alfvén instabilities dominate the linear growth phase. However, drift-Alfvén turbulence delays the onset of the grassy ELM and enhances the energy loss with the fluctuation extending to pedestal top region. Simulations further show that if the peeling drive is removed, the fluctuation amplitude drops by an order of magnitude and the ELM crashes disappear. The divertor heat flux width is ~2 times larger than the estimates based on the HD model and the Eich's ITPA multi-tokamak scaling (or empirical Eich scaling) due to the strong radial turbulence transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Comparison of MHD stability properties between QH-mode and ELMy H-mode plasmas by considering plasma rotation and ion diamagnetic drift effects

Magnetohydrodynamic (MHD) stability at tokamak edge pedestal in a quiescent H-mode (QH-mode) and type-I ELMy H-mode plasmas in DIII-D experiment was analyzed by considering plasma rotation and ion diamagnetic drift effects. QH-mode plasma is marginally stable to kink/peeling mode (K/PM), but ELMy H-mode one is almost unstable to peeling-ballooning mode (PBM). It was identified that there are three physics features responsible for the difference in the MHD stability properties between QH-mode plasma and ELMy H-mode one. These are the distance of pedestal foot from the last closed flux surface (LCFS), the amount of the ion diamagnetic drift frequency at pedestal, and impact of coupled rotation and ion diamagnetic drift effects. These features were confirmed through the numerical experiments that the stability properties of the QH-mode plasma can be changed to that of the ELMy H-mode one by shifting the plasma profiles inward in the radial direction and halving the ion diamagnetic drift frequency. The reasons of the change in the stability properties are thought as that K/PM is stabilized due to the inward shift of the bootstrap current profile, and PBM is destabilized due to the reduction of the coupled rotation and ion diamagnetic drift stabilizing effect. Importance of these features was validated through numerical experiments with experimental data of other QH-mode plasmas in DIII-D. All the results show that MHD stability properties of QH-mode plasma can be obtained in case that pedestal foot is close to LCFS, ion diamagnetic drift frequency is large due to high ion temperature, and strong rotation shear exists near pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

X-point effects on the ideal MHD modes in tokamaks in the description of dual-poloidal-region safety factor

The flux coordinates with dual-region safety factor (q) in the poloidal direction are developed in this work. The X-point effects on the ideal MHD modes in tokamaks are then analyzed using this coordinate system. Since the X-point effects mainly affect the edge region, the modes localized at the tokamak edge are particularly examined. Two types of modes are studied. The first is related to the conventional peeling or peeling-ballooning modes. The mode existence aligned with the local magnetic field in the poloidally core region as observed experimentally is confirmed. The X points are shown to contribute to a stabilizing effect for the conventionally treated modes with the surface-averaged q and with the tokamak edge portion truncated. The other is the axisymmetric modes localized in the vicinity of X points, which can affect the cross-field-line transport near the X points. The existence of axisymmetric modes points to the possibility of applying a toroidally axisymmetric resonant magnetic perturbation (RMP) in the X-point area for mitigating the edge localized modes, which can be an alternative to the current RMP design. The dual q description also has important implications for the existing non-axisymmetric RMP concept. It helps to understand why the RMP suppression of edge localized modes is difficult to achieve in the double-null tokamak configurations and points to the possibility of further improving the current RMP concept by considering the alignment to the local q.

Fourier analysis↗

DOE OSTI2020

H-mode pedestal improvements with neon injection in DIII-D

Improvements in the H-mode pedestal pressure and global energy confinement were determined on DIII-D with both 3 and 6 MW of neutral beam (NBI) heating with neon injection. The pedestal pressure increased, primarily associated with a significant increase in the pedestal density with a relatively smaller decrease in the pedestal temperature. At the same neon injection rate, a rapid decrease in edge localized mode (ELM) frequency and ~35% increase in plasma stored energy were observed in the 3 MW discharge, while in the 6 MW discharge, the ELM frequency was gradually decreased and the plasma stored energy gradually increased by up to 20%. The measured pedestal widths from both 3 and 6 MW discharges matched the EPED1.0 width scaling, ${\Delta _{{\psi _N}}} = 0.089\beta _{p,ped}^{1/2}$ ,data-, within 20% deviations before and after neon injection. The peeling-ballooning mode (PBM) stability analysis showed that after neon injection, the ballooning stability boundary was increased, while the peeling stability boundary only dropped once the neon level was relatively high. The ballooning stability boundary increase with neon injection is consistent with increase of the ion diamagnetic stabilization (i.e. including the contribution from carbon, neon as well as deuterium ions) which is used in the criterion to determine when the calculated linear growth rates of PBM are unstable. Consequently, the new ELITE computed PBM boundaries were more consistent with the operating points than when only deuterium was included. This in turn suggests that the impurity may help to improve pedestal pressure through affecting ion diamagnetic frequency in the pedestal region, which could positively lead to the ballooning stability boundary extension.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Elements of H-mode pedestal structure

Abstract This paper reviews current understanding of key physics elements that control the H-mode pedestal structure, which exists at the boundary of magnetically confined plasmas. The structure of interest is the width, height and gradient of temperature, density and pressure profiles in the pedestal. Emphasis is placed on understanding obtained from combined experimental, theoretical and simulation work and on results observed on multiple machines. Pedestal profiles are determined by the self-consistent interaction of sources, transport and magnetohydrodynamic limits. The heat source is primarily from heat deposited in the core and flowing to the pedestal. This source is computed from modeling of experimental data and is generally well understood. Neutrals at the periphery of the plasma provide the dominant particle source in current machines. This source has a complex spatial structure, is very difficult to measure and is poorly understood. For typical H-mode operation, the achievable pedestal pressure is limited by repetitive, transient magnetohydrodynamic instabilities. First principles models of peeling–ballooning modes are generally able to explain the observed limits. In some regimes, instability occurs below the predicted limits and these remain unexplained. Several mechanisms have been identified as plausible sources of heat transport. These include neoclassical processes for ion heat transport and several turbulent processes, driven by the steep pedestal gradients, as sources of electron and ion heat transport. Reduced models have successfully predicted the pedestal or density at the pedestal top. Firming up understanding of heat and particle transport remains a primary challenge for developing more complete predictive pedestal models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Multi-scale Interaction Mechanism for Edge-Localized-Mode Suppression in the Tokamak Edge

A central challenge in fusion energy is reconciling the high-confinement mode required for reactor performance with the intense intermittent relaxation events it produces, known as edge-localized modes. These instabilities arise in the steep pressure pedestal at the plasma edge when magnetohydrodynamic thresholds are crossed, inflicting damaging heat loads on reactor components. Here, we show that multiscale interactions between microscopic turbulence and macroscopic magnetohydrodynamic modes provide encouraging prospects for self-organized edge-localized modes regulation. Using direct quantitative measurements of multiscale modes, eddy dynamics, and turbulent flux, we show that small-scale electron drift wave turbulence actively scatters the large-scale peeling-ballooning modes. This scattering decorrelates the pressure and velocity fields of the instability, so arresting its growth. Our modeling and theoretical analysis confirm this suppression mechanism is effective even when conventional linear stability thresholds are exceeded. This work establishes a nonlinear principle for edge-localized modes stability, revealing how ambient micro-turbulence can be leveraged to maintain a macro-stable, high-performance pedestal for future fusion reactors.

experimental nuclear physics↗

DOE OSTI2025

Effects beyond ideal MHD on stability of wide and enhanced pedestal regimes in NSTX

Stability of edge-localized modes (ELMs) in spherical tokamaks is explored using the extended MHD model. Linear NIMROD simulations have been performed for three NSTX discharges 132543, 132588, and 141133, to investigate the role of resistivity, diffusivity, and shear flows on the onset of ELMs. The first discharge represents the wide pedestal regime and the later two discharges represent ELM-free enhanced pedestal H-mode. We first present the effect of toroidal rotation shear and find a flow shear destabilizing effect in these NSTX discharges. Simulations are also extended to include the two-fluid and ion gyroviscosity effects. Simulations show that the flow shear can shift the mode spectrum and alter the critical condition of ELM onset. We also uncover that ELM onset prediction in spherical tokamaks requires effects beyond MHD, in particular gyroviscosity and diamagnetic terms could stabilize Peeling-Ballooning modes consistent with the experimental observation of ELM-free regimes in NSTX. The findings give new insight into the nature of the interplay between resistivity, flows, and diamagnetic stabilization in ELM suppression and have potential applications to ELM control schemes in NSTX-U and next-generation spherical tokamaks. This study identifies the essential physical effects that must be included in future predictive and validation simulations.

NSTX↗