A Particle-in-cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

TL;DR

This paper introduces a novel particle-in-cell method for plasma simulation using a Hamiltonian formulation with generalized momentum, achieving high accuracy and efficiency, especially near the Debye length scale.

Contribution

The paper develops a new particle-in-cell approach employing a Hamiltonian formulation with generalized momentum and high-order spatial derivatives, improving accuracy and computational efficiency.

Findings

Mesh-independent numerical heating observed

Fewer particles needed per cell for accuracy

Method demonstrates high-order convergence in space

Abstract

This paper formulates a new particle-in-cell method for the Vlasov-Maxwell system. Under the Lorenz gauge condition, Maxwell's equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized momentum. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are shown to converge at the same rate as the fields in the both time and space. The field solver we consider in this work is first-order accurate in time and fifth-order accurate in space and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems, and the efficacy of the proposed formulation is demonstrated through a comparison with standard methods presented in the literature. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is close to the grid spacing. The use of high-order spatial approximations in the new method means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to standard explicit methods, which permits further computational savings.