An efficient energy-preserving algorithm for the Lorentz force system

摘要

In this paper, by utilizing the invariant energy quadratization method to transform the original Hamiltonian energy functional into a quadratic form, we propose a novel energy-preserving scheme to solve the Lorentz force system. The method is a linear-implicit scheme for canonical Hamiltonian system, and can efficiently simulate the motion of charged particles in constant magnetic field. With comparison to the well-used Boris method and a similar energy-preserving BDLI method [26], numerical experiments are presented to demonstrate the energy-preserving property and computational efficiency of the method.