Magnetohydrodynamic Stability Research Articles

Reduced electrostatic turbulence in the quasi-isodynamic stellarator configuration CIEMAT-QI4

Abstract CIEMAT-QI4 is a quasi-isodynamic stellarator configuration that simultaneously features very good fast-ion confinement in a broad range of β values, low neoclassical transport and bootstrap current, and ideal magnetohydrodynamic stability up to β=5%. In this paper it is shown that CIEMAT-QI4 also exhibits reduced turbulent transport. This is demonstrated through nonlinear electrostatic simulations with the gyrokinetic code stella, including kinetic ions and electrons. The relation between reduced turbulent transport and the fact that CIEMAT-QI4 very approximately satisfies the so-called maximum-J property is discussed.

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

Abstract 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.

The Effect of External Magnetic Field on Electron Scale Kelvin–Helmholtz Instability

We use particle-in-cell, fully electromagnetic, plasma kinetic simulation to study the effect of external magnetic field on electron scale Kelvin–Helmholtz instability (ESKHI). The results are applicable to collisionless plasmas when, e.g., solar wind interacts with planetary magnetospheres or a magnetic field is generated in AGN jets. We find that as in the case of magnetohydrodynamic (MHD) KHI, in the kinetic regime, the presence of an external magnetic field reduces the growth rate of the instability. In the MHD case, there is a known threshold magnetic field for KHI stabilization, while for ESKHI this is to be analytically determined. Without a kinetic analytical expression, we use several numerical simulation runs to establish an empirical dependence of ESKHI growth rate, Γ(B 0)ω pe, on the strength of the applied external magnetic field. We find the best fit is hyperbolic, Γ(B0)ωpe=Γ0ωpe/(A+BB¯0) , where Γ0 is the ESKHI growth rate without an external magnetic field and B¯0=B0/BMHD is the ratio of external and two-fluid MHD stability threshold magnetic field, derived here. An analytical theory to back up this growth rate dependence on the external magnetic field is needed. The results suggest that in astrophysical settings where a strong magnetic field pre-exists, the generation of an additional magnetic field by the ESKHI is suppressed, which implies that nature provides a “safety valve”—natural protection not to “over-generate” magnetic field by the ESKHI mechanism. Remarkably, we find that our two-fluid MHD threshold magnetic field is the same (up to a factor γ0 ) as the DC saturation magnetic field, previously predicted by fully kinetic theory.

Early internal detection of magnetic tearing and implications for tokamak magnetohydrodynamic stability

Internal fluctuation measurements with Faraday-effect polarimetry in the DIII-D tokamak reveal the onset of a tearing mode with toroidal mode number n = 3 well before it is detected by the sensing coils external to the plasma. This mode appears before the n = 2, 1 modes and is first detected with internal measurements at a lower value of the ideal-wall kink beta limit than is indicated at the time of first detection by the sensing coils. When the mode is first detected, the linear resistive stability parameter, Δ′, indicates marginal stability and continues to do so until later when the mode amplitude begins increasing linearly with time—together suggesting a neoclassical origin for this mode.

Effect of symmetry-breaking on the MHD edge stability limit of tokamak plasmas

We demonstrate for the first time that symmetry-breaking in tokamak plasmas reduces the window of stable edge operational space, leading to a lower achievable edge pressure gradient. Predictive simulations with the linear magnetohydrodynamic stability code CASTOR3D show a reduction of the critical pressure gradient by up to 30%, in agreement with experimental observations. The analysis has been extended to experimental plasmas considering, for the first time, ion diamagnetic drift effects in realistic non-axisymmetric (3D) tokamak geometry. The 3D geometry in plasmas with edge localised modes (ELMs) increases the growth rate of edge instabilities compared to the axisymmetric case. Further reducing the edge pressure gradient eliminates the instabilities corresponding to an experiment with suppressed ELMs, reproducing the empirically observed threshold for ELM suppression. Our findings highlight the importance of predictions of the ELM occurrence in full 3D geometry for future tokamak devices with intentionally and unintentionally broken axisymmetry.

Development of a neural network model for peeling–ballooning stability analysis in the KSTAR tokamak pedestals

