Energetic Particle Mode Dynamics in Tokamaks
F. Zonca,S. Briguglio,L. Chen,G. Fogaccia,G. Vlad,L.-J. Zheng
2000-01-01
Abstract:Energetic Particle Modes [1] (EPM) are strongly driven oscillationsexcited via wave-particle resonant interactions at the characteristic frequencies of the energetic ions [1, 2], ωtE, ωBE and/orωdE , i.e., respectively the transit frequency for circulating particles and the bounce and precessional drift frequencies for trapped ions. A sharp transition in the plasma stability at the critical EPM excitation threshold has been observed by nonperturbative gyrokinetic codes in terms of changes in normalized growth rate ( γ/ωA, with ωA = vA/qR0), real frequency (ωr/ωA) and parallel wave vector ( k‖qR0) both asα = −R0qβ [3, 4] of the thermal plasma and that, αE [3, 5, 6], of fast ions are varied. The present work further explores theoretical aspects of EPM excitations by spatially localized particle sources, possibly associated with frequency chirping, which can radially trap the EPM in the region where the free energy source is strongest. Results of a nonperturbative3D Hybrid MHD Gyrokinetic code [5] are also presented to emphasize that nonlinear behaviors of EPM’s are different from those of Toroidal Alfv ́ en Eigenmodes (TAE) [7] and Kinetic TAE (KTAE) [8] and that particle losses and mode saturation are consistent with themode-particle pumpingmodel [9] (particle radial convection ). Results of theoretical analyses of nonlinear EPM dynamics are also presented and the possible overlap with more general nonlinear dynamics problems is discussed. 1. EPM Excited by Ion Cyclotron Radio Frequency (ICRF) Heating. A one-dimensional fast particle simulation code (FAPS-1D) has been developed to study the trapped energetic particle effects on the short wavelength (high toroidal mode number n) Alfv énic modes in the tokamak configuration. The code employs the gyrokinetic (GKE) magnetohydrodynamic (MHD) hybrid simulation scheme. Thus, the core plasma is described by the MHD fluid equations, while the energetic particle species are described by the linearized GKE equations. This allows us to study in detail the kinetic effects of the energetic particles, such as wave-particle resonance and finite orbit size effects. The MHD fluid (or vorticity) equation is solved via a predictor corrector algorithm. The GKE equation for the energetic particles is solved by aδf particle-in-cell simulation method. Specially, simulations are performed to study the interactions between the Alfv ́ en modes and the trapped energetic ions produced by the ICRF heating [3]. It is found that for the trapped energetic ions the wave-particle resonances are mainly due to the precessional drift resonance, as indicated by the fact that the real part of the mode frequency Ωr = ωr/ωA increases withkθρA, ρA = vA/ωcE , for smallkθρA. This is further confirmed by examining the phase diagram of the averaged nonadiabatic distribution function, ∑ gi, which shows that there are two resonances (180-degree-phase shift): one peaks at a lower resonance energy (normalized to the effective ICRF tail ion thermal energy) E l rs ≈ 0.75 and the other has a higher value, E rs ≈ 3. It can be demonstrated that the lower one corresponds to the bounce resonance, while the higher one to the precessional drift resonance. Both the TAEs and energetic particle modes (EPMs) are found to be excited by the trapped energetic particles. Transitions between various modes are observed in the simulations, as shown 0.12 0.10 0.08 0.06 0.04 0.02 0.00 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.14