New applications for the Boris Spectral Deferred Correction algorithm for plasma simulations

Abstract

The paper investigates two new use cases for the Boris Spectral Deferred Corrections (Boris-SDC) time integrator for plasma simulations. First, we show that using Boris-SDC as a particle pusher in an electrostatic particle-in-cell (PIC) code can, at least in the linear regime, improve simulation accuracy compared with the standard second order Boris method. In some instances, the higher order of Boris-SDC even allows a much larger time step, leading to modest computational gains. Second, we propose a modification of Boris-SDC for the relativistic regime. Based on an implementation of Boris-SDC in the \textsc{runko} PIC code, we demonstrate for a relativistic Penning trap that Boris-SDC retains its high order of convergence for velocities ranging from $0.5c$ to $>0.99c$. We also show that for the force-free case where acceleration from electric and magnetic field cancel, Boris-SDC produces less numerical drift than Boris.

References in corpus (6)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Structure-preserving second-order integration of relativistic charged particle trajectories in electromagnetic fields
• Coherent Electromagnetic Emission from Relativistic Magnetized Shocks
• Explicit high-order symplectic integrators for charged particles in general electromagnetic fields
• 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