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.

The geodesic acoustic mode (GAM) is a coherently oscillating zonal flow that may regulate turbulence in toroidal plasmas. Uniquely, the complete poloidal and toroidal structure of the magnetic component of the turbulence-driven GAM has been mapped in the TCV tokamak. Radially localized measurements of the fluctuating density, ECE radiative temperature and poloidal flow show that the GAM is a fully coherent, radially propagating wave. These observations are consistent with electrostatic, gyrokinetic simulations.

The response of passing electrons in ion temperature gradient (ITG) and trapped electron mode (TEM) microturbulence regimes is investigated in tokamak geometry making use of the flux-tube version of the gyrokinetic code GENE. Results are obtained using two different electron models, fully kinetic and hybrid in which passing particles are forced to respond adiabatically while trapped are handled kinetically. Comparing linear eigenmodes obtained with these two models enables to systematically isolate fine radial structures located at corresponding mode rational surfaces, clearly resulting from the non-adiabatic passing electron response. Non-linear simulations show that these fine structures on the non-axisymmetric modes survive in the turbulent phase. Furthermore, through non-linear coupling to axisymmetric modes, they induce radial modulations in the effective profiles of density, ion/electron temperature and E × B shearing rate. Finally, the passing electron channel is shown to significantly contribute to the transport levels, at least in our ITG case. Also shown is that the passing electrons significantly influence the E × B saturation mechanism of turbulence fluxes.

This document is intended for publication in the open literature. It is made available on the clear understanding that it may not be further circulated and extracts or references may not be published prior to publication of the original when applicable, or without the consent of the Publications Officer, EUROfusion Programme Management Unit,

Reducing the uncertainty on physical input parameters derived from experimental measurements is essential towards improving the reliability of gyrokinetic turbulence simulations. This can be achieved by introducing physical constraints. Amongst them, the zero particle flux condition is considered here. A first attempt is also made to match as well the experimental ion/electron heat flux ratio. This procedure is applied to the analysis of a particular Tokamak à Configuration Variable discharge. A detailed reconstruction of the zero particle flux hyper-surface in the multi-dimensional physical parameter space at fixed time of the discharge is presented, including the effect of carbon as the main impurity. Both collisionless and collisional regimes are considered. Hyper-surface points within the experimental error bars are found. The analysis is done performing gyrokinetic simulations with the local version of the GENE code, computing the fluxes with a Quasi-Linear (QL) model and validating the QL results with non-linear simulations in a subset of cases.