Xavier Garbet

Gyrokinetic simulations of turbulent transport

This overview is an assessment of the gyrokinetic framework and simulations to compute turbulent transport in fusion plasmas. It covers an introduction to the gyrokinetic theory, the principal numerical techniques which are being used to solve the gyrokinetic equations, fundamentals in gyrokinetic turbulence and the main results which have been brought by simulations with regard to transport in fusion devices and fluctuation measurements.

A drift-kinetic semi-Lagrangian 4D code for ion turbulence simulation

A new code is presented here, named Gyrokinetic SEmi-LAgragian (GYSELA) code, which solves 4D drift-kinetic equations for ion temperature gradient driven turbulence in a cylinder (r,θ,z). The code validation is performed with the slab ITG mode that only depends on the parallel velocity. This code uses a semi-Lagrangian numerical scheme, which exhibits good properties of energy conservation in non-linear regime as well as an accurate description of fine spatial scales. The code has been validated in the linear and non-linear regimes. The GYSELA code is found to be stable over long simulation times (more than 20 times the linear growth rate of the most unstable mode), including for cases with a high resolution mesh (δr ∼ 0.1 Larmor radius, δz ∼ 10 Larmor radius).

Curvature effects on the dynamics of tearing modes in tokamaks

The curvature effects on the dynamics of magnetic island evolution in tokamaks are investigated both theoretically and numerically. By taking into account perpendicular and parallel heat diffusion, a new dispersion relation is derived for tearing modes that match the linear and nonlinear results. This evolution equation allows a quantitative description over the whole range of island sizes. It predicts a nonlinear instability, i.e., growing magnetic islands in linearly stable magnetic configurations. All these predictions are in excellent agreement with full tridimensional linear and nonlinear magnetohydrodynamic (MHD) computations with the latest version of XTOR [K. Lerbinger and J. F. Luciani, J. Comput. Phys. 97, 444 (1991)]. These results have important consequences on the onset of neoclassical tearing modes because they predict a resistive MHD threshold.

Large scale dynamics in flux driven gyrokinetic turbulence

The turbulent transport governed by the toroidal ion temperature gradient driven instability is analysed with the full-f global gyrokinetic code GYSELA (Grandgirard et al 2007 Plasma Phys. Control. Fusion 49 B173) when the system is driven by a prescribed heat source. Weak, yet finite, collisionality governs a neoclassical ion heat flux that can compete with the turbulent driven transport. In turn, the ratio of turbulent to neoclassical transport increases with the source magnitude, resulting in the degradation of confinement with additional power. The turbulent flux exhibits avalanche-like events, characterized by intermittent outbursts which propagate ballistically roughly at the diamagnetic velocity. Locally, the temperature gradient can drop well below the linear stability threshold. Large outbursts are found to correlate with streamer-like structures of the convection cells albeit their Fourier spectrum departs significantly from that of the most unstable linear modes. Last, the poloidal rotation of turbulent eddies is essentially governed by the radial electric field at moderate density gradient.

Trapped-ion driven turbulence in tokamak plasmas

Ion turbulence is expected to play an important role in anomalous transport. Thus an original approach is proposed here to study trapped-ion instability. A Vlasov code is used to determine the behaviour of the instability near the threshold and compared with analytical solutions of the Vlasov equation. Some interesting features which appear in the nonlinear regime are discussed.

TOKAM-3D: A 3D fluid code for transport and turbulence in the edge plasma of Tokamaks

We present a new code aiming at giving a global and coherent approach for transport and turbulence issues in the edge plasma of Tokamaks. The TOKAM-3D code solves 3D fluid drift equations in full-torus geometry including both closed field lines and SOL physics. No scale separation is assumed so that interactions between large scale flows and turbulence are coherently treated. Moreover, the code can be run in transport regimes ranging from purely anomalous diffusion to fully established turbulence. Specific numerical schemes have been developed which can solve the model equations whether the presence of a limiter in the plasma is taken into account or not. Example cases giving an overview of the field of application of the code as well as verification results are also presented.

Gyrokinetic modeling: A multi-water-bag approach

Predicting turbulent transport in nearly collisionless fusion plasmas requires one to solve kinetic (or, more precisely, gyrokinetic) equations. In spite of considerable progress, several pending issues remain; although more accurate, the kinetic calculation of turbulent transport is much more demanding in computer resources than fluid simulations. An alternative approach is based on a water-bag representation of the distribution function that is not an approximation but rather a special class of initial conditions, allowing one to reduce the full kinetic Vlasov equation into a set of hydrodynamic equations while keeping its kinetic character. The main result for the water-bag model is a lower cost in the parallel velocity direction since no differential operator associated with some approximate numerical scheme has to be carried out on this variable v∥. Indeed, a small bag number is sufficient to correctly describe the ion temperature gradient instability.

Nonlinear dynamics of magnetic islands imbedded in small-scale turbulence.

The nonlinear dynamics of magnetic tearing islands imbedded in a pressure gradient driven turbulence is investigated numerically in a reduced magnetohydrodynamic model. The study reveals regimes where the linear and nonlinear phases of the tearing instability are controlled by the properties of the pressure gradient. In these regimes, the interplay between the pressure and the magnetic flux determines the dynamics of the saturated state. A secondary instability can occur and strongly modify the magnetic island dynamics by triggering a poloidal rotation. It is shown that the complex nonlinear interaction between the islands and turbulence is nonlocal and involves small scales.

A 5D gyrokinetic full-f global semi-lagrangian code for flux-driven ion turbulence simulations

This paper addresses non-linear gyrokinetic simulations of ion temperature gradient (ITG) turbulence in tokamak plasmas. The electrostatic Gysela code is one of the few international 5D gyrokinetic codes able to perform global, full-f and flux-driven simulations. Its has also the numerical originality of being based on a semi-Lagrangian (SL) method. This reference paper for the Gysela code presents a complete description of its multi-ion species version including: (i) numerical scheme, (ii) high level of parallelism up to 500k cores and (iii) conservation law properties.

Neoclassical physics in full distribution function gyrokinetics

Treatment of binary Coulomb collisions when the full gyrokinetic distribution function is evolved is discussed here. A spectrum of different collision operators is presented, differing through both the physics that can be addressed and the numerics they are based on. Eulerian-like (semi-Lagrangian) and particle in cell (PIC) (Monte-Carlo) schemes are successfully cross-compared, and a detailed confrontation to neoclassical theory is shown.