A DUAL GRID GEOMETRIC ELECTROMAGNETIC PARTICLE IN CELL METHOD

Geometric particle-in-cell discretizations have been derived based on a discretization of the fields that is conforming with the de Rham structure of the Maxwell's equations and a standard particle-in-cell ansatz for the fields by deriving the equations of motion from a discrete action principle. While earlier work has focused on finite-element discretization of the fields based on the theory of finite-element exterior calculus, we propose in this article an alternative formulation of the field equations that is based on the ideas conveyed by mimetic finite differences, the needed duality being expressed by the use of staggered grids. We construct a finite-difference formulation based on degrees of freedom defined as point values, edge, face, and volume integrals on a primal and its dual grid. Compared to the finite-element formulation, no mass matrix inversion is involved in the formulation of the Maxwell solver. In numerical experiments, we verify the conservation properties of the novel method and study the influence of the various parameters in the discretization.