The neural network model, MISHKA-NN is developed to mitigate the computational burden associated with the linear ideal magnetohydrodynamic (MHD) stability analysis of the pedestal based on the peeling–ballooning (P–B) model. By utilizing both 1D plasma profiles (current density, pressure gradient, and safety factor) and 0D parameters (plasma geometry, total current, and toroidal mode number), the model predicts linear growth rate of edge-localized ideal MHD instability in a given equilibrium state. By enabling the prediction of each instability within a second, the model reduces the time required for plotting a pedestal P–B stability diagram ( j−α diagram) from approximately 100 CPU hours to a few CPU minutes. Notably, even with the utilization of parametric pressure and current profiles and plasma boundary shapes for the training dataset, the model shows a satisfactory level of performance in benchmarking the j−α diagram for the reconstructed equilibrium from a KSTAR tokamak experiment. We anticipate the model to serve as a versatile alternative to 2D linear MHD stability codes, alleviating numerical costs.

Reduction of electrostatic turbulence in a quasi-helically symmetric stellarator via critical gradient optimization

We present a stellarator configuration optimized for a large threshold (‘critical gradient’) for the onset of the ion temperature gradient (ITG) driven mode, which achieves the largest critical gradient we have seen in any stellarator. Above this threshold, gyrokinetic simulations show that the configuration has low turbulence levels over an experimentally relevant range of the drive strength. The applied optimization seeks to maximize the drift curvature, leading to enhanced local-shear stabilization of toroidal ITG modes, and the associated turbulence. These benefits are combined with excellent quasi-symmetry, yielding low neoclassical transport and vanishingly small alpha particle losses. Analysis of the resulting configuration suggests a trade-off between magnetohydrodynamic (MHD) and ITG stability, as the new configuration possesses a vacuum magnetic hill.

Deep learning approaches to recover the plasma current density profile from the safety factor based on Grad–Shafranov solutions across multiple tokamaks

Many magnetohydrodynamic stability analyses require generation of a set of equilibria with a fixed safety factor q-profile while varying other plasma parameters. A neural network (NN)-based approach is investigated that facilitates such a process. Both multilayer perceptron (MLP)-based NN and convolutional neural network (CNN) models are trained to map the q-profile to the plasma current density J-profile, and vice versa, while satisfying the Grad–Shafranov radial force balance constraint. When the initial target models are trained, using a database of semi-analytically constructed numerical equilibria, an initial CNN with one convolutional layer is found to perform better than an initial MLP model. In particular, a trained initial CNN model can also predict the q- or J-profile for experimental tokamak equilibria. The performance of both initial target models is further improved by fine-tuning the training database, i.e. by adding realistic experimental equilibria with Gaussian noise. The fine-tuned target models, referred to as fine-tuned MLP and fine-tuned CNN, well reproduce the target q- or J-profile across multiple tokamak devices. As an important application, these NN-based equilibrium profile convertors can be utilized to provide a good initial guess for iterative equilibrium solvers, where the desired input quantity is the safety factor instead of the plasma current density.

The maximum-J property in quasi-isodynamic stellarators

Some stellarators tend to benefit from favourable average magnetic curvature for trapped particles when the plasma pressure is sufficiently high. This so-called maximum- $J$ property has several positive implications, such as good fast-particle confinement, magnetohydrodynamic stability and suppression of certain trapped-particle instabilities. This property cannot be attained in quasisymmetric stellarators, in which deeply trapped particles experience average bad curvature and therefore precess in the diamagnetic direction close to the magnetic axis. However, quasi-isodynamic stellarators offer greater flexibility and allow the average curvature to be favourable and the precession to be reversed. We find that it is possible to design such stellarators so that the maximum- $J$ condition is satisfied for the great majority of all particles, even when the plasma pressure vanishes. The qualitative properties of such a stellarator field can be derived analytically by examining the most deeply and the most shallowly trapped particles, although some small fraction of the latter will inevitably not behave as desired. However, through numerical optimisation, we construct a vacuum field in which 99.6 % of all trapped particles satisfy the maximum- $J$ condition.

Progress of open systems at Budker Institute of Nuclear Physics

This paper is based on a report at the 2nd International Fusion Plasma Conference & 13th International Conference on Open Magnetic Systems for Plasma Confinement (iFPC & OS 2023), August 21–25, 2023, Busan, Korea and provides a brief overview of the status of work at the Budker Institute on the study of hot plasma confinement in open-type magnetic traps with a linear axisymmetric configuration. The main attention is paid to key problems: magnetohydrodynamics (MHD) stability in regimes with extremely high relative pressure, longitudinal electronic thermal conductivity, stability with respect to the development of kinetic modes and transverse transport. This paper provides an overview of the methods we are developing to address these problems, the experimental and theoretical results achieved and plans for future development. The last section of the article provides brief information about the preliminary design of the gas-dynamic multiple-mirror trap device, the development of which has been completed.

