Home / Research / Aerospace / MHD Modeling
Technical Deep Dive  ·  Plasma Physics  ·  USS Olympias

Magnetohydrodynamics
Modeling Hierarchy

From Ideal MHD to Two-Fluid — Complete Derivations for Plasma-Window Firewall Analysis  ·  Hypha Labs LLC

Ideal MHD Resistive MHD Hall MHD Two-Fluid KaTeX Rendered
Section 01
§1

Single-Fluid MHD Approximation

Derivation from Two-Fluid Equations

Single-fluid MHD emerges from the two-fluid (ion + electron) description under three approximations. Starting from separate ion and electron fluid equations, we apply: (1) quasi-neutrality — charge separation on scales much larger than the Debye length is negligible; (2) center-of-mass velocity — since \(m_i \gg m_e\), the bulk velocity \(\mathbf{v} \approx \mathbf{v}_i\); and (3) electron inertia drop — the electron momentum equation has \(m_e \to 0\), converting it into a constraint (generalized Ohm's law) rather than a dynamic equation.

1 Quasi-neutrality: \(n_e \approx n_i \equiv n\). Net charge density \(\rho_q = e(n_i - n_e) \approx 0\) on MHD scales \(L \gg \lambda_D\).
2 Center-of-mass frame: \(\mathbf{v} = (m_i \mathbf{v}_i + m_e \mathbf{v}_e)/(m_i + m_e) \approx \mathbf{v}_i\). Mass density \(\rho = n(m_i + m_e) \approx n m_i\).
3 Drop electron inertia: Set \(m_e \to 0\) in the electron momentum equation. This gives Ohm's law: \(\mathbf{E} + \mathbf{v} \times \mathbf{B} = 0\) (ideal case) or \(= \eta \mathbf{j}\) (resistive case).

Resulting Equation System

Continuity

Mass conservation
$$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0$$

Momentum — Full Derivation

Start from the Euler momentum equation for a conducting fluid element, adding the Lorentz body force \(\mathbf{f} = \mathbf{j} \times \mathbf{B}\):

$$\rho \left(\frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla)\mathbf{v}\right) = -\nabla p + \mathbf{j} \times \mathbf{B}$$

Substitute Ampère's law (quasi-static, dropping displacement current \(\partial_t \mathbf{E} \approx 0\) at low frequencies \(\omega \ll \omega_{pe}\)):

Ampère's law (quasi-static)
$$\mathbf{j} = \frac{1}{\mu_0} \nabla \times \mathbf{B}$$

The Lorentz force becomes:

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

Apply the vector identity \((\nabla \times \mathbf{B}) \times \mathbf{B} = (\mathbf{B} \cdot \nabla)\mathbf{B} - \nabla(B^2/2)\):

$$\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 back, the complete ideal MHD momentum equation is:

MHD Momentum Equation
$$\rho \frac{D\mathbf{v}}{Dt} = -\nabla\!\left(p + \frac{B^2}{2\mu_0}\right) + \frac{1}{\mu_0}(\mathbf{B} \cdot \nabla)\mathbf{B}$$
Physical Interpretation

The term \(-\nabla(B^2/2\mu_0)\) is magnetic pressure — it pushes fluid from high-\(B\) to low-\(B\) regions. The term \((1/\mu_0)(\mathbf{B} \cdot \nabla)\mathbf{B}\) is magnetic tension — it acts like a restoring force along curved field lines. In the plasma windows, magnetic pressure confines the FW-1/FW-2/FW-3 plasma columns against the detonation impulse; magnetic tension maintains the field topology during the microsecond blast cycle.

Adiabatic Energy

Adiabatic / isentropic closure
$$\frac{D}{Dt}\!\left(\frac{p}{\rho^\gamma}\right) = 0 \quad \Longleftrightarrow \quad p \propto \rho^\gamma$$

Where \(\gamma = 5/3\) for a monatomic ideal gas. This is appropriate for fast-timescale phenomena where heat conduction is negligible — exactly the detonation-cycle regime.

Induction Equation

