Finite difference form of the continuity condition at the polar axis, with application to the gyrokinetic simulation in a magnetic fusion torus
Abstract
Abstract A new computational method to solve the hyperbolic (gyrokinetic Vlasov) equation and the elliptic (Poisson-like) equation at the polar axis is proposed. It is shown that the value of a scalar function at the polar axis can be predicted by its neighbouring values based on the continuity condition. This continuity condition systematically solves the pole problems including the singular factor 1/ r in the hyperbolic equation and the inner boundary in the elliptic equation. The proposed method is applied to the global gyrokinetic simulation of the tokamak plasma with the magnetic axis included.
Fields of Study & Subjects
Magnetic confinement fusion research
Differential Equations and Boundary Problems
Nonlinear Waves and Solitons
Publication details
- Journal / source
- Plasma Science and Technology
- Publisher
- IOP Publishing
- Volume
- 28
- Issue
- 3
- Pages
- 035101–035101
- Publication date
- 9 January 2026