Magnetohydrodynamic hybrid simulation of Alfven eigenmodes in Heliotron J, a low shear helical axis stellarator/heliotron

Alfven eigenmodes destabilized by energetic particles in Heliotron J, a low magnetic shear helical axis stellarator/heliotron with four toroidal equilibrium field periods, are investigated with MEGA, a hybrid simulation code for energetic particles interacting with a magnetohydrodynamic fluid. The objectives of this study are to validate the MEGA simulation results on Alfven eigenmodes in Heliotron J, and to clarify the properties of Alfven eigenmodes. The experimentally observed n/m = 2/4 global Alfven eigenmode (GAE) was reproduced with MEGA simulations based on the experimental bulk plasma temperature and density profiles. Beside n/m = 2/4 GAE, n/m = 1/2 GAE was also observed with lower amplitude and linear growth rate. The frequency and spatial location of n/m = 1/2 GAE are close to those of the experimentally observed n/m = 1/2 energetic particle mode (EPM). For n/m = 2/4 GAE, it was found that the three-dimensional spatial profile of the GAE was primarily composed of n = 2 harmonics while the contribution from the other toroidal mode numbers such as n= 2 +/- iN(fp) was weak, where 'i' and 'N-fp' are arbitrary integer and toroidal field period, respectively. This indicates that the coupling of harmonics with the same toroidal mode number through toroidicity is dominant for the GAE spatial profile, and the coupling with the different toroidal mode numbers through helicity and bumpiness has minor effects. For n/m = 1/2 GAE, the coupling of harmonic through toroidicity is weaker than that of n/m = 2/4 GAE. The non-linear evolutions of the n/m = 2/4 GAE and n/m = 1/2 GAE destabilized by energetic particles with a bump-on-tail distribution and a slowing-down distribution were compared. Weak dependence of the linear growth rate for both n/m = 2/4 and n/m = 1/2 modes on charge-exchange time were implied.