A general theoretical framework for investigating the nonlinear dynamics of phase space zonal structures is presented in this work. It is then, more specifically, applied to the limit where the nonlinear evolution time scale is smaller or comparable to the wave-particle trapping period. In this limit, both theoretical and numerical simulation studies show that nonadiabatic frequency chirping and phase locking could lead to secular resonant particle transport on meso-or macro-scales. The interplay between mode structures and resonant particles then provides the crucial ingredient to properly understand and analyze the nonlinear dynamics of Alfvén wave instabilities excited by nonperturbative energetic particles in burning fusion plasmas. Analogies with autoresonance in nonlinear dynamics and with superradiance in free-electron lasers are also briefly discussed.

The role of energetic particles (EPs) in fusion plasmas is unique as they could act as mediators of cross-scale couplings. More specifically, EPs can drive instabilities on the macro-and meso-scales and intermediate between the microscopic thermal ion Larmor radius and the macroscopic plasma equilibrium scale lengths. On one hand, EP driven shear Alfvén waves (SAWs) could provide a nonlinear feedback onto the macro-scale system via the interplay of plasma equilibrium and fusion reactivity profiles. On the other hand, EP-driven instabilities could also excite singular radial mode structures at SAW continuum resonances, which, by mode conversion, yield microscopic fluctuations that may propagate and be absorbed elsewhere, inducing nonlocal behaviors. The above observations thus suggest that a theoretical approach based on advanced kinetic treatment of both EPs and thermal plasma is more appropriate for burning fusion plasmas. Energetic particles, furthermore, may linearly and nonlinearly (via SAWs) excite zonal structures, acting, thereby, as generators of nonlinear equilibria that generally evolve on the same time scale of the underlying fluctuations. These issues are presented within a general theoretical framework, discussing evidence from both numerical simulation results and experimental observations. Analogies of fusion plasmas dynamics with problems in condensed matter physics, nonlinear dynamics, and accelerator physics are also emphasized.

We present a series of numerical simulation experiments set up to illustrate the fundamental physics processes underlying the nonlinear dynamics of Alfv enic modes resonantly excited by energetic particles in tokamak plasmas and of the ensuing energetic particle transports. These phenomena are investigated by following the evolution of a test particle population in the electromagnetic fields computed in self-consistent MHD-particle simulation performed by the HMGC code. Hamiltonian mapping techniques are used to extract and illustrate several features of wave-particle dynamics. The universal structure of resonant particle phase space near an isolated resonance is recovered and analyzed, showing that bounded orbits and untrapped trajectories, divided by the instantaneous separatrix, form phase space zonal structures, whose characteristic non-adiabatic evolution time is the same as the nonlinear time of the underlying fluctuations. Bounded orbits correspond to a net outward resonant particle flux, which produces a flattening and/or gradient inversion of the fast ion density profile around the peak of the linear wave-particle resonance. The connection of this phenomenon to the mode saturation is analyzed with reference to two different cases: a Toroidal Alfv en eigenmode in a low shear magnetic equilibrium and a weakly unstable energetic particle mode for stronger magnetic shear. It is shown that, in the former case, saturation is reached because of radial decoupling (resonant particle redistribution matching the mode radial width) and is characterized by a weak dependence of the mode amplitude on the growth rate. In the latter case, saturation is due to resonance detuning (resonant particle redistribution matching the resonance width) with a stronger dependence of the mode amplitude on the growth rate. [

It is shown, both analytically and by numerical simulations, that, in the presence of thermal ion kinetic effects, the beta induced Alfvén eigenmode (BAE)shear Alfvén wave continuous spectrum can be discretized into radially trapped eigenstates known as kinetic BAE (KBAE). While thermal ion compressibility gives rise to finite BAE accumulation point frequency, the discretization occurs via the finite Larmor radius and finite orbit width effects. Simulations and analytical theories agree both qualitatively and quantitatively. Simulations also demonstrate that KBAE can be readily excited by the finite radial gradients of energetic particles.

Adopting the theoretical framework for the generalized fishbonelike dispersion relation, an extended hybrid magnetohydrodynamics gyrokinetic simulation model has been derived analytically by taking into account both thermal ion compressibility and diamagnetic effects in addition to energetic particle kinetic behaviors. The extended model has been used for implementing an eXtended version of Hybrid Magnetohydrodynamics Gyrokinetic Code (XHMGC) to study thermal ion kinetic effects on Alfvénic modes driven by energetic particles, such as kinetic beta induced Alfvén eigenmodes in tokamak fusion plasmas.

Global gyrokinetic particle simulations of reversed shear Alfvén eigenmode ͑RSAE͒ have been successfully performed and verified. We have excited the RSAE by initial perturbation, by external antenna, and by energetic ions. The RSAE excitation by antenna provides verifications of the mode structure, the frequency, and the damping rate. When the kinetic effects of the background plasma are artificially suppressed, the mode amplitude shows a near-linear growth. With kinetic thermal ions, the mode amplitude eventually saturates due to the thermal ion damping. The damping rates measured from the antenna excitation and from the initial perturbation simulation agree very well. The RSAE excited by fast ions shows an exponential growth. The finite Larmor radius effects of the fast ions are found to significantly reduce the growth rate. With kinetic thermal ions and electron pressure, the mode frequency increases due to the elevation of the Alfvén continuum by the geodesic compressibility. The nonperturbative contributions from the fast ions and kinetic thermal ions modify the mode structure relative to the ideal magnetohydrodynamic ͑MHD͒ theory. The gyrokinetic simulations have been benchmarked with extended hybrid MHD-gyrokinetic simulations.

The Divertor Tokamak Test (DTT) facility is proposed for studying power exhaust solutions as well as integrated physics and technology aspects for the demonstration power plant (DEMO). To illuminate the richness of new novel plasma physics that can be explored in this device, linear stability properties and shear Alfvén fluctuation spectra of a typical DTT reference scenario are investigated by self-consistent hybrid magnetohydrodynamic-gyrokinetic simulations. The DTT core plasma can be divided into two regions, characterized by reverse shear Alfvén eigenmode in the central core, and by toroidal Alfvén eigenmode in the outer core region. The non-perturabtive effect of energetic particles (EPs), the wave-EP resonance condition as well as power transfer are analyzed in great detail, demonstrating the peculiar role played by EPs in multi-scale dynamics. The most unstable mode number of dominant Alfvénic fluctuations are shown to be of the order of 10, consistent with the typical orbit widths of the EPs normalized to the plasma minor radius and the DTT target design.

The beta-induced Alfvén eigenmode ͑BAE͒ in toroidal plasmas is studied using global gyrokinetic particle simulations. The BAE real frequency and damping rate measured in the initial perturbation simulation and in the antenna excitation simulation agree well with each other. The real frequency is slightly higher than the ideal magnetohydrodynamic ͑MHD͒ accumulation point frequency due to the kinetic effects of thermal ions. Simulations with energetic particle density gradient show exponential growth of BAE with a growth rate sensitive to the energetic particle temperature and density. The nonperturbative contributions by energetic particles modify the mode structure and reduce the frequency relative to the MHD theory. The finite Larmor radius effects of energetic particles reduce the BAE growth rate. Benchmarks between gyrokinetic particle simulation and hybrid MHD-gyrokinetic simulation show good agreement in BAE real frequency and mode structure.