advanced35 min studyLesson 299 in Physics

Magnetohydrodynamic equilibria

Balance pressure against magnetic pressure and tension, calculate a slab current, and separate equilibrium from stability.

On this page

What you’ll understand

  • Derive static scalar-pressure force balance.
  • Compute pressure and current for a one-dimensional magnetic slab.
  • Explain why force balance does not establish stability or determine all current components.

A plasma equilibrium is a spatial arrangement in which pressure gradients and electromagnetic forces balance under specified external constraints. It is not simply a large magnetic field surrounding a gas. Magnetic stress depends on field geometry as well as magnitude, and the current required by that geometry must satisfy Ampère's law. Even a mathematically consistent force balance leaves a separate question: whether perturbations grow.

Static balance constrains directions

Set bulk acceleration to zero in the scalar-pressure momentum equation and neglect gravity, viscosity, and electrostatic charge-density force at the retained quasineutral order. The result is ∇p=j×B\nabla p=\mathbf j\times\mathbf B. Dotting with the magnetic field shows that pressure is constant along each field line in this model. It need not be constant between distinct lines. Dotting with current likewise constrains pressure variation along current paths where the fields are regular.

The force equation fixes only the component of current perpendicular to the field. A parallel current supplies no magnetic force locally and must instead be determined consistently by field geometry, current continuity, and boundaries. A force-free configuration has current parallel to magnetic field and can have nonzero current at constant scalar pressure. Thus finding a nonzero current is not enough to infer a pressure gradient, and finding constant pressure is not enough to infer a vacuum field.

∇p=j×B,j=∇×Bμ0,B⋅∇p=0\nabla p=\mathbf j\times\mathbf B,\qquad \mathbf j=\frac{\nabla\times\mathbf B}{\mu_0},\qquad \mathbf B\cdot\nabla p=0
Pressure p is Pa, current j is A/m², magnetic field B is T, and μ₀ is vacuum permeability. These equations describe static nonrelativistic scalar-pressure MHD without gravity or flow. Anisotropic pressure, rotation, and gravity change the equilibrium equation and its along-field conclusion.↔ Scroll wide equations horizontally.

Magnetic stress includes curvature

Using the magnetic divergence constraint, the Lorentz force separates into a magnetic-pressure gradient and a field-line tension term. The pressure part pushes away from larger field magnitude; the tension part depends on variation of the field vector along itself. A straight field varying transversely can balance gas pressure through magnetic pressure alone. Curved fields can supply another force component, so replacing magnetic stress everywhere by one scalar pressure discards relevant geometry.

For the special slab B=Bz(x)z^\mathbf B=B_z(x)\hat{\mathbf z}, field-line tension vanishes and Ampère's law gives jy=−dBz/dx /μ0j_y=-dB_z/dx\,/\mu_0. Force balance then integrates to constant gas plus magnetic pressure across x. An increasing positive axial field requires negative y current and decreasing gas pressure. These signs are worth checking directly with the cross product rather than trusting an unsigned magnetic-pressure estimate. The slab is a genuine divergence-free geometry despite its transverse variation.

j×B=−∇ ⁣(B22μ0)+(B⋅∇)Bμ0,p(x)+Bz(x)22μ0=Ptot\mathbf j\times\mathbf B=-\nabla\!\left(\frac{B^2}{2\mu_0}\right)+\frac{(\mathbf B\cdot\nabla)\mathbf B}{\mu_0},\qquad p(x)+\frac{B_z(x)^2}{2\mu_0}=P_{\rm tot}
The decomposition assumes ∇·B = 0. The final scalar relation is restricted to the straight one-dimensional axial slab, with constant total transverse pressure P_tot in Pa. Magnetic pressure has units Pa, while magnetic tension is represented by the directional-derivative term.↔ Scroll wide equations horizontally.

Boundaries and stability remain independent work

The constant total slab pressure is fixed by boundary data, not by the differential force equation alone. A proposed magnetic profile can demand negative gas pressure somewhere, revealing incompatible boundary values rather than an exotic physical equilibrium. At an interface a field jump can require a surface current, and the normal stress must match the external medium. A vacuum region has no gas pressure but still exerts magnetic stress through its field.

