Energy-conserving discontinuous Galerkin methods for the Vlasov–Maxwell system

• Published in: Journal of computational physics Vol. 279; pp. 145 - 173
• Main Authors: Cheng, Yingda, Christlieb, Andrew J., Zhong, Xinghui
• Format: Journal Article
• Language: English
• Published: Elsevier Inc 15-12-2014
• Subjects: Conservation; Discontinuous Galerkin methods; Discretization; Energy conservation; Galerkin methods; Magnetic fields; Mathematical models; Maxwell's equations; Plasmas; Streaming Weibel instability; Symplectic integrators; Temporal logic; Vlasov–Maxwell system; Vlasov–Maxwell system; Symplectic integrators; Discontinuous Galerkin methods; Streaming Weibel instability; Energy conservation
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2014.08.041

In this paper, we generalize the idea in our previous work for the Vlasov–Ampère (VA) system (Y. Cheng, A.J. Christlieb, and X. Zhong (2014) [10]) and develop energy-conserving discontinuous Galerkin (DG) methods for the Vlasov–Maxwell (VM) system. The VM system is a fundamental model in the simulation of collisionless magnetized plasmas. Compared to Y. Cheng, A.J. Christlieb, and X. Zhong (2014) [10], additional care needs to be taken for both the temporal and spatial discretizations to achieve similar type of conservation when the magnetic field is no longer negligible. Our proposed schemes conserve the total particle number and the total energy at the same time, therefore can obtain accurate and physically relevant numerical solutions. The main components of our methods include second order and above, explicit or implicit energy-conserving temporal discretizations, and DG methods for Vlasov and Maxwell's equations with carefully chosen numerical fluxes. Benchmark numerical tests such as the streaming Weibel instability are provided to validate the accuracy and conservation of the schemes.