A Particle-in-cell Method for Plasmas with A Generalized Momentum Formulation, Part III: A family of Gauge Conserving Methods

In this paper, we introduce a new family of spatially co-located field solvers for particle-in-cell applications which evolve the potential formulation of Maxwell’s equations under the Lorenz gauge. Our recent work [2] introduced the concept of time-consistency, which connects charge conservation to the preservation of the gauge at the semi-discrete level. It will be shown that there exists a large family of time discretizations which satisfy this property. Additionally, it will be further shown that for large classes of time marching methods, the satisfaction of the gauge condition automatically implies the satisfaction of Gauss’s law for electricity, with the potential formulation ensuring that that Gauss’s law for magnetism is satisfied by definition. We focus on popular time marching methods including centered differences, backward differences, and diagonally-implicit Runge-Kutta methods, which are coupled to a spectral discretization in space. We demonstrate the theory by testing the methods on a relativistic Weibel instability and a drifting cloud of electrons.. Authors: Christlieb, Andrew J. [Michigan State Univ., East Lansing, MI (United States)]; Sands, William A. [Univ. of Delaware, Newark, DE (United States)]; White, Stephen R. [Michigan State Univ., East Lansing, MI (United States)] (ORCID:0000000283088110). DOE Contract: SC0023164. Subjects: 97 MATHEMATICS AND COMPUTING; Vlasov-Maxwell system; generalized momentum; integral solution; method-of-lines-transpose; particle-in-cell