A dual grid geometric electromagnetic particle in cell method

Abstract

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 equation 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.

27 pages, 8 figures

References in corpus (9)

• Finite element exterior calculus: from Hodge theory to numerical stability
• GEMPIC: Geometric ElectroMagnetic Particle-In-Cell Methods
• A Hamiltonian splitting for the Vlasov-Maxwell system
• Physics-compatible discretization techniques on single and dual grids, with application to the Poisson equation of volume forms
• Detailed analysis of the effects of stencil spatial variations with arbitrary high-order finite-difference Maxwell solver
• Hamiltonian integration methods for Vlasov-Maxwell equations
• Energy-conserving time propagation for a geometric particle-in-cell Vlasov--Maxwell solver
• Comment on "Hamiltonian splitting for the Vlasov-Maxwell equations"
• Subcycling of particle orbits in variational, geometric electromagnetic particle-in-cell methods