Fast and spectrally accurate evaluation of gyroaverages in non-periodic gyrokinetic-poisson simulations

J. Guadagni, Antoine Cerfon

Research output: Contribution to journalArticle

Abstract

We present a fast and spectrally accurate numerical scheme for the evaluation of the gyroaveraged electrostatic potential in non-periodic gyrokinetic-Poisson simulations. Our method relies on a reformulation of the gyrokinetic-Poisson system in which the gyroaverage in Poisson’s equation is computed for the compactly supported charge density instead of the non-periodic, non-compactly supported potential itself. We calculate this gyroaverage with a combination of two Fourier transforms and a Hankel transform, which has the near optimal run-time complexity O(Nρ(P+ P)log(P+ P)), where P is the number of spatial grid points, P the number of grid points in Fourier space and Nρ the number of grid points in velocity space. We present numerical examples illustrating the performance of our code and demonstrating geometric convergence of the error.

Original languageEnglish (US)
Article number905830407
JournalJournal of Plasma Physics
Volume83
Issue number4
DOIs
StatePublished - Aug 1 2017

Fingerprint

grids
evaluation
simulation
Poisson equation
electrostatics

Keywords

  • Intense particle beams
  • Magnetized plasmas
  • Plasma simulation

ASJC Scopus subject areas

  • Condensed Matter Physics

Cite this

Fast and spectrally accurate evaluation of gyroaverages in non-periodic gyrokinetic-poisson simulations. / Guadagni, J.; Cerfon, Antoine.

In: Journal of Plasma Physics, Vol. 83, No. 4, 905830407, 01.08.2017.

Research output: Contribution to journalArticle

@article{08fa2b1ca86c4bfabc1c54262a4bbb5e,
title = "Fast and spectrally accurate evaluation of gyroaverages in non-periodic gyrokinetic-poisson simulations",
abstract = "We present a fast and spectrally accurate numerical scheme for the evaluation of the gyroaveraged electrostatic potential in non-periodic gyrokinetic-Poisson simulations. Our method relies on a reformulation of the gyrokinetic-Poisson system in which the gyroaverage in Poisson’s equation is computed for the compactly supported charge density instead of the non-periodic, non-compactly supported potential itself. We calculate this gyroaverage with a combination of two Fourier transforms and a Hankel transform, which has the near optimal run-time complexity O(Nρ(P+ P)log(P+ P)), where P is the number of spatial grid points, P the number of grid points in Fourier space and Nρ the number of grid points in velocity space. We present numerical examples illustrating the performance of our code and demonstrating geometric convergence of the error.",
keywords = "Intense particle beams, Magnetized plasmas, Plasma simulation",
author = "J. Guadagni and Antoine Cerfon",
year = "2017",
month = "8",
day = "1",
doi = "10.1017/S0022377817000599",
language = "English (US)",
volume = "83",
journal = "Journal of Plasma Physics",
issn = "0022-3778",
publisher = "Cambridge University Press",
number = "4",

}

TY - JOUR

T1 - Fast and spectrally accurate evaluation of gyroaverages in non-periodic gyrokinetic-poisson simulations

AU - Guadagni, J.

AU - Cerfon, Antoine

PY - 2017/8/1

Y1 - 2017/8/1

N2 - We present a fast and spectrally accurate numerical scheme for the evaluation of the gyroaveraged electrostatic potential in non-periodic gyrokinetic-Poisson simulations. Our method relies on a reformulation of the gyrokinetic-Poisson system in which the gyroaverage in Poisson’s equation is computed for the compactly supported charge density instead of the non-periodic, non-compactly supported potential itself. We calculate this gyroaverage with a combination of two Fourier transforms and a Hankel transform, which has the near optimal run-time complexity O(Nρ(P+ P)log(P+ P)), where P is the number of spatial grid points, P the number of grid points in Fourier space and Nρ the number of grid points in velocity space. We present numerical examples illustrating the performance of our code and demonstrating geometric convergence of the error.

AB - We present a fast and spectrally accurate numerical scheme for the evaluation of the gyroaveraged electrostatic potential in non-periodic gyrokinetic-Poisson simulations. Our method relies on a reformulation of the gyrokinetic-Poisson system in which the gyroaverage in Poisson’s equation is computed for the compactly supported charge density instead of the non-periodic, non-compactly supported potential itself. We calculate this gyroaverage with a combination of two Fourier transforms and a Hankel transform, which has the near optimal run-time complexity O(Nρ(P+ P)log(P+ P)), where P is the number of spatial grid points, P the number of grid points in Fourier space and Nρ the number of grid points in velocity space. We present numerical examples illustrating the performance of our code and demonstrating geometric convergence of the error.

KW - Intense particle beams

KW - Magnetized plasmas

KW - Plasma simulation

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

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

U2 - 10.1017/S0022377817000599

DO - 10.1017/S0022377817000599

M3 - Article

AN - SCOPUS:85070536648

VL - 83

JO - Journal of Plasma Physics

JF - Journal of Plasma Physics

SN - 0022-3778

IS - 4

M1 - 905830407

ER -