Full f gyrokinetic method for particle simulation of tokamak transport

Jukka A. Heikkinen (Corresponding Author), S. J. Janhunen, T. P. Kiviniemi, F. Ogando

Research output: Contribution to journalArticleScientificpeer-review

67 Citations (Scopus)

Abstract

A gyrokinetic particle-in-cell approach with direct implicit construction of the coefficient matrix of the Poisson equation from ion polarization and electron parallel nonlinearity is described and applied in global electrostatic toroidal plasma transport simulations. The method is applicable for calculation of the evolution of particle distribution function f including as special cases strong plasma pressure profile evolution by transport and formation of neoclassical flows. This is made feasible by full f formulation and by recording the charge density changes due to the ion polarization drift and electron acceleration along the local magnetic field while particles are advanced. The code has been validated against the linear predictions of the unstable ion temperature gradient mode growth rates and frequencies. Convergence and saturation in both turbulent and neoclassical limit of the ion heat conductivity is obtained with numerical noise well suppressed by a sufficiently large number of simulation particles. A first global full f validation of the neoclassical radial electric field in the presence of turbulence for a heated collisional tokamak plasma is obtained. At high Mach number (Mp∼1) of the poloidal flow, the radial electric field is significantly enhanced over the standard neoclassical prediction. The neoclassical radial electric field together with the related GAM oscillations is found to regulate the turbulent heat and particle diffusion levels particularly strongly in a large aspect ratio tokamak at low plasma current.
Original languageEnglish
Pages (from-to)5582-5609
JournalJournal of Computational Physics
Volume227
Issue number11
DOIs
Publication statusPublished - 2008
MoE publication typeA1 Journal article-refereed

Fingerprint

Plasmas
Electric fields
Ions
electric fields
simulation
Polarization
linear prediction
heat
particle diffusion
ions
toroidal plasmas
plasma pressure
electron acceleration
Electrons
plasma currents
Poisson equation
polarization
ion temperature
Charge density
Mach number

Keywords

  • particle simulation
  • plasma
  • turbulence
  • Tokamak
  • fusion energy

Cite this

Heikkinen, J. A., Janhunen, S. J., Kiviniemi, T. P., & Ogando, F. (2008). Full f gyrokinetic method for particle simulation of tokamak transport. Journal of Computational Physics, 227(11), 5582-5609. https://doi.org/10.1016/j.jcp.2008.02.013
Heikkinen, Jukka A. ; Janhunen, S. J. ; Kiviniemi, T. P. ; Ogando, F. / Full f gyrokinetic method for particle simulation of tokamak transport. In: Journal of Computational Physics. 2008 ; Vol. 227, No. 11. pp. 5582-5609.
@article{c8385e86b6494a51b5c247598c9a3221,
title = "Full f gyrokinetic method for particle simulation of tokamak transport",
abstract = "A gyrokinetic particle-in-cell approach with direct implicit construction of the coefficient matrix of the Poisson equation from ion polarization and electron parallel nonlinearity is described and applied in global electrostatic toroidal plasma transport simulations. The method is applicable for calculation of the evolution of particle distribution function f including as special cases strong plasma pressure profile evolution by transport and formation of neoclassical flows. This is made feasible by full f formulation and by recording the charge density changes due to the ion polarization drift and electron acceleration along the local magnetic field while particles are advanced. The code has been validated against the linear predictions of the unstable ion temperature gradient mode growth rates and frequencies. Convergence and saturation in both turbulent and neoclassical limit of the ion heat conductivity is obtained with numerical noise well suppressed by a sufficiently large number of simulation particles. A first global full f validation of the neoclassical radial electric field in the presence of turbulence for a heated collisional tokamak plasma is obtained. At high Mach number (Mp∼1) of the poloidal flow, the radial electric field is significantly enhanced over the standard neoclassical prediction. The neoclassical radial electric field together with the related GAM oscillations is found to regulate the turbulent heat and particle diffusion levels particularly strongly in a large aspect ratio tokamak at low plasma current.",
keywords = "particle simulation, plasma, turbulence, Tokamak, fusion energy",
author = "Heikkinen, {Jukka A.} and Janhunen, {S. J.} and Kiviniemi, {T. P.} and F. Ogando",
year = "2008",
doi = "10.1016/j.jcp.2008.02.013",
language = "English",
volume = "227",
pages = "5582--5609",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press",
number = "11",

}

