Abstract
Here we present details of an operator-split, implicit–explicit numerical scheme for the solution of the gyrokinetic-Poisson system of equations in the local limit. This scheme has been implemented in a new code called stella, which is capable of evolving electrostatic fluctuations with full kinetic electron effects and an arbitrary number of ion species in general magnetic geometry. We demonstrate the advantages of this mixed approach over a fully explicit treatment and provide linear and nonlinear benchmark comparisons for both axisymmetric and non-axisymmetric magnetic equilibria.
Original language | English (US) |
---|---|
Pages (from-to) | 365-380 |
Number of pages | 16 |
Journal | Journal of Computational Physics |
Volume | 391 |
DOIs | |
State | Published - Aug 15 2019 |
Externally published | Yes |
All Science Journal Classification (ASJC) codes
- Numerical Analysis
- Modeling and Simulation
- Physics and Astronomy (miscellaneous)
- General Physics and Astronomy
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics
Keywords
- Gyrokinetics
- Magnetic confinement fusion
- Stellarator
- Transport
- Turbulence