The GBS code for the self-consistent simulation of plasma turbulence and kinetic neutral dynamics in the tokamak boundary

TL;DR

The paper introduces an enhanced GBS simulation code that models plasma turbulence and neutral dynamics in the entire tokamak volume with flexible magnetic configurations, improved computational efficiency, and validated results including a TCV discharge simulation.

Contribution

The new GBS version extends the simulation domain, incorporates flexible magnetic configurations, implements a faster iterative solver, and enhances parallelization for better performance.

Findings

Extended simulation domain to include core-edge-SOL interplay.

Achieved significant speed-up with iterative solvers over direct methods.

Successfully simulated a TCV tokamak discharge.

Abstract

A new version of GBS (Ricci et al. Plasma Phys. Control. Fusion 54, 124047, 2012; Halpern et al. J. Comput. Phys. 315, 388-408, 2016; Paruta et al. Phys. Plasmas 25, 112301, 2018) is described. GBS is a three-dimensional, flux-driven, global, two-fluid turbulence code developed for the self-consistent simulation of plasma turbulence and kinetic neutral dynamics in the tokamak boundary. In the new version presented here, the simulation domain is extended to encompass the whole plasma volume, avoiding an artificial boundary with the core, hence retaining the core-edge-SOL interplay. A new toroidal coordinate system is introduced to increase the code flexibility, allowing for the simulation of arbitrary magnetic configurations (e.g. single-null, double-null and snowflake configurations), which can also be the result of the equilibrium reconstruction of an experimental discharge. The implementation of a new iterative solver for the Poisson and Amp\`ere equations is presented, leading to a remarkable speed-up of the code with respect to the use of direct solvers, therefore allowing for efficient electromagnetic simulations that avoid the use of the Boussinesq approximation. The self-consistent kinetic neutral model, initially developed for limited configurations, is ported to the magnetic configurations considered by the present version of GBS and carefully optimized. A new MPI parallelisation is implemented to evolve the plasma and neutral models in parallel, thus improving the code scalability. The numerical implementation of the plasma and neutral models is verified by means of the method of manufactured solutions. As an example of the simulation capabilities of the new version of GBS, a simulation of a TCV tokamak discharge is presented.