Simulations of internal kink modes and sawtooth crashes for SPARC baseline-like scenarios using the M3D-C1 code

Using the high-fidelity M3D-C1 code, this study investigates low-n MHD instabilities in SPARC baseline scenarios, revealing that internal kink modes drive sawtooth crashes through magnetic reconnection and pressure mixing, with crash timing highly sensitive to the on-axis safety factor and temperature profiles.

The Gist

The Big Picture: Building a Better Fusion Reactor

Imagine scientists are trying to build a miniature sun inside a machine called a tokamak (specifically, one named SPARC) to generate limitless clean energy. To make this work, they need to trap super-hot plasma (a soup of charged particles) inside a magnetic cage.

The problem? This plasma is temperamental. It likes to wiggle, twist, and sometimes suddenly collapse, spilling its energy and cooling down. One of the most common and dangerous "wiggles" is called a sawtooth crash.

This paper is like a weather forecast and crash test simulation for the SPARC reactor. The researchers used a super-computer code (M3D-C1) to simulate what happens inside the machine when things get unstable, so they can predict how to keep the reactor running smoothly.

---

The Main Characters: The "Kink" and the "Sawtooth"

1. The Internal Kink (The Twisted Rubber Band)

Imagine the magnetic field holding the plasma is like a rubber band wrapped around a bundle of spaghetti.

• The Setup: In the center of the reactor, the magnetic field lines twist. If they twist too much, the "rubber band" wants to snap back or twist itself into a knot.
• The Instability: This is the Internal Kink Mode. It's like a rubber band that has been twisted so tightly it suddenly snaps into a new shape. In the paper, the scientists found that the most dangerous twist happens when the magnetic field lines make exactly one full loop for every turn around the machine (called the $n=1$ mode).
• The Trigger: This snap happens when two things go wrong at the same time:
  1. Too much current: The electric current in the center gets too strong (like twisting the rubber band too hard).
  2. Too much pressure: The plasma gets too hot and puffy (like the rubber band is filled with high-pressure air).

2. The Sawtooth Crash (The Heart Attack)

When that "kink" snaps, it causes a Sawtooth Crash.

• The Analogy: Think of a pot of soup that is boiling vigorously. Suddenly, the heat source gets turned off, and the soup settles down, but in a messy way.
• What happens: The hot, energetic plasma in the very center gets mixed up with the cooler plasma on the outside. The temperature in the center drops sharply (a "crash"), and the pressure profile gets "hollowed out" (like a donut where the middle is empty).
• Why it matters: If the center gets too cold, the fusion reaction slows down or stops. We need to know exactly when and why this happens to prevent it.

---

What the Scientists Did (The Simulation)

The researchers didn't just guess; they ran a digital experiment using the SPARC design.

Step 1: The Linear Scan (Finding the Weak Spot)

First, they ran a "stress test" on the machine's design. They asked: "If we change the temperature or the magnetic twist slightly, how fast does the instability grow?"

• The Finding: They found that the machine is extremely sensitive. If the central magnetic twist ($q_0$) is just a tiny bit below 1, the instability grows fast. If it's just a tiny bit above 1, it might stop. It's like balancing a pencil on its tip; a tiny nudge makes it fall.
• The Pressure Factor: They discovered that the heat (pressure) of the plasma makes the instability much worse. A hot, puffy plasma is much more likely to crash than a cool, tight one.

Step 2: The Nonlinear Simulation (The Crash Test)

Next, they let the simulation run forward in time to watch the crash happen in 3D.

• Scenario A (Just Current): If they only had the current issue (but low heat), the crash happened, but it was a "clean" snap. The plasma flattened out, but the center didn't get weirdly hollow.
• Scenario B (Just Heat): If they only had the heat issue (but the magnetic twist was safe), the plasma wiggled a bit but didn't crash hard.
• Scenario C (The Real SPARC Case): This is the big one. They combined High Current AND High Heat.
  • The Result: A massive, violent crash.
  • The "Hollow" Effect: The simulation showed something unique. The plasma didn't just flatten; it formed a hollow bubble in the middle.
  • The Metaphor: Imagine a balloon that is being squeezed from the sides (current) while the air inside is expanding (heat). Instead of just popping, the balloon twists, the air swirls around, and the center gets sucked out, creating a donut shape.

The Two Theories: Kadomtsev vs. Wesson

