Log In Sign Up

The anatomy of Boris type solvers and large time-step plasma simulations

by   Siu A. Chin, et al.

This work gives a Lie operator derivation of various Boris solvers via a detailed study of first and second-order trajectory errors in a constant magnetic field. These errors in the gyro-circle center and gyro-radius are the foundational basis for why Boris solvers existed, independent of any finite-difference schemes. The elimination of these errors then forces the second-order solver's trajectory to be exactly on the gyro-circle. By revisiting some historical calculations, it is found that many publications do not distinguish the poorly behaved first-order leap-frog solver with the correct second-order Boris algorithm. This work shows that this second-order Boris solver is much more accurate then previously thought and that its trajectory remains close to the exact orbit in a combined nonuniform electric and magnetic field at time-steps greater than the cyclotron period.


page 1

page 2

page 3

page 4


Structure preserving schemes for Fokker-Planck equations of irreversible processes

In this paper, we introduce second order and fourth order space discreti...

A Fully Fourth Order Accurate Energy Stable FDTD Method for Maxwell's Equations in Metamaterials

We present a novel fully fourth order in time and space finite differenc...

An arbitrary order time-stepping algorithm for tracking particles in inhomogeneous magnetic fields

The Lorentz equations describe the motion of electrically charged partic...

Second-order finite difference approximations of the upper-convected time derivative

In this work, new finite difference schemes are presented for dealing wi...

An efficient jet marcher for computing the quasipotential for 2D SDEs

We present a new algorithm, the efficient jet marching method (EJM), for...