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. doi:10.1016/j.cpc.2009.07.001

Research output not available from this repository, contact author.

Request Changes to record.

Abstract

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 summary

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 [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. Runningtime: (On recent x86-64 hardware) -10 minutes for a short linear problem; 48 hours for typical nonlinear kinetic run. Reference: [1] 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
ISSN:	0010-4655
Official Date:	December 2009
Dates:	Date Event December 2009 Published
Volume:	Vol.180
Number:	No.12
Number of Pages:	23
Page Range:	pp. 2650-2672
DOI:	10.1016/j.cpc.2009.07.001
Status:	Peer Reviewed
Publication Status:	Published
Access rights to Published version:	Restricted or Subscription Access

Data sourced from Thomson Reuters' Web of Knowledge

Request changes or add full text files to a record

Repository staff actions (login required)

View Item