Heikkinen, JA, Janhunen, SJ, Kiviniemi, TP & Ogando, F 2008, 'Full f gyrokinetic method for particle simulation of tokamak transport', Journal of Computational Physics, vol. 227, no. 11, pp. 5582-5609. https://doi.org/10.1016/j.jcp.2008.02.013

Full f gyrokinetic method for particle simulation of tokamak transport. / Heikkinen, Jukka A. (Corresponding Author); Janhunen, S. J.; Kiviniemi, T. P.; Ogando, F.

In: Journal of Computational Physics, Vol. 227, No. 11, 2008, p. 5582-5609.

Research output: Contribution to journalArticleScientificpeer-review

TY - JOUR

T1 - Full f gyrokinetic method for particle simulation of tokamak transport

AU - Heikkinen, Jukka A.

AU - Janhunen, S. J.

AU - Kiviniemi, T. P.

AU - Ogando, F.

PY - 2008

Y1 - 2008

N2 - A gyrokinetic particle-in-cell approach with direct implicit construction of the coefficient matrix of the Poisson equation from ion polarization and electron parallel nonlinearity is described and applied in global electrostatic toroidal plasma transport simulations. The method is applicable for calculation of the evolution of particle distribution function f including as special cases strong plasma pressure profile evolution by transport and formation of neoclassical flows. This is made feasible by full f formulation and by recording the charge density changes due to the ion polarization drift and electron acceleration along the local magnetic field while particles are advanced. The code has been validated against the linear predictions of the unstable ion temperature gradient mode growth rates and frequencies. Convergence and saturation in both turbulent and neoclassical limit of the ion heat conductivity is obtained with numerical noise well suppressed by a sufficiently large number of simulation particles. A first global full f validation of the neoclassical radial electric field in the presence of turbulence for a heated collisional tokamak plasma is obtained. At high Mach number (Mp∼1) of the poloidal flow, the radial electric field is significantly enhanced over the standard neoclassical prediction. The neoclassical radial electric field together with the related GAM oscillations is found to regulate the turbulent heat and particle diffusion levels particularly strongly in a large aspect ratio tokamak at low plasma current.

AB - A gyrokinetic particle-in-cell approach with direct implicit construction of the coefficient matrix of the Poisson equation from ion polarization and electron parallel nonlinearity is described and applied in global electrostatic toroidal plasma transport simulations. The method is applicable for calculation of the evolution of particle distribution function f including as special cases strong plasma pressure profile evolution by transport and formation of neoclassical flows. This is made feasible by full f formulation and by recording the charge density changes due to the ion polarization drift and electron acceleration along the local magnetic field while particles are advanced. The code has been validated against the linear predictions of the unstable ion temperature gradient mode growth rates and frequencies. Convergence and saturation in both turbulent and neoclassical limit of the ion heat conductivity is obtained with numerical noise well suppressed by a sufficiently large number of simulation particles. A first global full f validation of the neoclassical radial electric field in the presence of turbulence for a heated collisional tokamak plasma is obtained. At high Mach number (Mp∼1) of the poloidal flow, the radial electric field is significantly enhanced over the standard neoclassical prediction. The neoclassical radial electric field together with the related GAM oscillations is found to regulate the turbulent heat and particle diffusion levels particularly strongly in a large aspect ratio tokamak at low plasma current.

KW - particle simulation

KW - plasma

KW - turbulence

KW - Tokamak

KW - fusion energy

U2 - 10.1016/j.jcp.2008.02.013

DO - 10.1016/j.jcp.2008.02.013

M3 - Article

VL - 227

SP - 5582

EP - 5609

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 11

ER -