Explosive eruption cycles in a rotating Z-pinch

A transonic shear flow directed along magnetic field lines can linearly stabilize a steep pressure gradient in a confined magnetohydrodynamic (MHD) plasma. In Z-pinch geometry, we show that, like the edge pedestal in tokamak devices, this transport barrier – which we call the ``MHD pedestal’’ – is metastable, i.e., unstable to finite-amplitude displacements of flux tubes. We simulate the slow formation of an MHD pedestal in a heated and sheared Z-pinch, which collapses on reaching a critical height, expelling an order-unity fraction of the confined thermal energy. The MHD pedestal then rebuilds and the process repeats, in a manner analogous to the ELM cycle seen in fusion experiments. We show that the available energy of the metastable equilibrium, and the most energetically favorable amount of ejected plasma, can be calculated from first principles via combinatorial optimization of flux-tube interchanges.

💡 Research Summary

The authors investigate a novel transport barrier that forms in a low‑β (β≪1) rotating Z‑pinch plasma when a transonic shear flow is aligned with the magnetic field. By deriving the radial force balance for a flux tube in terms of three Lagrangian invariants—specific entropy (s), specific magnetic flux (χ), and specific angular momentum (ℓ)—they show that a sufficiently steep increase of ℓ²/χ can linearly stabilize an equally steep decrease of s/χ, thereby creating a sharp pressure gradient region they term the “MHD pedestal”. The linear stability criterion (L>0) is expressed in Eq. (5) and depends on the sound speed and the azimuthal Mach number.

Although the pedestal is linearly stable, the authors prove that for adiabatic indices γ<2 the system is only metastable: large‑amplitude radial displacements eventually make the s/χ term dominate the ℓ²/χ term in the acceleration expression (Eq. 4), leading to a loss of stability. The critical pedestal height is limited by the ram pressure of the confining azimuthal flow (p < ½ ρ u_ϕ², Eq. 6).

To explore the nonlinear dynamics, they perform two‑dimensional axisymmetric simulations with the Athena++ MHD code. Two regimes are examined:

1. L‑mode – no azimuthal drag (τ_drag,ϕ → ∞). Heating at a prescribed radius drives an exchange‑turbulence that mixes the plasma until s/χ becomes spatially uniform; no pedestal forms.

2. H‑mode – finite azimuthal drag (τ_drag,ϕ ≈ τ_out). A drag term applied inside the would‑be pedestal radius enforces a transonic shear (Ma_ϕ ≈ 1) that builds a steep ℓ²/χ gradient. The pedestal grows until the pressure profile reaches the marginal‑stability condition (Eq. 12).

When the marginal state is reached, a saddle‑node bifurcation occurs: two equilibrium positions for a given flux tube coalesce and annihilate. The flux tube then erupts outward with an explosive acceleration that follows the analytic trajectory (Eqs. 7–8), forming a mushroom‑shaped plume characteristic of Rayleigh–Taylor‑like nonlinear growth. Hot plasma is expelled, cold peripheral plasma rushes inward, and the pressure profile collapses to a level well below the H‑mode prediction.

After the eruption, interchange turbulence in the core re‑mixes the remaining hot plasma with the incoming cold plasma, restoring a uniform s/χ distribution. The pedestal subsequently reforms, and the cycle repeats. The authors document multiple such cycles, showing the temporal evolution of the exchange acceleration (Fig. 4a) and the rebuilding of the pedestal between eruptions.

A key contribution is the quantification of the energy budget using the concept of Available Potential Energy (APE). APE is defined as the difference between the total magnetic, thermal, and azimuthal kinetic energy of the current state and the minimum energy attainable through any combination of (i) nonlinear flux‑tube interchanges and (ii) draining to the outer sink. In the β → 0 limit, the otherwise NP‑hard nonlinear assignment problem reduces to a linear sum‑assignment problem, enabling the authors to compute the stored energy of the metastable pedestal and the maximum energy that can be released during an eruption from first principles.

The overall picture mirrors the Edge‑Localized Mode (ELM) cycle observed in tokamaks: a linearly stable edge transport barrier (the pedestal) that is metastable, a rapid, quasi‑ballistic eruption when a stability limit is crossed, and a subsequent recovery phase. By demonstrating that the same physics arises in a simple Z‑pinch geometry, the work provides a clean testbed for studying pedestal metastability, saddle‑node collapse, and energy release mechanisms. It also suggests practical control knobs—heating power, shear flow strength, and drag timescales—that could be tuned to modify the eruption frequency and amplitude, offering potential strategies for ELM mitigation in larger fusion devices.