Numerical Stability Improvements for the Pseudo-Spectral EM PIC Algorithm
Abstract
The pseudo-spectral analytical time-domain (PSATD) particle-in-cell (PIC) algorithm solves the vacuum Maxwell's equations exactly, has no Courant time-step limit (as conventionally defined), and offers substantial flexibility in plasma and particle beam simulations. It is, however, not free of the usual numerical instabilities, including the numerical Cherenkov instability, when applied to relativistic beam simulations. This paper presents several approaches that, when combined with digital filtering, almost completely eliminate the numerical Cherenkov instability. The paper also investigates the numerical stability of the PSATD algorithm at low beam energies.
Keywords
Cite
@article{arxiv.1309.0116,
title = {Numerical Stability Improvements for the Pseudo-Spectral EM PIC Algorithm},
author = {Brendan B. Godfrey and Jean-Luc Vay and Irving Haber},
journal= {arXiv preprint arXiv:1309.0116},
year = {2014}
}
Comments
Based on presentation at the 2013 IEEE Pulsed Power and Plasma Science Conference; 8 pages, 8 figures