Derivation of Magnetohydrodynamic Equilibrium Equations

In the field of plasma physics, the Magnetohydrodynamics (MHD) model serves as the fundamental framework for analyzing the macroscopic behavior of high-temperature plasmas. When a plasma system evolves on a timescale much slower than its internal microscopic processes, it can be treated as being in a quasi-static state. Analyzing the MHD equilibrium is a critical prerequisite for understanding plasma confinement performance in magnetic confinement fusion (MCF) devices, such as Tokamaks and Stellarators. This article provides a systematic derivation of the MHD equilibrium equations and explores their profound physical implications.

The Ideal MHD Framework

The MHD model treats the plasma as a single, electrically conducting fluid. Its dynamics are governed by the coupling of fluid mechanics equations with Maxwell's equations. To derive the equilibrium state, we must first establish the governing equations of ideal MHD:

  • The Momentum Equation: Describes the motion of the plasma under various forces:
    $$ \rho \frac{d\mathbf{v}}{dt} = -\nabla p + \mathbf{J} \times \mathbf{B} + \rho \mathbf{g} $$
  • Ampère’s Law (neglecting the displacement current):
    $$ \mu_0 \mathbf{J} = \nabla \times \mathbf{B} $$
  • Gauss’s Law for Magnetism:
    $$ \nabla \cdot \mathbf{B} = 0 $$

In these expressions, $\rho$ represents the fluid mass density, $\mathbf{v}$ is the velocity field, $p$ is the kinetic (thermal) pressure, $\mathbf{J}$ is the current density, $\mathbf{B}$ is the magnetic induction, and $\mu_0$ is the permeability of free space. In the context of magnetic confinement, the gravitational term $\rho \mathbf{g}$ is typically negligible compared to the electromagnetic forces and is therefore omitted from the equilibrium analysis.

Assumptions for Equilibrium

A state of equilibrium implies that the macroscopic physical properties of the system are stationary over time, meaning $\partial / \partial t = 0$. Furthermore, we assume a static equilibrium, where the macroscopic plasma velocity $\mathbf{v}$ is zero or sufficiently small to be neglected.

By applying these assumptions to the momentum equation, the inertial term $\rho \frac{d\mathbf{v}}{dt}$ vanishes. The equation then simplifies to a force balance condition, indicating that the internal pressure gradient is exactly offset by the Lorentz force.

Derivation of the Equilibrium Equation

Under the static equilibrium assumption, the momentum equation reduces to:

$$ -\nabla p + \mathbf{J} \times \mathbf{B} = 0 $$

To gain deeper insight into the relationship between the magnetic field and the plasma pressure, we express the Lorentz force $\mathbf{J} \times \mathbf{B}$ explicitly in terms of the magnetic field $\mathbf{B}$. By substituting Ampère’s Law into the Lorentz force term, we obtain:

$$ \mathbf{J} \times \mathbf{B} = \frac{1}{\mu_0} (\nabla \times \mathbf{B}) \times \mathbf{B} $$

Using the vector identity $(\nabla \times \mathbf{B}) \times \mathbf{B} = (\mathbf{B} \cdot \nabla)\mathbf{B} - \nabla \left( \frac{B^2}{2} \right)$, the Lorentz force can be decomposed into two distinct components:

$$ \mathbf{J} \times \mathbf{B} = \frac{1}{\mu_0} (\mathbf{B} \cdot \nabla)\mathbf{B} - \nabla \left( \frac{B^2}{2\mu_0} \right) $$

Substituting this back into the force balance equation and rearranging the terms yields the general vector form of the MHD equilibrium equation:

$$ \nabla p + \nabla \left( \frac{B^2}{2\mu_0} \right) = \frac{1}{\mu_0} (\mathbf{B} \cdot \nabla)\mathbf{B} $$

In this form, the left-hand side represents the sum of the thermal pressure gradient and the magnetic pressure gradient, while the right-hand side represents the magnetic tension.

Physical Interpretation

The derived equation reveals that the magnetic field exerts two fundamentally different types of forces on the plasma:

  1. Magnetic Pressure: Represented by the term $-\nabla \left( \frac{B^2}{2\mu_0} \right)$. The quantity $p_B = \frac{B^2}{2\mu_0}$ has the same dimensions as fluid pressure and acts to push the plasma. It represents the tendency of the magnetic field to expand and resist being compressed.
  2. Magnetic Tension: Represented by the term $\frac{1}{\mu_0} (\mathbf{B} \cdot \nabla)\mathbf{B}$. Much like the restoring force of a stretched rubber band, magnetic tension acts along the magnetic field lines, resisting any curvature or bending of the lines.

For a stable equilibrium to exist, the plasma's thermal pressure gradient must be balanced by a combination of magnetic pressure and magnetic tension. To quantify the efficiency of this confinement, physicists use the plasma beta ($\beta$) parameter:

$$ \beta = \frac{p}{B^2/(2\mu_0)} $$

A higher $\beta$ indicates that the plasma can sustain a higher thermal pressure relative to the magnetic pressure, signifying more efficient magnetic use. However, excessively high $\beta$ values often trigger MHD instabilities, which can degrade confinement.

Case Study: One-Dimensional Cylindrical Equilibrium

To illustrate these principles in a concrete scenario, consider a simple cylindrical plasma column (resembling a Z-pinch). We assume cylindrical symmetry, where all physical quantities depend only on the radial coordinate $r$ ($\partial / \partial \theta = 0$ and $\partial / \partial z = 0$).

Let the magnetic field consist of an azimuthal component $B_\theta(r)$ and an axial component $B_z(r)$, such that $\mathbf{B} = B_\theta \hat{\theta} + B_z \hat{z}$. Due to the symmetry, the magnetic field lines are not curved in a way that produces a radial tension component, meaning the radial projection of $(\mathbf{B} \cdot \nabla)\mathbf{B}$ is zero.

In this case, the equilibrium equation simplifies to a radial balance between the pressure gradients:

$$ \frac{dp}{dr} + \frac{d}{dr}\left( \frac{B_\theta^2 + B_z^2}{2\mu_0} \right) = 0 $$

By integrating this equation from the center to the plasma boundary $a$ (assuming $p(a) = 0$), and considering the external magnetic field, we find that the average plasma pressure $\langle p \rangle$ is supported by two mechanisms:

  • The azimuthal magnetic pressure (the "pinch effect") caused by the plasma current.
  • The difference between the internal and external axial magnetic pressure (the "magnetic plug effect").

Conclusion

The derivation of the MHD equilibrium equations unveils the mechanical essence of magnetic confinement: the thermal pressure gradient of the plasma must be perfectly counteracted by the gradients of magnetic pressure and the force of magnetic tension. By decomposing the Lorentz force, we gain a clear physical understanding of how magnetic fields can manipulate and contain high-energy fluids. Whether applied to simple cylindrical models or complex toroidal geometries like Tokamaks, these equations remain the indispensable theoretical cornerstone for plasma design, parameter optimization, and stability analysis in fusion research.