Magnetohydrodynamic Equilibrium Research Articles

Three-dimensional magnetohydrodynamic equilibrium of quiescent H-modes in tokamak systems

Three dimensional free boundary magnetohydrodynamic equilibria that recover saturated ideal kink/peeling structures are obtained numerically. Simulations that model the JET tokamak at fixed with a large edge bootstrap current that flattens the q-profile near the plasma boundary demonstrate that a radial parallel current density ribbon with a dominant m /n = 5/1 Fourier component at MA develops into a broadband spectrum when the toroidal current It is increased to 2.5 MA.

Saturated ideal kink/peeling formations described as three-dimensional magnetohydrodynamic tokamak equilibrium states

Free boundary magnetohydrodynamic equilibrium states with spontaneous three dimensional deformations of the plasma-vacuum interface are computed for the first time. The structures obtained have the appearance of saturated ideal external kink/peeling modes. High edge pressure gradients yield toroidal mode number n = 1 corrugations for a high edge bootstrap current and larger n distortions when this current is small. Deformations in the plasma boundary region induce a nonaxisymmetric Pfirsch-Schlüter current driving a field-aligned current ribbon consistent with reported experimental observations. A variation in the 3D equilibrium confirms that the n = 1 mode is a kink/peeling structure. We surmise that our calculated equilibrium structures constitute a viable model for the edge harmonic oscillations and outer modes associated with a quiescent H-mode operation in shaped tokamak plasmas.

On the Validation of Gyrokinetic L-Mode Simulations

Over the last decade, plasma turbulence simulations based on gyrokinetic theory have reached an amazing degree of physical comprehensiveness and realism. In contrast to early gyrokinetic studies, which were restricted to qualitative statements, state-of-the-art investigations may now be compared quantitatively, therefore enabling validation and detailed analysis of their predictive capabilities. Here, particular attention is paid to outer-core L-mode discharges for which some previous gyrokinetic studies have found an underprediction of ion heat transport by almost one order of magnitude, the so-called shortfall. Carrying out radially local and nonlocal GENE simulations using actual plasma profiles and parameters and magnetohydrodynamic equilibria, and employing as much physics as available, only a mild underprediction is found, which can, furthermore, be overcome by varying the ion temperature gradient within the error bars associated with the experimental measurement. The significance and reliability of these simulations is furthermore demonstrated by extensive comparison with experimental measurements. The latter involve sophisticated synthetic beam emission spectroscopy and correlation electron cyclotron emission data analysis. The agreement found between the measurements and the state-of-the-art postprocessed simulation data confirms the high degree of realism.

Numerical solution of ideal MHD equilibrium via radial basis functions collocation and moving least squares approximation methods

In this study, two different meshfree methods consisting of the Radial Basis Functions (RBFs) and the Moving Least Square Method (MLS) are applied to solve the Grad–Shafranov (GS) equation for the axisymmetric equilibrium of plasma in the tokamak. The validity and the effectiveness of the proposed schemes are studied by several test problems through absolute and Root Mean Squared (RMS) errors. Although, during the past few years, a meshfree method is normally applied in magnetohydrodynamic (MHD) studies to the numerical solution of partial differential equations (PDEs) but to the best of our knowledge, its application in MHD equilibrium of the tokamak plasma investigations is rare. The future more extensive studies regarding this numerical method would definitely have a significant impact on improving tokamak numerical tools.

Investigation of island formation due to RMPs in DIII-D plasmas with the SIESTA resistive MHD equilibrium code

The SIESTA magnetohydrodynamic (MHD) equilibrium code has been used to compute a sequence of ideally stable equilibria resulting from numerical variation of the helical resonant magnetic perturbation (RMP) applied to an axisymmetric DIII-D plasma equilibrium. Increasing the perturbation strength at the dominant $m=2$, $n=-1$ resonant surface leads to lower MHD energies and increases in the equilibrium island widths at the $m=2$ (and sidebands) surfaces, in agreement with theoretical expectations. Island overlap at large perturbation strengths leads to stochastic magnetic fields which correlate well with the experimentally inferred field structure. The magnitude and spatial phase (around the dominant rational surfaces) of the resonant (shielding) component of the parallel current are shown to change qualitatively with the magnetic island topology.

Radial magnetic compression in the expelled jet of a plasma deflagration accelerator

