Hyper Boris integrators for kinetic plasma simulations
Abstract
We propose a family of numerical solvers for the nonrelativistic Newton--Lorentz equation in kinetic plasma simulations. The new solvers extend the standard 4-step Boris procedure, which has second-order accuracy in time, in three ways. First, we repeat the 4-step procedure multiple times, using an n-times smaller timestep ( t/n). We derive a formula for the arbitrary subcycling number n, so that we obtain the result without repeating the same calculations. Second, prior to the 4-step procedure, we apply Boris-type gyrophase corrections to the electromagnetic field. In addition to a well-known correction to the magnetic field, we correct the electric field in an anisotropic manner to achieve higher-order (N=2,4,6 …th order) accuracy. Third, combining these two methods, we propose a family of high-accuracy particle solvers, the hyper Boris solvers, which have two hyperparameters of the subcycling number n and the order of accuracy, N. The n-cycle Nth-order solver gives a numerical error of ( t/n)N at affordable computational cost.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.