The paper explains this crash using two famous theories, like two different ways to explain why a car crashed:

1. The Kadomtsev Model (The Magnetic Reconnection): This theory says the magnetic field lines break and reconnect, like a tangled headphone cord that you untangle by cutting and re-tying. This flattens the plasma.
2. The Wesson Model (The Pressure Bubble): This theory says the pressure is so high it creates a swirling flow (like a whirlpool) that pushes the hot stuff out of the center, creating a hollow shape.

The Paper's Conclusion: The SPARC crash is a hybrid. It starts with the magnetic lines breaking (Kadomtsev), but the intense heat creates a swirling flow that hollows out the center (Wesson). It's a "perfect storm" of both effects.

Why Does This Matter?

If we don't understand these crashes, the SPARC reactor might fail to reach its goals.

• Alpha Particles: In a real fusion reactor, the reaction creates super-fast "alpha particles" that act as the fuel to keep the fire burning. If a sawtooth crash happens too violently, it kicks these particles out of the center, and the reactor cools down.
• The Fix: By knowing exactly how the crash happens, engineers can design the machine (or the heating systems) to prevent the "perfect storm" of high current and high pressure from happening at the same time.

The Takeaway

This paper is a warning label and a blueprint. It tells us that the SPARC reactor is powerful, but it has a "weak spot" where the magnetic twist and the heat combine to cause a violent crash. However, by simulating it, the scientists now understand the mechanics of that crash (the hollowing, the reconnection) and can start designing ways to keep the plasma stable, ensuring the future of clean fusion energy.

Technical Summary

1. Problem Statement

The SPARC tokamak, a high-field, high-current compact device designed to achieve $Q > 10$ (fusion power exceeding input power), operates under extreme plasma conditions. A critical concern for its operation is the stability of Magnetohydrodynamic (MHD) modes, specifically sawtooth oscillations.

• The Issue: Sawtooth crashes involve the periodic collapse of the central plasma temperature and density profiles. In SPARC, these events can lead to the rapid redistribution or loss of high-energy alpha particles (the primary heating source for self-sustained fusion), potentially preventing the device from reaching breakeven or high-performance regimes.
• The Gap: While previous studies (e.g., using the TRANSP code with the Porcelli model) have predicted sawtooth behavior, they rely on semi-empirical models that may not fully capture the complex 3D nonlinear physics, particularly the interplay between current-driven and pressure-driven instabilities in SPARC's high-$\beta$ (plasma pressure to magnetic pressure ratio) regime. There is a need for high-fidelity 3D simulations to understand the triggering mechanisms and the specific nature of the crash (e.g., magnetic reconnection vs. interchange instability).

2. Methodology

The authors utilized the M3D-C1 code, a high-fidelity, nonlinear, 3D extended-MHD code, to simulate SPARC scenarios.

• Equilibrium Setup:
  • The study used a "relaxed" baseline equilibrium based on the SPARC Primary Reference Discharge (PRD) design point.
  • Key Parameters: On-axis temperature $T_{e,0} = T_{i,0} = 20$ keV, density $n_0 = 4.5 \times 10^{20}$ m$^{-3}$, toroidal field $B_0 = 12.2$ T, and plasma current $I_p = 8.7$ MA.
  • The safety factor profile has $q_0 = 0.93$ (on-axis) and $q=1$ at normalized radius $\rho \approx 0.53$.
  • Note: The pedestal region was removed to focus on core physics and avoid edge-localized mode (ELM) complexities, which are beyond the scope of this study.
• Simulation Strategy:
  • Linear Phase: Scanned resistivity (Lundquist number $S$), on-axis safety factor ($q_0$), and plasma beta ($\beta_0$) to identify the dominant unstable mode and its scaling laws. A simplified 1D eigenvalue solver was also used for qualitative comparison.
  • Nonlinear Phase: Simulated the saturation and crash of the dominant mode. Three distinct scenarios were tested to isolate driving mechanisms:
    1. Baseline Case: High $\beta$ ($2.41\%$) and low $q_0$ ($0.93$).
    2. Current-Drive Dominant: High $\beta$ but $q_0$ shifted closer to 1 ($0.991$ and $1.008$) to reduce current drive.
    3. Pressure-Drive Dominant: Low $\beta$ (relaxed temperature profiles) with low $q_0$ ($\approx 0.92$).
  • Physics Models: Included Spitzer resistivity, Ohmic heating, and two-fluid effects (via extended-MHD equations). Alpha particle heating and external heating (ICRH) were not included in the baseline simulation, serving as a "worst-case" instability scenario.

