On Energy and Momentum Conservation in Particle-in-Cell Simulation

Abstract

Particle-in-cell (PIC) plasma simulations are a productive and valued tool for the study of nonlinear plasma phenomena, yet there are basic questions about the simulation methods themselves that remain unanswered. Here we study one such question: energy and momentum conservation by PIC. We employ both analysis and simulations of one-dimensional, electrostatic plasmas to understand why PIC simulations are either energy or momentum conserving but not both, what the role of numerical stability is in non-conservation, and how do errors in conservation scale with the numerical parameters. Conserving both momentum and energy make it possible to model problems such as Jeans' -type equilibria. Avoiding numerical instability is useful, but so is being able to identify when its effect on the results may be important. Designing simulations to achieve the best possible accuracy with the least expenditure of effort requires results on the scaling of error with the numerical parameters.. Our results identify the central role of Gauss' law in conservation of both momentum and energy, and the significant differences in numerical stability and error scaling between energy-conserving and momentum-conserving simulations.

Cited by in corpus (14)

• Exactly Energy Conserving Semi-Implicit Particle in Cell Formulation
• SHARP: A Spatially Higher-order, Relativistic Particle-in-Cell Code
• A dynamical adaptive tensor method for the Vlasov-Poisson system
• A fully implicit, conservative, non-linear, electromagnetic hybrid particle-ion/fluid-electron algorithm
• A new hybrid code (CHIEF) implementing the inertial electron fluid equation without approximation
• Modeling of very high frequency large-electrode capacitively coupled plasmas with a fully electromagnetic particle-in-cell code
• Time-step dependent force interpolation scheme for suppressing numerical Cherenkov instability in relativistic particle-in-cell simulations
• On the validity of quasilinear theory applied to the electron bump-on-tail instability
• Arbitrary-order time-accurate semi-Lagrangian spectral approximations of the Vlasov-Poisson system
• A Semi-Lagrangian Spectral Method for the Vlasov-Poisson System based on Fourier, Legendre and Hermite Polynomials
• Moment tracking and their coordinate transformations for macroparticles with an application to plasmas around black holes
• A Particle-in-cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition
• Grid instability growth rates for explicit, electrostatic momentum- and energy-conserving particle-in-cell algorithms
• Finite spatial-grid effects in energy-conserving particle-in-cell algorithms