Runaway electron generation in ITER mitigated disruptions with improved physics models

This study utilizes an enhanced 1D simulation framework to demonstrate that while staggered or low-neon shattered pellet injections can effectively mitigate runaway electron generation in ITER L-mode scenarios, avoiding multi-megaampere beams in DT H-mode disruptions requires specific conditions like long pre-thermal quench durations and high deuterium assimilation to counteract nuclear seed effects.

The Gist

Imagine a fusion reactor like ITER as a giant, super-hot pot of soup (the plasma) that we are trying to keep boiling without it exploding. Sometimes, the soup gets too unstable and threatens to spill over or boil over violently. This is called a disruption.

When this happens, two bad things occur:

1. The Heat Spill: The soup dumps all its heat onto the walls of the pot, potentially melting them.
2. The Lightning Storm: A massive electric field forms inside the soup, accelerating a few stray electrons to near the speed of light. These become Runaway Electrons (REs). If they hit the wall, they act like a high-powered laser beam, drilling holes through the reactor.

The goal of this paper is to figure out how to stop this "lightning storm" from forming when we try to save the reactor.

The Rescue Plan: The "Shattered Pellet"

To stop the heat spill, scientists use a technique called Shattered Pellet Injection (SPI). Imagine throwing a handful of frozen snowballs (pellets made of hydrogen and neon gas) into the hot soup.

• The Goal: The snowballs melt, cool the soup down, and make it thick enough to stop the runaway electrons from speeding up.
• The Problem: In the past, our computer models for this were a bit like using a 2D map to navigate a 3D city. They missed some crucial details, leading to overly optimistic or pessimistic predictions.

The New "Super-Model"

The authors of this paper upgraded their computer simulation (called Dream) with four new "physics engines" to make it more realistic. Here is what they added, using simple analogies:

1. The Leaky Boat (Vertical Motion):

   • Old Model: Assumed the soup stayed perfectly still in the middle of the pot.
   • New Model: Realizes that when things go wrong, the soup often tilts and moves up or down (like a boat listing). As it moves, the edges of the soup scrape against the pot walls.
   • Why it helps: This "scraping" acts like a drain, letting the dangerous runaway electrons escape before they can build up into a giant beam.

2. The Drifting Snowball (Plasmoid Drift):

   • Old Model: Assumed the melted snowball stayed exactly where it was thrown.
   • New Model: Realizes that in the curved magnetic field of the reactor, the cold cloud of gas drifts sideways (like a leaf floating down a curved river).
   • Why it matters: If the gas drifts too far to the edge, it never cools the center of the soup effectively. The model now accounts for this drift to see if we actually get enough gas in the right place.

3. The Safety Valve (Hyper-Resistivity):

   • Old Model: Sometimes the computer simulation created weird, thin spikes of electricity that didn't make physical sense, causing the simulation to crash or give wrong answers.
   • New Model: Added a "safety valve" that mimics how real magnetic fields naturally smooth out these spikes. It prevents the simulation from getting stuck on fake, unphysical problems.

4. The New Wall (Updated Compton Seed):

   • Old Model: Used an old design for the reactor's inner wall (made of Beryllium).
   • New Model: Updated to the new design (made of Tungsten).
   • Why it matters: The wall material affects how much radiation bounces back into the soup and creates new runaway electrons. The new wall is slightly different, so the model had to be recalibrated.

What Did They Find?

After running thousands of simulations with this new, smarter model, they found some surprising truths:

• Timing is Everything: To stop the runaway electrons, you need to cool the soup slowly at first. If you cool it too fast, you create a "hot tail" of electrons that escape the cooling and become runaways. You need a "pre-cooling" phase to let the electrons settle down before the big crash.
• The "Goldilocks" Mixture: You need just the right amount of Neon (the cooling gas). Too much Neon cools the soup too fast, creating runaways. Too little, and the soup doesn't cool enough.
• The Nuclear Problem: In the most powerful scenarios (using Deuterium and Tritium fuel), there is a background "seed" of runaway electrons created by nuclear radiation (gamma rays).
  • The Bad News: In these nuclear scenarios, it is very hard to stop the runaways completely because this nuclear seed is so strong.
  • The Good News: If we can combine the new "leaky boat" effect (electrons escaping as the soup tilts) with the "safety valve" (smoothing out the current), we might be able to keep the runaway beam small enough to survive, even with the nuclear seed.

