Variability of MHD Instabilities in Benign Termination of High-Current Runaway Electron Beams in the JET and DIII-D Tokamaks

This study analyzes high-current runaway electron beam terminations in JET and DIII-D tokamaks to demonstrate that the success of benign termination is determined by the interplay between current profile peaking and MHD perturbation amplitudes, rather than ideal MHD growth rates alone.

The Gist

The Big Picture: Taming the "Runaway" Electron Storm

Imagine a fusion reactor (like the JET or DIII-D machines mentioned in the paper) as a giant, super-hot pot of soup. Inside this soup, charged particles (plasma) swirl around at incredible speeds. Sometimes, things go wrong, and the soup gets too hot or unstable. This is called a disruption.

When a disruption happens, a tiny fraction of the particles can get "super-charged" and turn into Runaway Electrons (REs). Think of these as a swarm of angry, super-fast bees that have broken out of their hive. If they hit the walls of the reactor, they can melt the metal, causing massive damage.

The Goal: Scientists want to stop these bees safely. They want to "deconfine" them—spread them out gently so they hit a large area of the wall with low energy, rather than focusing all their power on one tiny spot. This is called a "Benign Termination."

The Strategy: The "Hydrogen Bomb" (But Gentle)

To stop the runaway electrons, scientists inject a massive amount of hydrogen gas into the reactor.

• The Good Outcome (Benign): The hydrogen mixes with the plasma, cools it down, and the runaway electrons lose their energy and spread out harmlessly. It's like throwing a net over the angry bees, calming them down, and letting them land softly.
• The Bad Outcome (Non-Benign): The hydrogen fails to mix properly. The bees stay angry and focused. They slam into the wall with a concentrated beam of energy, potentially damaging the reactor.

What This Paper Discovered

The researchers looked at data from two major fusion machines: JET (in the UK) and DIII-D (in the US). They wanted to figure out why some attempts to stop the electrons worked (Benign) and others failed (Non-Benign), especially when the reactor was running at very high power.

Here are the key findings, explained simply:

1. The Shape of the "Bees" Matters (Current Peaking)

Imagine the swarm of runaway electrons as a river of water.

• Benign Case: The river is wide and spread out. It's easy to calm down.
• Non-Benign Case: The river is squeezed into a tiny, high-pressure hose. It's very focused and powerful.

The paper found that when the reactor is running at high power (High Current), the "river" tends to get squeezed into a tight, focused beam. This happens because the electrons are generated in a way that makes them cluster in the center.

• The Analogy: If you try to put out a fire with a wide spray of water, it works. If the water is forced through a tiny nozzle, it just shoots a high-pressure jet that might miss the fire or cause damage. The "tight beam" makes it hard for the hydrogen to mix in and calm the electrons down.

2. The "Safety Factor" (The Edge of the Cliff)

Scientists use a number called the Edge Safety Factor ($q_{edge}$) to measure how "twisted" the magnetic field is at the edge of the plasma. Think of this as the steepness of a cliff.

• On DIII-D: They found they could stop the electrons safely even when the cliff was very steep (low safety factor).
• On JET: When the reactor was running at very high power, they could only stop the electrons safely if the cliff was less steep (higher safety factor). If they tried to stop it when the cliff was steep, the "bees" would just crash into the wall.

3. The "Shake" (MHD Instabilities)

To stop the electrons, scientists intentionally make the magnetic field wobble (an instability). This wobble is supposed to shake the electrons loose from their tight path so they spread out.

• The Finding: In the "Benign" (successful) cases, the magnetic field shook hard and loud. It was a big, chaotic wobble that scattered the electrons everywhere.
• The Failure: In the "Non-Benign" (failed) cases, the magnetic field barely shook at all. It was a weak, quiet wobble. Because the shake was too weak, the electrons stayed in their tight beam and crashed into the wall.

Why did the shake fail? The paper suggests that because the electron beam was so tightly focused (high "peaking"), the magnetic field couldn't grab hold of it effectively to start the big wobble. It's like trying to shake a heavy, compacted snowball; it doesn't move much. But a fluffy pile of snow (a spread-out beam) shakes apart easily.

4. The "Re-ionization" Problem

There is a tricky twist. When the hydrogen is injected, it's supposed to cool the plasma. But if the electron beam is too focused and energetic, it acts like a laser. Instead of cooling down, the electrons re-ionize the hydrogen gas (they rip the electrons off the hydrogen atoms again).

• The Result: The cooling agent turns back into hot plasma, and the runaway electrons keep running. This is why high-power attempts on JET often failed—the beam was too strong for the hydrogen to handle.

The Bottom Line

The paper concludes that how the electron beam is shaped is the most important factor.

• If the beam is spread out, the magnetic "wobble" works, the hydrogen cools it down, and everything is safe.
• If the beam is too focused (which happens at high power), the wobble is too weak, the hydrogen gets re-ionized, and the beam crashes into the wall.

Why does this matter for the future?
The next big fusion reactor, ITER, will run at even higher powers than JET. If we don't figure out how to keep the electron beam spread out (or how to make the magnetic wobble stronger), we risk damaging the reactor during an emergency. This research helps scientists design better "safety nets" for future fusion power plants.

Summary in One Sentence

To safely stop a runaway electron storm in a fusion reactor, you need the storm to be spread out enough so that a magnetic "shake" can scatter it; if the storm is too tightly focused, the shake fails, and the reactor gets hit with a dangerous beam of energy.

Technical Summary

1. Problem Statement

Mitigating Runaway Electron (RE) beams is critical for future high-current fusion reactors (e.g., ITER) to prevent catastrophic thermal and mechanical loads on vessel components. The "benign termination" strategy involves injecting hydrogenic gas into the RE beam to recombine the companion plasma, triggering a Magnetohydrodynamic (MHD) instability that deconfines REs uniformly over a large surface area.

