2016/09/16 by Stephen D. Webb, Webb, Stephen D., Dan T. Abell +7
Mathematics · Physics and Astronomy · #Magnetic confinement fusion research #Numerical methods for differential equations #Physics of Superconductivity and Magnetism #physics.acc-ph #physics.comp-ph #physics.plasm-ph
paper · pdf · doi:10.48550/arxiv.1609.05095
22 pages -- Text contains an error in the basis functions, and requires some work to fix. An updated version will be posted when this is done
arxiv created 2017/05/09 · arxiv updated 2017/05/10
Symplectic integrators for Hamiltonian systems have been quite successful for studying few-body dynamical systems. These integrators are frequently derived using a formalism built on symplectic maps. There have been recent efforts to extend the symplectic approach to plasmas, which have focused primarily on discrete Lagrangian mechanics. In this paper, we derive a a symplectic electromagnetic macroparticle algorithm using the map formalism. The resulting algorithm is designed to prevent numerical instabilities such as numerical Čerenkov, which result from incorrect dispersion relations for the fields, as well as the artificial heating of plasmas, which arise from the non-symplectic nature of conventional particle-in-cell algorithms. This is the first self-consistent electromagnetic algorithm derived using a map-based approach.