Linear and nonlinear benchmark of gyrokinetic simulation of energetic particle driven toroidal Alfven eigenmodes in ITPA TAE benchmark case

A new gyrokinetic code, TEK, was benchmarked in simulating energetic particle (EP) driven toroidal Alfven eigenmodes (TAEs) in the simple tokamak configuration chosen by the ITPA-EP group for code benchmarking purposes. Linear benchmark has been well established by other codes, whereas nonlinear benchmark for this case is lacking. This paper presents, besides the linear benchmark, nonlinear results for both single-n and multiple-n simulations (n is the toroidal mode number). The nonlinear results are in good agreement with an analytical theory on zonal field beat-driven by Alfven eigenmodes, partially verifying the correctness of the nonlinear simulations. The saturation level and the resulting EP transport are examined. This provides data for future inter-code nonlinear benchmarking. In TEK, all species (electrons, thermal ions, EPs) are treated on the same footing using the gyrokinetic model (with electrons in the zero Larmor radius limit). The electromagnetic cancellation problem is mitigated by using the mixed-variable pull- back method. Numerical details related to electromagnetic gyrokinetic simulation are discussed.