Confining a thermonuclear plasma at temperatures exceeding one hundred million Kelvin requires eliminating material wall contact using closed magnetic geometries. In an axisymmetric torus, such as a tokamak, force balance between thermal kinetic pressure and electromagnetic Lorentz forces is governed by a singular, elegant nonlinear partial differential equation. In this treatise, we derive the Grad-Shafranov Equation from first principles and analyze the structural equilibrium of toroidal fusion devices.
1. The Magnetostatic Equilibrium Condition
In a static, steady-state plasma ((mathbf{v} = 0, partial_t = 0)), the momentum conservation equation reduces to the magnetostatic force balance:
[
abla P = mathbf{J} imes mathbf{B}
]
where (P) is scalar isotropic plasma pressure, (mathbf{J}) is current density, and (mathbf{B}) is magnetic field. Ampère’s law and the solenoidal constraint give:
[
abla imes mathbf{B} = mu_0 mathbf{J}, quad
abla cdot mathbf{B} = 0
]
From (
abla P = mathbf{J} imes mathbf{B}), two immediate topological consequences emerge:
[
mathbf{B} cdot
abla P = 0, quad mathbf{J} cdot
abla P = 0
]
Magnetic field lines and electric current lines lie entirely on constant pressure isobaric surfaces known as Magnetic Flux Surfaces.
2. Derivation of the Grad-Shafranov Equation
In cylindrical coordinates ((R, phi, Z)) with axisymmetry ((partial/partialphi = 0)), the condition (
abla cdot mathbf{B} = 0) guarantees that the poloidal magnetic field (mathbf{B}_p = B_R hat{mathbf{R}} + B_Z hat{mathbf{Z}}) can be expressed in terms of a scalar Poloidal Flux Function (psi(R, Z)):
[
mathbf{B}_p = rac{1}{R}
abla psi imes hat{boldsymbol{phi}} = -rac{1}{R}rac{partial psi}{partial Z} hat{mathbf{R}} + rac{1}{R}rac{partial psi}{partial R} hat{mathbf{Z}}
]
The total magnetic field is (mathbf{B} = mathbf{B}_p + B_phi hat{boldsymbol{phi}}). Defining the poloidal current stream function (F(R, Z) quiv R B_phi), Ampère’s law dictates that current density (mathbf{J} = rac{1}{mu_0}
abla imes mathbf{B}) decomposes into poloidal and toroidal components:
[
mathbf{J} = rac{1}{mu_0 R}
abla F imes hat{boldsymbol{phi}} + J_phi hat{boldsymbol{phi}}
]
[
J_phi = -rac{1}{mu_0 R} Delta^* psi
]
where (Delta^*) is the elliptic Shafranov Differential Operator:
[
Delta^* psi quiv R rac{partial}{partial R}left( rac{1}{R} rac{partial psi}{partial R}
ight) + rac{partial^2 psi}{partial Z^2} = R^2
abla cdot left( rac{
abla psi}{R^2}
ight)
]
Substituting (mathbf{B}) and (mathbf{J}) into (mathbf{J} imes mathbf{B} =
abla P) and resolving components yields (P = P(psi)) and (F = F(psi)): pressure and poloidal current are pure flux functions. The radial component yields the celebrated Grad-Shafranov Equation:
[
Delta^* psi = -mu_0 R^2 rac{dP}{dpsi} – F rac{dF}{dpsi}
]
3. Rotational Transform and the Safety Factor (q)
Pure toroidal fields cannot confine a plasma because grad-(B) and curvature drifts cause vertical charge separation, inducing an electric field (mathbf{E}) that ejects the plasma radially outwards via (mathbf{E} imes mathbf{B}) drift. To cancel this drift, a poloidal magnetic field is superimposed to twist field lines into nested helical spirals. The pitch of this helix is quantified by the Safety Factor ((q)):
[
q(psi) = rac{1}{2pi} oint rac{B_phi}{R B_p} ds = rac{F(psi)}{2pi} oint rac{ds}{R^2 B_p}
]
In a large aspect ratio circular tokamak ((R_0 gg a)), (q(r) approx rac{r B_phi}{R_0 B_ heta}). The safety factor determines whether field lines close on themselves after rational numbers of toroidal and poloidal circuits, dictating both magnetohydrodynamic stability and the maximum plasma beta (beta = rac{2mu_0 langle P
angle}{B^2}) achievable before catastrophic disruption.
