Conservative Projection Between Finite Element and Particle Bases

Abstract

Particle-in-Cell (PIC) methods employ particle representations of unknown fields, but also employ continuum fields for other parts of the problem. Thus projection between particle and continuum bases is required. Moreover, we often need to enforce conservation constraints on this projection. We derive a mechanism for enforcement based on weak equality, and implement it in the PETSc libraries. Scalability is demonstrated to more than 1B particles.

References in corpus (7)

• An Energy- and Charge-conserving, Implicit, Electrostatic Particle-in-Cell Algorithm
• The Energy Conserving Particle-in-Cell Method
• Structure and structure-preserving algorithms for plasma physics
• Efficient mesh management in Firedrake using PETSc-DMPlex
• Metriplectic Integrators for the Landau Collision Operator
• Conservative discretization of the Landau collision integral
• Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC