Global gyrokinetic nonlinear simulations of kinetic infernal modes in reversed shear tokamaks

Abstract

The nonlinear evolution of electromagnetic instabilities in reversed shear plasmas is investigated by means of global gyrokinetic simulations. It is found that the kinetic infernal mode (KIM), which is a pressure-driven instability with low to intermediate toroidal mode number excited in a region of low magnetic shear, is unstable at high β, while the ion temperature gradient mode is unstable at low β, where β is the ratio of the plasma kinetic pressure to the magnetic pressure. The β threshold of the KIM is much lower than that of the kinetic ballooning mode (KBM) appearing in a normal shear plasma, while both the KIM and KBM are strong at the unfavorable curvature region, and the KIM has the same parity as the KBM. Nonlinear simulations show that the KIM gets saturated by exciting strong zonal flows and fluctuations of low toroidal mode number. The amplitude of the KIM turbulence is similar to that of the KBM turbulence in spite of the fact that the linear growth rate of the KIM is much higher than that of the KBM. This is because the excitation of zonal flows and fluctuations at low toroidal mode number is stronger in the reversed shear plasma than that of the normal shear plasma. On the other hand, the energy flux and particle flux due to the KIM turbulence are about two or three times larger than those by the KBM turbulence.