The nonlinear gyro-kinetic flux tube code GKW

A.G. Peeters, Y. Camenen, F.J. Casson, W.A. Hornsby, A.P. Snodin, D. Strintzi and G. Szepesi

Computer Physics Communications, 180, 2650 (2009)

Copyright Elsevier (2009), doi:10.1016/j.cpc.2009.07.001
Download PDF or go to the journal link

Scaling graph

A new nonlinear gyro-kinetic flux tube code (GKW) for the simulation of micro instabilities and turbulence in magnetic confinement plasmas is presented in this paper. The code incorporates all physics effects that can be expected from a state of the art gyro-kinetic simulation code in the local limit: kinetic electrons, electromagnetic effects, collisions, full general geometry with a coupling to a MHD equilibrium code, and E×B shearing. In addition the physics of plasma rotation has been implemented through a formulation of the gyro-kinetic equation in the co-moving system. The gyro-kinetic model is five-dimensional and requires a massive parallel approach. GKW has been parallelised using MPI and scales well up to 8192+ cores. The paper presents the set of equations solved, the numerical methods, the code structure, and the essential benchmarks.


This paper is the main reference for the GKW code. It details the model equations used, and describes the implementation and benchmarks. Note that there is an even more detailed manual available in the documentation directory of the distributed code.

The figure shows the result of a strong scaling on a physically relevant problem.