3. Key Contributions

• First High-Fidelity 3D Study: This is the first comprehensive 3D nonlinear MHD study of low-$n$ instabilities specifically for SPARC baseline-like equilibria.
• Mechanism Identification: The paper distinguishes between pure current-driven (Kadomtsev) and pressure-driven (Wesson) sawtooth mechanisms and demonstrates that SPARC's baseline scenario is a hybrid of both.
• Validation of Scaling Laws: Confirmed theoretical scaling laws for resistive kink and tearing modes in the SPARC parameter regime.
• Insight into Hollowed Profiles: Provided a physical explanation for the "hollowed" pressure profiles observed in high-performance tokamaks, linking them to pressure-driven interchange motions triggered by the kink mode.

4. Key Results

Linear Stability Analysis

• Dominant Mode: The most unstable mode is the $n=1, m=1$ internal kink mode.
• Scaling: The growth rate ($\gamma$) scales with the Lundquist number ($S$) consistent with resistive kink ($\gamma \propto S^{-1/3}$) and tearing-like ($\gamma \propto S^{-3/5}$) regimes, depending on the resistivity level.
• Sensitivity: The instability is highly sensitive to $q_0$. As $q_0$ approaches 1 (from 0.991 to 1.008), the growth rate drops by an order of magnitude.
• Beta Effect: High plasma beta ($\beta_0 > 1\%$) significantly enhances the growth rate (by 1–2 orders of magnitude) and broadens the radial mode structure. The 1D eigenvalue solver confirmed that the pressure gradient term ($g(r)$) is crucial for destabilizing the mode and smoothing the eigenfunction at the $q=1$ surface.

Nonlinear Simulations

• Baseline Case (High $\beta$, Low $q_0$):
  • A strong sawtooth crash occurs shortly after $q_0$ drops below unity.
  • Profile Flattening: The central pressure drops by 40–50%.
  • Hollowed Structure: Unlike standard Kadomtsev crashes which flatten the profile, this case produces a hollowed pressure profile (a "bubble" of low pressure in the center).
  • Mechanism: The crash involves a clear magnetic reconnection event (Kadomtsev-like) where a new magnetic island replaces the old one. However, simultaneously, a strong circulatory flow driven by the high-$\beta$ pressure gradient induces an interchange-type motion (Wesson-like), creating the hollowed structure.
• Current-Drive Only (High $\beta$, $q_0 \approx 1$):
  • When $q_0$ is shifted to 0.991, the crash still occurs but is weaker.
  • When $q_0$ reaches 1.008, the mode becomes marginally stable, and no significant crash occurs, highlighting the sensitivity of the instability to $q_0$ near unity.
• Pressure-Drive Only (Low $\beta$, Low $q_0$):
  • Sawtooth-like oscillations occur with periodic crashes.
  • The pressure flattening is moderate (20–30%) and does not produce a hollowed structure.
  • The crashes are driven primarily by the current profile, resembling the classic Kadomtsev model.

5. Significance and Future Outlook

• Performance Assessment: The results suggest that without sufficient stabilizing effects (such as alpha particle stabilization or external heating), SPARC could experience severe sawtooth crashes that drastically reduce central temperature and potentially degrade alpha particle confinement.
• Hybrid Mechanism: The study establishes that SPARC's sawtooth trigger is likely a combined effect of current peaking (driving the kink) and high pressure (driving the interchange/hollowing). This challenges simple models that treat these mechanisms separately.
• Timescales: The Ohmic heating recovery time in the simulation is estimated at ~1.6 seconds, suggesting long periods between crashes if no other heating is present. Future work must include ICRH and alpha heating to determine if faster, keV-level sawtooth cycles are possible.
• Next Steps: The authors plan to incorporate kinetic effects (hybrid kinetic-MHD models) to account for alpha particle stabilization, include realistic heating sources, and explore a wider range of equilibria to optimize SPARC's operational window and prevent disruptive MHD events.

In conclusion, this work provides a foundational understanding of MHD stability in SPARC, demonstrating that the device's high-performance parameters make it uniquely susceptible to a complex, hybrid sawtooth instability that requires advanced modeling to predict and control.