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 nonlinear dynamics of beta-induced Alfvén eigenmodes (BAEs) driven by energetic particles (EPs) in the presence of ion-temperature-gradient turbulence is investigated, by means of selfconsistent global gyrokinetic simulations and analytical theory. A tokamak magnetic equilibrium with large aspect ratio and reversed shear is considered. A previous study of this configuration has shown that the electron species plays an important role in determining the nonlinear saturation level of a BAE in the absence of turbulence (Biancalani et al 2020 J. Plasma Phys.). Here, we extend the study to a turbulent plasma. The EPs are found modify the heat fluxes by introducing energy at the large spatial scales, mainly at the toroidal mode number of the dominant BAE and its harmonics. In this regime, BAEs are found to carry a strong electron heat flux. The feed-back of the global relaxation of the temperature profiles induced by the BAE, and on the turbulence dynamics, is also discussed.

The nonlinear dynamics of energetic-particle (EP) driven geodesic acoustic modes (EGAM) is investigated here. A numerical analysis with the global gyrokinetic particle-in-cell code ORB5 is performed, and the results are interpreted with the analytical theory, in close comparison with the theory of the beam-plasma instability. Only axisymmetric modes are considered, with a nonlinear dynamics determined by wave-particle interaction. Quadratic scalings of the saturated electric field with respect to the linear growth rate are found for the case of interest. As a main result, the formula for the saturation level is provided. Near the saturation, we observe a transition from adiabatic to non-adiabatic dynamics, i.e. the frequency chirping rate becomes comparable to the resonant EP bounce frequency. The numerical analysis is performed here with electrostatic simulations with circular flux surfaces, and kinetic effects of the electrons are neglected.

The linear destabilization and nonlinear saturation of energetic-particle driven Alfvénic instabilities in tokamaks strongly depend on the damping channels. In this work, the collisionless damping mechanisms of Alfvénic modes are investigated within a gyrokinetic framework, by means of global simulations with the particle-in-cell code ORB5, and compared with the eigenvalue code LIGKA and reduced models. In particular, the continuum damping and the Landau damping (of ions and electrons) are considered. The electron Landau damping is found to be dominant on the ion Landau damping for experimentally relevant cases. As an application, the linear and nonlinear dynamics of toroidicity induced Alfvén eigenmodes and energetic-particle driven modes in ASDEX Upgrade is investigated theoretically and compared with experimental measurements.

The linear properties of the geodesic acoustic modes (GAM) in tokamaks are investigated by means of the comparison of analytical theory and gyrokinetic numerical simulations. The dependence on the value of the safety factor, finite-orbit-width of the ions in relation to the radial mode width, magnetic-flux-surface shaping, and electron/ion mass ratio are considered. Nonuniformities in the plasma profiles (such as density, temperature, or safety factor), electro-magnetic effects, collisions and presence of minority species are neglected. Also, only linear simulations are considered, focusing on the local dynamics. We use three different gyrokinetic codes: the lagrangian (particle-in-cell) code ORB5, the eulerian code GENE and semi-lagrangian code GYSELA. One of the main aims of this paper is to provide a detailed comparison of the numerical results and analytical theory, in the regimes where this is possible. This helps understanding better the behavior of the linear GAM dynamics in these different regimes, the behavior of the codes, which is crucial in the view of a future work where more physics is present, and the regimes of validity of each specific analytical dispersion relation. 1 theories derived so far treat the m = 0 component of the electrons as adiabatic (whereas the m=0 component of the electron density perturbation is imposed to zero). We will refer to this model for treating the electrons as "adiabatic". The importance of having an analytical description is twofold. On the one hand, it allows a direct understanding of the physical mechanisms leading to the each different effect under investigation. On the other hand, it allows a detailed linear verification process of the numerical tools, which is at the basis of the development of gyrokinetic codes aimed at a rigorous turbulence investigation.Many numerical investigations of the linear GAM dynamics and comparison with analytical theory or benchmark among codes have been carried out in the last few decades, most of which treating the electrons as adiabatic. As a non-extensive list of example, we mention here simulations performed with the gyrokinetic codes GTC [20,21], ORB5 [22] (where the effect of the elongation was studied), TEMPEST [16] (where the effect of high-order terms of the finite ion orbit width was studied), GYRO [23] (where the effect of finite orbit width and the application to the radial velocity in the large-q limit was studied), GYSELA [24], ELMFIRE [25], and GENE with GKW [26]. In particular, a first verification of ORB5 against analytical theories, for circular geometries and low values of k r ρ i , was started in Ref. [27]. A numerical study of the effect of kinetic electrons in circular plasmas has been described in Ref. [28].In this paper, we aim at performing a comprehensive cross-code verification and benchmark of several gyrokinetic codes, in different regimes. We perform numerical simulations with three different gyrokinetic codes, adopting equivalent physical models for the dynamics of the ions (which is the most basic ...