A filtered two-step variational integrator for charged-particle dynamics in a moderate or strong magnetic field

This article is concerned with a new filtered two-step variational integrator for solving the charged-particle dynamics in a mildly non-homogeneous moderate or strong magnetic field with a dimensionless parameter $\epsilon$ inversely proportional to the strength of the magnetic field. In the case of a moderate magnetic field ($\epsilon=1$), second-order error bounds and long time energy and momentum conservations are obtained. Moreover, the proof of the long-term analysis is accomplished by the backward error analysis. For the strong magnetic field ($0<\epsilon \ll1$), this paper clarifies the behaviour of the filtered variational integrator for both a large stepsize $h^2 \geq C\epsilon$ and a smaller stepsize $ h \leq c\epsilon$. The approach to analysing the error bounds for these two stepsizes is based on comparing the modulated Fourier expansions of the exact and the numerical solutions. It is shown that the proposed integrator achieves a second-order accuracy $\mathcal{O}(h^2)$ in the position and in the parallel velocity for a large step size and an $\mathcal{O}(\epsilon)$ accuracy for a smaller stepsize. This paper also yields the long time energy and magnetic moment conservations for the strong magnetic field by developing the modulated Fourier expansion of the proposed scheme. All the theoretical results of the error behaviour and long-term conservations are numerically demonstrated by four numerical experiments.