Plasma beta, β=2μ0p/B2\beta=2\mu_0p/B^2, compares gas pressure with magnetic energy density locally. It does not by itself quantify curvature forces, solve the boundary problem, or prove stability. A pinch can satisfy static force balance while admitting kink or interchange perturbations. Stability requires an additional dynamical or energy analysis with specified admissible displacements. Resistive evolution is another separate timescale: a nearly stationary force balance can slowly change even when it is stable to ideal perturbations.

Find the current and pressure of a slab

Across a 0.20 m slab, a positive axial field rises linearly from 0.10 T to 0.20 T. At the weak-field face gas pressure is 20,000 Pa. Use μ0=1.25×10−6 H/m\mu_0=1.25\times10^{-6}\,\mathrm{H/m} for easy arithmetic. Find current density and gas pressure at the strong-field face in static slab balance.

  1. The field slope is (0.20−0.10)/0.20=0.50 T/m(0.20-0.10)/0.20=0.50\,\mathrm{T/m}. Ampère's law gives jy=−0.50/(1.25×10−6)=−4.0×105 A/m2j_y=-0.50/(1.25\times10^{-6})=-4.0\times10^5\,\mathrm{A/m^2}, with the sign fixed by the curl of an axial field varying along x.

  2. At the first face magnetic pressure is 0.102/(2μ0)=4000 Pa0.10^2/(2\mu_0)=4000\,\mathrm{Pa}. Adding its gas pressure fixes Ptot=24,000 PaP_{\rm tot}=24{,}000\,\mathrm{Pa}. This boundary value supplies the integration constant.

  3. At the second face magnetic pressure is 0.202/(2μ0)=16,000 Pa0.20^2/(2\mu_0)=16{,}000\,\mathrm{Pa}. Subtraction gives gas pressure 24,000−16,000=8000 Pa24{,}000-16{,}000=8000\,\mathrm{Pa}. The nonlinear squared-field term means pressure is not linear in x even though the field is.

The result

The slab requires a uniform −400,000 A/m² y current and pressure 8000 Pa at the stronger-field face. This satisfies static transverse force balance; it does not establish stability or specify how the current is maintained.

Explore the model

Interactive lab

Interactive lab · function

Read a pressure gradient from a quadratic

Treat y = ax² as a normalized magnetic-pressure contribution on x ≥ 0. Reset and set both controls to a = 1 and x = 0.5; predict slope 1. Reset and set both controls to a = 1 and x = 1; predict slope 2. The required gas-pressure slope has the opposite sign. The lab calculates a quadratic and tangent, not a magnetostatic equilibrium or current density.

Read a pressure gradient from a quadraticThis lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Current readings: y at chosen x: 0.608; Tangent slope · 2ax: 1.35; Signed area · ax³/3: 0.182. All parameters are adjustable with the labeled controls below.−2.4−1.201.22.4−9−4.504.59xy(x, ax²)
y = ax²TangentSigned integral
y at chosen x
0.608
Tangent slope · 2ax
1.35
Signed area · ax³/3
0.182
t = 0 s

This lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forA field magnitude linear in position has a magnetic pressure quadratic in position; matching total pressure requires a compensating gas-pressure gradient.

Check your understanding

What follows from static scalar-pressure MHD force balance when gravity and flow are neglected?

Try it yourself

A proposed slab has the same total transverse pressure 24,000 Pa24{,}000\,\mathrm{Pa} as the worked example, but its axial field reaches 0.25 T. Using the same permeability, calculate the required gas pressure and the largest field compatible with nonnegative gas pressure. Explain what the result diagnoses.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Set gas pressure equal to total pressure minus magnetic pressure; its lower physical limit is zero.

Compare with the solution

At 0.25 T, magnetic pressure is 0.252/(2.50×10−6)=25,000 Pa0.25^2/(2.50\times10^{-6})=25{,}000\,\mathrm{Pa}, requiring gas pressure −1000 Pa. This is incompatible with the stated scalar gas model and boundary total pressure. The largest allowed field is 2μ0Ptot=0.060≃0.245 T\sqrt{2\mu_0P_{\rm tot}}=\sqrt{0.060}\simeq0.245\,\mathrm T. The contradiction diagnoses an inadmissible proposed profile with those boundary data; it does not prove instability of a valid equilibrium or a universal upper bound on magnetic field in other geometries.

Take it with you

  • Static force balance constrains pressure variation along field lines but permits parallel current.
  • Magnetic pressure and magnetic tension must be distinguished.
  • Boundary compatibility, stability, and resistive evolution are separate equilibrium questions.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.