On the mechanism of Pedestal Relaxation Events -- Insights gained by turbulence simulations with GRILLIX

This study utilizes global trans-collisional fluid simulations with GRILLIX to demonstrate that Pedestal Relaxation Events (PREs) in ASDEX Upgrade I-mode discharges are triggered by Micro-Tearing Modes (MTMs), successfully reproducing experimental characteristics and validating a qualitative cycle model through agreement with linear growth-rate theory.

The Gist

The Big Picture: Taming the Fusion Fire

Imagine a fusion reactor as a giant, super-hot pot of soup (plasma) that we are trying to keep boiling without spilling over the sides. To get enough energy out of it, we need to keep the soup very hot and dense right at the edge of the pot. This hot, dense layer is called the "pedestal."

Sometimes, this pedestal gets unstable and suddenly dumps a little bit of energy out. In the world of fusion, we have two types of these "spills":

1. The Big Spills (ELMs): These are like massive tsunamis crashing over the wall. They are dangerous and can damage the reactor.
2. The Small Burps (PREs): These are the focus of this paper. They are tiny, periodic "burps" of energy. They are much smaller than the big spills (only about 1% of the energy), but they still happen frequently, especially when the reactor is operating in a special, efficient mode called "I-mode."

Scientists have known these "burps" happen, but they didn't know exactly why or how they started. This paper uses a super-computer simulation to figure it out.

The Detective Work: Finding the Culprit

The researchers used a software tool called GRILLIX (think of it as a high-tech weather forecast for plasma) to simulate a specific fusion experiment. They watched the simulation run for a few milliseconds and saw three of these "burps" (PREs) happen.

They asked: What is causing these burps?

They looked for clues, much like a detective looking for fingerprints at a crime scene. They found three main clues that pointed to a specific suspect: Micro-Tearing Modes (MTMs).

• Clue 1: The Heat Pattern. When the burp happened, the heat (electron temperature) flattened out, but the density didn't change much. This is exactly what you'd expect if the "tearing" was happening.
• Clue 2: The Magnetic Shape. They looked at the magnetic fields inside the plasma. The pattern looked like a "tear" in the fabric of the magnetic field. In physics, this specific shape is called "tearing parity," and it is the signature of MTMs.
• Clue 3: The Speed. They measured how fast the waves were moving. The speed matched the theoretical prediction for MTMs perfectly.

The Verdict: The "burps" are caused by tiny, electromagnetic rips (MTMs) in the magnetic field that let heat escape quickly.

The Cycle: How a "Burp" Happens

The paper sketches a cycle of how these events repeat, like a rubber band being stretched and snapped:

1. The Stretch: The temperature gradient (how fast the heat changes from the center to the edge) gets steeper and steeper. Think of this as stretching a rubber band.
2. The Snap: Eventually, the rubber band gets too tight. The Micro-Tearing Mode (MTM) suddenly wakes up and starts growing.
3. The Release: The MTM creates a "stochastic" (chaotic) magnetic field, acting like a shortcut for heat to escape. The temperature gradient flattens out instantly.
4. The Calm: Because the gradient is now flat, the MTM loses its fuel (the steep temperature difference) and dies out.
5. Repeat: The system starts stretching the rubber band again, and the cycle begins anew.

The Secret Ingredient: The "Landau" Recipe

One of the most important findings in this paper is about the math used to run the simulation.

To simulate plasma, scientists have to make choices about how to calculate heat flow.

• The Old Recipe (Braginskii): This is like using a simple rule of thumb. When the researchers used this, the simulation was calm. No burps happened.
• The New Recipe (Landau-fluid): This is a more complex, "non-local" method. It accounts for the fact that particles can travel far without bumping into each other (low collisionality). When they used this recipe, the "burps" appeared!

The Takeaway: The "burps" only happen when you use the advanced math that accounts for these long-distance particle movements. This suggests that in the real, low-collisionality edge of a fusion reactor, these burps are real and driven by this specific physics.

A Note of Caution: The Simulation vs. Reality

The authors are very honest about one difference between their simulation and the real experiment:

• In the Experiment: The "burp" happens, and the stored energy goes down (the pot cools slightly).
• In the Simulation: The "burp" happens, but the stored energy goes up.

Why? It's a quirk of how they set up the simulation. When the heat escapes, the computer automatically pumps in more power to keep the temperature steady, which accidentally adds more energy than was lost. However, the authors argue that the mechanism (the tearing mode causing the heat to escape) is still correct, even if the energy balance is slightly off due to this setup.

The "Why" Behind the "When"

Finally, the paper asks: "If the real experiment (ASDEX Upgrade) didn't have these burps at this specific moment, why did our simulation show them?"

They suspect it's because of resistivity (how much the plasma resists electric current). The math they used (Spitzer resistivity) might underestimate how much resistance there is at very high temperatures. If the resistance were actually higher, it would dampen (stop) the "tearing" modes, preventing the burps. Since their math underestimated the resistance, the "tearing" modes grew too easily in the simulation.

Summary

This paper uses advanced computer simulations to show that small, periodic energy releases (PREs) in fusion reactors are caused by tiny magnetic "tears" (MTMs). These tears grow when the temperature gradient gets too steep, snap open to let heat escape, and then die out, only to repeat the cycle. The study highlights that using the correct, advanced math (Landau-fluid closure) is essential to seeing these phenomena, and it suggests that improving how we calculate electrical resistance in our models will help us predict exactly when and where these events will happen in real fusion reactors.

Technical Summary

