Targeting Realistic Geometry in Tokamak Code Gysela^★

¹ Inria, F-54600 Villers-lèes-Nancy, France, hosted by CEA/IRFM, F-13108 Saint-Paul-lez-Durance, France
² Max Planck Institute for Plasma Physics and Technische Universität München, D-85748, Germany
³ CEA/IRFM, F-13108 Saint-Paul-lez-Durance, France
⁴ IRMA, Université de Strasbourg, 7, rue René Descartes, F-67084 Strasbourg & INRIA-Nancy Grand-Est, projet TONUS

Abstract

In magnetically confined plasmas used in Tokamak, turbulence is respon-sible for specific transport that limits the performance of this kind of reactors. Gyroki-netic simulations are able to capture ion and electron turbulence that give rise to heat losses, but require also state-of-the-art HPC techniques to handle computation costs. Such simulations are a major tool to establish good operating regime in Tokamak such as ITER, which is currently being built. Some of the key issues to address more re- alistic gyrokinetic simulations are: efficient and robust numerical schemes, accurate geometric description, good parallelization algorithms. The framework of this work is the Semi-Lagrangian setting for solving the gyrokinetic Vlasov equation and the Gy-sela code. In this paper, a new variant for the interpolation method is proposed that can handle the mesh singularity in the poloidal plane at r = 0 (polar system is used for the moment in Gysela). A non-uniform meshing of the poloidal plane is proposed instead of uniform one in order to save memory and computations. The interpolation method, the gyroaverage operator, and the Poisson solver are revised in order to cope with non-uniform meshes. A mapping that establish a bijection from polar coordinates to more realistic plasma shape is used to improve realism. Convergence studies are provided to establish the validity and robustness of our new approach.