An energy-conserving and asymptotic-preserving charged-particle orbit implicit time integrator for arbitrary electromagnetic fields

Abstract

We present a new implicit asymptotic preserving time integration scheme for charged-particle orbit computation in arbitrary electromagnetic fields. The scheme is built on the Crank-Nicolson integrator and continues to recover full-orbit motion in the small time-step limit, but also recovers all the first-order guiding center drifts as well as the correct gyroradius when stepping over the gyration time-scale. In contrast to previous efforts in this direction, the new scheme also features exact energy conservation. In the derivation of the scheme, we find that a new numerical time-scale is introduced. This scale is analyzed and the resulting restrictions on time-step are derived. Based on this analysis, we develop an adaptive time-stepping strategy the respects these constraints while stepping over the gyration scale when physically justified. It is shown through numerical tests on single-particle motion that the scheme's energy conservation property results in tremendous improvements in accuracy, and that the scheme is able to transition smoothly between magnetized and unmagnetized regimes as a result of the adaptive time-stepping.

References in corpus (6)

• An Energy- and Charge-conserving, Implicit, Electrostatic Particle-in-Cell Algorithm
• A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov-Darwin particle-in-cell algorithm
• An energy- and charge-conserving, nonlinearly implicit, electromagnetic 1D-3V Vlasov-Darwin particle-in-cell algorithm
• An efficient mixed-precision, hybrid CPU-GPU implementation of a fully implicit particle-in-cell algorithm
• Fluid preconditioning for Newton-Krylov-based, fully implicit, electrostatic particle-in-cell simulations
• A fast, high-order solver for the Grad-Shafranov equation

Cited by in corpus (7)

• Enabling particle applications for exascale computing platforms
• An implicit, conservative and asymptotic-preserving electrostatic particle-in-cell algorithm for arbitrarily magnetized plasmas in uniform magnetic fields
• Quasi-Helmholtz Decomposition, Gauss' Laws and Charge Conservation for Finite Element Particle-in-Cell
• The anatomy of Boris type solvers and the Lie operator formalism for deriving large time-step magnetic field integrators
• Exploring exponential time integration for strongly magnetized charged particle motion
• Large-stepsize integrators for charged-particle dynamics over multiple time scales
• Nyström Type Exponential Integrators for Strongly Magnetized Charged Particle Dynamics