While successful on devices like DIII-D, recent experiments on the Joint European Torus (JET) at high pre-disruptive currents ($I_p \ge 2.5$ MA) have revealed significant challenges. Many high-current JET discharges result in non-benign terminations, where REs strike the wall locally, causing damage. The paper addresses the lack of a conclusive explanation for why high-current JET discharges fail to achieve benign termination while DIII-D discharges often succeed, and seeks to unify the phenomenology across different tokamaks.

2. Methodology

The study employs a comparative empirical analysis of approximately 40 JET discharges (2019–2023) and 20 DIII-D discharges (2018, 2021, 2022).

• Selection Criteria: Only discharges with pure hydrogenic injections (protium/deuterium) and impurity fractions $<5\%$ were analyzed to isolate the effects of recombination and MHD dynamics.
• Diagnostics:
  • Fast Magnetic Sensors: Mirnov coils were used to measure time-varying magnetic fields ($\dot{B}$) for mode decomposition ($n=1, 2$), growth rate extraction ($\gamma$), and perturbation amplitude ($\delta B$).
  • Equilibrium Reconstruction: EFIT++ (JET) and EFIT (DIII-D) codes solved the Grad-Shafranov equation to reconstruct plasma boundaries, safety factor profiles ($q$), and internal inductance ($l_i$) as a proxy for current profile peaking.
• Modeling: Linear resistive MHD simulations were performed using the CASTOR3D code to map stability boundaries in the $l_i$ vs. $q_{edge}$ parameter space.

3. Key Contributions

• High-Current JET Analysis: First systematic analysis of benign termination limits at $I_p \approx 3$ MA on JET, identifying a distinct failure mode compared to lower-current regimes.
• Unification of Phenomenology: A comparative framework linking JET and DIII-D behaviors, resolving discrepancies regarding edge safety factor ($q_{edge}$) thresholds for benign termination.
• Role of Current Peaking: Identification of internal inductance ($l_i$) and current profile peaking as the dominant ordering parameters determining the MHD stability boundary, rather than just the edge safety factor alone.
• Distinction of Mechanisms: Differentiation between ideal MHD timescales (Alfvénic) and resistive dynamics in governing the termination process.

4. Key Results

A. Divergent Behavior at High Currents

• JET: Non-benign terminations at high $I_p$ ($\ge 2.5$ MA) occur at low edge safety factors ($q_{edge} \approx 2$). These events are preceded by intermittent, non-terminating MHD activity at higher rational $q$ values (e.g., $q=3, 4$).
• DIII-D: Displays a broader range of terminating $q_{edge}$. Benign terminations often occur at low $q_{edge} \approx 2$ (via center-post compression), while non-benign cases tend to occur at slightly higher $q_{edge}$.

B. The Role of Current Profile Peaking ($l_i$)

• JET: Non-benign cases exhibit higher internal inductance ($l_i$), indicating more peaked RE current profiles. This peaking is linked to the high pre-disruptive current and strong avalanche generation.
  • Mechanism: Highly peaked profiles may hinder the recombination of the companion plasma or cause re-ionization of the background plasma due to high local current densities. This prevents the formation of a stable, low-density plasma required for ideal MHD deconfinement.
• DIII-D: Shows a wider distribution of $l_i$, with benign cases generally having lower $l_i$ (flatter profiles).

C. MHD Perturbation Amplitude ($\delta B$) as the Discriminator

• Growth Rates: The measured exponential growth rates ($\gamma$) of the terminating $n=1$ modes are similar for both benign and non-benign cases. This suggests that ideal MHD timescales (dependent on electron density) alone cannot explain the difference in termination outcomes.
• Perturbation Amplitude: The decisive factor is the amplitude of the magnetic perturbation ($\delta B$).
  • Benign cases: Exhibit large $\delta B$ amplitudes, sufficient to drive stochasticity and widespread RE deconfinement.
  • Non-benign cases: Exhibit suppressed $\delta B$ amplitudes, failing to trigger a catastrophic, global deconfinement event.
• Interpretation: The suppression of $\delta B$ in non-benign cases is likely due to the interplay between ideal and resistive dynamics. Peaked current profiles (high $l_i$) may favor internal kink modes or resistive effects that do not efficiently stochastize the field, whereas flatter profiles allow for stronger external kink modes that drive global stochasticity.

D. Modeling Validation

• CASTOR3D modeling confirmed that the observed stability boundaries align with kink mode theory.
  • JET (Low $l_i$): Stability is governed by external kink modes.
  • DIII-D (High $l_i$): Stability is influenced by resistive internal kink modes.
• The modeling supports the hypothesis that the specific RE current profile formation (driven by machine-specific avalanche physics) dictates which stability boundary is encountered.

5. Significance and Implications

• Reactor Extrapolation: The findings suggest that simply scaling up current is insufficient for benign termination in future reactors like ITER. The current profile evolution during the RE plateau phase is critical.
• Operational Strategy: To achieve benign termination at high currents, operators must actively steer the RE current profile to be less peaked (lower $l_i$). This could prevent re-ionization and ensure the plasma remains in a regime where large-amplitude MHD instabilities can trigger uniform deconfinement.
• Physics Understanding: The study clarifies that the transition from benign to non-benign is not solely a function of density or growth rate, but a complex interplay between current profile peaking, re-ionization physics, and the amplitude of magnetic perturbations required to bridge ideal and resistive MHD regimes.

Conclusion: The paper establishes that the success of benign RE termination depends on maintaining a current profile that allows for large-amplitude MHD perturbations. High-current JET failures are attributed to peaked profiles causing re-ionization and suppressing these necessary perturbations, a challenge that must be addressed through profile control in future fusion devices.