Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aees_v1_0.tar.gz(206 Kbytes)
Manuscript Title: The nonlinear gyro-kinetic flux tube code GKW
Authors: A.G. Peeters, Y. Camenen, F.J. Casson, W.A. Hornsby, A.P. Snodin, D. Strintzi
Program title: GKW
Catalogue identifier: AEES_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 180(2009)2650
Programming language: Fortran 95.
Computer: Not computer specific.
Operating system: Any for which a Fortran 95 compiler is available.
Has the code been vectorised or parallelized?: Yes. The program can efficiently utilise 4096+ processors, depending on problem and available computer. 128 processors is reasonable for a typical nonlinear kinetic run on the latest x86-64 machines.
RAM: ~128MB-1GB for a linear run; 25GB for typical nonlinear kinetic run (30 million grid points).
Keywords: gyro-kinetic, flux tube, drift wave, tokamak, plasma turbulence.
PACS: 52.25.Fi, 52.25.Xz, 52.30.Gz, 52.35.Qz, 52.55.Fa, 52.65.Tt.
Classification: 19.8, 19.9, 19.11.

External routines: None required, although the functionality of the program is somewhat limited without a MPI implementation (preferably MPI-2) and the FFTW3 library.

Nature of problem:
Five dimensional gyro-kinetic Vlasov equation in general flux tube tokamak geometry with kinetic electrons, electro-magnetic effects and collisions.

Solution method:
Pseudo-spectral and finite difference with explicit time integration.

Additional comments:
The MHD equilibrium code CHEASE [1] is used for the general geometry calculations. This code has been developed in CRPP Lausanne and is not distributed together with GKW, but can be downloaded separately. The geometry module of GKW is based on the version 7.1 of CHEASE, which includes the output for Hamada coordinates.

Running time:
(On recent x86-64 hardware) ~10 minutes for a short linear problem; 48 hours for typical nonlinear kinetic run.

References:
[1] H. Lütjens, A. Bondeson, O. Sauter, Comput. Phys. Commun 97 (1996) 219 (http://cpc.cs.qub.ac.uk/summaries/ADDH_v1_0.html)