Technical Summary: Mechanism of Pedestal Relaxation Events in I-Mode Discharges

Problem Statement
Future fusion reactors require operational regimes that offer high energy confinement without exceeding material limits at the vessel wall. The Improved Confinement (I-mode) regime is a candidate for such conditions, characterized by elevated confinement, low impurity content, and the absence of large Type-I Edge Localised Modes (ELMs). However, I-mode discharges, particularly those near the I-H transition, can exhibit Pedestal Relaxation Events (PREs). While PREs involve periodic energy ejections similar to ELMs, they occur at significantly lower energy levels (~1% of stored energy) and are believed to be driven by a different mechanism. Unlike ELMs, which are triggered by Peeling-Ballooning modes, the underlying physics of PREs remains less understood. Previous attempts to simulate PREs using gyrofluid models in simple geometries suggested Micro-Tearing Modes (MTMs) as candidates but observed a transition to interchange-like turbulence without significant electromagnetic transport. Furthermore, simulating the low-collisionality edge of fusion devices (where $\rho_{pol} > 0.9$) with fluid models is challenging due to the breakdown of standard collisional assumptions.

Methodology
The authors employed the global trans-collisional fluid turbulence code GRILLIX to simulate an I-mode discharge from the ASDEX Upgrade tokamak (Shot #37980). The simulation utilized a drift-reduced electromagnetic Braginskii model extended for trans-collisional regimes, specifically incorporating:

• A Landau-fluid closure for parallel conductive heat fluxes to address low collisionality.
• A neoclassical extension for ion viscosity.
• A three-moment neutral gas model.
• Adaptive source terms to maintain core boundary values ($n, T_e, T_i$) at $\rho_{pol} = 0.91$.

The simulation spanned over 5 ms of physical time. To validate the findings, the authors:

1. Compared simulation outputs (profiles, magnetic activity, event duration) with experimental data.
2. Performed a detailed analysis of mode parity and dispersion relations to identify the microinstability.
3. Traced the system's evolution in gradient length space ($a/L_{Te}$ vs. $a/L_n$) against a linear growth rate estimate derived from a simplified slab geometry analysis.
4. Conducted a comparative simulation using the standard Braginskii closure with heat-flux limiters to isolate the effect of the Landau-fluid closure.
5. Performed linear stability analysis for both closures to determine the conditions for MTM growth.

Key Results

• Observation of PREs: The simulation reproduced three intermittent bursts identified as PREs. These events matched experimental characteristics, including strong magnetic activity, a duration of approximately 0.3 ms, and a primary flattening of the electron temperature profile rather than the density profile.
• Identification of MTMs: Detailed analysis confirmed that the PREs are driven by Micro-Tearing Modes. Evidence includes:
  • Tearing Parity: The parallel vector potential ($\tilde{A}_{\parallel}$) exhibited even parity, while the electrostatic potential ($\tilde{\phi}$) showed odd parity.
  • Transport Characteristics: The events were dominated by radial electromagnetic electron heat transport (magnetic flutter).
  • Dispersion Relation: The computed dispersion relation matched linear kinetic estimates for MTMs, showing propagation in the electron diamagnetic direction with a peak at $k_y\rho_s \approx 0.2$.
• Global Cycle Mechanism: The authors sketched a qualitative PRE cycle:
  1. The electron temperature gradient ($a/L_{Te}$) steepens over time.
  2. Once the gradient exceeds a critical threshold, the MTM grows (linear growth rate $\hat{\gamma} > 0$).
  3. The growing MTM induces strong radial electron heat transport, rapidly flattening the temperature gradient.
  4. The system enters a linearly stable region ($\hat{\gamma} < 0$), causing the mode to vanish.
  5. The gradient slowly steepens again, restarting the cycle.
• Role of Fluid Closure: The Landau-fluid closure was found to be essential. A comparative simulation using the Braginskii closure (with flux limiters) resulted in a quiescent state with no PREs. Linear analysis revealed that the Landau-fluid closure introduces non-local effects (via the $k_{\parallel}/|k_{\parallel}|$ term) that yield a positive growth rate for MTMs, whereas the Braginskii closure predicts linear stability.
• Discrepancy in Stored Energy: Unlike experiments where stored energy decreases after a PRE, the simulation showed an increase. The authors attribute this to the adaptive boundary conditions: as radial transport increases, the adaptive source injects more heating power to maintain the core target values, artificially boosting stored energy.
• Resistivity Limitations: The authors note that PREs were observed in the simulation for a discharge where they did not occur experimentally. They hypothesize this is due to the Spitzer resistivity model used, which underestimates resistivity at high temperatures. This underestimation reduces the damping term in the growth rate equation, artificially destabilizing the MTMs.

Significance and Claims
The paper claims to provide the first global fluid simulation of I-mode discharges that successfully reproduces the qualitative features of PREs and identifies MTMs as the triggering mechanism. The work highlights the critical importance of the Landau-fluid closure for modeling low-collisionality edge turbulence, demonstrating that standard fluid closures (Braginskii) may fail to capture these instabilities.

The authors modestly state that while the simulation predicts PREs earlier than observed in the specific experimental discharge (likely due to resistivity modeling limitations), the qualitative picture of the PRE cycle—driven by the interplay between electron temperature gradients and MTM growth—is robust. They conclude that understanding this mechanism, specifically the reliance on the electron temperature gradient, could aid in experimental strategies to mitigate PREs. The paper suggests that future work should improve resistivity models within fluid codes or utilize higher-fidelity gyrokinetic codes (like GENE-X) to address low-collisionality regimes more accurately.