Gyrokinetic Simulation of Energetic Particles
The SciDAC-3 project Gyrokinetic Simulation of Energetic Particles (GSEP) was started in 2008 as cooperative agreement research with General Atomic (GA), University of California at Irvine (UCI lead institution), Oak Ridge National Laboratory, and others. This Technical Report covers only the GA work. Understanding and predicting the confinement of energetic particles (EPs) in fusion plasma is an important research area vital to the successful operation of fusion reactors. Exterior Neutral Beam Injected (NBI) EPs are a primary source of plasma heating in current tokamaks like DIII-D. Internally generated fusion alpha particles in burning plasmas like ITER must be nearly all confined to heat the plasma to fusion conditions without much exterior heating. The thermal hydrogen plasma suffer energy confinement loses from high-n micro-turbulent driven transport. Electrostatic microturbulence has a small but not negligible impart on EP confinement loss. However is has been known since the early 90’ that NBI and particularly fusion alphas EPs can drive low-n MHD-like low-n electromagnetic Alfven Eigenmodes (AEs) unstable with the possibility significant loss of plasma heating. Even without significant heating loss, low levels of high energy escaping alphas can damage and erode the tokamak wall. The main purpose and goal of the GSEP project was to apply GA’s GYRO and UCI’s GTC gyrokinetic codes to delineate the linear stability boundaries of the EP driven AE modes in the most physically accurate way going beyond the kinetic EP perturbed MHD codes of the 90’s. Most importantly the nonlinear gyrokinetic simulations were to treat the EP driven low-n AE modes embedded in the high-n microturbulence to actually quantify the EP transport loss of confinement. The gyrokinetic codes GYRO, GTC (and GEM) where to be verified and cross-compared with simpler gyro-fluid (TAEFL) and kinetic EP perturbed MHD codes (NOVA-K) to project ITER performance, as well as validated against DIII-D experiments on AE driven NBI confinement losses. The following documents and briefly summarizes the published GA-GSEP contributions to the success of this 9-year project. The work was done by the R.E. Waltz (GA-GSEP PI) in collaboration with UCSD Research Scientist E. M. Bass also supported under this grant. Some of the later work was done in collaboration with a Peking University (Beijing, China) graduate student He Sheng (not funded by this grant).