The Winning Strategy: The Two-Stage Injection

The paper proposes a specific recipe for saving the reactor, especially in the most dangerous nuclear scenarios:

1. Stage 1 (The Soft Landing): Throw in a pellet that is mostly Hydrogen with a tiny bit of Neon. This cools the soup gently and thickens it without triggering a sudden crash. This stops the "hot tail" electrons from forming.
2. Stage 2 (The Big Chill): Wait a few milliseconds, then throw in a pellet full of Neon. Now that the soup is thick and cool, this triggers a controlled, safe thermal quench (the "TQ") that stops the reaction safely.

The Bottom Line

The paper concludes that while stopping runaway electrons in a massive reactor like ITER is extremely difficult and requires perfect conditions, it is theoretically possible.

It's like trying to stop a tsunami with a sandbag. If the sandbag is too small, or the water moves too fast, you fail. But if you use the right mix of sand (Hydrogen/Neon), throw it at the right time (Staggered injection), and rely on the water naturally spilling over the side (Vertical motion losses), you might just be able to save the day.

The authors warn that we are walking a tightrope. If the reactor behaves slightly differently than our models predict (e.g., if the magnetic fields don't smooth out as expected), the runaway electrons could still win. But with these new, more realistic models, we now have a much better map for how to keep ITER safe.

Technical Summary

1. Problem Statement

Tokamak disruptions pose a critical threat to the reliable operation of ITER and future fusion reactors. During a disruption, the rapid release of thermal and magnetic energy induces a strong toroidal electric field, which can accelerate electrons to relativistic speeds, creating Runaway Electrons (REs). In a 15 MA ITER plasma, this can result in multi-megampere (MA) RE beams. If such a beam strikes the vessel wall, it can cause catastrophic damage.

The primary mitigation strategy for ITER is Shattered Pellet Injection (SPI), which injects cryogenic material (a mixture of neon and hydrogen) to cool the plasma and increase collisional drag before REs can form. However, predicting the efficacy of SPI in reactor-scale conditions is challenging due to complex, non-linear physics. Previous simulations often relied on simplified models that neglected key phenomena, leading to uncertainties in whether SPI can reliably suppress REs to safe levels (<150 kA).

2. Methodology

The authors utilized the Dream code, a 1D disruption simulation framework, to assess RE generation in ITER. The core innovation of this study was the integration of four improved physics models to address limitations in previous simulations:

1. Vertical Displacement Event (VDE) Scrape-off Model: A reduced model accounting for RE losses when the plasma moves vertically during the Current Quench (CQ), causing magnetic flux surfaces to intersect the vessel wall (open field lines).
2. Plasmoid Drift Model: A semi-analytical model describing the drift of ablated material (plasmoids) across magnetic field lines toward the low-field side due to charge separation in curved tokamak geometry. This corrects previous overestimations of core material deposition.
3. Adaptive Hyper-Resistive Transport: A model to suppress unphysical, narrow current channels that form in 1D fluid simulations during the CQ. This mimics the effect of 3D MHD instabilities that flatten current profiles and enhance transport.
4. Updated Compton Seed Source: An updated calculation for the RE seed generated by Compton scattering of gamma rays from the Tungsten (W) first wall (replacing the previous Beryllium model), including a time-dependent flux reduction as fusion reactions cease.

Scenarios Simulated:

• 15 MA Full Current:
  • H26: Hydrogen L-mode (non-nuclear, no nuclear seeds).
  • DTHmode24: Deuterium-Tritium H-mode (nuclear, includes tritium beta decay and Compton seeds).
• 7.5 MA Intermediate Current:
  • DD7.5MA: Deuterium H-mode (non-nuclear).

3. Key Contributions

• Physics Integration: The first comprehensive 1D study to simultaneously model vertical motion losses, plasmoid drift, hyper-resistivity, and updated wall materials for ITER SPI scenarios.
• Identification of Critical Thresholds: The study establishes that avoiding multi-MA RE beams requires a specific combination of conditions: a long pre-thermal quench (TQ) duration, high deuterium assimilation with limited neon, and a seed current comparable to a single relativistic electron.
• Scenario-Specific Insights:
  • L-mode vs. H-mode: Demonstrated that H-mode plasmas suffer from an "H-mode penalty" where plasmoid drift is stronger, hindering core deposition and making mitigation more difficult than in L-mode.
  • Nuclear Seeds: Quantified the dominance of Compton scattering seeds in DT scenarios, which often overwhelms mitigation strategies effective in non-nuclear cases.
• Proposed Mitigation Strategy: Developed a theoretical two-stage SPI strategy for DT H-mode plasmas that successfully limits RE currents even in the presence of nuclear seeds.

4. Key Results

A. Conditions for RE Suppression

Complete avoidance of multi-MA RE beams is only possible in a confined parameter space:

1. Long Pre-TQ: The time before the thermal quench must be long enough to thermalize "hot-tail" electrons.
2. Material Balance: High assimilation of hydrogen/deuterium is needed to increase collisional drag, but neon must be limited to prevent premature radiative collapse and excessive target electrons for avalanche multiplication.
3. Seed Suppression: The representative seed current must be extremely low (ideally <1 A).

B. Impact of New Models

• Plasmoid Drift: In H-mode (hot plasma), pure hydrogen pellets drift significantly outward, reducing core assimilation and triggering early TQs. This leads to large hot-tail seeds. Doping with small amounts of neon anchors the plasmoid, improving core deposition.
• Scrape-off Losses (VDE): In scenarios where REs are nearly suppressed, the inclusion of VDE scrape-off losses is often the deciding factor between a negligible residual current and a multi-MA beam. The vertical motion opens field lines, allowing REs to escape and reducing the poloidal flux available for avalanche multiplication.
• Hyper-Resistivity: Including hyper-resistive transport during the TQ flattens the current profile, reducing the poloidal flux variation ($\Delta \psi_p$). This can reduce the potential avalanche gain by up to 4 orders of magnitude ($10^{-4}$), significantly lowering the final RE current.
• CQ Radial Transport: If magnetic perturbations ($\delta B/B$) exceed a critical threshold ($\sim 4 \times 10^{-4}$), radial transport during the CQ can suppress RE currents to negligible levels, even for large seeds.

C. Scenario Outcomes

• L-mode (H26): Staggered injections (protium followed by neon) and low-neon single injections successfully suppress REs by extending the pre-TQ and ensuring deep core deposition.
• H-mode DT (DTHmode24): Standard SPI strategies generally fail to suppress REs due to strong plasmoid drift and dominant nuclear (Compton) seeds. However, a two-stage injection (Stage 1: low-Ne protium to dilute/cool; Stage 2: high-Ne to trigger TQ) combined with hyper-resistivity and VDE losses can reduce the final RE current to 0.24 MA, a theoretically viable level.
• Reduced Current (7.5 MA): RE suppression is more achievable at lower currents due to reduced avalanche gain. Staggered injection with a long TQ can fully suppress REs, tolerating seed currents up to ~0.4 A.

5. Significance

This paper provides a more realistic and conservative assessment of ITER's ability to mitigate runaway electrons.

• Feasibility: It confirms that while RE avoidance is challenging, it is theoretically possible in ITER, particularly in non-nuclear scenarios or with optimized injection strategies.
• Operational Constraints: It highlights that successful mitigation is highly sensitive to injection timing, pellet composition, and the specific plasma regime (L-mode vs. H-mode).
• Design Implications: The results suggest that relying solely on SPI may be insufficient for DT H-mode scenarios without the synergistic effects of VDEs, current profile flattening (hyper-resistivity), and potentially external magnetic perturbations.
• Future Pathways: The proposed two-stage SPI strategy offers a concrete theoretical pathway for ITER operations, emphasizing the need for precise control over pellet composition and injection timing to manage plasmoid drift and seed generation.

In conclusion, the study underscores that ITER lies close to a threshold where RE avoidance is marginally accessible, requiring the simultaneous optimization of multiple physical mechanisms to ensure machine safety.