Skip to main content

Lund University Publications

LUND UNIVERSITY LIBRARIES

CENTORI: A global toroidal electromagnetic two-fluid plasma turbulence code

Knight, P.J. ; Thyagaraja, A. ; Edwards, T.D. ; Hein, Joachim LU ; Romanelli, M. and McClements, K.G. (2012) In Computer Physics Communications 183(11). p.2346-2363
Abstract
A new global two-fluid electromagnetic turbulence code, CENTORI, has been developed for the purpose of studying magnetically-confined fusion plasmas on energy confinement timescales. This code is used to evolve the combined system of electron and ion fluid equations and Maxwell equations in toroidal configurations with axisymmetric equilibria. Uniquely, the equilibrium is co-evolved with the turbulence, and is thus modified by it. CENTORI is applicable to tokamaks of arbitrary aspect ratio and high plasma beta. A predictor–corrector, semi-implicit finite difference scheme is used to compute the time evolution of fluid quantities and fields. Vector operations and the evaluation of flux surface averages are speeded up by choosing the... (More)
A new global two-fluid electromagnetic turbulence code, CENTORI, has been developed for the purpose of studying magnetically-confined fusion plasmas on energy confinement timescales. This code is used to evolve the combined system of electron and ion fluid equations and Maxwell equations in toroidal configurations with axisymmetric equilibria. Uniquely, the equilibrium is co-evolved with the turbulence, and is thus modified by it. CENTORI is applicable to tokamaks of arbitrary aspect ratio and high plasma beta. A predictor–corrector, semi-implicit finite difference scheme is used to compute the time evolution of fluid quantities and fields. Vector operations and the evaluation of flux surface averages are speeded up by choosing the Jacobian of the transformation from laboratory to plasma coordinates to be a function of the equilibrium poloidal magnetic flux. A subroutine, GRASS, is used to co-evolve the plasma equilibrium by computing the steady-state solutions of a diffusion equation with a pseudo-time derivative. The code is written in Fortran 95 and is efficiently parallelised using Message Passing Interface (MPI). Illustrative examples of output from simulations of a tearing mode in a large aspect ratio tokamak plasma and of turbulence in an elongated conventional aspect ratio tokamak plasma are provided. (Less)
Please use this url to cite or link to this publication:
author
; ; ; ; and
organization
publishing date
type
Contribution to journal
publication status
published
subject
keywords
Two-fluid and multi-fluid plasmas, Drift waves, Tokamaks, Spherical tokamaks, Plasma turbulence, Magnetohydrodynamic and fluid equations
in
Computer Physics Communications
volume
183
issue
11
pages
2346 - 2363
publisher
Elsevier
external identifiers
  • wos:000308122000004
  • scopus:84864444421
ISSN
0010-4655
DOI
10.1016/j.cpc.2012.06.002
language
English
LU publication?
yes
id
bb54ecdd-8cd2-4586-9b96-eb338b0c3aff (old id 3195298)
date added to LUP
2016-04-01 10:52:27
date last changed
2022-02-17 22:07:31
@article{bb54ecdd-8cd2-4586-9b96-eb338b0c3aff,
  abstract     = {{A new global two-fluid electromagnetic turbulence code, CENTORI, has been developed for the purpose of studying magnetically-confined fusion plasmas on energy confinement timescales. This code is used to evolve the combined system of electron and ion fluid equations and Maxwell equations in toroidal configurations with axisymmetric equilibria. Uniquely, the equilibrium is co-evolved with the turbulence, and is thus modified by it. CENTORI is applicable to tokamaks of arbitrary aspect ratio and high plasma beta. A predictor–corrector, semi-implicit finite difference scheme is used to compute the time evolution of fluid quantities and fields. Vector operations and the evaluation of flux surface averages are speeded up by choosing the Jacobian of the transformation from laboratory to plasma coordinates to be a function of the equilibrium poloidal magnetic flux. A subroutine, GRASS, is used to co-evolve the plasma equilibrium by computing the steady-state solutions of a diffusion equation with a pseudo-time derivative. The code is written in Fortran 95 and is efficiently parallelised using Message Passing Interface (MPI). Illustrative examples of output from simulations of a tearing mode in a large aspect ratio tokamak plasma and of turbulence in an elongated conventional aspect ratio tokamak plasma are provided.}},
  author       = {{Knight, P.J. and Thyagaraja, A. and Edwards, T.D. and Hein, Joachim and Romanelli, M. and McClements, K.G.}},
  issn         = {{0010-4655}},
  keywords     = {{Two-fluid and multi-fluid plasmas; Drift waves; Tokamaks; Spherical tokamaks; Plasma turbulence; Magnetohydrodynamic and fluid equations}},
  language     = {{eng}},
  number       = {{11}},
  pages        = {{2346--2363}},
  publisher    = {{Elsevier}},
  series       = {{Computer Physics Communications}},
  title        = {{CENTORI: A global toroidal electromagnetic two-fluid plasma turbulence code}},
  url          = {{http://dx.doi.org/10.1016/j.cpc.2012.06.002}},
  doi          = {{10.1016/j.cpc.2012.06.002}},
  volume       = {{183}},
  year         = {{2012}},
}