A spectroscopic study of a pulsed plasma deflagration accelerator is carried out that confirms the existence of a strong compression in the emerging jet at the exit plane of the device. An imaging spectrometer is used to collect broadened Hα emission from a transaxial slice of the emerging jet at high spatial resolution, and the radial plasma density profile is computed from Voigt fits of the Abel inverted emissivity profiles. The plasma temperature, determined via Doppler broadening of impurity line emission, is compared against the temperature predictions of a radial magnetohydrodynamic equilibrium model applied to the measured density profiles. Empirical scaling laws developed for the plasma density, combined with the measured and predicted temperatures, indicate that a radially equilibrated Z-pinch is formed within the expelled plasma jet at the exit plane during the deflagration process.

Gyrokinetic particle simulation of a field reversed configuration

Gyrokinetic particle simulation of the field-reversed configuration (FRC) has been developed using the gyrokinetic toroidal code (GTC). The magnetohydrodynamic equilibrium is mapped from cylindrical coordinates to Boozer coordinates for the FRC core and scrape-off layer (SOL), respectively. A field-aligned mesh is constructed for solving self-consistent electric fields using a semi-spectral solver in a partial torus FRC geometry. This new simulation capability has been successfully verified and driftwave instability in the FRC has been studied using the gyrokinetic simulation for the first time. Initial GTC simulations find that in the FRC core, the ion-scale driftwave is stabilized by the large ion gyroradius. In the SOL, the driftwave is unstable on both ion and electron scales.

A 3-D MHD equilibrium description of nonlinearly saturated ideal external kink/peeling structures in tokamaks

Novel free boundary magnetohydrodynamic equilibrium states with spontaneous three-dimensional (3-D) deformations of the plasma–vacuum interface are computed. The structures obtained look like saturated ideal external kink/peeling modes. Large edge pressure gradients yield toroidal mode number $n=1$ distortions when the edge bootstrap current is large and higher $n$ corrugations when this current is small. Linear ideal MHD stability analyses confirm the nonlinear saturated ideal kink equilibrium states produced and we can identify the Pfirsch–Schlüter current as the main linear instability driving mechanism when the edge pressure gradient is large. The dominant non-axisymmetric component of this Pfirsch–Schlüter current drives a near resonant helical parallel current density ribbon that aligns with the near vanishing magnetic shear region caused by the edge bootstrap current. This current ribbon is a manifestation of the outer mode previously found on JET (Solano 2010). We claim that the equilibrium corrugations describe structures that are commonly observed in quiescent H-mode tokamak discharges.

Structure and Stability of Filamentary Clouds Supported by Lateral Magnetic Field

AbstractWe have constructed two types of analytical models for an isothermal filamentary cloud supported mainly by magnetic tension. The first one describes an isolated cloud while the second considers filamentary clouds spaced periodically. The filamentary clouds are assumed to be highly flattened in both the models. The former is proved to be the asymptotic limit of the latter in which each filamentary cloud is much thinner than the distance to the neighboring filaments. These models show that the mass to flux ratio is crucial for the magnetohydrodynamical equilibrium. The upper bound for the line density, i.e., the mass per unit length, is proportional to the magnetic flux. The mass to flux ratio is slightly larger than the critical value, ($2 \pi \sqrt{G}$)−1, in the first model and lower in the second model. The first model is unstable against fragmentation and the wavelength of the fastest growing mode is several times longer than the cloud diameter. The second model is likely to be unstable only when the mass to flux ratio is supercritical.

Validation of equilibrium tools on the COMPASS tokamak

Various MHD (magnetohydrodynamic) equilibrium tools, some of which being recently developed or considerably updated, are used on the COMPASS tokamak at IPP Prague. MHD equilibrium is a fundamental property of the tokamak plasma, whose knowledge is required for many diagnostics and modelling tools. Proper benchmarking and validation of equilibrium tools is thus key for interpreting and planning tokamak experiments. We present here benchmarks and comparisons to experimental data of the EFIT++ reconstruction code (Appel et al., 2006), the free-boundary equilibrium code FREEBIE (Artaud and Kim, 2012), and a rapid plasma boundary reconstruction code VacTH (Faugeras et al., 2014). We demonstrate that FREEBIE can calculate the equilibrium and corresponding poloidal field (PF) coils currents consistently with EFIT++ reconstructions from experimental data. Both EFIT++ and VacTH can reconstruct equilibria generated by FREEBIE from synthetic, optionally noisy diagnostic data. Hence, VacTH is suitable for real-time control. Optimum reconstruction parameters are estimated.

