A linear momentum conserving interpolation for the electric field in a three-dimensional toroidal particle-in-cell gyrokinetic plasma simulation in a tokamak configuration is found unstable due to the false divergence of Ε X Β flow of particles. A paradigm 8-point finite difference for interpolation of the radial electric field on the poloidal plane is proposed which stabilizes the simulation.
- Fusion plasma physics
- gyrokinetic particle simulation
- momentum conserving interpolation
Heikkinen, J. A., Korpilo, T., Janhunen, S. J., Kiviniemi, T. P., Leerink, S., & Ogando, F. (2012). Interpolation for momentum conservation in 3D toroidal gyrokinetic particle simulation of plasmas. Computer Physics Communications, 183(8), 1719-1727. https://doi.org/10.1016/j.cpc.2012.03.013