Magnetohydrodynamics (MHD)

Why MHD matters

• Solar physics. The corona, sunspots, prominences, flares, and coronal mass ejections are all MHD phenomena. The 11-year solar cycle is a global dynamo; the slow and fast solar wind originate from open and closed coronal field structures.
• Fusion plasmas. Tokamak and stellarator confinement physics is dominated by MHD equilibrium and stability analysis. Designing stable, high-pressure plasmas means navigating a phase diagram of MHD instabilities (kink, ballooning, tearing, neoclassical tearing, ELMs).
• Geomagnetism. Earth's magnetic field is a self-excited dynamo in the molten iron outer core. MHD also governs the magnetosphere — the field's interaction with solar wind — and drives the aurora through reconnection at the magnetopause.
• Astrophysics. Accretion disks, jets, neutron-star magnetospheres, supernova remnants, and galactic-scale magnetic fields all live in the MHD regime, often relativistic (RMHD).
• Liquid-metal engineering. Fast-reactor coolant pumps, continuous casting in steelmaking, electromagnetic stirring, and proposed ITER blanket modules use MHD effects on liquid sodium, lithium, and lead-bismuth.
• MHD propulsion. Ship propulsion via ionic-current Lorentz force in seawater (low efficiency but silent — used in submarine experiments). Hall-effect thrusters and other plasma propulsion in space draw on MHD theory.
• Industrial processes. Aluminum smelting (Hall-Héroult cell electrolyte motion is MHD), arc furnace dynamics, and electromagnetic levitation for ultra-pure casting.

The MHD equation set

1. Continuity: ∂ρ/∂t + ∇·(ρu) = 0.
2. Momentum (Navier-Stokes + Lorentz): ρ(∂u/∂t + u·∇u) = −∇p + J×B + μ∇²u + ρg.
3. Induction (from Maxwell + Ohm's law): ∂B/∂t = ∇×(u×B) + (1/μ₀σ)∇²B.
4. Solenoidal: ∇·B = 0 (no monopoles, exact at all times).
5. Energy: evolution of internal energy or pressure with Joule heating J²/σ and resistive plus viscous dissipation.
6. Equation of state: p(ρ, T) — ideal gas, polytropic, or more elaborate.

Alfvén waves and other MHD modes

• Alfvén wave. Transverse, incompressible, propagating along B at v_A = B/√(μ₀ρ). Pure magnetic-tension restoring force.
• Slow magnetosonic wave. Compressional and partly thermal-pressure-driven; speed below both sound and Alfvén speeds.
• Fast magnetosonic wave. Compressional with combined magnetic-pressure and thermal-pressure restoring force; speed above both v_s and v_A.
• Shocks. MHD admits fast, slow, and intermediate shocks; the bow shock at Earth's magnetopause is a fast MHD shock.

Frozen-in flux in three sentences

Take the curl of Ohm's law in the perfect-conductor limit (E + u×B = 0). Faraday's law then becomes ∂B/∂t = ∇×(u×B), which is exactly the Lie-derivative statement that magnetic flux through any surface co-moving with the fluid is conserved. Equivalently, the magnetic field lines move with the fluid as if attached — stretch the fluid and the field lines stretch; compress it and they bunch.

Magnetic reconnection — when frozen-in fails

Strict frozen-in requires σ → ∞. Real plasmas have finite resistivity (or kinetic effects) that allow field lines to slip through fluid in thin current sheets. When two oppositely directed field lines meet, they can break and reconnect — releasing magnetic energy as kinetic and thermal energy in plasma jets. Reconnection drives:

• Solar flares (energy release ~10²⁵ J in minutes).
• Substorms in Earth's magnetotail (auroral brightening).
• Disruptions in tokamaks (sudden loss of confinement).
• Sawtooth oscillations in tokamak cores.

Sweet-Parker (resistive) reconnection is too slow for solar flares; observations require Petschek or plasmoid-instability modes that produce reconnection rates ~0.1 v_A — a still-active research frontier.

A taxonomy of MHD instabilities

• Kink (m=1). A current-carrying flux tube becomes helically deformed when the safety factor q drops below 1.
• Sausage (m=0). Axisymmetric pinching of a flux tube; controlled by axial field stiffening.
• Tearing. Resistive instability that creates magnetic islands by reconnection at rational surfaces.
• Ballooning. Pressure-driven instability that grows on the outer (low-field) side of a tokamak.
• Edge-localized modes (ELMs). Periodic bursts at the H-mode pedestal; potentially damaging to ITER first wall.
• Magnetorotational instability (MRI). In differentially rotating disks with a weak vertical field — drives angular-momentum transport in accretion disks.

Common misconceptions

• "MHD needs full ionization." Any electrically conducting fluid works — liquid metals at 300 K, partially ionized chromospheric plasma, even sufficiently salty water. The conductivity σ controls how strongly fluid and field are coupled.
• "Instabilities are rare." The opposite: ideal MHD has a rich spectrum of instabilities (kink, sausage, tearing, ballooning, MRI), and most laboratory plasmas operate close to one or more thresholds.
• "Field and fluid are decoupled." J×B is a force on the fluid; u×B is an EMF that drives currents. Field and fluid are tightly two-way coupled — that is the entire content of MHD.
• "Frozen-in is exact." Frozen-in holds in the σ → ∞ limit only. Resistivity, Hall effect, and electron inertia all break it on small scales — and that is why magnetic reconnection happens.
• "Alfvén waves are slow." Alfvén speed depends on B/√ρ. In the corona (B ~ 10 G, ρ ~ 10⁻¹⁵ kg/m³), v_A ~ 10³ km/s; in Earth's outer core, ~10⁻⁴ m/s. Each regime has its own characteristic timescale.
• "MHD is just Maxwell + Navier-Stokes added." The tight coupling through Ohm's law and the Lorentz force makes MHD qualitatively distinct: it has different waves, different conservation laws (cross-helicity, magnetic helicity), and instabilities absent from either parent.