A high-order Boris integrator

Abstract

This work introduces the high-order Boris-SDC method for integrating the equations of motion for electrically charged particles in an electric and magnetic field. Boris-SDC relies on a combination of the Boris-integrator with spectral deferred corrections (SDC). SDC can be considered as preconditioned Picard iteration to compute the stages of a collocation method. In this interpretation, inverting the preconditioner corresponds to a sweep with a low-order method. In Boris-SDC, the Boris method, a second-order Lorentz force integrator based on velocity-Verlet, is used as a sweeper/preconditioner. The presented method provides a generic way to extend the classical Boris integrator, which is widely used in essentially all particle-based plasma physics simulations involving magnetic fields, to a high-order method. Stability, convergence order and conservation properties of the method are demonstrated for different simulation setups. Boris-SDC reproduces the expected high order of convergence for a single particle and for the center-of-mass of a particle cloud in a Penning trap and shows good long-term energy stability.

References in corpus (5)

• A multi-level spectral deferred correction method
• Interweaving PFASST and Parallel Multigrid
• Symplectic Integration of Magnetic Systems
• Inexact spectral deferred corrections
• Spectral deferred correction methods for second-order problems

Cited by in corpus (13)

• PIC methods in astrophysics: Simulations of relativistic jets and kinetic physics in astrophysical systems
• Multi-Level Spectral Deferred Corrections Scheme for the Shallow Water Equations on the Rotating Sphere
• Advanced fuel fusion, phase space engineering, and structure-preserving geometric algorithms
• Convergence analysis of multi-level spectral deferred corrections
• An Analytic Boris Pusher for Plasma Simulation
• An arbitrary order time-stepping algorithm for tracking particles in inhomogeneous magnetic fields
• Performance of the BGSDC integrator for computing fast ion trajectories in nuclear fusion reactors
• Spectral deferred correction methods for second-order problems
• New applications for the Boris Spectral Deferred Correction algorithm for plasma simulations
• A multigrid perspective on the parallel full approximation scheme in space and time
• Parallelizing spectral deferred corrections across the method
• Hyper Boris integrators for kinetic plasma simulations
• A Comparative Analysis of Relativistic Particle Pushers vis-à-vis Computation Time & Accuracy