Asymptotic quasisymmetric high-beta three-dimensional magnetohydrodynamic equilibria near axisymmetry

Quasisymmetry (QS), a hidden symmetry of the magnetic field strength, is known to support nested flux surfaces and provide superior particle confinement in stellarators. In this work, we study the ideal magnetohydrodynamic (MHD) equilibrium and stability of high-beta plasma in a large-aspect-ratio stellarator. In particular, we show that the lowest-order description of a near-axisymmetric equilibrium vastly simplifies the problem of three-dimensional quasisymmetric MHD equilibria, which can be reduced to a standard elliptic Grad–Shafranov equation for the flux function. We show that any large-aspect-ratio tokamak, deformed periodically in the vertical direction, is a stellarator with approximate volumetric QS. We discuss exact analytical solutions and numerical benchmarks. Finally, we discuss the ideal ballooning and interchange stability of some of our equilibrium configurations.

On the magnetohydrodynamic stability of channel flow of Navier–Stokes–Voigt fluid

The linear stability of a pressure-driven channel flow of an electrically conducting Navier–Stokes–Voigt type of the viscoelastic fluid subject to a transverse magnetic field is investigated. The validity of Squire's theorem is proved, and the generalized eigenvalue problem for two-dimensional modes is obtained by adopting the Galerkin method, which is subsequently solved using the QZ-algorithm. Although the base flow retains its Newtonian fluid characteristics, the noticeable influence of the Kelvin–Voigt parameter in conjunction with the Hartmann number on the stability of fluid flow is perceived. Instability is exclusively identified within a specific range of the Kelvin–Voigt parameter, markedly affected by the Hartmann number. Furthermore, closed neutral stability curves arise, indicating the requirement of two values of the Reynolds number to completely assess the linear stability criteria, in contrast to the typical single value observed in the case of Newtonian fluids.

Experimental issues of energy balance in open magnetic trap

The paper presents an overview of experimental results of an investigation of different energy loss channels in the gas dynamic trap (GDT), which is a magnetic mirror plasma confinement device in the Budker Institute of Nuclear Physics. Energy losses along magnetic field lines are considered as well as losses onto radial limiters, which restrict the plasma column radius and provide its magnetohydrodynamic stability via the ‘vortex confinement’ mechanism. The losses along the field lines were measured using a set of pyroelectric bolometers on the plasma absorber and the losses onto the limiters were determined with thermistors from their temperature rise. Additionally, the losses due to charge exchange of fast plasma ions on the residual neutral gas in the GDT were measured using a longitudinal array of pyroelectric bolometers mounted on the wall of the central cell. An attempt was made to draw up the energy balance in the GDT in order to identify the predominant loss channels and reduce those losses in the future.

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).

Prospects for a high-field, compact break-even axisymmetric mirror (BEAM) and applications

This paper explores the feasibility of a break-even-class mirror referred to as BEAM (break-even axisymmetric mirror): a neutral-beam-heated simple mirror capable of thermonuclear-grade parameters and $Q\sim 1$ conditions. Compared with earlier mirror experiments in the 1980s, BEAM would have: higher-energy neutral beams, a larger and denser plasma at higher magnetic field, both an edge and a core and capabilities to address both magnetohydrodynamic and kinetic stability of the simple mirror in higher-temperature plasmas. Axisymmetry and high-field magnets make this possible at a modest scale enabling a short development time and lower capital cost. Such a $Q\sim 1$ configuration will be useful as a fusion technology development platform, in which tritium handling, materials and blankets can be tested in a real fusion environment, and as a base for development of higher- $Q$ mirrors.

Enhancing Magnetohydrodynamic Stability in Channel Flows through Porous Media Incorporating Darcy-Forchheimer Modifications into the Orr-Sommerfeld Analysis

This research paper investigates the stability characteristics of magnetohydrodynamic (MHD) flow within channels filled with porous media. Specifically, the study explores modifications to the classical Orr-Sommerfeld equation using the Darcy-Forchheimer model to account for the presence of porous media. The analysis aims to elucidate the combined influence of porous media and magnetic fields on flow stability in such systems. Numerical simulations and analytical techniques are employed to assess the impact of varying parameters on the stability behavior of the flow. The study revealed that the velocity profile increases with increase in Schmidt and Prandtl’s number and stabilizes with the increase in Grashof number and the buoyancy. The temperature profile increases with Reynold’s Number and Hartmann number. It stabilizes with Schmidt, Prandtl and Grashof number. The profile of species concentration stabilizes with thermo-phoretic parameter, Prandtl number and Darcy-Forchheimer’s parameter augumentation.

