Finite spatial-grid effects in energy-conserving particle-in-cell algorithms
From MaRDI portal
Publication:6041005
Abstract: Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space-time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to resolve Debye lengths spatially. We also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas.
Cites work
- A charge- and energy-conserving implicit, electrostatic particle-in-cell algorithm on mapped computational meshes
- A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov-Darwin particle-in-cell algorithm
- An energy- and charge-conserving, implicit, electrostatic particle-in-cell algorithm
- An energy- and charge-conserving, nonlinearly implicit, electromagnetic 1D-3V Vlasov-Darwin particle-in-cell algorithm
- An implicit method for electromagnetic plasma simulation in two dimensions
- Anderson acceleration for fixed-point iterations
- Development of a consistent and stable fully implicit moment method for Vlasov-Ampère particle in cell (PIC) system
- Energy-conserving numerical approximations for Vlasov plasmas
- Exactly energy conserving semi-implicit particle in cell formulation
- Finite grid instability and spectral fidelity of the electrostatic particle-in-cell algorithm
- Jacobian-free Newton-Krylov methods: a survey of approaches and applications.
- Numerical computation of polynomial zeros by means of Aberth's method
- On energy and momentum conservation in particle-in-cell plasma simulation
- On the numerical dispersion of electromagnetic particle-in-cell code: finite grid instability
- Plasma self-heating and saturation due to numerical instabilities
- The energy conserving particle-in-cell method
Cited in
(5)- A pseudospectral implicit particle-in-cell method with exact energy and charge conservation
- Semi-implicit particle-in-cell methods embedding sparse grid reconstructions
- An energy- and charge-conserving electrostatic implicit particle-in-cell algorithm for simulations of collisional bounded plasmas
- On the Effect of Angular and Spatial Discretization on Perturbation Calculations
- Eliminating finite-grid instabilities in gyrokinetic particle-in-cell simulations
This page was built for publication: Finite spatial-grid effects in energy-conserving particle-in-cell algorithms
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6041005)
