Search

Global simulations of electromagnetic turbulence, collisionless tearing modes, and Alfven Eigenmodes in the presence of fast particles are carried out using the gyrokinetic particle-in-cell codes ORB5 (E. Lanti et al, Comp. Phys. Comm., ${\bf 251}$, 107072 (2020)) and EUTERPE (V. Kornilov et al, Phys. Plasmas, ${\bf 11}$, 3196 (2004)) in tokamak and stellarator geometries. Computational feasibility of simulating such complex coupled systems is demonstrated.

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.

Turbulence in tokamaks generates radially sheared zonal flows. Their oscillatory counterparts, geodesic acoustic modes (GAMs), appear due to the action of the magnetic field curvature. The GAMs can be driven unstable by an anisotropic energetic particle (EP) population leading to the formation of global radial structures, called EGAMs. The EGAMs can redistribute EP energy to the bulk plasma through collisionless wave-particle interaction. In such a way, the EGAMs might contribute to the plasma heating. Thus, investigation of EGAM properties, especially in the velocity space, is necessary for precise understanding of the transport phenomena in tokamak plasmas. In this work, the nonlinear dynamics of EGAMs without considering the mode interaction with the turbulence is investigated with the help of a Mode-Particle-Resonance (MPR) diagnostic implemented in the global gyrokinetic particle-in-cell code ORB5. An ASDEX Upgrade discharge is chosen as a reference case for this investigation due to its rich EP nonlinear dynamics. An experimentally relevant magnetic field configuration, thermal species profiles and an EP density profile are taken for EGAM chirping modelling and its comparison with available empirical data. The same magnetic configuration is used to explore energy transfer by the mode from the energetic particles to the thermal plasma including kinetic electron effects. For a given EGAM level the plasma heating by the mode can be significantly enhanced by varying the EP parameters. Electron dynamics decreases the EGAM saturation amplitude and consequently reduces the plasma heating, even though the mode transfers its energy to thermal ions much more than to electrons.

This paper presents the current state of the global gyrokinetic code ORB5 as an update of the previous reference [Jolliet et al., Comp. Phys. Commun. 177 409 (2007)]. The ORB5 code solves the electromagnetic Vlasov-Maxwell system of equations using a PIC scheme and also includes collisions and strong flows. The code assumes multiple gyrokinetic ion species at all wavelengths for the polarization density and drift-kinetic electrons. Variants of the physical model can be selected for electrons such as assuming an adiabatic response or a ``hybrid'' model in which passing electrons are assumed adiabatic and trapped electrons are drift-kinetic. A Fourier filter as well as various control variates and noise reduction techniques enable simulations with good signal-to-noise ratios at a limited numerical cost. They are completed with different momentum and zonal flow-conserving heat sources allowing for temperature-gradient and flux-driven simulations. The code, which runs on both CPUs and GPUs, is well benchmarked against other similar codes and analytical predictions, and shows good scalability up to thousands of nodes.

We present a portable platform, called PIC_ENGINE, for accelerating Particle-In-Cell (PIC) codes on heterogeneous many-core architectures such as Graphic Processing Units (GPUs). The aim of this development is efficient simulations on future exascale systems by allowing different parallelization strategies depending on the application problem and the specific architecture. To this end, this platform contains the basic steps of the PIC algorithm and has been designed as a test bed for different algorithmic options and data structures. Among the architectures that this engine can explore, particular attention is given here to systems equipped with GPUs. The study demonstrates that our portable PIC implementation based on the OpenACC programming model can achieve performance closely matching theoretical predictions. Using the Cray XC30 system, Piz Daint, at the Swiss National Supercomputing Centre (CSCS), we show that PIC_ENGINE running on an NVIDIA Kepler K20X GPU can outperform the one on an Intel Sandybridge 8-core CPU by a factor of 3.4.