Tearing Stability Prediction Combining Toroidal Calculations With a Two-Fluid Slab Layer

This paper presents a rapid and robust simulation workflow that combines the toroidal STRIDE code with the two-fluid slab layer solver SLAYER to accurately predict classical tearing mode stability and growth rates across a wide range of reactor-relevant plasma conditions.

The Gist

Imagine you are trying to keep a giant, swirling pot of super-hot soup (the plasma) stable inside a magnetic bowl (the tokamak reactor). The goal is to keep this soup boiling so hot that it creates energy, but without it spilling over or collapsing.

One of the biggest threats to this soup is a "tearing mode." Think of this like a magnetic zipper on the surface of the soup. If the zipper gets stuck or starts to rip, it opens up a hole (a magnetic island). Once that hole opens, the hot soup leaks out, the temperature drops, and the whole experiment fails. In the worst case, the whole pot crashes, damaging the machine.

This paper introduces a new, super-fast weather forecasting tool for these magnetic zippers. Here is how the authors built it, explained simply:

1. The Problem: Too Much Math, Too Slow

To predict if a zipper will rip, scientists usually have to run massive, slow computer simulations that try to model every single particle in the entire pot of soup. It's like trying to predict the weather by simulating every single raindrop in the atmosphere. It's accurate, but it takes too long to be useful for real-time safety checks.

2. The Solution: The "Split-Brain" Approach

The authors realized they didn't need to simulate the whole pot to know if a zipper is about to rip. They only needed to look at two specific parts:

• The Outer Shell (The Big Picture): This is the overall shape of the magnetic bowl. Is it round? Is it squashed? This part is handled by a code called STRIDE. Think of STRIDE as a satellite view of the storm, telling you the general wind patterns and pressure systems.
• The Inner Layer (The Tiny Crack): This is the tiny, thin slice of soup right where the zipper is trying to rip. This is where the real physics of the "tear" happens. This part is handled by a code called SLAYER. Think of SLAYER as a high-powered microscope looking at the exact fibers of the zipper to see if they are fraying.

3. How They Work Together (The "Handshake")

The magic of this new workflow is how these two codes talk to each other:

1. STRIDE looks at the big picture and says, "Hey, the magnetic pressure here is high; the zipper is under a lot of stress."
2. SLAYER zooms in on that specific stress point. It uses a simplified "slab" model (like looking at a flat piece of fabric instead of a curved balloon) to calculate exactly how the zipper fibers react. It asks: "Will this stress cause a tear, or will the soup's own rotation and heat stop it?"
3. They combine their answers to give a single "Yes/No" prediction on whether the mode will grow and how fast.

4. The "Glasser" Safety Net

The paper also adds a new safety feature called Glasser stabilization. Imagine the magnetic bowl isn't just a static container; it has a natural "curvature" that acts like a safety net. If the soup tries to tear, the curvature of the bowl can sometimes push it back together, like a spring.
The authors figured out how to calculate exactly how strong this "spring" is, even when heat flows through the soup. If the spring is strong enough, it can stop the tear even if the stress is high.

5. Why This Matters

The authors tested their new tool against:

• Mathematical theories: It matched perfectly.
• Other complex codes: It agreed with them but ran much faster.
• Realistic scenarios: They simulated a "shaped" plasma (like a DIII-D reactor) with complex currents and temperatures.

The Result: This new tool is like a fast-forward button for safety analysis.

• Old way: Wait hours to run a simulation to see if a specific plasma shape is safe.
• New way: Get the answer in seconds.

This allows engineers to quickly map out "safe zones" for future fusion reactors. They can design the startup and shutdown phases of the reactor to avoid those dangerous "tearing" zippers entirely, ensuring the machine runs smoothly and safely without crashing.

In a nutshell: They built a fast, accurate "zipper detector" that combines a satellite view of the magnetic bowl with a microscopic view of the tear, allowing us to design safer, more stable fusion reactors.

Technical Summary

Here is a detailed technical summary of the paper "Tearing Stability Prediction Combining Toroidal Calculations With a Two-Fluid Slab Layer" by Burgess et al.

1. Problem Statement

Tearing Modes (TMs) are resistive magnetohydrodynamic (MHD) instabilities in tokamaks that can lead to the formation of magnetic islands. If these islands grow large enough, they can "lock" to the vessel wall, causing plasma rotation braking and potentially leading to major disruptions.

Predicting the onset and growth of TMs is critical for designing safe discharge trajectories (ramp-up, flat-top, ramp-down) for future fusion reactors. However, existing modeling approaches face trade-offs:

• "Straight-through" codes (e.g., NIMROD, M3D-C1) solve resistive MHD across the entire plasma volume. While physically comprehensive, they are computationally expensive and often too slow for rapid scenario mapping.
• Asymptotic matching approaches (e.g., RDCON) are faster but often rely on simplified inner-layer models that may lack necessary physics (such as two-fluid effects, drifts, and thermal transport) or fail to capture complex toroidal shaping effects accurately.

There is a need for a workflow that is computationally rapid, numerically robust, and includes high-fidelity physics (two-fluid effects, toroidal shaping, and thermal conduction) to predict classical tearing stability across reactor-relevant conditions.

2. Methodology: The STRIDE+SLAYER Workflow

The authors developed a new hybrid workflow that combines a detailed inner-layer solver with a full-toroidal outer-region solver using asymptotic matching.

A. Inner-Layer Physics (SLAYER)

