Explicit Runge–Kutta integrator with Hamiltonian correction for long-time simulations of guiding-center orbit in tokamak configurations

Abstract

Hamiltonian correction method is proposed to improve the variable time-step fourth-order Runge–Kutta methods in computing guiding-center orbits in a tokamak. It is found that the new method can significantly improve the computation efficiency of the conventional Runge–Kutta method in simulation of the long-time behavior of the guiding-center orbits.