High order Gauss-Seidel schemes for charged particle dynamics

Abstract

Gauss-Seidel projection methods are designed for achieving desirable long-term computational efficiency and reliability in micromagnetics simulations. While conventional Gauss-Seidel schemes are explicit, easy to use and furnish a better stability as compared to Euler's method, their order of accuracy is only one. This paper proposes an improved Gauss-Seidel methodology for particle simulations of magnetized plasmas. A novel new class of high order schemes are implemented via composition strategies. The new algorithms acquired are not only explicit and symmetric, but also volume-preserving together with their adjoint schemes. They are highly favorable for long-term computations. The new high order schemes are then utilized for simulating charged particle motions under the Lorentz force. Our experiments indicate a remarkable satisfaction of the energy preservation and angular momentum conservation of the numerical methods in multi-scale plasma dynamics computations.