Explicit integrators for the magnetized equations of motion in Particle in Cell codes
Name
ExplicitMover.pdf
Size
669.57 KB
Format
Adobe PDF
Checksum (MD5)
50adf0913cf1f32f69ab7f95dac0bea8
Author(s)
Patacchini, Leonardo
Hutchinson, Ian Horner
Alternative Title
Explicit time-reversible orbit integration in Particle In Cell codes with static homogeneous magnetic field
Date Issued
December 2008
Journal
Journal of Computational Physics
Publisher
Academic Press
Citation
Patacchini, L., and I.H. Hutchinson. “Explicit time-reversible orbit integration in Particle In Cell codes with static homogeneous magnetic field.” Journal of Computational Physics 228.7 (2009): 2604-2615. © 2009 Elsevier Inc.
Version
Original manuscript
Abstract
A new explicit time-reversible orbit integrator for the equations of motion in a static homogeneous magnetic field – called Cyclotronic integrator – is presented. Like Spreiter and Walter’s Taylor expansion algorithm, for sufficiently weak electric field gradients this second order method does not require a fine resolution of the Larmor motion; it has however the essential advantage of being symplectic, hence time-reversible. The Cyclotronic integrator is only subject to a linear stability constraint ([OmegaDelta t] < pi, [Omega] being the Larmor angular frequency), and is therefore particularly suitable to electrostatic Particle In Cell codes with uniform magnetic field where [Omega]is larger than any other characteristic frequency, yet a resolution of the particles’ gyromotion is required. Application examples and a detailed comparison with the well-known (time-reversible) Boris algorithm are presented; it is in particular shown that implementation of the Cyclotronic integrator in the kinetic codes SCEPTIC and Democritus can reduce the cost of orbit integration by up to a factor of ten.
MIT Department
Massachusetts Institute of Technology. Department of Nuclear Science and Engineering
Massachusetts Institute of Technology. Plasma Science and Fusion Center
Terms of Use
Article is made available in accordance with the publisher's policy and may be subject to US copyright law. Please refer to the publisher's site for terms of use.
Persistent DSpace Link
DOI of Published Version
https://doi.org/10.1016/j.jcp.2008.12.021
