Empirical impact of near-separatrix plasma and neutral transport on the pedestal in the transition between EDA and ELMy H-modes on Alcator C-Mod

This study investigates the transition between ELMy and EDA H-modes on Alcator C-Mod by validating pedestal density models, identifying a distinct sensitivity to neutral sources in ELMy regimes, and demonstrating that incorporating a resistive ballooning mode transport channel significantly improves predictions for EDA regimes while reducing predicted pedestal density by approximately 20%.

The Gist

Imagine a fusion reactor as a giant, super-hot pot of soup (the plasma) that we are trying to keep boiling without it splashing out and burning the pot (the reactor walls).

To keep this soup hot enough to create energy, we need to build a "wall" of heat right at the edge of the pot. This wall is called the pedestal. The taller and steeper this wall is, the hotter the soup gets in the middle. However, if the wall gets too steep, it becomes unstable and collapses, sending a massive wave of hot soup crashing into the walls. In fusion terms, this is called an ELM (Edge Localized Mode), and it's like a giant splash that could damage the reactor.

Scientists want to find a way to build a tall heat wall that doesn't collapse. This paper studies two different ways the "soup" behaves in a specific machine called Alcator C-Mod to figure out how to build that perfect, stable wall.

Here is the story of their findings, broken down into simple concepts:

1. The Two Types of "Soup" Behavior

The researchers looked at two different states of the plasma:

• The "Splashing" Mode (ELMy H-mode): Think of this like a pot of water boiling violently. The heat wall builds up until it gets too high, then it suddenly snaps and releases a big burst of energy (a big ELM). In this mode, the amount of "soup" (density) is very sensitive to how much water you add. If you turn up the faucet (fueling), the soup level rises quickly.
• The "Steady Stream" Mode (EDA H-mode): This is the "holy grail" the scientists are looking for. Imagine a pot where the water is swirling so smoothly that it never splashes. Instead of big bursts, it releases a tiny, continuous trickle of energy. In this mode, the soup level is stiff. Even if you turn the faucet up or down, the level of the soup inside the wall doesn't change much. It seems to have its own internal thermostat.

2. The Secret Ingredient: Neutral Particles

The researchers discovered why these two modes behave differently. It comes down to "ghosts" in the machine—neutral atoms that haven't been turned into plasma yet.

• In the Splashing Mode, these ghosts are easy to push around. If you add more fuel, the ghosts push the plasma level up easily.
• In the Steady Stream Mode, the ghosts get stuck. The plasma becomes so dense that the ghosts can't penetrate deep inside. The plasma level becomes "stiff" because the transport of particles is being controlled by a different, more chaotic mechanism (turbulence) rather than just how much fuel you add.

3. The "Quasi-Coherent Mode" (The Heartbeat)

In the Steady Stream mode, the plasma develops a specific rhythm, like a heartbeat, called the Quasi-Coherent Mode (QCM).

• Think of this as a gentle, rhythmic breathing that keeps the pot from overflowing.
• The researchers found that this "heartbeat" gets stronger and more organized as the soup gets denser, but only up to a point. If it gets too dense, the heartbeat starts to get weak again, and the soup might become unstable in a different way.
• This heartbeat is the key mechanism that prevents the big, damaging splashes (Type-I ELMs).

4. Building a Better Map (The Models)

Scientists use computer models to predict how high the heat wall will get.

• The Old Map: They tested a model called "Saarelma-Connor." It worked great for the "Splashing Mode" but failed to predict the "Steady Stream" mode correctly. It thought the wall would get too high.
• The New Map: They realized that in the Steady Stream mode, there is an extra "leak" in the system caused by a specific type of turbulence (called Resistive Ballooning Modes). When they added this "leak" to their map, the predictions finally matched reality. The model now correctly predicts that the wall height will be limited, preventing the big splashes.

5. Looking to the Future: SPARC

Finally, the team used their new understanding to predict what will happen in SPARC, a new, bigger fusion machine currently being built.

