Discretization form of the continuity condition at the polar axis, with application to the gyrokinetic simulation in a magnetic fusion torus
Abstract
A new computational method to solve the hyperbolic (Vlasov) equation coupled to 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.
- Publication:
-
arXiv e-prints
- Pub Date:
- February 2024
- DOI:
- 10.48550/arXiv.2402.03657
- arXiv:
- arXiv:2402.03657
- Bibcode:
- 2024arXiv240203657W
- Keywords:
-
- Physics - Computational Physics;
- Physics - Plasma Physics
- E-Print:
- 26 pages, 8 figures, 1 table