Helical Core Formation and MHD Stability in ITER-Scale Plasmas with Fusion-born Alpha Particles
Abstract
The effect of fusion-born alpha particles on the helical core (HC), a long-lived ideal saturation state of the $m/n=1/1$ kink/quasi-interchange mode, is studied in the ITER-scale hybrid scenario where a core plasma has a low magnetic shear $q\gtrsim1$. The HC state is determined by 3-D MHD force balance and all factors that contribute to it, such as plasma shaping, the safety factor profile, and the pressure profiles of all particle species. An incomplete but useful measure of the HC is the displacement of the magnetic axis, $\delta_\mathrm{HC}$. Using MHD-PIC simulations, we find that $\delta_\mathrm{HC}$ is enhanced by increasing alpha particle pressure $\beta_\mathrm{\alpha}$. Within the ITER operating alpha pressure $\beta_\mathrm{\alpha}(0) \lesssim 1\%$, $\beta_\mathrm{\alpha}$ can be approximately treated as part of the total MHD pressure. In this regime, there is no notable flattening of the pressure profile, indicating that the HC preserves the omnigenity of the plasma. If one increases $\beta_\mathrm{\alpha}(0)$ beyond $1\%$, $\delta_\mathrm{HC}$ continues to increase with $\beta_\mathrm{\alpha}$ until it reaches an upper limit at $\beta_\mathrm{\alpha}(0)=3\%$ for our reference case. At this limit, both the bulk and alpha pressure profiles are partially flattened, indicating a reduction in omnigenity. After HC formation, a resistive pressure-driven MHD mode can become unstable, which is localized along the compressed magnetic flux region of the HC. This secondary mode consists of a broad spectrum of short-wavelength Fourier components that grow at same rates and are thus part of a single coherent entity. Our present simulation model is insufficient to adequately represent such a secondary mode; however, preliminary results suggest that it can facilitate magnetic chaos, which affects plasma confinement.