Even when an MHD equilibrium satisfies the Grad-Shafranov equation exactly, the configuration is not necessarily stable. Vast reservoirs of free energy stored in plasma pressure gradients and internal magnetic currents can trigger violent macroscopic instabilities that rupture magnetic confinement within microseconds. In this article, we formulate Bernstein’s Energy Principle ((delta W)), derive the Kruskal-Shafranov kink limit, and analyze interchange and ballooning modes.
1. The Ideal MHD Energy Principle ((delta W))
Consider an equilibrium plasma subjected to an arbitrary virtual displacement field (boldsymbol{xi}(mathbf{r}, t)). The linearized equation of motion is (
ho_0 ddot{boldsymbol{xi}} = mathbf{F}(boldsymbol{xi})), where (mathbf{F}) is the self-adjoint force operator. Multiplying by (dot{boldsymbol{xi}}) and integrating over the volume yields kinetic energy and the potential energy perturbation (delta W):
[
delta W(boldsymbol{xi}, boldsymbol{xi}) = -rac{1}{2} int_V boldsymbol{xi} cdot mathbf{F}(boldsymbol{xi}) , d^3x
]
According to Bernstein et al. (1958), a plasma equilibrium is linearly stable if and only if (delta W > 0) for all admissible displacements satisfying boundary conditions. The potential energy decomposes into fluid, surface, and vacuum contributions (delta W = delta W_F + delta W_S + delta W_V), where the fluid potential energy is:
[
delta W_F = rac{1}{2} int_V left[ rac{|mathbf{Q}_perp|^2}{mu_0} + rac{B^2}{mu_0}left|
abla cdot boldsymbol{xi}_perp + 2boldsymbol{xi}_perp cdot boldsymbol{kappa}
ight|^2 + gamma P (
abla cdot boldsymbol{xi})^2 – 2(boldsymbol{xi}_perp cdot
abla P)(boldsymbol{xi}_perp cdot boldsymbol{kappa}) – J_parallel (boldsymbol{xi}_perp imes mathbf{B}) cdot mathbf{Q}_perp
ight] d^3x
]
where (mathbf{Q} quiv
abla imes (boldsymbol{xi} imes mathbf{B})) is the perturbed magnetic field and (boldsymbol{kappa} quiv (hat{mathbf{b}} cdot
abla)hat{mathbf{b}}) is the magnetic field line curvature vector.
2. Current-Driven Instabilities: Sausage and Kink Modes
Instabilities powered by the parallel current density (J_parallel) are Current-Driven (Kink) Modes:
- Sausage Mode ((m = 0)): In a cylindrical (Z)-pinch, axisymmetric radial pinching creates constrictions. At narrowed necks, the surface azimuthal field (B_ heta propto 1/r) intensifies, increasing the magnetic compression pressure (B_ heta^2/2mu_0) and driving a runaway neck collapse.
- Kink Mode ((m = 1)): Helical bending of the plasma column compresses magnetic field lines on the concave side of the bend while rarefying them on the convex side. The resulting magnetic pressure differential amplifies the bend, throwing the column into the vessel wall.
Stability against the external kink mode requires that the magnetic field lines twist less than one full toroidal turn over the column length (L), establishing the fundamental Kruskal-Shafranov Criterion:
[
q(a) = rac{2pi a B_z}{L B_ heta(a)} > 1
]
If (q(a) < 1), the (m=1, n=1) kink mode becomes violently unstable on the Alfvén timescale ( au_A sim a / v_A).
3. Pressure-Driven Instabilities: Interchange and Ballooning Modes
Instabilities powered by the pressure gradient (
abla P) represent the plasma analog of the hydrodynamic Rayleigh-Taylor instability:
- Interchange (Flute) Instabilities: Adjacent magnetic flux tubes swap places. If plasma moves into a region of weaker magnetic field ((
abla P cdot boldsymbol{kappa} > 0), “bad curvature”), potential energy is released ((delta W < 0)). Stabilization is achieved via Magnetic Shear (s quiv rac{r}{q}rac{dq}{dr}) and favorable average curvature, quantified by the Mercier Criterion:
[
-D_I = rac{1}{4}left(rac{r}{q}rac{dq}{dr}
ight)^2 + rac{2mu_0 r}{B_z^2}rac{dP}{dr}left(1 – q^2
ight) > 0
] - Ballooning Modes: At high plasma pressure ((beta gg 1)), perturbations localize predominantly in regions of unfavorable curvature on the outboard midplane of a torus, ballooning outwards and setting the strict Troyon Beta Limit (beta_{ ext{max}} approx beta_N rac{I_p}{a B_t}) for magnetic confinement.
