Second order nonlinear gyrokinetic theory : From the particle to the gyrocenter

A gyrokinetic reduction is based on a specific ordering of the different small parameters characterizing the background magnetic field and the fluctuating electromagnetic fields. In this tutorial, we consider the following ordering of the small parameters: $\epsilon\_B=\epsilon\_\delta^2$ where $\epsilon\_B$ is the small parameter associated with spatial inhomogeneities of the background magnetic field and $\epsilon\_\delta$ characterizes the small amplitude of the fluctuating fields. In particular, we do not make any assumption on the amplitude of the background magnetic field. Given this choice of ordering, we describe a self-contained and systematic derivation which is particularly well suited for the gyrokinetic reduction, following a two-step procedure. We follow the approach developed in [Sugama, Physics of Plasmas 7, 466 (2000)]:In a first step, using a translation in velocity, we embed the transformation performed on the symplectic part of the gyrocentre reduction in the guiding-centre one. In a second step, using a canonical Lie transform, we eliminate the gyroangle dependence from the Hamiltonian. As a consequence, we explicitly derive the fully electromagnetic gyrokinetic equations at the second order in $\epsilon\_\delta$

Lagrangian and Hamiltonian constraints for guiding-center Hamiltonian theories

A consistent guiding-center Hamiltonian theory is derived by Lie-transform perturbation method, with terms up to second order in magnetic-field nonuniformity. Consistency is demonstrated by showing that the guiding-center transformation presented here satisfies separate Jacobian and Lagrangian constraints that have not been explored before. A new first-order term appearing in the guiding-center phase-space Lagrangian is identified through a calculation of the guiding-center polarization. It is shown that this new polarization term also yields a simpler expression of the guiding-center toroidal canonical momentum, which satisfies an exact conservation law in axisymmetric magnetic geometries. Lastly, an application of the guiding-center Lagrangian constraint on the guiding-center Hamiltonian yields a natural interpretation for its higher-order corrections.Comment: 11 page

Creation of a Transport Barrier for the E x B drift in magnetized plasmas

International audienceWe modelize the chaotic dynamics of charged test-particles in a turbulent electric field, across the confining magnetic field in controlled thermonuclear fusion devices by a 1.5 degrees of freedom Hamiltonian dynamical system. The external electric field E is given by a some potential V and the magnetic field B is considered uniform. We prove that, by introducing a small additive control term to the external electric field, it is possible to create a transport barrier. The robustness of this control method is also numerically investigated

Verification of Gyrokinetic codes: theoretical background and applications

In fusion plasmas the strong magnetic field allows the fast gyro-motion to be systematically removed from the description of the dynamics, resulting in a considerable model simplification and gain of computational time. Nowadays, the gyrokinetic (GK) codes play a major role in the understanding of the development and the saturation of turbulence and in the prediction of the subsequent transport. Naturally, these codes require thorough verification and validation. Here we present a new and generic theoretical framework and specific numerical applications to test the faithfulness of the implemented models to theory and to verify the domain of applicability of existing GK codes. For a sound verification process, the underlying theoretical GK model and the numerical scheme must be considered at the same time, which has rarely been done and therefore makes this approach pioneering. At the analytical level, the main novelty consists in using advanced mathematical tools such as variational formulation of dynamics for systematization of basic GK code's equations to access the limits of their applicability. The verification of numerical scheme is proposed via the benchmark effort. In this work, specific examples of code verification are presented for two GK codes: the multi-species electromagnetic ORB5 (PIC) and the radially global version of GENE (Eulerian). The proposed methodology can be applied to any existing GK code. We establish a hierarchy of reduced GK Vlasov-Maxwell equations implemented in the ORB5 and GENE codes using the Lagrangian variational formulation. At the computational level, detailed verifications of global electromagnetic test cases developed from the CYCLONE Base Case are considered, including a parametric $\beta$-scan covering the transition from ITG to KBM and the spectral properties at the nominal $\beta$ value.Comment: 16 pages, 2 Figures, APS DPP 2016 invited pape

Transport barrier for the radial diffusion due to the ExB drift motion of guiding centers in cylindrical confinement geometry

13 pages, 9 figures, 1 columnInternational audienceWe consider the radial transport of test particles due to the ExB drift motion in the guiding center approximation. Using an explicit expression to modify the electrostatic potential, we show that it is possible to construct a transport barrier which suppresses radial transport. We propose an algorithm for the implementation of this local modification computed from an electrostatic potential known on a spatio-temporal grid. The number of particles which escape the inner region defined by the barrier measures the efficiency of the control. We show that the control is robust by showing a significant reduction of radial transport, when applied with a reduced number of probes aligned on a circle