Search NASA⌕ Search

SEARCH · Search NASA

Results for “MHD 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Ideal MHD stability and characteristics of edge localized modes on CFETR

Investigation on the equilibrium operation regime, its ideal magnetohydrodynamics (MHD) stability and edge localized modes (ELM) characteristics is performed for the China Fusion Engineering Test Reactor (CFETR). The CFETR operation regime study starts with a baseline scenario (R = 5.7 m, B T = 5 T) derived from multi-code integrated modeling, with key parameters ${{\beta }_{{\rm N}}},{{\beta }_{{\rm T}}},{{\beta }_{{\rm p}}}$ varied to build a systematic database. These parameters, under profile and pedestal constraints, provide the foundation for the engineering design. The long wavelength low-n global ideal MHD stability of the CFETR baseline scenario, including the wall stabilization effect, is evaluated by GATO. It is found that the low-n core modes are stable with a wall at r/a = 1.2. An investigation of intermediate wavelength ideal MHD modes (peeling ballooning modes) is also carried out by multi-code benchmarking, including GATO, ELITE, BOUT++ and NIMROD. A good agreement is achieved in predicting edge-localized instabilities. Nonlinear behavior of ELMs for the baseline scenario is simulated using BOUT++. A mix of grassy and type I ELMs is identified. When the size and magnetic field of CFETR are increased (R = 6.6 m, B T = 6 T), collisionality correspondingly increases and the instability is expected to shift to grassy ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

MHD stability of negative triangularity DIII-D plasmas

Abstract Negative triangularity (NT) experiments in DIII-D point to an emergent reactor scenario free of sawteeth, endowed with benign, nondisruptive n = 2 tearing modes, which experience q min ⩾ 1 similar to the positive triangularity hybrid scenario. Plasmas exhibiting this behavior attain β N > 3 , high enough to reconsider long held views of NT stability. Ideal MHD and tearing stability analysis of well-diagnosed equilibrium reconstructions of experimental hybrid-like plasmas predict that among shape parameters, MHD stability limits are only sensitive to average triangularity. Operation is predicted to be possible at β N relevant to upcoming tokamaks and commercial NT reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Special Issue on the 20th Workshop on MHD Stability Control

The 20th workshop on magnetohydrodynamic (MHD) stability control took place November 22–24, 2015, in Princeton Plasma Physics Laboratory (PPPL), following the American Physical Society—Division of Plasma Physics annual meeting on November 16–20 in Savannah, GA. The purpose of this workshop is to stimulate in depth discussion and motivate future research in the areas of MHD stability physics and control of magnetically confined plasmas. Furthermore, the workshop was organized jointly by Auburn University, Columbia University, General Atomics, Princeton Plasma Physics Laboratory, University of Wisconsin-Madison, and the Los Alamos National Laboratory, and under the auspices of the US/Japan Collaboration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Identification of a window for quiescent H-mode operation in MHD stability diagram of DIII-D plasmas

A window for quiescent H-mode (QH-mode) operation in DIII-D plasmas in the diagram of pedestal magnetohydrodynamic (MHD) stability was identified through linear MHD stability analysis, accounting for the effects of plasma rotation and ion diamagnetic drift. The operation window lies between the stability boundaries of the kink/peeling mode (K/PM) identified with and without the inclusion of plasma rotation effects. In this region, the mode remains unstable unless rotation effects are considered alongside the ion diamagnetic drift, which is consistently included in the analysis. The stabilization of the MHD mode, facilitated by the coupled effects of plasma rotation and ion diamagnetic drift, plays a crucial role in widening the window, enabling the attainment of the QH-mode state. Furthermore, the suppression of edge localized modes (ELM) can be achieved by controlling the pedestal structure to ensure the plasma state remains within the operation window. The location of the operation window in the stability diagram depends on the K/PM stability properties. Therefore, optimizing conditions for QH-mode requires adjustments based on changes in stability characteristics. A pressure pedestal and its associated bootstrap current density near the last closed flux surface are advantageous for situating the plasma state within the window. However, excessive current density can trigger ELMs. This trend was confirmed through comparisons of MHD stability diagrams between QH-mode and ELMy H-mode plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 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 OSTI2020

Special issue on the 24th workshop on MHD stability control