• They ran two scenarios: one that looks like the old "Splashing" mode and one that tries to achieve the "Steady Stream" mode.
• The Good News: If SPARC can operate in the high-density "Steady Stream" mode, the extra "leak" (turbulence) they discovered will act as a safety valve. It will naturally limit how high the heat wall gets, keeping it stable and preventing the dangerous splashes that could damage the machine.
• The Result: This suggests that running fusion reactors at very high densities might be the safest and most efficient way to go, as it naturally avoids the destructive "splashes."

The Big Takeaway

This paper is like a detective story where scientists figured out how to stop a hot pot from splashing. They found that by cranking up the density and letting a specific "rhythmic breathing" (the QCM) take over, the plasma creates its own safety valve. This discovery gives them confidence that future fusion reactors can run hot and stable without breaking their walls.

Technical Summary

1. Problem Statement

The primary challenge in fusion energy research is integrating a hot core plasma with a cold edge plasma while avoiding large, damaging Edge Localized Modes (ELMs). Type-I ELMy H-modes, while high-performing, expel large fractions of plasma energy onto divertor components, a scaling issue that threatens future reactors like SPARC and ITER.
Alternative regimes, such as the Enhanced D-alpha (EDA) H-mode and Quasi-Continuous Exhaust (QCE), offer high-density operation with small or no ELMs, improving divertor survivability. However, the physical mechanisms governing the transition between ELMy and EDA regimes, particularly regarding pedestal density limits, neutral transport, and turbulence, are not fully understood. This lack of understanding hinders the predictive capability required to design high-performance, ELM-free scenarios for next-step devices.

2. Methodology

The authors utilized a high-resolution experimental database from Alcator C-Mod, focusing on a specific run day where plasmas transitioned from EDA to ELMy H-modes by varying the line-averaged density. The study employed a multi-faceted approach:

• Experimental Diagnostics:
  • Thomson Scattering (ETS): Used to measure electron density ($n_e$) and temperature ($T_e$) profiles. Two fitting functions were applied: a modified hyperbolic tangent (mtanh) for the pedestal top and an exponential decay fit for the separatrix and near-Scrape-Off Layer (SOL) to determine separatrix density ($n_{sep}$) and gradient scale lengths.
  • Phase Contrast Imaging (PCI): Used to analyze density fluctuation spectra, specifically tracking the Quasi-Coherent Mode (QCM) characteristic of EDA and comparing it to inter-ELM fluctuations in ELMy regimes.
  • Neutral Pressure: Outer midplane neutral pressure ($p_0^{OMP}$) was measured to correlate with fueling sources.
• Modeling and Validation:
  • Saarelma-Connor Model: A reduced model combining neutral penetration and plasma transport was validated against the C-Mod dataset. The model was extended to include a Resistive Ballooning Mode (RBM) driven transport channel ($D_{RBM}$) to account for high-density effects.
  • EPED Model: The EPED (Pedestal) code was used to perform scans of pedestal density ($n_{ped}$) at varying separatrix-to-pedestal density ratios ($n_{sep}/n_{ped}$) to test the validity of the Kinetic Ballooning Mode (KBM) and Peeling-Ballooning Mode (PBM) stability constraints for these regimes.
  • KN1D Simulations: A 1D kinetic neutral code was used to parameterize the separatrix neutral density boundary condition ($n_{sep}^0$) based on experimental wall pressure measurements.

3. Key Contributions

