Speaker
Description
In the pedestal, large bootstrap and Pfirsch-Schluter currents, arising from the steep pressure gradient, drive kink instabilities and enhance the drive for the peeling-ballooning modes. The latter triggers periodic plasma eruptions at the edge of a tokamak that significantly influence the pedestal properties, limiting the pedestal pressure gradient and degrading core confinement. The conventional gyrokinetic theory usually employed to describe the physics of the pedestal is not sufficient to recover the full spectrum of ideal MHD results, e.g. it cannot capture the physics of kink modes. Indeed, it adopts a Maxwellian equilibrium and therefore does not retain the current density gradient along the magnetic field lines. These current density gradient effects can typically be ignored in a core plasma of a conventional tokamak, but must be retained in the core of a spherical tokamak plasma or in the vicinity of the pedestal of a conventional tokamak where the pressure gradient is the strongest. We have developed a new electromagnetic gyrokinetic theory that captures the strong gradients in the pedestal. The particle distribution function is written in the form: $f=F + \delta f$, where $F$ describes the equilibrium piece (in the gyrokinetic ordering), i.e. it varies on transport time scales, and $\delta f$ is the fluctuating piece (the particle species indices are dropped for simplicity). To recover ideal MHD, we must retain neoclassical effects in the equilibrium distribution function and thus solve perturbatively the drift kinetic equation that describes $F$, expanding it in powers of $\delta = \rho_{\vartheta}/L \ll 1$, where $\rho_{\vartheta}$ is the poloidal Larmor radius and $L$ is the characteristic size of the system. Here the leading order corresponds to the Maxwellian, employed in conventional gyrokinetics. The following order introduces the neoclassical flows (e.g. bootstrap current flows), required to recover kink physics. The gyro-angle dependent piece of the second order correction in a $\delta$ expansion must be retained as well to ensure a consistent ordering. This $F$ is then employed to obtain the extended gyrokinetic equation for fluctuations. To solve the gyrokinetic equation for $\delta f$, we employ an expansion in $\Delta = \rho/L$, where $\rho$ is the particle Larmor radius. We provide two sets of equations: (a) global, integro-differential gyrokinetic equations that allow one to capture the global properties of modes and (b) local gyrokinetic equations which employ an eikonal ansatz for the fluctuating electromagnetic field. The latter provides the theoretical formalism in a form ready for implementing in existing local gyrokinetic codes.