Navigation and service


GYrokinetic SEmi-LAgrangian code for plasma tubulence simulations

Scientific area

Simulations of the electrostatic branch of the ion temperature gradient turbulence in tokamak plasmas.

Short description

The GYSELA code is a non-linear 5D global gyrokinetic full-f code which performs flux-driven simulations of ion temperature gradient driven turbulence (ITG) in the electrostatic limit with adiabatic electrons. It solves the standard gyrokinetic equation for the full-f distribution function, i.e. no assumption on scale separation between equilibrium and perturbations is done. This 5D equation is self-consistently coupled to a 3D quasineutrality equation. The code also includes ion-ion collisions in global simplified magnetic geometry (concentric toroidal magnetic surfaces with circular cross-sections). The system is driven out of thermodynamical equilibrium by versatile sources which are capable of injecting separately parallel momentum and heat. Alternate sources have also been studied, in order to excite supra-thermal particles. Finally, the code has the originality to be based on a semi-Lagrangian scheme and it is parallelized into an hybrid OpenMP/MPI paradigm. Recent simulations have investigated: (i) the impact on Energetic-Particle-Driven Geodesic Acoustic Modes on Turbulence, (ii) ion transport barriers triggered by plasma polarization and (iii) turbulent momentum transport in core tokamak plasmas.


  • 458,752 cores (1,835,008 parallel threads) on BlueGene/Q (JUQUEEN)
  • 65,536 cores on x86 (Curie and Helios)

Relative efficiency (weak scaling) of GYSELA on JUQUEENRelative efficiency (weak scaling) of GYSELA on JUQUEEN

Programming language and model

  • Fortran 90, C
  • MPI, OpenMP, pthreads

Tested on platforms

  • BlueGene/Q
  • x86

Application developers

F-13108 Saint-Paul-Lez-Durance cedex


Virginie Grandgirard
IRFM-DSM-CEA Cadarache