The nonlinear gyro-kinetic flux tube code GKW
Peeters, A. G., Camenen, Y., Casson, F. J. (Francis James), Hornsby, W. A., Snodin, A. P., Strintzi, D. and Szepesi, Gabor. (2009) The nonlinear gyro-kinetic flux tube code GKW. Computer Physics Communications, Vol.180 (No.12). pp. 2650-2672. ISSN 0010-4655Full text not available from this repository.
Official URL: http://dx.doi.org/10.1016/j.cpc.2009.07.001
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 x 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 fivedimensional 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.
Program title: GKW Catalogue identifier: AEES_vI_0 Program summary URL: http://cpc.cs.qLib.ac.uk/summaries/AEES-vl-O.html
Program obtainablefrom: CK Program Library, Queen's University, Belfast, N. Ireland
Licensing provisions: GNU GPL v3
No. of lines in distributed program, including test data, etc.: 29998
No. of bytes in distributed program, including test data, etc.: 206943 Distribution format. tar.gz
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 parallelised?: Yes. The program can efficiently utilise 8192+ processors, depending on problem and available computer. 128 processors is reasonable for a typical nonlinear kinetic run on the latest x86-64 machines. RAM: similar to 128 MB-1 GB for a linear run; 25 GB for typical nonlinear kinetic run (30 million grid points)
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  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. Runningtime: (On recent x86-64 hardware) -10 minutes for a short linear problem; 48 hours for typical nonlinear kinetic run. Reference:  H. Lutjens, A. Bondeson, O. Sauter, Comput. Phys. Comm. 97 (1996) 219, http://cpc.cs.qub.ac.uk/ surnrnaries/ADDH_v1_0.html Crown Copyright (C) 2009 Published by Elsevier B.V. All rights reserved.
|Item Type:||Journal Article|
|Subjects:||Q Science > QA Mathematics > QA76 Electronic computers. Computer science. Computer software
Q Science > QC Physics
|Divisions:||Faculty of Science > Physics|
|Journal or Publication Title:||Computer Physics Communications|
|Publisher:||Elsevier Science BV|
|Official Date:||December 2009|
|Number of Pages:||23|
|Page Range:||pp. 2650-2672|
|Access rights to Published version:||Restricted or Subscription Access|
Actions (login required)
Downloads per month over past year