Finite grid instability and spectral fidelity of the electrostatic Particle-In-Cell algorithm
arXiv:1508.03360 · doi:10.1016/j.cpc.2016.05.021
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 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.
References in corpus (5)
- Variational formulation of particle algorithms for kinetic plasma simulations
- Explicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems
- Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm
- Numerical instability due to relativistic plasma drift in EM-PIC simulations
- On the Numerical Dispersion of Electromagnetic Particle-In-Cell Code : Finite Grid Instability
Cited by in corpus (13)
- A fully implicit, conservative, non-linear, electromagnetic hybrid particle-ion/fluid-electron algorithm
- Explicit Structure-Preserving Geometric Particle-in-Cell Algorithm in Curvilinear Orthogonal Coordinate Systems and Its Applications to Whole-Device 6D Kinetic Simulations of Tokamak Physics
- Accurately simulating nine-dimensional phase space of relativistic particles in strong fields
- Efficient Fourier Basis Particle Simulation
- Time-step dependent force interpolation scheme for suppressing numerical Cherenkov instability in relativistic particle-in-cell simulations
- Local Energy Conservation Law for Spatially-Discretized Hamiltonian Vlasov-Maxwell System
- Fundamental form of the electrostatic -PIC algorithm and discovery of a converged numerical instability
- A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm
- Symplectic particle-in-cell methods for hybrid plasma models with Boltzmann electrons and space-charge effects
- Particle Acceleration in Relativistic Electron-Ion Outlfows
- Finite spatial-grid effects in energy-conserving particle-in-cell algorithms
- An \mathcal{O}(N) Maxwell solver with improved numerical dispersion properties
- On Geometric Fourier Particle In Cell Methods