From Faraday's law \(\partial_t \mathbf{B} = -\nabla \times \mathbf{E}\) and the ideal Ohm's law \(\mathbf{E} = -\mathbf{v} \times \mathbf{B}\):

Ideal MHD Induction
$$\frac{\partial \mathbf{B}}{\partial t} = \nabla \times (\mathbf{v} \times \mathbf{B})$$

This encodes flux-freezing: magnetic field lines are advected with the fluid. In the plasma window, it means the magnetic bottle topology is maintained as long as resistivity is negligible \((R_m \gg 1)\).

Solenoidal Constraint

$$\nabla \cdot \mathbf{B} = 0$$

Not a dynamic equation — an initial condition that is preserved by the induction equation. In numerical MHD codes, divergence-cleaning schemes (e.g., Powell's 8-wave, constrained transport) enforce this at machine precision.

Validity and Limitations

Single-fluid ideal MHD is valid when: the plasma is collision-dominated (\(\nu_{ei} \gg \omega\)); length scales satisfy \(L \gg d_i, r_{Li}\) (ion inertial length, ion Larmor radius); and time scales satisfy \(t \gg \Omega_i^{-1}\) (inverse ion cyclotron frequency). For the USS Olympias plasma windows at steady-state operation, \(R_m \sim 10^3\text{–}10^4\) and single-fluid MHD is appropriate. During microsecond detonation transients, Hall and two-fluid effects become significant — that is where higher-fidelity models are required.

Section 02
§2

Resistive MHD Extension

Ideal MHD assumes perfect conductivity. In practice, plasma has finite resistivity \(\eta\) (in \(\Omega \cdot\text{m}\)) that allows field lines to diffuse through the fluid. This is not a flaw in the plasma window design — it is the regeneration mechanism.

Derivation of the Resistive Induction Equation

Start from Faraday's law:

$$\frac{\partial \mathbf{B}}{\partial t} = -\nabla \times \mathbf{E}$$

Generalized Ohm's law with finite resistivity \(\eta\):

Resistive Ohm's Law
$$\mathbf{E} + \mathbf{v} \times \mathbf{B} = \eta \mathbf{j}$$

Express \(\mathbf{j}\) via Ampère and substitute \(\mathbf{E}\):

$$\mathbf{E} = \eta \mathbf{j} - \mathbf{v} \times \mathbf{B} = \frac{\eta}{\mu_0}\nabla \times \mathbf{B} - \mathbf{v} \times \mathbf{B}$$

Insert into Faraday's law:

$$\frac{\partial \mathbf{B}}{\partial t} = -\nabla \times \!\left(\frac{\eta}{\mu_0}\nabla \times \mathbf{B} - \mathbf{v} \times \mathbf{B}\right) = \nabla \times (\mathbf{v} \times \mathbf{B}) - \frac{\eta}{\mu_0}\nabla \times (\nabla \times \mathbf{B})$$

Use the double-curl identity \(\nabla \times (\nabla \times \mathbf{B}) = \nabla(\nabla \cdot \mathbf{B}) - \nabla^2 \mathbf{B} = -\nabla^2 \mathbf{B}\) (since \(\nabla \cdot \mathbf{B} = 0\)):

Resistive MHD Induction
$$\boxed{\frac{\partial \mathbf{B}}{\partial t} = \nabla \times (\mathbf{v} \times \mathbf{B}) + \frac{\eta}{\mu_0}\nabla^2 \mathbf{B}}$$

Defining the magnetic diffusivity \(\eta_m = \eta/\mu_0\) (in \(\text{m}^2/\text{s}\)):

$$\frac{\partial \mathbf{B}}{\partial t} = \underbrace{\nabla \times (\mathbf{v} \times \mathbf{B})}_{\text{advection}} + \underbrace{\eta_m \nabla^2 \mathbf{B}}_{\text{diffusion}}$$

Magnetic Reynolds Number

The ratio of advection to diffusion defines the magnetic Reynolds number:

Magnetic Reynolds Number
$$R_m = \frac{Lv}{\eta_m} = \frac{\mu_0 L v}{\eta}$$

where \(L\) is the characteristic length scale and \(v\) the characteristic velocity. When \(R_m \gg 1\): advection dominates, flux-freezing holds approximately, and the field topology is maintained. When \(R_m \sim 1\): diffusion is significant, reconnection can occur, and the field can leak through the plasma. For the plasma windows during steady-state magnetized operation, \(R_m \sim 10^3\), confirming robust magnetic confinement. During detonation transients when plasma is momentarily disrupted, \(R_m\) drops toward unity — this is precisely when resistive effects drive field reconnection and re-threading.

Spitzer Resistivity

The classical (Spitzer) resistivity of a fully ionized plasma is:

Spitzer Resistivity
$$\eta_S = \frac{\pi Z e^2 m_e^{1/2} \ln\Lambda}{(4\pi\epsilon_0)^2 (k_B T_e)^{3/2}}$$

The key scaling is \(\eta_S \propto T_e^{-3/2}\): hotter plasma is a better conductor. A plasma at \(T_e = 10\text{ eV}\) has roughly \(10^3\times\) higher resistivity than one at \(T_e = 1\text{ keV}\). The Coulomb logarithm \(\ln\Lambda \approx 10\text{–}20\) for typical plasma parameters.

Critical Insight — Regeneration Mechanism

Resistivity is not a failure mode for the plasma windows — it is the self-healing mechanism. After each detonation pulse disrupts the plasma column, Joule heating \(\mathbf{j} \cdot \mathbf{E} = \eta j^2\) rapidly re-heats the plasma. As \(T_e\) rises, \(\eta_S \propto T_e^{-3/2}\) falls, the plasma becomes nearly ideal, flux-freezing re-establishes, and the magnetic bottle reforms — all within the microsecond recovery window between detonation events.

Key Physical Effects

▸ Magnetic diffusion: Field diffuses through the plasma on the resistive timescale \(\tau_\eta = L^2/\eta_m = \mu_0 L^2/\eta\). For a plasma window of diameter \(L \sim 0.5\text{ m}\) and Spitzer resistivity at \(T_e = 10\text{ eV}\), \(\tau_\eta \sim 1\text{ ms}\) — long compared to the \(\sim 100\text{ μs}\) detonation cycle.
▸ Magnetic reconnection: Where field lines of opposite polarity are forced together (e.g., at plasma boundary), resistivity enables topology change. Reconnection releases stored magnetic energy as heat and bulk kinetic energy — a secondary plasma heating contribution.
▸ Joule heating: Ohmic dissipation at rate \(Q_J = \mathbf{j} \cdot \mathbf{E} = \eta j^2 = \eta |\nabla \times \mathbf{B}|^2 / \mu_0^2\) continuously feeds energy back into the plasma. This sustains \(T_e\) against radiation and conduction losses.
Section 03
§3

Hall MHD Extension

Hall MHD adds one more term to the generalized Ohm's law: the Hall term \((\mathbf{j} \times \mathbf{B})/(ne)\). This term becomes important when the current layer width \(\delta\) approaches the ion inertial length \(d_i\), causing ion and electron motion to decouple.

Generalized Ohm's Law with Hall Term

Hall MHD Ohm's Law
$$\mathbf{E} + \mathbf{v} \times \mathbf{B} = \eta\mathbf{j} + \frac{\mathbf{j} \times \mathbf{B}}{ne}$$

The Hall term \(\mathbf{j} \times \mathbf{B} / (ne)\) has units of electric field and represents the difference between ion and electron drift velocities in the presence of a magnetic force.

Hall MHD Induction Equation

Substituting the generalized Ohm's law into Faraday's law:

$$\frac{\partial \mathbf{B}}{\partial t} = -\nabla \times \mathbf{E} = \nabla \times \!\left(\mathbf{v} \times \mathbf{B} - \eta\mathbf{j} - \frac{\mathbf{j} \times \mathbf{B}}{ne}\right)$$

Using \(\mathbf{j} = (\nabla \times \mathbf{B})/\mu_0\):

Hall MHD Induction
$$\frac{\partial \mathbf{B}}{\partial t} = \nabla \times (\mathbf{v} \times \mathbf{B}) - \nabla \times \!\left(\frac{(\nabla \times \mathbf{B}) \times \mathbf{B}}{\mu_0 ne}\right) + \eta_m \nabla^2 \mathbf{B}$$

The middle term is the Hall correction. It introduces whistler-wave dispersion and allows fast reconnection at scales \(\sim d_i\).

Alternative: Electron Velocity Form

The Hall term can be written more elegantly by noting that the magnetic field is frozen to the electron fluid (not the bulk), with electron velocity:

$$\mathbf{v}_e = \mathbf{v} - \frac{\mathbf{j}}{ne} = \mathbf{v} - \frac{\nabla \times \mathbf{B}}{\mu_0 ne}$$

The induction equation then reads:

$$\frac{\partial \mathbf{B}}{\partial t} = \nabla \times (\mathbf{v}_e \times \mathbf{B}) + \eta_m \nabla^2 \mathbf{B}$$

This is physically transparent: the field is advected by electrons, not ions. At scales below \(d_i\), electrons and ions decouple — ions are effectively unmagnetized while electrons remain tied to field lines.

Ion Inertial Length and Ordering

The key Hall parameter is the ion inertial length:

Ion Inertial Length
$$d_i = \frac{c}{\omega_{pi}} = c\sqrt{\frac{\epsilon_0 m_i}{n e^2}}$$

where \(\omega_{pi}\) is the ion plasma frequency. Hall effects are significant when the current layer width \(\delta \lesssim d_i\). For a hydrogen plasma at number density \(n = 10^{20}\text{ m}^{-3}\), \(d_i \approx 0.7\text{ m}\). For the compressed plasma in the FW-1/2/3 windows at typical operating density, \(d_i \sim 0.1\text{–}0.5\text{ m}\) — comparable to the plasma window diameter.

Relevance to Detonation Transients

During the microsecond blast cycle, plasma density drops sharply as the detonation products expand. As \(n\) decreases, \(d_i \propto n^{-1/2}\) grows — pushing the system from the resistive-MHD regime into Hall MHD territory. Whistler waves propagating at \(v_w \sim (k d_i) v_A\) (with Alfvén speed \(v_A = B/\sqrt{\mu_0 \rho}\)) provide a fast-reconnection pathway that helps re-establish the magnetic topology before the next detonation event.

New Wave Physics — Whistler Dispersion

Whistler wave dispersion
$$\omega_W = \frac{k^2 d_i v_A}{1 + k^2 d_i^2} \approx k^2 d_i v_A \quad (k d_i \ll 1)$$

Whistler waves carry magnetic helicity and can drive fast reconnection, explain electron jet observations, and establish a Hall-dominated current sheet structure that connects to the larger-scale MHD solution.

Section 04
§4

Two-Fluid MHD (Complete System)

The most complete fluid description treats ions and electrons as two interpenetrating, charge-carrying fluids coupled through electromagnetic fields and collisions. All lower-fidelity models (Hall, resistive, ideal MHD) are limiting cases of this system.

§4a — Continuity Equations

Control Volume Derivation

Consider an arbitrary control volume \(V\) bounded by surface \(S\) with outward normal \(\hat{n}\). The rate of change of particle number in \(V\) equals the flux through \(S\):

$$\frac{d}{dt}\int_V n_s \, dV = -\oint_S n_s \mathbf{v}_s \cdot d\mathbf{S}$$

Apply the divergence theorem to the right-hand side:

$$\int_V \frac{\partial n_s}{\partial t} \, dV = -\int_V \nabla \cdot (n_s \mathbf{v}_s) \, dV$$

Since \(V\) is arbitrary, the integrands must be equal pointwise, yielding the continuity equations for species \(s \in \{i, e\}\):

Ion Continuity
$$\frac{\partial n_i}{\partial t} + \nabla \cdot (n_i \mathbf{v}_i) = 0$$
Electron Continuity
$$\frac{\partial n_e}{\partial t} + \nabla \cdot (n_e \mathbf{v}_e) = 0$$

Under quasi-neutrality \(n_e \approx n_i \equiv n\), subtracting the two equations gives the current constraint:

Current continuity
$$\nabla \cdot \mathbf{j} = 0, \qquad \mathbf{j} = en(\mathbf{v}_i - \mathbf{v}_e)$$

In Lagrangian (material derivative) form, following an ion fluid parcel:

$$\frac{Dn_i}{Dt} + n_i \nabla \cdot \mathbf{v}_i = 0, \qquad \frac{D}{Dt} \equiv \frac{\partial}{\partial t} + \mathbf{v}_i \cdot \nabla$$

§4b — Momentum Equations

Derivation from Newton's Second Law

For a fluid element of volume \(\delta V\) containing \(n_s \delta V\) particles of species \(s\), Newton's second law gives:

$$m_s n_s \delta V \frac{D\mathbf{v}_s}{Dt} = \underbrace{q_s n_s \delta V (\mathbf{E} + \mathbf{v}_s \times \mathbf{B})}_{\text{Lorentz force}} - \underbrace{\nabla p_s \, \delta V}_{\text{pressure gradient}} + \underbrace{\mathbf{R}_{s} \, \delta V}_{\text{collisions}}$$

Taking the limit \(\delta V \to 0\) and dividing by \(\delta V\):

Ion Momentum
$$m_i n\!\left(\frac{\partial \mathbf{v}_i}{\partial t} + (\mathbf{v}_i \cdot \nabla)\mathbf{v}_i\right) = en(\mathbf{E} + \mathbf{v}_i \times \mathbf{B}) - \nabla p_i + \mathbf{R}_{ie}$$
Electron Momentum
$$m_e n\!\left(\frac{\partial \mathbf{v}_e}{\partial t} + (\mathbf{v}_e \cdot \nabla)\mathbf{v}_e\right) = -en(\mathbf{E} + \mathbf{v}_e \times \mathbf{B}) - \nabla p_e + \mathbf{R}_{ei}$$

Collision Friction Term

Electron–ion momentum exchange
$$\mathbf{R}_{ei} = -m_e n \nu_{ei}(\mathbf{v}_e - \mathbf{v}_i) = \frac{m_e n \nu_{ei}}{en}\mathbf{j} \equiv -\eta n e \mathbf{j}$$

By Newton's third law, \(\mathbf{R}_{ie} = -\mathbf{R}_{ei}\). The collision frequency \(\nu_{ei} \propto n T_e^{-3/2}\) (Spitzer).

Recovery of Generalized Ohm's Law

This is the key derivation linking two-fluid to Hall MHD. Take the electron momentum equation and set \(m_e \to 0\) (drop electron inertia):

$$0 = -en(\mathbf{E} + \mathbf{v}_e \times \mathbf{B}) - \nabla p_e + \mathbf{R}_{ei}$$

Divide by \(-en\) and substitute \(\mathbf{v}_e = \mathbf{v} - \mathbf{j}/(ne)\) and \(\mathbf{R}_{ei} = -\eta n e \mathbf{j}\):

$$\mathbf{E} + \mathbf{v}_e \times \mathbf{B} = -\frac{\nabla p_e}{en} + \frac{\mathbf{R}_{ei}}{-en} = \frac{\nabla p_e}{-en} \cdot(-1) + \eta\mathbf{j}$$

Writing \(\mathbf{v}_e \times \mathbf{B} = (\mathbf{v} - \mathbf{j}/(ne)) \times \mathbf{B} = \mathbf{v} \times \mathbf{B} - (\mathbf{j} \times \mathbf{B})/(ne)\):

Full Generalized Ohm's Law
$$\mathbf{E} + \mathbf{v} \times \mathbf{B} = \eta \mathbf{j} + \frac{\mathbf{j} \times \mathbf{B}}{ne} - \frac{\nabla p_e}{ne}$$
Hierarchy Proof

This single derivation closes the hierarchy. Dropping the electron pressure gradient \(\nabla p_e/(ne)\) recovers Hall MHD. Further dropping the Hall term \((\mathbf{j} \times \mathbf{B})/(ne)\) recovers resistive MHD. Setting \(\eta = 0\) recovers ideal MHD. Hall MHD is therefore a systematic limiting case of two-fluid, and ideal MHD is a further limit of Hall.

§4c — Energy Equations

From the first law of thermodynamics for a fluid element, tracking thermal energy per unit volume \(\mathcal{E}_s = \tfrac{3}{2} n k_B T_s\):

Ion Energy

Ion energy (Lagrangian)
$$\frac{3}{2} n k_B \frac{DT_i}{Dt} + p_i \nabla \cdot \mathbf{v}_i = Q_{ei}$$

The left-hand side is adiabatic compression heating; the right-hand side is collisional energy exchange from electrons to ions. Ions heat slowly — they gain energy only through collisional equilibration with the (hotter) electrons.

Electron Energy

Electron energy (Lagrangian)
$$\frac{3}{2} n k_B \frac{DT_e}{Dt} + p_e \nabla \cdot \mathbf{v}_e = \underbrace{\mathbf{j} \cdot \mathbf{E}}_{\text{Joule heating}} + \underbrace{Q_{ie}}_{\text{collisional}} + \nabla \cdot (\kappa_e \nabla T_e)$$

Electrons absorb Joule heating \(\mathbf{j} \cdot \mathbf{E} = \eta j^2\) directly and conduct energy along field lines (heat conductivity \(\kappa_e \propto T_e^{5/2}\)).

Collisional Energy Exchange

Collisional heat exchange
$$Q_{ei} = \frac{3 m_e n k_B (T_e - T_i)}{m_i \tau_{ei}}, \qquad Q_{ie} = -Q_{ei}$$

The equilibration timescale \(\tau_{eq} \approx (m_i/m_e)\tau_{ei}\) is much longer than the electron collision time by the mass ratio. For a hydrogen plasma at \(T_e = 100\text{ eV}\) and \(n = 10^{19}\text{ m}^{-3}\), \(\tau_{eq} \sim 1\text{ μs}\) — comparable to the detonation cycle period. This means \(T_e\) and \(T_i\) can decouple significantly during the blast transient: electrons spike hot as they absorb Joule and shock heating, while ions remain cooler. The two-temperature model is essential here; single-fluid MHD with a single temperature would miss this physics entirely.

§4d — Maxwell's Equations (System Closure)

Faraday's Law
$$\frac{\partial \mathbf{B}}{\partial t} = -\nabla \times \mathbf{E}$$
Ampère's Law (low-freq.)
$$\nabla \times \mathbf{B} = \mu_0 \mathbf{j}$$

With displacement current (high-frequency / full Maxwell):

$$\nabla \times \mathbf{B} = \mu_0 \mathbf{j} + \mu_0 \epsilon_0 \frac{\partial \mathbf{E}}{\partial t}$$

The displacement current term \(\mu_0\epsilon_0\,\partial_t\mathbf{E}\) is negligible at frequencies well below the plasma frequency \(\omega \ll \omega_{pe}\), which holds throughout the MHD regime. Gauss's law for charge is satisfied automatically by quasi-neutrality: \(\nabla \cdot \mathbf{E} \approx 0\). The system is closed by specifying initial conditions and boundary conditions on all six fluid variables \((n_i, \mathbf{v}_i, T_i, n_e, \mathbf{v}_e, T_e)\) plus the six electromagnetic field components \((\mathbf{E}, \mathbf{B})\).

Section 05
§5

Modeling Hierarchy Summary

Each level of the hierarchy is obtained by a systematic set of approximations. Dropping electron inertia from the two-fluid system yields Hall MHD; dropping the Hall term yields resistive MHD; setting resistivity to zero yields ideal MHD. The chain of approximations is:

Two-Fluid
→ \(m_e \to 0\) →
Hall MHD
→ drop \(d_i\) term →
Resistive MHD
→ \(\eta \to 0\) →
Ideal MHD
Model When to Use What It Captures Plasma Window Relevance Comp. Cost
Ideal MHDLow \(R_m \gg 1\), \(L \gg d_i\), low-freq. Alfvén waves, flux-freezing, magnetosonic modes, bulk plasma equilibrium Steady-state magnetic confinement geometry; structural equilibrium of FW-1/2/3 columns Lowest. Explicit PDE solver, \(\sim 10^6\) cells feasible
Resistive MHDLow \(R_m \sim 1\text{–}10^3\), reconnection important Magnetic diffusion, Joule heating, resistive instabilities (tearing, kink), reconnection at MHD scales Post-detonation plasma regeneration; field reconnection during low-density transient window; Joule self-heating ~2× ideal MHD. Requires resistivity profile \(\eta(T)\)
Hall MHDMedium \(\delta \lesssim d_i\), fast reconnection, \(\mu\)s transients Whistler waves, fast reconnection (\(\sim 0.1\, v_A\)), electron-ion decoupling, helicity transport Microsecond blast transient when \(n\) drops and \(d_i\) grows; fast reconnection enabling magnetic topology recovery ~5–10× resistive MHD. Stiff: requires small \(\Delta t\)
Two-FluidHigh Full kinetic-MHD coupling, separate \(T_e \neq T_i\), edge physics All of the above + \(T_e/T_i\) decoupling, electron heat conduction, full kinetic wave spectrum, ion cyclotron effects Electron spike heating during detonation; separate \(T_e/T_i\) evolution; electron heat loss to chamber walls ~20–50× resistive MHD. Requires \(\Delta t \ll \Omega_e^{-1}\)
Kinetic (PIC/Vlasov)Extreme Sub-Larmor scale physics, wave-particle resonance, non-Maxwellian distributions Landau damping, stochastic heating, ion beam instabilities, full distribution function Beyond current scope. Relevant for VASIMR-heritage RF heating of plasma window seed plasma Intractable at MHD scales. Limited to local patch simulations
Section 06
§6

Applications in Fusion Plasmas

The MHD hierarchy developed here is the same framework used in ITER, SPARC, and NIF diagnostics. This is not academic — it is directly transferable institutional knowledge. The plasma window firewall physics borrows from decades of fusion MHD research.

Magnetic Reconnection
Tokamak sawteeth (internal kink mode \(q=1\) surface), ELMs (edge localized modes), and disruptions all involve Hall-MHD reconnection. The same two-fluid equations govern the current sheet physics in plasma windows during detonation transients.
Edge / SOL Physics
Scrape-off layer blob transport, filamentary structures, and turbulent cross-field transport in tokamaks are modeled with resistive MHD + two-fluid extensions. The same physics governs plasma-gas boundary conditions at the plasma window interface.
Pellet Injection Analogy
Pellet injection into tokamak plasmas involves ablation, plasmoid formation, and rapid \(T_e\) drop around the pellet — directly analogous to the Bump Accelerator pellet entering the plasma window firewall. Two-fluid models with ablation source terms are state-of-the-art here.
RF Heating
VASIMR's RF heating (ICRH regime) heats ions while electrons remain cold — a two-temperature, two-fluid effect that ideal MHD cannot capture. The NASA VASIMR heritage in the USS Olympias design means this \(T_i \neq T_e\) physics is directly applicable.
Fast Ion Physics
Alpha particle confinement in burning plasmas, neutral beam injection, and energetic-particle-driven modes all require two-fluid or kinetic descriptions. Detonation products (fast ions, alphas from fusion if applicable) in the Bump Accelerator exhaust channel present analogous fast-ion MHD problems.
Disruption Mitigation
Massive gas injection and shattered pellet injection for disruption mitigation in ITER are studied with resistive and Hall MHD codes (e.g., M3D-C1, JOREK). The controlled plasma disruption and reformation every detonation cycle is the plasma window's analogous controlled-disruption regime.
Bottom Line

Every numerical tool built for ITER and SPARC — JOREK, M3D-C1, NIMROD, BOUT++, GENE — solves variants of the equations in Sections 1–4. The USS Olympias plasma window firewall operates in the same physical parameter space as the tokamak edge and SOL. The derivations above are not exotic: they are the standard industrial toolkit of magnetic confinement fusion, now applied to a pulsed propulsion architecture.

← Back to USS Olympias Bump Accelerator Pellet Biosynthesis All Research