Some numerical aspects of the conservative PSM scheme in a 4D drift-kinetic code

Jean-Philippe Braeunig (1, 2) , Nicolas Crouseilles (1) , Virginie Grandgirard (3) , Guillaume Latu (1, 3) , Michel Mehrenberger (1, 4) , Eric Sonnendrücker (1, 4)

1 CALVI - Scientific computation and visualization
2 DAM/DIF - DAM Île-de-France
3 IRFM - Institut de Recherche sur la Fusion par confinement Magnétique
4 IRMA - Institut de Recherche Mathématique Avancée

The purpose of this work is simulation of magnetised plasmas in the ITER pro ject framework. In this context, Vlasov-Poisson like models are used to simulate core turbulence in the tokamak in a toroidal geometry. This leads to heavy simulation because a 6D dimensional problem has to be solved, 3D in space and 3D in velocity. The model is reduced to a 5D gyrokinetic model, taking advantage of the particular motion of particles due to the presence of a strong magnetic field. However, accurate schemes, parallel algorithms need to be designed to bear these simulations. This paper describes the numerical studies to improve robustness of the conservative PSM scheme in the context of its development in the GYSELA code. In this paper, we only consider the 4D drift-kinetic model which is the backbone of the 5D gyrokinetic models and relevant to build a robust and accurate numerical method.

Mots clés: numerical simulation; conservative scheme; ITER; plasma turbulence

Domaines: Analyse numérique [cs.NA]

Jean-Philippe Braeunig, Nicolas Crouseilles, Virginie Grandgirard, Guillaume Latu, Michel Mehrenberger, et al.. Some numerical aspects of the conservative PSM scheme in a 4D drift-kinetic code. [Research Report] RR-7109, INRIA. 2009, pp.37. ⟨inria-00435203v1⟩