EUTERPE gyrokinetic code

EUTERPE is a gyrokinetic particle in cell and delta-f code specifically designed for three-dimensional geometríes [Jost PoP 2001; Kornilov NF 2004]. It was created as a global linear code at CRPP Lausanne, and afterward, it was further developed at Max Planck IPP Greifswald, where it is presently developed and maintained. Many new features have been added since its creation in 2001.

Presently, the code implements a variety of models of increasing complexity and can perform either electrostatic or electromagnetic simulations. It can manage up to three fully kinetic species: ions, electrons, and a third (fast or impurity) species with arbitrary charge/mass. It implements a pitch angle scattering collision operator [Kauffman et al. JoPCS 2010], and more recently, an energy and momentum conserving operator was implemented [Regaña et al PPCF 2013]. Both linear and nonlinear simulations can now be carried out, and several tools for numerical noise control are implemented in it [Sánchez et al. JPP 2020].

The code EUTERPE has been benchmarked in linear settings against several codes, such as TORB [Sánchez et al. IEEE TPP 2010], GENE [Helander et al. NF 2015; Görler et al. PoP 2016, Monreal et al. PPCF 2016, Monreal et al. PPCF 2017, Sanchez et al. NF 2021], GYSELA, GKW and ORB5 (Görler et al. PoP 2016, Sánchez et al. JPP 2020), XGC [Cole et al. PoP 2019], GTC [Wang et al. POP 2020], and stella [Sanchez et al. NF 2021]. In a turbulence nonlinear setting, it has been compared with TORB for electrostatic simulations in a screw pinch geometry [Sánchez et al. IEEE TPP 2010], with ORB5 [Sánchez et al JPP 2020], and with GENE-3D in [Sánchez et al. NF 2023].

For more details about the code, the reader is referred to [Jost et al. PoP 2001, Kornilov et al. PoP 2004, Slaby et al. NF 2018, Kleiber2023 (preparation)].