Defining an equilibrium state in global full-f gyrokinetic models

Guilhem Dif-Pradalier (1) , Virginie Grandgirard (1) , Yanick Sarazin (1) , Xavier Garbet (1) , Philippe Ghendrih (1)

1 IRFM - Institut de Recherche sur la Fusion par confinement Magnétique

This paper tackles the delicate choice of the initial distribution function in full-f gyrokinetic codes such as GYSELA 5D, aiming at predicting the turbulent transport level in low collisional tokamak plasmas. It is found, both analytically and numerically, that a Maxwellian distribution function with constant profiles on magnetic flux surfaces leads to the fast generation of a large scale electric field. Such a field opposes the up–down charge separation governed by the inhomogeneity of the equilibrium magnetic field. If large enough, the shearing rate induced by the resulting poloidal E Â B velocity could efficiently reduce the plasma micro-instabilities which account for the development of the turbulence. Starting in the ab initio code GYSELA 5D from an equilibrium distribution function depending on motion invariants only is shown to cure such a problem. In this case, charge separation is counterbalanced by parallel flow, and the standard fluid force balance is recovered.

Keywords: Tokamak plasma; Nonlinear gyrokinetic simulations; Collisionless equilibrium definition

Domains: Nonlinear Sciences [physics]; Plasma Physics [physics.plasm-ph]

Guilhem Dif-Pradalier, Virginie Grandgirard, Yanick Sarazin, Xavier Garbet, Philippe Ghendrih. Defining an equilibrium state in global full-f gyrokinetic models. Communications in Nonlinear Science and Numerical Simulation, 2008, 13, pp.65 - 71. ⟨10.1016/j.cnsns.2007.05.004⟩. ⟨cea-01468358⟩