Integrated physics analysis of plasma start-up scenario of helical reactor FFHR-d1

1D physics analysis of the plasma start-up scenario of the large helical device (LHD)-type helical reactor FFHR-d1 was conducted. The time evolution of the plasma profile is calculated using a simple model based on the LHD experimental observations. A detailed assessment of the magnetohydrodynamic equilibrium and neo-classical energy loss was conducted using the integrated transport analysis code TASK3D. The robust controllability of the fusion power was confirmed by feedback control of the pellet fuelling and a simple staged variation of the external heating power with a small number of simple diagnostics (line-averaged electron density, edge electron density and fusion power). A baseline operation control scenario (plasma start-up and steady-state sustainment) of the FFHR-d1 reactor for both self-ignition and sub-ignition operation modes was demonstrated.

Conceptual design study of the moderate size superconducting spherical tokamak power plant

A new conceptual design of the superconducting spherical tokamak (ST) power plant was proposed as an attractive choice for tokamak fusion reactors. We reassessed a possibility of the ST as a power plant using the conservative reactor engineering constraints often used for the conventional tokamak reactor design. An extensive parameters scan which covers all ranges of feasible superconducting ST reactors was completed, and five constraints which include already achieved plasma magnetohydrodynamic (MHD) and confinement parameters in ST experiments were established for the purpose of choosing the optimum operation point. Based on comparison with the estimated future energy costs of electricity (COEs) in Japan, cost-effective ST reactors can be designed if their COEs are smaller than 120 mills kW−1 h−1 ($2013). We selected the optimized design point: A = 2.0 and Rp = 5.4 m after considering the maintenance scheme and TF ripple. A self-consistent free-boundary MHD equilibrium and poloidal field coil configuration of the ST reactor were designed by modifying the neutral beam injection system and plasma profiles. The MHD stability of the equilibrium was analysed and a ramp-up scenario was considered for ensuring the new ST design. The optimized moderate-size ST power plant conceptual design realizes realistic plasma and fusion engineering parameters keeping its economic competitiveness against existing energy sources in Japan.

Adaptive Grids in Simulations of Toroidal Plasma Starting from Magneto-Hydrodynamic Equilibrium

This paper introduces the notion of Tokamak Magneto-Hydrodynamics (TMHD), which explicitly reflects the anisotropy of a high temperature tokamak plasma. The set of TMHD equations is formulated for simulation of macroscopic plasma dynamics and disruptions in tokamaks. Free from the Courant restriction on the time step, this set of equations is adequate to plasma dynamics with realistic parameters of high performance plasmas and does not require any extension of the MHD plasma model. At the same time, TMHD requires the use of magnetic field aligned numerical grids. Examples of their use in 2-dimensional cases of tokamak equilibria and dynamics of the wall touching kink mode are presented. For the 3-dimensional case of an ergodic magnetic field, this paper introduces the reference magnetic coordinates as a practical algorithm for generating adaptive grids for TMHD.

Magnetohydrodynamic equilibria with incompressible flows: Symmetry approach

We identify and discuss a family of azimuthally symmetric, incompressible, magnetohydrodynamic plasma equilibria with poloidal and toroidal flows in terms of solutions of the Generalized Grad Shafranov (GGS) equation. These solutions are derived by exploiting the incompressibility assumption, in order to rewrite the GGS equation in terms of a different dependent variable, and the continuous Lie symmetry properties of the resulting equation and, in particular, a special type of “weak” symmetries.

Trigonometric and hyperbolic functions method for constructing analytic solutions to nonlinear plane magnetohydrodynamics equilibrium equations

In this paper, we present a solution method for constructing exact analytic solutions to magnetohydrodynamics (MHD) equations. The method is constructed via all the trigonometric and hyperbolic functions. The method is applied to MHD equilibria with mass flow. Applications to a solar system concerned with the properties of coronal mass ejections that affect the heliosphere are presented. Some examples of the constructed solutions which describe magnetic structures of solar eruptions are investigated. Moreover, the constructed method can be applied to a variety classes of elliptic partial differential equations which arise in plasma physics.

Magnetic islands and singular currents at rational surfaces in three-dimensional magnetohydrodynamic equilibria

