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 pullback scheme is implemented in the global gyrokinetic particle-incell code ORB5 [S. Jolliet et al, Comp. Phys. Comm., 177, 409 (2007)] to mitigate the cancellation problem in electromagnetic simulations. The equations and the discretisation used by the code are described. Numerical simulations of the Toroidal Alfvén Eigenmodes are performed in linear and nonlinear regimes to verify the scheme. A considerable improvement in the code efficiency is observed. For the internal kink mode, it is shown that the pullback mitigation efficiently cures a numerical instability which would make the simulation more costly otherwise.

This paper presents a study of the interaction between Alfvén modes and zonal structures, considering a realistic ASDEX Upgrade equilibrium. The results of gyrokinetic simulations with the global, electromagnetic, particle-in-cell code ORB5 are presented, where the modes are driven unstable by energetic particles with a bump-on-tail equilibrium distribution function, with radial density gradient. Two regimes have been observed. At low energetic particle concentration, the zonal structure (identified as an energetic particle-driven geodesic acoustic mode) is more unstable than the Alfvén mode. In the regime at high energetic particle concentration, the Alfvén mode is more unstable than the zonal structure. The interplay between the modes leads to a modification of their growth rates as well as to a modification of their saturation levels. The theoretical explanation of the mode interaction is given in terms of three-wave coupling of the energetic particle-driven geodesic acoustic mode and Alfvén mode, mediated by the curvature–pressure coupling term of the energetic particles.

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.