• Empirical Characterization of the Transition: The study provides a comprehensive map of the operational space between ELMy and EDA H-modes, identifying critical thresholds in collisionality ($\alpha_t$) and electron pressure ($\beta_e$) that govern the transition.
• Neutral vs. Transport Control: It demonstrates a fundamental difference in density control mechanisms: ELMy H-modes are sensitive to neutral fueling sources, whereas EDA H-modes exhibit "stiff" density profiles where pedestal density is regulated by internal particle transport rather than external fueling.
• Fluctuation Physics: The paper links the emergence and strength of the QCM to dimensionless parameters ($\nu^* \alpha_{MHD}$ and the ratio of RBM to Electromagnetic turbulence wavenumbers, $k_{RBM}/k_{EM}$), suggesting the QCM arises from the coupling of drift-Alfvén waves and resistive ballooning modes.
• Model Extension: The Saarelma-Connor model was successfully extended to predict EDA H-mode pedestals by incorporating an RBM-driven transport channel, significantly improving prediction accuracy at high densities ($n_{ped}^e > 2.0 \times 10^{20} \text{ m}^{-3}$).
• SPARC Predictions: The study applies these findings to predict the edge pedestal for the SPARC tokamak, contrasting a standard Primary Reference Discharge (PRD) with a proposed high-density, ELM-free scenario.

4. Key Results

• Density Stiffness and Neutrals:
  • In ELMy H-modes ($p_0^{OMP} < 0.1$ mTorr), $n_{ped}^e$ scales linearly with neutral pressure.
  • In EDA H-modes ($p_0^{OMP} > 0.1$ mTorr), $n_{ped}^e$ saturates at $\approx 2.5 \times 10^{20} \text{ m}^{-3}$ regardless of increasing neutral pressure. The separatrix density ($n_{sep}^e$) continues to rise, flattening the pedestal gradient ($n_{sep}^e / n_{ped}^e \to 0.8$).
• Fluctuation Evolution:
  • The QCM amplitude ($B$) grows with density and dimensionless parameters ($\nu^* \alpha_{MHD}$), saturating at high densities.
  • Background turbulence levels ($A$) increase linearly with $\alpha_t$ and $\beta_e$ in the EDA regime, suggesting enhanced filamentary transport.
  • The QCM strength correlates inversely with the ratio $k_{RBM}/k_{EM}$, peaking when these wavenumbers are comparable ($k_{RBM} \sim k_{EM}$), supporting the theory of DAW-RBM coupling.
• Model Validation:
  • The standard Saarelma-Connor model (without RBM) underpredicts $n_{ped}^e$ for EDA regimes.
  • Adding an RBM transport term ($D_{RBM} \propto \alpha_t$ or $1/k_{RBM}^2 \hat{q}_{cyl}$) successfully reproduces the experimental saturation of $n_{ped}^e$ up to $3.0 \times 10^{20} \text{ m}^{-3}$.
• EPED Comparison:
  • EPED predictions (based on KBM-PBM stability) agree well with large ELMy discharges but overpredict pressure for many EDA cases.
  • Increasing the ratio $n_{sep}^e / n_{ped}^e$ shifts the peeling-ballooning transition to lower densities and reduces the maximum attainable pressure in the ballooning branch.
  • High-density EDA pedestals often deviate from KBM scaling, exhibiting wider profiles than predicted, likely due to resistive instabilities.
• SPARC Predictions:
  • For the standard PRD scenario, the model predicts a pedestal density consistent with EPED expectations.
  • For the proposed high-density scenario, the inclusion of RBM transport significantly limits the pedestal density gradient near the separatrix. This reduces the predicted $n_{ped}^e$ by approximately 20% compared to models without RBM, potentially avoiding Type-I ELMs by keeping the pressure gradient below the PBM limit.

5. Significance

This work bridges the gap between empirical observation and predictive modeling for ELM-free high-density regimes. By identifying the specific role of resistive turbulence (RBM) and neutral transport in limiting pedestal density, the authors provide a physical basis for accessing stable, high-performance regimes in future reactors.
The findings suggest that operating at high densities (EDA-like) may be a viable strategy for SPARC and ITER to mitigate ELM damage. The validated models offer a pathway to predict edge profiles in regimes where traditional MHD stability limits (EPED) are insufficient, highlighting the necessity of coupling local turbulence physics with global stability constraints. This is a critical step toward designing fusion power plants that can sustain high performance without catastrophic heat loads on plasma-facing components.