Linear gyrokinetic simulation of high-n toroidal Alfven eigenmodes in a burning plasma

A hybrid gyrokinetic ions/massless fluid electron model is used to study the stability of high-n toroidal Alfven eigenmodes (TAEs) in ITER [M. Shimada et al., Nucl. Fusion 47, S1 (2007)]. The hybrid model has been implemented in the particle-in-cell turbulence simulation code GEM[Y. Chen and S. E. Parker, J. Comput. Phys. 220, 839 (2007)]. The adequacy of the hybrid model for simulating TAEs has been previously demonstrated [J. Lang et al., Phys. Plasmas 16, 102101 (2009)] by comparing the simulated TAE mode frequency and structure with an eigenmode analysis, and the thermal ion kinetic damping effect with analytic theory. By using a global particle-in-cell code the effects of large orbit width and nonlocal mode structures can be accurately included. Damping rate due to numerical filtering is carefully monitored, and convergence with respect to particle number, grid resolution, etc., is thoroughly tested. The simulations show that the most unstable modes in ITER lie in the rage of 10<n<20. Thermal ion pressure effect and alpha particle nonperturbative effect are important in determining the mode radial location and stability threshold. The thermal ion Landau damping rate and radiative damping rate from the simulations are compared with analytical estimates. The thermal ion Landau damping is the dominant damping mechanism. Plasma elongation has a strong stabilizing effect on the alpha driven TAEs. The central alpha particle pressure threshold for the most unstable n=15 mode is about βα(0)=0.7% for the fully shaped ITER equilibrium.