Abstract | ||
---|---|---|
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. Introduces Boris-SDC, an arbitrary order extension of the Boris integrator.Boris-SDC is based on merging the Boris method with spectral deferred corrections.SDC for a 2nd-order problem using interpretation as preconditioned iteration.Analyzes properties of high-order Boris-SDC for particles in a Penning trap. |
Year | DOI | Venue |
---|---|---|
2014 | 10.1016/j.jcp.2015.04.022 | Journal of Computational Physics |
Keywords | Field | DocType |
collocation method,magnetic field | Mathematical optimization,Preconditioner,Mathematical analysis,Fixed-point iteration,Penning trap,Integrator,Rate of convergence,Equations of motion,Collocation method,Mathematics,Lorentz force | Journal |
Volume | Issue | ISSN |
abs/1409.5677 | C | Journal of Computational Physics 295, pp. 456-474, 2015 |
Citations | PageRank | References |
3 | 0.45 | 8 |
Authors | ||
3 |
Name | Order | Citations | PageRank |
---|---|---|---|
Mathias Winkel | 1 | 23 | 2.12 |
R. Speck | 2 | 11 | 2.03 |
Daniel Ruprecht | 3 | 71 | 10.02 |