A finite element Poisson solver for gyrokinetic particle simulations in a global field aligned mesh

Yasutaro Nishimura, Z. Lin, J. L.V. Lewandowski, S. Ethier

Research output: Contribution to journalArticle

24 Citations (Scopus)

Abstract

A new finite element Poisson solver is developed and applied to a global gyrokinetic toroidal code (GTC) which employs the field aligned mesh and thus a logically non-rectangular grid in a general geometry. Employing test cases where the analytical solutions are known, the finite element solver has been verified. The CPU time scaling versus the matrix size employing portable, extensible toolkit for scientific computation (PETSc) to solve the sparse matrix is promising. Taking the ion temperature gradient modes (ITG) as an example, the solution from the new finite element solver has been compared to the solution from the original GTC's iterative solver which is only efficient for adiabatic electrons. Linear and nonlinear simulation results from the two different forms of the gyrokinetic Poisson equation (integral form and the differential form) coincide each other. The new finite element solver enables the implementation of advanced kinetic electron models for global electromagnetic simulations.

Original languageEnglish
Pages (from-to)657-671
Number of pages15
JournalJournal of Computational Physics
Volume214
Issue number2
DOIs
Publication statusPublished - 2006 May 20

Fingerprint

mesh
Electrons
Poisson equation
ion temperature
Thermal gradients
Program processors
temperature gradients
electrons
simulation
grids
electromagnetism
scaling
Kinetics
Geometry
Ions
kinetics
matrices
geometry

All Science Journal Classification (ASJC) codes

  • Numerical Analysis
  • Modelling and Simulation
  • Physics and Astronomy (miscellaneous)
  • Physics and Astronomy(all)
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics

Cite this

@article{c7fb3e53ba88449c8286f0d1fd25e8e3,
title = "A finite element Poisson solver for gyrokinetic particle simulations in a global field aligned mesh",
abstract = "A new finite element Poisson solver is developed and applied to a global gyrokinetic toroidal code (GTC) which employs the field aligned mesh and thus a logically non-rectangular grid in a general geometry. Employing test cases where the analytical solutions are known, the finite element solver has been verified. The CPU time scaling versus the matrix size employing portable, extensible toolkit for scientific computation (PETSc) to solve the sparse matrix is promising. Taking the ion temperature gradient modes (ITG) as an example, the solution from the new finite element solver has been compared to the solution from the original GTC's iterative solver which is only efficient for adiabatic electrons. Linear and nonlinear simulation results from the two different forms of the gyrokinetic Poisson equation (integral form and the differential form) coincide each other. The new finite element solver enables the implementation of advanced kinetic electron models for global electromagnetic simulations.",
author = "Yasutaro Nishimura and Z. Lin and Lewandowski, {J. L.V.} and S. Ethier",
year = "2006",
month = "5",
day = "20",
doi = "10.1016/j.jcp.2005.10.011",
language = "English",
volume = "214",
pages = "657--671",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",
number = "2",

}

A finite element Poisson solver for gyrokinetic particle simulations in a global field aligned mesh. / Nishimura, Yasutaro; Lin, Z.; Lewandowski, J. L.V.; Ethier, S.

In: Journal of Computational Physics, Vol. 214, No. 2, 20.05.2006, p. 657-671.

Research output: Contribution to journalArticle

TY - JOUR

T1 - A finite element Poisson solver for gyrokinetic particle simulations in a global field aligned mesh

AU - Nishimura, Yasutaro

AU - Lin, Z.

AU - Lewandowski, J. L.V.

AU - Ethier, S.

PY - 2006/5/20

Y1 - 2006/5/20

N2 - A new finite element Poisson solver is developed and applied to a global gyrokinetic toroidal code (GTC) which employs the field aligned mesh and thus a logically non-rectangular grid in a general geometry. Employing test cases where the analytical solutions are known, the finite element solver has been verified. The CPU time scaling versus the matrix size employing portable, extensible toolkit for scientific computation (PETSc) to solve the sparse matrix is promising. Taking the ion temperature gradient modes (ITG) as an example, the solution from the new finite element solver has been compared to the solution from the original GTC's iterative solver which is only efficient for adiabatic electrons. Linear and nonlinear simulation results from the two different forms of the gyrokinetic Poisson equation (integral form and the differential form) coincide each other. The new finite element solver enables the implementation of advanced kinetic electron models for global electromagnetic simulations.

AB - A new finite element Poisson solver is developed and applied to a global gyrokinetic toroidal code (GTC) which employs the field aligned mesh and thus a logically non-rectangular grid in a general geometry. Employing test cases where the analytical solutions are known, the finite element solver has been verified. The CPU time scaling versus the matrix size employing portable, extensible toolkit for scientific computation (PETSc) to solve the sparse matrix is promising. Taking the ion temperature gradient modes (ITG) as an example, the solution from the new finite element solver has been compared to the solution from the original GTC's iterative solver which is only efficient for adiabatic electrons. Linear and nonlinear simulation results from the two different forms of the gyrokinetic Poisson equation (integral form and the differential form) coincide each other. The new finite element solver enables the implementation of advanced kinetic electron models for global electromagnetic simulations.

UR - http://www.scopus.com/inward/record.url?scp=33644957068&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=33644957068&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2005.10.011

DO - 10.1016/j.jcp.2005.10.011

M3 - Article

VL - 214

SP - 657

EP - 671

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 2

ER -