Determination of confinement regime boundaries via separatrix parameters on Alcator C-Mod based on a model for interchange-drift-Alfv\'en turbulence

TL;DR

This study validates a turbulence-based model predicting confinement regime boundaries in Alcator C-Mod, using experimental data to improve understanding of transitions like L-H and density limits, with implications for future fusion devices.

Contribution

The paper demonstrates the applicability of the SepOS model to Alcator C-Mod, extending its predictive capability across various operational boundaries and conditions, including the L-H transition and density limits.

Findings

The SepOS model accurately predicts key confinement boundaries in C-Mod.

Empirical scaling laws relate turbulence parameters to operational regimes.

The model's predictions align with experimental observations across diverse plasma conditions.

Abstract

The separatrix operational space (SepOS) model [Eich \& Manz, \emph{Nuclear Fusion} (2021)] is shown to predict the L-H transition, the L-mode density limit, and the ideal MHD ballooning limit in terms of separatrix parameters for a wide range of Alcator C-Mod plasmas. The model is tested using Thomson scattering measurements across a wide range of operating conditions on C-Mod, spanning $\overline{n}_{e} = 0.3 - 5.5 \times 10^{20}$m$^{-3}$, $B_{t} = 2.5 - 8.0$ T, and $B_{p} = 0.1 - 1.2$ T. An empirical regression for the electron pressure gradient scale length, $\lambda_{p_{e}}$, against a turbulence control parameter, $\alpha_{t}$, and the poloidal fluid gyroradius, $\rho_{s,p}$, for H-modes is constructed and found to require positive exponents for both regression parameters, indicating turbulence widening of near-SOL widths at high $\alpha_{t}$ and an inverse scaling with $B_{p}$, consistent with results on AUG. The SepOS model is also tested in the unfavorable drift direction and found to apply well to all three boundaries, including the L-H transition as long as a correction to the Reynolds energy transfer term, $\alpha_\mathrm{RS} < 1$ is applied. I-modes typically exist in the unfavorable drift direction for values of $\alpha_{t} \lesssim 0.35$. Finally, an experiment studying the transition between the type-I ELMy and EDA H-mode is analyzed using the same framework. It is found that a recently identified boundary at $\alpha_{t} = 0.55$ excludes most EDA H-modes but that the balance of wavenumbers responsible for the L-mode density limit, namely $k_\mathrm{EM} = k_\mathrm{RBM}$, may better describe the transition on C-Mod. The ensemble of boundaries validated and explored is then applied to project regime access and limit avoidance for the SPARC primary reference discharge parameters.