Rubrics for Charge Conserving Current Mapping in Finite Element Particle in Cell Methods
arXiv:2101.12128 · doi:10.1109/TPS.2021.3122410
Abstract
Modeling of kinetic plasmas using electromagnetic particle in cell methods (EM-PIC) is a problem that is well worn, in that methods developed have been used extensively both understanding physics and exploiting them for device design. EM-PIC tools have largely relied on finite difference methods coupled with particle representations of the distribution function. Refinements to ensure consistency and charge conservation have largely been an ad-hoc efforts specific to finite difference methods. Meanwhile, solution methods for field solver have grown by leaps and bounds with significant performance metrics compared to finite difference methods. Developing new EM-PIC computational schemes that leverage modern field solver technology means re-examining analysis framework necessary for self-consistent EM-PIC solution. In this paper, we prescribe general rubrics for charge conservation, demonstrate how these are satisfied in conventional finite difference PIC as well as finite element PIC, and prescribe a novel charge conserving finite element PIC. Our effort leverages proper mappings on to de-Rham sequences and lays a groundwork for understanding conditions that must be satisfied for consistency. Several numerical results demonstrate the applicability of these rubrics.
References in corpus (4)
- Explicit Structure-Preserving Geometric Particle-in-Cell Algorithm in Curvilinear Orthogonal Coordinate Systems and Its Applications to Whole-Device 6D Kinetic Simulations of Tokamak Physics
- A Set of Benchmark Tests for Validation of 3D Particle In Cell Methods
- Time Integrator Agnostic Charge Conserving Finite Element PIC
- Quasi-Helmholtz Decomposition, Gauss' Laws and Charge Conservation for Finite Element Particle-in-Cell