Interplay of plasma resistivity and rotation on β limits in free boundary diverted tokamaks

A plasma resistivity-β driving mechanism aimed at explaining the appearance of long wavelength global instabilities in free boundary high-β tokamak plasmas with a divertor is presented. These perturbations resemble very closely the resistive wall mode phenomenon. Performing a proper toroidal analysis, we show that the magnetohydrodynamic stability is worsened by the interplay of plasma β and resistivity. By modelling the effect of a magnetic separatrix through a careful positioning of the resonant surfaces, we find that in an ideal plasma wall effects are effectively screened, so that the ideal β limit becomes independent of the wall position/physics. A lower wall dependent critical β is found if plasma resistivity is allowed. We find that global stability can be improved with a toroidal flow, small enough not to induce equilibrium modification. The rotation stabilisation effectiveness depends upon the proximity of the plasma equilibrium parameters to the resistive marginal boundary.

Trapped-particle precession and modes in quasisymmetric stellarators and tokamaks: a near-axis perspective

This paper presents the calculation of the bounce-averaged drift of trapped particles in a near-axis framework for axisymmetric and quasisymmetric magnetic fields that possess up-down and stellarator symmetry, respectively. This analytic consideration provides important insight on the dependence of the bounce-averaged drift on the geometry and stability properties of the field. In particular, we show that although the maximum- $\mathcal {J}$ property is unattainable in quasisymmetric stellarators, one may approach it through increased plasma $\beta$ and triangular shaping, albeit going through a reduced precession scenario with potentially higher particle losses. The description of trapped particles allows us to calculate the available energy of trapped electrons analytically in two asymptotic regimes, providing insight into the behaviour of this measure of turbulence. It is shown that the available energy is intimately related to magnetohydrodynamics (MHD) stability, providing a potential synergy between this measure of gyrokinetic turbulence and MHD stability.

The effect of pressure anisotropy on 3-D MHD stability for low magnetic field LHD equilibria

The magnetohydrodynamic stability of plasmas with an anisotropic pressure component is analysed for a low magnetic field configuration of the large helical device. Magnetic equilibria are calculated by the anisotropic Neumann inverse moments equilibrium code, an extension of the three-dimensional variational moments equilibrium code. A modified version of the bi-Maxwellian is used to model the anisotropic particle velocity distribution. Magnetohydrodynamic stability calculations for the $n=1$ mode family are carried out by TERPSICHORE, which has been expanded by the Kruskal–Oberman energy principle. For on-axis particle deposition, the growth rate and plasma displacement show that the parallel dominant plasmas are significantly more stable than isotropic or perpendicular dominant plasmas. For off-axis particle deposition, the growth rate and the Mercier criterion in the peripheral region $\rho =0.9$ , show that low field (LF) deposition perpendicular dominant plasmas are most unstable. For the most realistic off-axis deposition profile, it is found that parallel dominant plasmas are most stable for LF deposition, while perpendicular dominant plasmas are most stable for high field deposition. We conclude that, under low magnetic field conditions in the large helical device, tangential neutral beam injection heating has a stabilising influence on the plasma.

Equilibriumβ-limits dependence on bootstrap current in classical stellarators

While it is important to design stellarators with high magnetohydrodynamic stability$\beta$-limit, it is also crucial to ensure that good magnetic surfaces exist in a large range of$\beta$values. As$\beta$increases, pressure-driven currents perturb the vacuum magnetic field and often lead to the emergence of magnetic field line chaos, which can worsen the confinement and is the cause of another kind of$\beta$-limit, the so-called equilibrium$\beta$-limit. In this paper, we explore numerically the dependence of the equilibrium$\beta$-limit on the bootstrap current strength in a classical stellarator geometry using the stepped pressure equilibrium code. We develop a diagnostic to determine whether or not magnetic islands are expected to participate significantly to radial transport, and we build an analytical model to predict the expected equilibrium$\beta$-limit, which recovers the main features of the numerical results. This research opens the possibility to include additional targets in stellarator optimization functions, provides additional understanding on the existence of magnetic surfaces at large$\beta$, and is a step forward in the understanding of the equilibrium$\beta$-limit.