We report the 24th workshop on magnetohydrodynamic (MHD) stability control was held on the Columbia University campus 28–30th October, 2019. The purpose of this annual workshop is to discuss cutting edge results and form future plans for active control of MHD instabilities in magnetic confinement fusion reactor concepts. The program was organized by an international committee of representatives from Auburn University (USA), Columbia University (USA), Massachusetts Institute of Technology (USA), Princeton Plasma Physics Laboratory (USA), Princeton University (USA), Tottori University (Japan), Kyoto Institute of Technology (Japan), National Institute of Radiological Sciences (Japan), National Institute for Fusion Science (Japan), Institutes of Physical Science—Chinese Academy of Sciences (People's Republic of China), and Consiglio Nazionale delle Ricerche (Italian Republic). Dr Gerald Navratil and Columbia University, in particular, were responsible for organizing and hosting the meeting on their beautiful Morningside Heights campus.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

MHD stability of spherical tokamak equilibria with non-monotonic q-profiles

We use the 3D magnetohydrodynamic (MHD) code M3D-C1 [Jardin et al., Comput. Sci. Discovery 5, 014002 (2012)] to examine the MHD stability and subsequent evolution of NSTX shot 129169. This discharge had a period with a non-monotonic safety factor profile, q (reversed shear), which was terminated by a MHD event that abruptly lowered the central safety factor, q0, and greatly reduced the peakedness of the pressure profile. We show that the equilibrium just before the MHD event occurred was linearly unstable to many pressure-driven infernal modes. Modes with toroidal mode number n≥3 all had rational surfaces very close to the minimum value of q. However, a non-resonant pressure-driven (1, 1) mode was also present, and this dominated the nonlinear evolution. The final state in the simulation, after the MHD activity subsided, had a reduced and flattened pressure profile and a nearly monotonic q-profile, in qualitative agreement with experimental results. The initial state was also unstable to the resistive interchange criteria in the reversed-shear region, but the final state was stable everywhere. The “double tearing mode” (DTM) does not appear to play a role in the MHD activity of this discharge. In Appendix A, we show that in a torus, the DTM is strongly stabilized by pressure, but it is destabilized in cylindrical geometry (which has been the most extensively analyzed in the literature).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Editorial for the 28th Workshop on Magnetohydrodynamic (MHD) Stability Control

This special issue contains four papers which were presented at the 28th Workshop on Magnetohydrodymanic (MHD) Stability Control. The workshop was held at Auburn University and on Zoom from 3–5 October 2024. Furthermore, this annual workshop covers the active and passive control of MHD instabilities such as kink, tearing and edge-localized modes in tokamaks, reversed field pinches, spherical tori or other magnetic confinement devices.

Yang, J. [Princeton Plasma Physics Laboratory (PPP↗
DOE OSTI2019

Analysis of the MHD stability and energetic particles effects on EIC events in LHD plasma using a Landau-closure model

The aim of this study is to perform a theoretical analysis of the magnetohydrodynamic (MHD) stability and energetic particle effects on a LHD equilibria, calculated during a discharge where energetic-ion-driven resistive interchange mode (EIC) events were triggered. Here, we use 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 moments for the energetic particles species, including the effect of the acoustic modes, multiple energetic particles (EP) species, helical couplings and helically trapped EP. We add the Landau damping and resonant destabilization effects using a closure relation. The simulations suggest that the helically trapped EP driven by the perpendicular neutral beam injector (NBI) further destabilizes the 1/1 MHD-like mode located at the plasma periphery ( r / a = 0.88). If the β of the EP driven by the perpendicular NBI is larger than 0.0025 a 1/1 EIC with a frequency around 3 kHz is destabilized. If the effect of the passing EP driven by the tangential NBI is included on the model, any enhancement of the injection intensity of the tangential NBI below β leads to a decrease of the instability growth rate. The simulations indicate that the perpendicular NBI EP is the main driver of the EIC events, as it was observed in the experiment. If the effect of the helical couplings are added in the model, an 11/13 EIC is destabilized with a frequency around 9 kHz, inward shifted ( r / a = 0.81) compared to the 1/1 EIC. Thus, one possible explanation for the EIC frequency chirping down from 9 to 3 kHz is a transition between the 11/13 to the 1/1 EIC due to a weakening of the destabilizing effect of the high n modes, caused by a decrease of the EP drive due to a loss of helically trapped EP or a change in the EP distribution function after the EIC burst. The experimental data during the EIC bursting phase shows a complex mode structure and an inward shift of the instability, although no direct evidence of the proposed transition has been observed yet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Impact of ideal MHD stability limits on high-beta hybrid operation

Here, the hybrid scenario is a candidate for stationary high-fusion gain tokamak operation in ITER and DEMO. To obtain such performance, the energy confinement and the normalized pressure ${{\beta}_{N}}$ must be maximized, which requires operating near or above ideal MHD no-wall limits. New experimental findings show how these limits can affect hybrid operation. Even if hybrids are mainly limited by tearing modes, proximity to the no-wall limit leads to 3D field amplification that affects plasma profiles, e.g. rotation braking is observed in ASDEX Upgrade throughout the plasma and peaks in the core. As a result, even the small ASDEX Upgrade error fields are amplified and their effects become visible. To quantify such effects, ASDEX Upgrade measured the response to 3D fields applied by $8\times 2$ non-axisymmetric coils as ${{\beta}_{N}}$ approaches the no-wall limit. The full n = 1 response profile and poloidal structure were measured by a suite of diagnostics and compared with linear MHD simulations, revealing a characteristic feature of hybrids: the n = 1 response is due to a global, marginally-stable n = 1 kink characterized by a large m = 1, n = 1 core harmonic due to q min being just above 1. A helical core distortion of a few cm forms and affects various core quantities, including plasma rotation, electron and ion temperature, and intrinsic W density. In similar experiments, DIII-D also measured the effect of this helical core on the internal current profile, providing information useful to understanding of the physics of magnetic flux pumping, i.e. anomalous current redistribution by MHD modes that keeps ${{q}_{\text{min}}}>1$ . Thanks to flux pumping, a broad current profile is maintained in DIII-D even with large on-axis current drive, enabling fully non-inductive operation at high ${{\beta}_{N}}$ up to 3.5–4.

3D equilibrium↗
DOE OSTI2020

Analysis of MHD stability and active mode control on KSTAR for high confinement, disruption-free plasma

Long-pulse plasma operation at high normalized beta, $β_N$, above the $\textit{n}$= 1 ideal MHD no-wall stability limit in KSTAR is presently limited by tearing instabilities rather than resistive wall modes. H-mode plasma operation during the recent KSTAR device campaign produced discharges having strong $\textit{m/n}$= 2/1 tearing instabilities at $β_N$ lower than the ideal MHD no-wall beta limit. The unstable tearing mode consequently reduced plasma confinement and toroidal plasma rotation significantly. We report the experiment confirmed that an extended duration of electron cyclotron heating (ECH) at the initial phase of the discharge plays a critical role in mode destabilization. To study destabilizing mechanisms that affect the mode growth, the stability of the observed tearing modes from plasmas with significantly different $β_N$ is computed by using the resistive DCON code and the M3D-C 1 code. Equilibrium reconstructions that include constraints from internal profile diagnostics, and computed fast particle pressure are used as input for reliable computation of stability. The classical tearing stability index, Δ', from resistive DCON is computed to be unstable when the island is fully saturated with large amplitude, while the unstable mode is not computed by M3D-C 1 . The modified Rutherford equation (MRE) describing the evolution of the island width has been constructed for KSTAR plasmas by using plasma parameters computed by the TRANSP code. The MRE model estimates a saturated island width corresponding to ~10% of the plasma minor radius for equilibrium at high $β_N$ 3 having a stable Δ' from the resistive DCON. In preparation for long-pulse plasma operation at higher beta utilizing increased plasma heating power, a resistive wall mode (RWM) active feedback control algorithm has been completed and enabled on KSTAR. To accurately determine the $\textit{n}$= 1component produced by RWMs, an algorithm has been developed that includes magnetic sensor compensation of the prompt applied field and the field from the induced current on the passive conductors. Use of multiple toroidal sensor arrays is enabled by modifying the sensor toroidal angles assumed in mode decomposition to include the effect of varied mode helicities in the outboard region where the mode measurement is made. This analysis on stability, transport, and control provides the required foundation for disruption prediction and avoidance research on KSTAR.

resistive wall mode (RWM)↗
DOE OSTI2021

Mode spectrum characteristics and onset of the low-shear MHD stability regime

Equilibria with extended regions of weak magnetic shear, including some tokamak scenarios and stellarators, can be susceptible to pressure-driven internal MHD instabilities even though there is no mode rational surface in the plasma. Nonresonant modes, in particular, can have properties that are unattractive for confinement, including displacing substantial volumes of the plasma and leading to more efficient pressure gradient flattening in the nonlinear regime. The onset and linear properties of the low shear stability regime are examined using the initial-value, extended-MHD code M3D-C. For monotonic q-profiles, we demonstrate a clear correlation between the convergents associated with the continued fraction representation of q 0 and the spectrum of unstable modes. Nonresonant modes are observed to be destabilized preferentially to any other resonant instability with the same toroidal mode number when n > 1. Using the observed connection between the spectrum associated with q 0 and the overall equilibrium stability properties, we suggest a technique for reducing the uncertainty on both q 0 and magnetic shear in the core region, obtained either via measurement or through the analysis and reconstruction of experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Development and Application of Predictive Tools for MHD Stability Limits in Tokamaks

This is a project to develop and apply analytic and computational tools to answer physics questions relevant to the onset of non-ideal magnetohydrodynamic (MHD) instabilities in toroidal magnetic confinement plasmas. The focused goal of the research is to develop predictive tools for these instabilities, including an inner layer solution algorithm, a resistive wall with control coils, and energetic particle effects. The production phase compares studies of instabilities in such systems using analytic techniques, PEST- III and NIMROD. Two important physics puzzles are targeted as guiding thrusts for the analyses. The first is to form an accurate description of the physics determining whether the resistive wall mode or a tearing mode will appear first as β is increased at low rotation and low error fields in DIII-D. The second is to understand the physical mechanism behind recent NIMROD results indicating strong damping and stabilization from energetic particle effects on linear resistive modes. The work seeks to develop a highly relevant predictive tool for ITER, advance the theoretical description of this physics in general, and analyze these instabilities in experiments such as ASDEX Upgrade, DIII-D, JET, JT-60U and NTSX. The awardee on this grant is the University of Tulsa. The research efforts are supervised principally by Dr. Brennan. Support is included for two graduate students, and a strong collaboration with Dr. John M. Finn of LANL. The work includes several ongoing collaborations with General Atomics, PPPL, and the NIMROD team, among others.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Study of MHD Stability, Active Mode Control, and Disruption Avoidance in NSTX-U

Learning how to achieve sustained operation at high plasma pressure is an important goal of the National Spherical Torus Upgrade (NSTX-U) program and is critically important to any future fusion power plant. This research addresses high priority stability and control issues outlined by the NSTX-U Program. Reduced models of the stabilization physics learned continue to be created and tested against actual plasma operation in NSTX, NSTX-U, and other tokamak devices. This research is forming physics-based models for disruption forecasting that aims to unify results across different fusion machines for more confident extrapolation of the techniques studied to future devices including ITER. The Disruption Event Characterization and Forecasting (DECAF) code was expanded to analyze data from the TCV tokamak, including real-time MHD analysis signals. The presence of a mode and likelihood of its m/n numbers can be tracked, as well as the mode frequency, although a limitation at low frequency does not yet allow tracking modes to the locking point in TCV. Initial analyses of NSTX-U operation were conducted using the DECAF code. The timing of sawtooth oscillations was investigated taking advantage of the more extensive L-mode operation of the device in its first run campaign. An initial investigation of thermal collapses in NSTX-U was also conducted with specific comparison to the Greenwald density limit. A local island power balance model to explain the density limit in tokamaks has been tested with NSTX data and is being evaluated as a potential disruption forecaster in the DECAF code, presently showing similarities to the global Greenwald density evolution and limit. A module was written for the DECAF code that generalizes magnetic spectrogram analysis of MHD modes. This new portable code, written in Python, discriminates the toroidal mode number of instabilities and can track each mode as a DECAF event, monitoring the evolution of the mode through bifurcation and locking points. This is a critical new capability to allow disruption event characterization for general tokamak databases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