Screening effect of plasma flow on the resonant magnetic perturbation penetration in tokamak based on two-fluid model

Using the updated MDC two-fluid code, this study reveals that bootstrap current enables finite mode penetration at zero rotation and demonstrates that sufficiently large diamagnetic drift flow stabilizes neoclassical tearing modes while inducing island width oscillations driven by negative pressure feedback.

The Gist

Imagine a Tokamak as a giant, high-tech donut-shaped oven designed to cook nuclear fuel (plasma) to temperatures hotter than the sun. To keep this super-hot soup contained, scientists use powerful magnetic fields, like invisible walls holding the liquid in place.

However, sometimes these magnetic walls get a little wobbly. They can develop "kinks" or ripples called magnetic islands. Think of these islands like bubbles forming in a pot of boiling water. If a bubble gets too big, it can break the pot (the plasma containment), causing the whole experiment to fail.

This paper is about a specific tool scientists use to try to fix or control these bubbles: Resonant Magnetic Perturbations (RMPs). You can think of RMPs as a "magnetic tuning fork" that scientists tap against the plasma to try to smooth out the ripples or lock the bubbles in a safe spot.

Here is what the researchers discovered, explained simply:

1. The "Seed" Problem

Sometimes, a tiny bubble (a "seed island") appears naturally. If the plasma is just sitting there, a small tap from the tuning fork (RMP) might just make the bubble wiggle a bit. But if the plasma has a special internal current (called bootstrap current, which acts like a self-sustaining engine), that same small tap can suddenly cause the bubble to explode in size.

• The Analogy: Imagine pushing a swing. If the swing is empty, you have to push hard to make it go high. But if the swing is already moving in rhythm with your push (the bootstrap current), even a tiny nudge can send it flying. The researchers found that without plasma flow, there is a "tipping point" where a small push suddenly creates a huge problem.

2. The "Wind" Effect (Plasma Flow)

The plasma inside the donut isn't still; it's spinning and flowing like a river. The researchers wanted to see how this "wind" affects the magnetic bubbles. They looked at two types of wind:

• Electric Drift: Like a wind blowing because of an electric charge.
• Diamagnetic Drift: Like a wind blowing because of pressure differences (like air rushing out of a tire).

The Discovery:
They found that if the plasma is spinning fast enough, it acts like a shield.

• The Analogy: Imagine trying to push a heavy door open. If the door is locked (no flow), a small push might just jiggle it. But if the door is on a fast-moving conveyor belt (plasma flow), the wind blowing past it actually pushes the door back, making it much harder for your "tuning fork" (RMP) to get inside and disturb the bubble. This is called the screening effect. The faster the plasma spins, the better it hides the bubble from the external magnetic taps.

3. The "Bouncing" Bubble (Oscillation)

Here is the most surprising part. When the plasma flow was very strong (specifically the pressure-driven "diamagnetic" wind), the magnetic bubble didn't just grow or shrink; it started pulsing or bouncing up and down in size.

• The Analogy: Imagine a balloon being squeezed. As you squeeze it, the air pressure inside builds up and pushes back, making the balloon expand again. Then it gets squeezed again.
• What happened in the paper: The magnetic bubble grew, which flattened the pressure inside it. This change in pressure altered the "wind" (diamagnetic flow), which then pushed back on the bubble, making it shrink. As it shrank, the pressure changed again, and the cycle repeated. It was a negative feedback loop: the bubble's own growth created the conditions to stop its growth, leading to a rhythmic dance of expanding and contracting.

4. Why This Matters for the Study

The researchers used a supercomputer simulation (their "MDC" code) to test these ideas. They found that:

• If you ignore the plasma flow, you might think a small magnetic tap will always cause a big problem.
• But if you include the flow, the plasma can actually protect itself (screening).
• However, if the flow is too strong and specific conditions are met, the bubble starts oscillating (bouncing) instead of staying still.

In a Nutshell:
This paper explains that the plasma in a fusion reactor isn't just a passive target; it's an active participant. It can spin fast enough to block outside magnetic disturbances, but under certain high-pressure conditions, it can also start "breathing" (oscillating) in a complex dance between pressure and magnetic fields. Understanding this dance helps scientists figure out how to keep the fusion reactor stable and prevent those dangerous magnetic bubbles from breaking the containment.

Technical Summary

1. Problem Statement

The research addresses the critical issue of Neoclassical Tearing Modes (NTMs) in tokamak plasmas, which are a primary limitation on core plasma temperature and a cause of major disruptions. NTMs are driven by the loss of bootstrap current within magnetic islands. A key challenge in controlling NTMs is understanding the interaction between Resonant Magnetic Perturbations (RMPs) and magnetic islands.

Specifically, the paper investigates the phenomenon of mode penetration, where an external RMP overcomes plasma flow screening to drive magnetic reconnection and form a magnetic island. While single-fluid models have been used to study this, they often fail to capture complex physics related to plasma rotation and finite Larmor radius effects. The authors aim to clarify:

