Charge-conserving, variational particle-in-cell method for the drift-kinetic Vlasov-Maxwell system

This paper proposes a charge-conserving, variational, spatio-temporal discretization for the drift-kinetic Vlasov-Maxwell system, utilizing finite-elements for the electromagnetic fields and the particle-in-cell approach for the Vlasov distribution. The proposed scheme is fully electromagnetic, dealing with fields instead of potentials, and includes the effects of polarization and magnetization in the Gauss and Ampere-Maxwell laws, a consequence of reducing the full particle dynamics to drift-center dynamics. There is, however, no need to invert the Gauss law: it is satisfied automatically at every time-step as a result of a discrete Noether symmetry, and the electric field is updated directly from the Ampere-Maxwell equation. The method provides an update for the magnetic field that is fully explicit, involving only local operations. The update for particles is implicit for each particle individually, also leading to local operations only. The update for the electric field is linearly implicit due to the presence of a finite-element mass matrix and polarization and magnetization effects in the Ampere-Maxwell equation, hence involving a sparse matrix inversion once at every time step. Because the scheme deals with the electromagnetic fields and not the potentials, it also provides the first serious attempt at constructing a structure-preserving numerical scheme for the mixed kinetic-ion--drift-kinetic-electron Vlasov-Maxwell model. Consequently, the proposed method could be used to simulate electromagnetic turbulence in fusion experiments or space plasmas that exhibit a strong background magnetic field while retaining all of the ion physics, most of the necessary electron physics, yet eliminating perhaps the biggest obstacle in reaching macroscopic transport time scales in kinetic simulations, namely the electron cyclotron time scale.