Using the recently developed multiregion, relaxed MHD (MRxMHD) theory, which bridges the gap between Taylor's relaxation theory and ideal MHD, we provide a thorough analytical and numerical proof of the formation of singular currents at rational surfaces in non-axisymmetric ideal MHD equilibria. These include the force-free singular current density represented by a Dirac δ-function, which presumably prevents the formation of islands, and the Pfirsch-Schlüter 1/x singular current, which arises as a result of finite pressure gradient. An analytical model based on linearized MRxMHD is derived that can accurately (1) describe the formation of magnetic islands at resonant rational surfaces, (2) retrieve the ideal MHD limit where magnetic islands are shielded, and (3) compute the subsequent formation of singular currents. The analytical results are benchmarked against numerical simulations carried out with a fully nonlinear implementation of MRxMHD.

Geodesic acoustic mode in anisotropic plasmas using double adiabatic model and gyro-kinetic equation

Geodesic acoustic mode in anisotropic tokamak plasmas is theoretically analyzed by using double adiabatic model and gyro-kinetic equation. The bi-Maxwellian distribution function for guiding-center ions is assumed to obtain a self-consistent form, yielding pressures satisfying the magnetohydrodynamic (MHD) anisotropic equilibrium condition. The double adiabatic model gives the dispersion relation of geodesic acoustic mode (GAM), which agrees well with the one derived from gyro-kinetic equation. The GAM frequency increases with the ratio of pressures, p⊥/p∥, and the Landau damping rate is dramatically decreased by p⊥/p∥. MHD result shows a low-frequency zonal flow existing for all p⊥/p∥, while according to the kinetic dispersion relation, no low-frequency branch exists for p⊥/p∥≳ 2.

Equilibrium solutions of relativistic rotating stars with mixed poloidal and toroidal magnetic fields

Stationary and axisymmetric solutions of relativistic rotating stars with strong mixed poloidal and toroidal magnetic fields are obtained numerically. Because of the mixed components of the magnetic field, the underlying stationary and axisymmetric spacetimes are no longer circular. These configurations are computed from the full set of the Einstein-Maxwell equations, Maxwell's equations, and from first integrals and integrability conditions of the magnetohydrodynamic equilibrium equations. After a brief introduction of the formulation of the problem, we present the first results for highly deformed magnetized rotating compact stars.

Stability analysis of the polytropic solar wind

A linear stability analysis of a simple polytropic model for the solar wind dynamics within the framework of a magnetohydrodynamic equilibrium configuration is theoretically proposed. The simplistic analysis is based on the model developed based on the data available from the Advanced Composition Explorer (ACE) spacecraft mission. A unique form of dispersion relation is derived by coupling the adiabatic and polytropic processes in the limit of ideal gas approximation for the solar wind gas in accordance with the standard Fourier technique. Applying usual variable-separation methodology on the dispersion relation, we obtain the linear growth rate of the fluctuations. It is seen that the growth rate is an explicitly nonlinear function of the variable polytropic index (α) and radial position (r) with respect to the considered center of the Sun. Numerical analyses are carried out to understand the physical insight of the growth profiles of the fluctuations. It is shown that the growth is maximized near the solar corona, where α ∼ 1, relative to that observed elsewhere in the entire solar plasma system. The source for this growth may be attributed to the free flow of energy coming from the dynamic equilibrium of the solar plasma itself. As compared with existing model predictions, our results are qualitatively capable of reproducing the average behavior of the solar wind fluctuation and stability behaviors on the astrophysical scales of space and time.

Self-consistent stationary MHD shear flows in the solar atmosphere as electric field generators

Magnetic fields and flows in coronal structures, for example, in gradual phases in flares, can be described by 2D and 3D magnetohydrostatic (MHS) and steady magnetohydrodynamic (MHD) equilibria. Within a physically simplified, but exact mathematical model, we study the electric currents and corresponding electric fields generated by shear flows. Starting from exact and analytically calculated magnetic potential fields, we solveid the nonlinear MHD equations self-consistently. By applying a magnetic shear flow and assuming a nonideal MHD environment, we calculated an electric field via Faraday's law. The formal solution for the electromagnetic field allowed us to compute an expression of an effective resistivity similar to the collisionless Speiser resistivity. We find that the electric field can be highly spatially structured, or in other words, filamented. The electric field component parallel to the magnetic field is the dominant component and is high where the resistivity has a maximum. The electric field is a potential field, therefore, the highest energy gain of the particles can be directly derived from the corresponding voltage. In our example of a coronal post-flare scenario we obtain electron energies of tens of keV, which are on the same order of magnitude as found observationally. This energy serves as a source for heating and acceleration of particles.