• The distinction between "driven reconnection" and true "mode penetration" in the presence of bootstrap current.
• The specific screening roles of $E \times B$ drift flow versus diamagnetic drift flow.
• The physical mechanisms behind observed oscillations in island width under high-flow conditions.

2. Methodology

The study utilizes a newly upgraded initial value code, MDC (MHD@Dalian Code), based on a two-fluid, four-field MHD model.

• Governing Equations: The model solves four coupled nonlinear equations for:
  1. Vorticity ($U$)
  2. Poloidal magnetic flux ($\psi$)
  3. Plasma pressure ($p$)
  4. Parallel ion velocity ($v$)
• Key Physics Included:
  • Bootstrap Current ($j_b$): Crucial for NTM destabilization.
  • Transport Effects: Both parallel ($\chi_{\parallel}$) and perpendicular ($\chi_{\perp}$) heat transport.
  • Two-Fluid Effects: Retains electron diamagnetic drift and $E \times B$ flow, accounting for finite Larmor radius (FLR) effects via the parameter $\delta$.
  • Geometry: Cylindrical geometry $(r, \theta, z)$ with periodic boundary conditions.
• Numerical Setup:
  • Simulations focus on a low-density, ohmically heated discharge scenario ($n_e \approx 2 \times 10^{19} \text{ m}^{-3}$, $B_0 = 2 \text{ T}$).
  • The resonant surface is set at $q = 3/2$ ($r \approx 0.4$).
  • RMPs are applied as boundary conditions to simulate external perturbation.
  • The study systematically varies RMP amplitude, plasma rotation frequency, bootstrap current fraction ($f_b$), resistivity ($\eta$), and transport coefficients.

3. Key Contributions and Results

A. Re-evaluation of Mode Penetration Thresholds

• Classical Tearing Mode ($f_b = 0$): Without bootstrap current, there is no threshold for mode penetration. The magnetic island width scales linearly with the RMP amplitude (driven reconnection).
• Neoclassical Tearing Mode ($f_b \neq 0$): When bootstrap current is included, a "mode penetration-like" phenomenon is observed.
  • At low RMP amplitudes, the island remains small.
  • Once the RMP exceeds a critical threshold, the island grows rapidly to a large size.
  • Distinction: The authors clarify that this is not a single threshold but a transition from a driven reconnection regime to an NTM growth regime. This explains why experiments often detect a finite penetration threshold even at zero rotation, which single-fluid models might misinterpret.

B. Screening Effects: $E \times B$ vs. Diamagnetic Drift

The study compares the screening capabilities of $E \times B$ flow and diamagnetic drift flow:

• Penetration Threshold: Both flows provide similar screening effects regarding the threshold required to penetrate the plasma. The threshold depends on the rotation difference between the resonant surface and the RMP, regardless of the flow type.
• Stabilization: A significant difference emerges in the saturated island width.
  • $E \times B$ Flow: Leads to larger saturated island widths.
  • Diamagnetic Drift Flow: Exerts a strong stabilizing effect, resulting in significantly smaller saturated island widths compared to $E \times B$ flow at the same rotation frequency. This is attributed to the $\delta \nabla_{\parallel} p$ term in the generalized Ohm's law.

C. Discovery of Island Width Oscillation

In regimes with high Lundquist numbers ($S$) and high parallel-to-perpendicular transport ratios ($\chi_{\parallel}/\chi_{\perp}$), the authors discovered a periodic oscillation of the magnetic island width after penetration.

• Mechanism: The oscillation is driven by a negative feedback loop between the magnetic island and plasma pressure.
  1. As the island grows, the pressure gradient inside the island flattens, but the gradient at the resonant surface changes.
  2. This alters the diamagnetic flow frequency ($\omega_{dia}$).
  3. The change in flow affects the island width via the two-fluid terms.
  4. The system oscillates as the pressure gradient and island width regulate each other.
• Verification: The oscillation disappears when:
  • Resistivity ($\eta$) is increased (diffusion dampens the effect).
  • Bootstrap current fraction ($f_b$) is reduced.
  • Parallel transport ($\chi_{\parallel}$) is reduced or perpendicular transport ($\chi_{\perp}$) is increased.

4. Significance and Conclusion

This work provides critical insights for the operation of future advanced tokamaks (e.g., ITER, CFETR), which rely heavily on high bootstrap current fractions.

1. Control Strategies: The findings suggest that while RMPs can be used to lock and locate NTMs for Electron Cyclotron Current Drive (ECCD) suppression, the diamagnetic drift plays a vital stabilizing role that should be leveraged.
2. Physics Modeling: The study highlights the necessity of using two-fluid models over single-fluid models to accurately predict NTM thresholds and island dynamics, particularly regarding the stabilizing influence of diamagnetic flows.
3. Operational Windows: The identification of oscillation regimes helps define operational boundaries where magnetic islands might behave unpredictably due to pressure-flow coupling, aiding in the design of robust disruption mitigation systems.

In summary, the paper demonstrates that plasma flow screening is not uniform; diamagnetic drift is more effective at stabilizing magnetic islands than $E \times B$ flow, and the interplay between pressure gradients and flow can lead to complex oscillatory behaviors in high-performance tokamak scenarios.