A high-order Boris integrator

JOURNAL OF COMPUTATIONAL PHYSICS

卷:295

Article

Winkel, Mathias¹  Speck, Robert^1,2  Ruprecht, Daniel¹ 

[1] Univ Lugano, Inst Computat Sci, Lugano, Switzerland

[2] Forschungszentrum Julich, Julich Supercomp Ctr, D-52425 Julich, Germany

关键词: Boris integrator;    Time integration;    Magnetic field;    High-order;    Spectral deferred corrections (SDC);    Collocation method;   

DOI  :  10.1016/j.jcp.2015.04.022

来源: Elsevier

【 摘 要 】

This work introduces the high-order Boris-SDC method for integrating the equations of motion for electrically charged particles in electric and magnetic fields. 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. (C) 2015 Elsevier Inc. All rights reserved.

【 授权许可】

Free