General Gyrokinetic Equations for Edge Plasmas

Abstract

During the pedestal cycle of H-mode edge plasmas in tokamak experiments, large-amplitude pedestal build-up and destruction coexist with small-amplitude drift wave turbulence. The pedestal dynamics simultaneously includes fast time-scale electromagnetic instabilities, long time-scale turbulence-induced transport processes, and more interestingly the interaction between them. To numerically simulate the pedestal dynamics from first principles, it is desirable to develop an effective algorithm based on the gyrokinetic theory. However, existing gyrokinetic theories cannot treat fully nonlinear electromagnetic perturbations with multi-scale-length structures in spacetime, and therefore do not apply to edge plasmas. A set of generalized gyrokinetic equations valid for the edge plasmas has been derived. This formalism allows large-amplitude, time-dependent background electromagnetic fields to be developed fully nonlinearly in addition to small-amplitude, short-wavelength electromagnetic perturbations. It turns out that the most general gyrokinetic theory can be geometrically formulated. The Poincaré-Cartan-Einstein 1-form on the 7D phase space determines particles' worldlines in the phase space, and realizes the momentum integrals in kinetic theory as fiber integrals. The infinitesimal generator of the gyro-symmetry is then asymptotically constructed as the base for the gyrophase coordinate of the gyrocenter coordinate system. This is accomplished by applying the Lie coordinate perturbation method to the Poincaré-Cartan-Einstein 1-form. General gyrokinetic Vlasov-Maxwell equations are then developed as the Vlasov-Maxwell equations in the gyrocenter coordinate system, rather than a set of new equations. Because the general gyrokinetic system developed is geometrically the same as the Vlasov-Maxwell equations, all the coordinate-independent properties of the Vlasov-Maxwell equations, such as energy conservation, momentum conservation, and phase space volume conservation, are automatically carried over to the general gyrokinetic system. The pullback transformation associated with the coordinate transformation is shown to be an indispensable part of the general gyrokinetic Vlasov-Maxwell equations. As an example, the pullback transformation in the gyrokinetic Poisson equation is explicitly expressed in terms of moments of the gyrocenter distribution function, with the important gyro-orbit squeezing effect due to the large electric field shearing in the edge and the full finite Larmour radius effect for short wavelength fluctuations. The familiar “polarization drift density” in the gyrocenter Poisson equation is replaced by a more general expression. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)