• Code: Uses SLAYER, which implements the Fitzpatrick-Cole-Waelbroeck four-field model in a slab geometry.
• Physics: Solves the two-fluid drift MHD equations to calculate the inner-layer response, $\Delta(\gamma, \omega)$.
• Key Effects Included:
  • Two-fluid effects (electron and ion diamagnetic drifts).
  • Finite rotation.
  • Semi-collisional and Hall effects.
  • Viscosity and thermal transport.
• Approximations: Assumes a narrow resistive layer (high Lundquist number), neglects pressure gradients and parallel flow within the slab. These are justified for $\sim$keV tokamak plasmas with low-to-moderate $\beta$.
• Solution: The model reduces to a single Ordinary Differential Equation (ODE) in Fourier space, solved rapidly via a Riccati transform to yield the normalized inner-layer stability index $\hat{\Delta}$.

B. Outer-Region Physics (STRIDE)

• Code: Uses STRIDE to calculate the toroidal ideal MHD tearing stability index, $\Delta'$.
• Physics: Computes the jump in the logarithmic derivative of the radial magnetic field perturbation across rational surfaces ($q = m/n$).
• Key Features: Includes full toroidal geometry, poloidal mode coupling, and generalized plasma shaping (elongation, triangularity).
• Output: Provides the driving force ($\Delta'$) based on the toroidal current profile gradient.

C. Glasser Stabilization and Dispersion Relation

• Glasser Stabilization: The workflow incorporates a critical threshold, $\Delta_{crit}$, which accounts for stabilization due to finite pressure and favorable magnetic curvature (Glasser-Greene-Johnson effect). This is extended to include thermal transport effects.
• Effective Stability Index: The net driving force is defined as $\Delta_{eff} = \Delta' - \Delta_{crit}$.
• Dispersion Relation: The growth rate ($\gamma$) and rotation frequency ($\omega$) are found by solving the matching condition:
  $$\Delta' - \Delta_{crit} = \Delta(\omega, \gamma)$$
  The solution involves scanning the $(\hat{\omega}, \hat{\gamma})$ space to find the intersection of real and imaginary contours.

3. Key Contributions

1. Novel Hybrid Workflow: Integration of the high-fidelity two-fluid slab model (SLAYER) with a full-toroidal ideal MHD solver (STRIDE) to create a rapid, robust tearing mode stability predictor.
2. Inclusion of Thermal Transport: Extension of the Glasser stabilization criterion ($\Delta_{crit}$) to include thermal conduction effects in a toroidal geometry, refining the stability threshold.
3. Comprehensive Benchmarking:
   • Validated against analytic linear growth rate regimes (Diffusive-Resistive, Viscous-Resistive, Semi-Collisional, Resistive-Inertial).
   • Benchmarked against the TJ code (a well-established linear stability code) across scans of plasma $\beta$ and inverse aspect ratio ($\epsilon$).
4. Application to Shaped Equilibria: Demonstrated the workflow's capability to handle realistic, shaped H-mode-like plasmas with self-consistent bootstrap currents and reverse-shear profiles.

4. Results

• Analytic Agreement: The SLAYER growth rates transition smoothly between the four distinct analytic linear regimes (DR, VR, SC, RI) and align closely with analytic limits, showing only minor offsets in magnitude.
• Code-to-Code Benchmark ($\beta$ Scan): In circular cross-section equilibria, the STRIDE+SLAYER workflow showed excellent agreement with the TJ code for both 2/1 and 3/1 tearing modes across a wide range of normalized beta ($\beta_N$).
• Inverse Aspect Ratio ($\epsilon$) Scan:
  • At low $\epsilon$ (large aspect ratio), STRIDE+SLAYER and TJ agreed well.
  • At high $\epsilon$ (low aspect ratio, $\epsilon > 0.45$), the codes diverged. STRIDE predicted higher growth rates as it approached an ideal limit, suggesting STRIDE captures higher-order geometric shaping effects that the $\epsilon$-expansion used in TJ misses.
• Glasser Stabilization: The toroidal $\Delta_{crit}$ calculated by the new workflow generally predicts less stabilization than large-aspect-ratio approximations, particularly at typical tokamak aspect ratios.
• Shaped Plasma Application: Applied to a synthetic DIII-D-like H-mode equilibrium. The workflow successfully calculated growth rates for 2/1 and 3/1 modes as the toroidal current profile gradient was varied. It demonstrated that while 2/1 modes remained unstable, 3/1 modes could be driven from stable to unstable as the current gradient increased.

5. Significance

• Operational Safety: The workflow provides a fast and reliable tool for mapping stable operational regimes and designing safe discharge trajectories (ramp-up/down) for future devices, where computational speed is essential for real-time or near-real-time analysis.
• Physics Fidelity vs. Speed: It successfully bridges the gap between slow, full-resistive MHD simulations and simplified analytic models, offering a "sweet spot" of high physics fidelity (two-fluid, toroidal shaping) with computational efficiency.
• Future Development: The authors note that while the current workflow handles decoupled rational surfaces well, future work will focus on incorporating surface-to-surface coupling and resistive wall effects to further align linear theory with experimental observations, particularly for Neoclassical Tearing Modes (NTMs).
• Open Source: The tool is open-source, facilitating broader adoption and validation within the fusion community.

In conclusion, the STRIDE+SLAYER workflow represents a significant advancement in classical tearing mode prediction, offering a robust, toroidally self-consistent method that is ready for application to both synthetic and experimental tokamak data.