A semi-implicit, energy- and charge-conserving particle-in-cell algorithm for the relativistic Vlasov-Maxwell equations

Abstract

Conventional explicit electromagnetic particle-in-cell (PIC) algorithms do not conserve discrete energy exactly. Time-centered fully implicit PIC algorithms can conserve discrete energy exactly, but may introduce large dispersion errors in the light-wave modes. This can lead to intolerable simulation errors where accurate light propagation is needed (e.g. in laser-plasma interactions). In this study, we selectively combine the leap-frog and Crank-Nicolson methods to produce an exactly energy- and charge-conserving relativistic electromagnetic PIC algorithm. Specifically, we employ the leap-frog method for Maxwell's equations, and the Crank-Nicolson method for the particle equations. The semi-implicit algorithm admits exact global energy conservation, exact local charge conservation, and preserves the dispersion properties of the leap-frog method for the light wave. The algorithm employs a new particle pusher designed to maximize efficiency and minimize wall-clock-time impact vs. the explicit alternative. It has been implemented in a code named iVPIC, based on the Los Alamos National Laboratory VPIC code (\url{https://github.com/losalamos/vpic}). We present numerical results that demonstrate the properties of the scheme with sample test problems: relativistic two-stream instability, Weibel instability, and laser-plasma instabilities.

References in corpus (8)

• An Energy- and Charge-conserving, Implicit, Electrostatic Particle-in-Cell Algorithm
• The Energy Conserving Particle-in-Cell Method
• Exactly Energy Conserving Semi-Implicit Particle in Cell Formulation
• Structure-preserving second-order integration of relativistic charged particle trajectories in electromagnetic fields
• A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov-Darwin particle-in-cell algorithm
• An energy- and charge-conserving, nonlinearly implicit, electromagnetic 1D-3V Vlasov-Darwin particle-in-cell algorithm
• Gauss's Law Satisfying Energy-Conserving Semi-Implicit Particle-in-Cell Method
• A fully implicit, conservative, non-linear, electromagnetic hybrid particle-ion/fluid-electron algorithm

Cited by in corpus (10)

• Energy-conserving time propagation for a geometric particle-in-cell Vlasov--Maxwell solver
• Enabling particle applications for exascale computing platforms
• Accurately simulating nine-dimensional phase space of relativistic particles in strong fields
• An implicit, conservative and asymptotic-preserving electrostatic particle-in-cell algorithm for arbitrarily magnetized plasmas in uniform magnetic fields
• Rubrics for Charge Conserving Current Mapping in Finite Element Particle in Cell Methods
• An Asymptotic-Preserving and Energy-Conserving Particle-In-Cell Method for Vlasov-Maxwell Equations
• Subcycling of particle orbits in variational, geometric electromagnetic particle-in-cell methods
• Energy conserving particle-in-cell methods for relativistic Vlasov--Maxwell equations of laser-plasma interaction
• Explicit energy-conserving modification of relativistic PIC method
• Charge-conserving, variational particle-in-cell method for the drift-kinetic Vlasov-Maxwell system