Finite grid instability and spectral fidelity of the electrostatic Particle-In-Cell algorithm

Abstract The origin of the Finite Grid Instability (FGI) is studied by resolving the dynamics in the 1D electrostatic Particle-In-Cell (PIC) model in the spectral domain at the single particle level and at the collective motion level. The spectral fidelity of the PIC model is contrasted with the underlying physical system or the gridless model. The systematic spectral phase and amplitude errors from the charge deposition and field interpolation are quantified for common particle shapes used in the PIC models. It is shown through such analysis and in simulations that the lack of spectral fidelity relative to the physical system due to the existence of aliased spatial modes is the major cause of the FGI in the PIC model.

[1]  Jean-Luc Vay,et al.  Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm , 2012, J. Comput. Phys..

[2]  Charles K. Birdsall,et al.  Plasma self-heating and saturation due to numerical instabilities☆ , 1980 .

[3]  C. Birdsall,et al.  Plasma Physics via Computer Simulation , 2018 .

[4]  Jianyuan Xiao,et al.  A variational multi-symplectic particle-in-cell algorithm with smoothing functions for the Vlasov-Maxwell system , 2013 .

[5]  Wei Lu,et al.  Numerical instability due to relativistic plasma drift in EM-PIC simulations , 2012, Comput. Phys. Commun..

[6]  A. Stamm,et al.  Variational formulation of particle algorithms for kinetic electromagnetic plasma simulations , 2013, 2013 Abstracts IEEE International Conference on Plasma Science (ICOPS).

[7]  H. Lewis,et al.  Energy-conserving numerical approximations for Vlasov plasmas , 1970 .

[8]  Liu Chen,et al.  I. Heating of Magnetized Plasmas by Large-Amplitude Electric Field. I. Reduction of the Grid Effects in Simulation Plasmas. , 1974 .

[9]  Jianyuan Xiao,et al.  Explicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems , 2015, 1510.06972.

[10]  Brendan B. Godfrey,et al.  Numerical Cherenkov instabilities in electromagnetic particle codes , 1974 .

[11]  Brian J. Albright,et al.  On the numerical dispersion of electromagnetic particle-in-cell code: Finite grid instability , 2015, J. Comput. Phys..

[12]  Hideo Okuda,et al.  Nonphysical noises and instabilities in plasma simulation due to a spatial grid , 1972 .

[13]  Paolo Ricci,et al.  Kinetic approach to microscopic-macroscopic coupling in space and laboratory plasmas , 2006 .

[14]  H. Lewis,et al.  A comparison of some particle-in-cell plasma simulation methods , 1972 .

[15]  Giovanni Lapenta,et al.  A method to suppress the finite-grid instability in plasma simulations , 1994 .

[16]  A. Bruce Langdon,et al.  EFFECTS OF THE SPATIAL GRID IN SIMULATION PLASMAS. , 1970 .

[17]  Luis Chacón,et al.  An energy- and charge-conserving, implicit, electrostatic particle-in-cell algorithm , 2011, J. Comput. Phys..

[18]  Stefano Markidis,et al.  The energy conserving particle-in-cell method , 2011, J. Comput. Phys..

[19]  Luis Chacón,et al.  Development of a Consistent and Stable Fully Implicit Moment Method for Vlasov-Ampère Particle in Cell (PIC) System , 2013, SIAM J. Sci. Comput..

[20]  B. D. Martino,et al.  Gridless finite-size-particle plasma simulation , 2001 .

[21]  E. Lindman,et al.  DISPERSION RELATION FOR COMPUTER-SIMULATED PLASMAS. , 1970 .