Use of reconstructed 3D equilibria to determine onset conditions of helical cores in tokamaks for extrapolation to ITER
The threshold for spontaneous symmetry breaking, which results in growth of m/n = 1/1 helical cores, in tokamaks is determined using VMEC. This model is based on a DIII-D hybrid discharge with helical core, and it predicts ITER (15MA scenario) to operate far in the helical core formation regime. Helical cores can only exist in tokamak discharges with monotonic but low, or reversed q-shear and qmin 1 in the core. The helical core is a saturated internal kink mode; its onset threshold is proportional to (dp/d)/B 2 t around q = 1. Below the threshold, applied 3-D fields can drive a helical core to finite size, as in DIIID. Above it, a random 3-D kick causes a bifurcation from axisymmetry and excites a spontaneous helical core. Helical cores occur frequently in C-Mod due to impurities; modeling shows the axisymmetric equilibrium becomes unstable at the end of the current ramp up due to the slow current penetration and a resultant reversed shear q-profile, which lowers the threshold. Lastly, a comparison of the helical core onset threshold for discharges from DIII-D, C-Mod and ITER confirms that while DIII-D is marginally stable, C-Mod and especially ITER are highly susceptible to helical core formation.