Symplectic Integration of Magnetic Systems
arXiv:1309.3723 · doi:10.1016/j.jcp.2014.03.049
Abstract
Dependable numerical results from long-time simulations require stable numerical integration schemes. For Hamiltonian systems, this is achieved with symplectic integrators, which conserve the symplectic condition and exactly solve for the dynamics for an approximate Hamiltonian. The conventional generating function formalism used for Hamiltonian systems is problematic for magnetic systems, as the Hamiltonian is not easily separable. This paper presents a derivation of symplectic integrators from a discretized action. Using this method, it is shown that the integrator due to Boris is a symplectic integrator in the non-relativistic limit, and a relativistic integrator is derived.
13 pages, 1 figure
Cited by in corpus (13)
- Novel aspects of cosmic ray diffusion in synthetic magnetic turbulence
- Explicit high-order symplectic integrators for charged particles in general electromagnetic fields
- A high-order Boris integrator
- Comment on "Symplectic integration of magnetic systems": a proof that the Boris algorithm is not variational
- Energetic particle dynamics in a simplified model of a solar wind magnetic switchback
- Invariant compact finite difference schemes
- Explicit volume-preserving numerical schemes for relativistic trajectories and spin dynamics
- Large-stepsize integrators for charged-particle dynamics over multiple time scales
- A two-step symmetric method for charged-particle dynamics in a normal or strong magnetic feld
- Error estimates of some splitting schemes for charged-particle dynamics under strong magnetic field
- Computational Frameworks for Patterned Two-Dimensional Magnetism
- A sixth order averaged vector field method
- Exponential energy-preserving methods for charged-particle dynamics in a strong and constant magnetic field