Time Integrator Agnostic Charge Conserving Finite Element PIC

Abstract

Developing particle-in-cell (PIC) methods using finite element basis sets, and without auxiliary divergence cleaning methods, was a long standing problem until recently. It was shown that if consistent spatial basis functions are used, one can indeed create a methodology that was charge conserving, albeit using a leap-frog time stepping method. While this is a significant advance, leap frog schemes are only conditionally stable and time step sizes are closely tied to the underlying mesh. Ideally, to take full advantage of advances in finite element methods (FEMs), one needs a charge conserving PIC methodology that is agnostic to the time stepping method. This is the principal contribution of this paper. In what follows, we shall develop this methodology, prove that both charge and Gauss' laws are discretely satisfied at every time step, provide the necessary details to implement this methodology for both the wave equation FEM and Maxwell Solver FEM, and finally demonstrate its efficacy on a suite of test problems. The method will be demonstrated by single particle evolution, non-neutral beams with space-charge, and adiabatic expansion of a neutral plasma, where the debye length has been resolved, and real mass ratios are used.

arXiv commentThe following article has been submitted to AIP Physics of Plasmas

References in corpus (8)

• An Energy- and Charge-conserving, Implicit, Electrostatic Particle-in-Cell Algorithm
• Geometric integration of the Vlasov-Maxwell system with a variational particle-in-cell scheme
• Exact Charge-Conserving Scatter-Gather Algorithm for Particle-in-Cell Simulations on Unstructured Grids: A Geometric Perspective
• Energy-conserving time propagation for a geometric particle-in-cell Vlasov--Maxwell solver
• 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
• The geometric theory of charge conservation in particle-in-cell simulations
• Rubrics for Charge Conserving Current Mapping in Finite Element Particle in Cell Methods

Cited by in corpus (3)

• Rubrics for Charge Conserving Current Mapping in Finite Element Particle in Cell Methods
• A Charge Conserving Exponential Predictor Corrector FEMPIC Formulation for Relativistic Particle Simulations
• A gauge-compatible Hamiltonian splitting algorithm for particle-in-cell simulations using finite element exterior calculus