Effect of Magnetic Field Configuration on Stability

In magnetically confined plasmas, the configuration of the magnetic field dictates far more than just particle orbits and confinement efficiency—it directly governs macroscopic stability. Magnetohydrodynamics (MHD) treats the plasma as a conducting fluid, where its stability can be evaluated through energy principles: while magnetic tension generally provides a stabilizing effect, magnetic pressure, current gradients, and field curvature supply the free energy necessary for instabilities to grow. When a magnetic geometry is poorly optimized, minute perturbations can rapidly evolve along field lines or in the radial direction into interchange, kink, tearing, or ballooning modes, ultimately causing severe confinement degradation or catastrophic disruptions.
Ideal magnetic confinement configurations are typically characterized by nested magnetic flux surfaces. The helical pitch of magnetic field lines winding around these surfaces is quantified by the safety factor, denoted as $q$:

  • In tokamaks, $q$ is defined as the ratio of toroidal transits to poloidal transits;
  • Rational surfaces occur where $q = m/n$ (with $m$ and $n$ being integers);
  • Magnetic shear, $s = \frac{r}{q}\frac{dq}{dr}$, measures the differential twist between adjacent flux surfaces.

Magnetic shear serves as a cornerstone parameter for stability. Strong shear creates phase mismatch for radial perturbations across different flux surfaces, effectively suppressing ideal kink modes and certain ballooning modes. However, excessively high shear or peaked current profiles can simultaneously drive resistive tearing modes near rational surfaces, leading to the formation of magnetic islands.

Curvature-Driven Instabilities: Interchange and Ballooning Modes

Magnetic field curvature plays a dual role in plasma stability. If the curvature vector and the pressure gradient point in mutually unfavorable directions—a condition known as "bad curvature"—the plasma behaves analogously to a heavy fluid resting atop a light fluid, rendering it susceptible to interchange instabilities. Typically, the outboard side of a tokamak resides in a bad curvature region, whereas the inboard side enjoys good curvature.

Primary mechanisms to stabilize bad curvature include:

  • Magnetic Wells: Designing a configuration where the magnetic field strength features a local minimum along a field line, creating an average minimum-$B$ geometry that suppresses interchange modes;
  • Magnetic Shear: Restricting the radial expansion of localized perturbations;
  • Finite Poloidal Fields: Utilizing poloidal magnetic fields in tokamaks to establish a localized magnetic well.

Ballooning modes represent localized modes driven jointly by pressure gradients and unfavorable curvature, with their stability thresholds conventionally scaled by the critical $\beta$ parameter. Pushing $\beta$ higher enhances fusion power density, but doing so inevitably drives the system closer to the ballooning limit. Consequently, magnetic configuration design constantly navigates a trade-off between energy confinement and operational stability.

Current Profiles and MHD Modes

The internal current density profile fundamentally dictates the $q$ profile, which in turn governs kink and tearing instabilities. The classic benchmark is the Kruskal-Shafranov limit: in standard tokamaks, if the edge safety factor drops below unity ($q_a < 1$), the ideal $m=1, n=1$ external/internal kink mode typically becomes unstable. Therefore, a foundational design criterion is maintaining $q_a > 1$.

Internal kink modes are closely tied to sawtooth oscillations. When the core $q$ value dips below unity, periodic magnetic reconnection can trigger within the $q=1$ rational surface, causing abrupt drops in core electron temperature. Tearing modes, by contrast, materialize at rational surfaces such as $q=2$ or $q=3/2$, where finite resistivity allows field lines to break and reconnect, producing magnetic islands. Neoclassical tearing modes (NTMs), driven by perturbed bootstrap currents, can further degrade confinement and trigger major disruptions.

Boundary Control and Three-Dimensional Configurations

The configuration of the plasma boundary is equally vital. In diverted configurations, the magnetic X-point guides open field lines into a divertor chamber, effectively managing heat loads and impurity influx. Edge localized modes (ELMs) are intimately linked to steep edge pressure gradients and localized bootstrap currents. Resonant magnetic perturbations (RMPs)—applied non-axisymmetric fields that generate controlled stochastic boundary layers—have proven effective at mitigating ELMs, though careful optimization is required to prevent excessive core confinement loss.

Stellarators abandon net toroidal plasma current in favor of fully three-dimensional magnetic field designs, which inherently avoid current-driven kink instabilities while enabling intrinsic magnetic wells, quasi-symmetry, and optimized particle drift orbits. Nevertheless, the engineering complexity of 3D coil fabrication and neoclassical transport optimization remains formidable. Alternative concepts, such as field-reversed configurations (FRCs) and compact toroids, explore alternative regions of the curvature-shear parameter space.

Numerical Modeling Workflow

Numerical investigation of plasma stability typically follows a two-step paradigm: equilibrium calculation followed by linear or non-linear stability analysis.

  1. Compute axisymmetric equilibria using the Grad-Shafranov equation, or 3D MHD equilibria utilizing specialized codes like VMEC;
  2. Extract critical profiles including the safety factor $q(r)$, pressure $p(r)$, magnetic shear $s$, and local field curvature;
  3. Solve the ideal or resistive MHD eigenvalue equations to determine growth rates and 2D/3D eigenmode structures;
  4. Perform parameter scans over $\beta$, $q_a$, current profile shapes, and boundary magnetic fields to map out operational stability windows.

For instance, adopting a parabolic cylindrical approximation $q(r) = q_0 + (q_a - q_0)(r/a)^2$, if $q_a < 1$, the $m=1, n=1$ kink mode is likely to grow. Raising $q_a$ above unity stabilizes this ideal kink, though resistive tearing modes may still emerge at the $q=2$ surface. While augmenting magnetic shear generally raises the ballooning $\beta$ threshold, it requires precise coordination with rational surface locations to avoid deleterious island growth.

Summary

Ultimately, magnetic field configuration dictates plasma stability through a complex interplay of magnetic curvature, shear, magnetic wells, current profiles, and boundary shaping. Effective design principles can be summarized as follows: avoiding extensive bad curvature regions, maintaining a sufficient edge safety factor $q_a$, tailoring current and shear profiles, leveraging magnetic wells along with 3D shaping, and exercising rigorous control over boundary resonances. Mastering these geometric controls remains the fundamental bedrock for advanced fusion device design and predictive plasma modeling.