The observation that fast ions stabilize ion-temperature-gradient-driven microturbulence has profound implications for future fusion reactors. It is also important in optimizing the performance of present-day devices. In this work, we examine in detail the phenomenology of fast ion stabilization and present a reduced model which describes this effect. This model is derived from the high-energy limit of the gyrokinetic equation and extends the existing "dilution" model to account for nontrivial fast ion kinetics. Our model provides a physically-transparent explanation for the observed stabilization and makes several key qualitative predictions. Firstly, that different classes of fast ions, depending on their radial density or temperature variation, have different stabilizing properties. Secondly, that zonal flows are an important ingredient in this effect precisely because the fast ion zonal response is negligible. Finally, that in the limit of highly-energetic fast ions, their response approaches that of the "dilution" model; in particular, alpha particles are expected to have little, if any, stabilizing effect on plasma turbulence. We support these conclusions through detailed linear and nonlinear gyrokinetic simulations.
arXiv:1801.00664v3 [physics.plasm-ph] 5 Jun 2018Effect of fast ions on microturbulence 2
Fast ion stabilization of ITG microturbulenceMicroturbulence fundamentally limits the confinement time in current and future tokamak experiments [2]. The very gradients that are required to achieve high central densities and temperatures also provide a source of free energy. This free energy drives transport of particles, momentum, and energy, usually far in excess of collisional transport. The ion temperature gradient (ITG) mode, for instance, limits the core temperature of tokamaks [3,4,5]. This turbulence occurs on the scale of the thermal ion Larmor radius ρ i .Gyrokinetics is the reduction of the Fokker-Planck kinetic equation that rigorously handles electromagnetic fields whose fluctuations vary on spatial scales similar to ρ i , but on timescales much slower than the gyro-frequency Ω i [6,7,8]. A multitude of computational tools have been developed to solve the nonlinear gyrokinetic equation [9,10,11,12]. These tools take as inputs the equilibrium magnetic field, density, and temperature profiles, and predict the ensuing turbulent fluctuations and the associated transport fluxes. By using experimentally-determined profiles and comparing the computed fluxes to experimentally-inferred fluxes (usually from power-balance calculations), and through more detailed comparisons of turbulence characteristics, the validity of the gyrokinetic approach can be verified [13]. This has been widely demonstrated [14,15,16,17].However, due to the scarcity of computational resources, these matching exercises of necessity entail the use of simulations that neglect various parts of the complex physics of the experimental setup. It was discovered during one of these exercises, in which experiments on the Joint European Tor...