Exact Charge-Conserving Scatter-Gather Algorithm for Particle-in-Cell Simulations on Unstructured Grids: A Geometric Perspective

Abstract

We describe a charge-conserving scatter-gather algorithm for particle-in-cell simulations on unstructured grids. Charge conservation is obtained from first principles, i.e., without the need for any post-processing or correction steps. This algorithm recovers, at a fundamental level, the scatter-gather algorithms presented recently by Campos-Pinto et al. [1] (to first-order) and by Squire et al. [2], but it is derived here in a streamlined fashion from a geometric viewpoint. Some ingredients reflecting this viewpoint are (1) the use of (discrete) differential forms of various degrees to represent fields, currents, and charged particles and provide localization rules for the degrees of freedom thereof on the various grid elements (nodes, edges, facets), (2) use of Whitney forms as basic interpolants from discrete differential forms to continuum space, and (3) use of a Galerkin formula for the discrete Hodge star operators (i.e., "mass matrices" incorporating the metric datum of the grid) applicable to generally irregular, unstructured grids. The expressions obtained for the scatter charges and scatter currents are very concise and do not involve numerical quadrature rules. Appropriate fractional areas within each grid element are identified that represent scatter charges and scatter currents within the element, and a simple geometric representation for the (exact) charge conservation mechanism is obtained by such identification. The field update is based on the coupled first-order Maxwell's curl equations to avoid spurious modes with secular growth (otherwise present in formulations that discretize the second-order wave equation). Examples are provided to verify preservation of discrete Gauss' law for all times.

References in corpus (2)

• Fast-ignition design transport studies: realistic electron source, integrated PIC-hydrodynamics, imposed magnetic fields
• Alternating-Order Interpolation in a Charge-Conserving Scheme for Particle-In-Cell Simulations

Cited by in corpus (16)

• GEMPIC: Geometric ElectroMagnetic Particle-In-Cell Methods
• Structure-Preserving Geometric Particle-in-Cell Methods for Vlasov-Maxwell Systems
• Axisymmetric Charge-Conservative Electromagnetic Particle Simulation Algorithm on Unstructured Grids: Application to Microwave Vacuum Electronic Devices
• Relativistic Extension of a Charge-Conservative Finite Element Solver for Time-Dependent Maxwell-Vlasov Equations
• 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
• Detection and Prediction of Equilibrium States in Kinetic Plasma Simulations via Mode Tracking using Reduced-Order Dynamic Mode Decomposition
• A Set of Benchmark Tests for Validation of 3D Particle In Cell Methods
• Numerical Electromagnetic Frequency Domain Analysis with Discrete Exterior Calculus
• Rubrics for Charge Conserving Current Mapping in Finite Element Particle in Cell Methods
• Finite Element Time-Domain Body-of-Revolution Maxwell Solver based on Discrete Exterior Calculus
• Time Integrator Agnostic Charge Conserving Finite Element PIC
• Quasi-Helmholtz Decomposition, Gauss' Laws and Charge Conservation for Finite Element Particle-in-Cell
• Geometric Electrostatic Particle-In-Cell Algorithm on Unstructured Meshes
• Accelerating Particle-in-Cell Kinetic Plasma Simulations via Reduced-Order Modeling of Space-Charge Dynamics using Dynamic Mode Decomposition
• A gauge-compatible Hamiltonian splitting algorithm for particle-in-cell simulations using finite element exterior calculus
• A Charge Conserving Exponential Predictor Corrector FEMPIC Formulation for Relativistic Particle Simulations