Stability of energy-conserving PIC algorithms with smoother interpolations and more accurate Poisson solves

ORAL

Abstract

We have performed periodic, one-dimensional particle-in-cell (PIC) simulations to test explicit energy-conserving PIC algorithms with a variety of properties: continuous interpolated electric fields, continuous first derivatives of the interpolated electric field, Lagrangian-derived field solves, and second- and forth-order accurate field solves. We present numerical measurements of the stability, accuracy, and noise for each method, as well as comment on the generalizability to complex boundary conditions. In all cases, we find that the algorithm suffers from a cold-beam instability when vthermal < vdrift and when vdrift < vcritical where vcritical is some algorithm-dependent critical velocity. Notably, all algorithms are stable to grid-heating (vdrift = 0) instabilities. This work confirms the semi-analytic work of Barnes and Chacon (2021), and extends their results by considering algorithms that have continuous first derivatives of the interpolated electric field and algorithms with a variety of field solves. We present an argument for the general failure of field solves to prevent cold-beam instabilities and derive vcritical analytically for some algorithms.

Presenters

  • Luke C Adams

    University of Colorado, Boulder

Authors

  • Luke C Adams

    University of Colorado, Boulder

  • Gregory R Werner

    University of Colorado, Boulder

  • John Robert Cary

    University of Colorado, Boulder, Tech-X Corporation