Andrew J. Christlieb,William A. Sands,Stephen White
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.
What problem does this paper attempt to address?
The main problems that this paper attempts to solve are the limitations in numerical heating, stability, and geometric adaptability of the traditional Particle - in - Cell (PIC) method in plasma simulations. Specifically:
1. **Numerical heating problem**: The traditional PIC method will have grid - related numerical heating phenomena when the grid spacing is close to the Debye length, which will affect the accuracy of the simulation.
2. **Stability problem**: When the grid spacing exceeds the Debye length, the traditional method may become unstable, especially when dealing with complex geometric structures.
3. **Geometric adaptability problem**: The traditional method needs to use cut - cells or staircase approximations when dealing with complex geometric shapes, which will introduce additional errors and computational complexity.
To solve these problems, the paper proposes a new Particle - in - Cell method based on the generalized momentum formulation, which has the following characteristics:
- **Grid - independent numerical heating**: The new method can maintain stable numerical performance when the grid spacing is close to the Debye length, avoiding the numerical heating problem.
- **High - order spatial approximation**: A high - order spatial discretization method is adopted, which reduces the number of required grid points, thereby improving computational efficiency.
- **Non - staggered grid**: A non - staggered grid is used, so that the field and particle data are co - located on the same grid point, simplifying the data mapping process.
- **Hamiltonian formulation**: Through the Hamiltonian formulation of the generalized momentum (i.e., conjugate momentum), the need for time derivatives under the Lorentz gauge is eliminated, ensuring the time accuracy of the field solver.
- **Analytical calculation of spatial derivatives**: The integral solution method is used to analytically calculate the spatial derivatives of the field, ensuring that the convergence rate of these derivatives is the same as that of the field itself.
### Main improvement points
1. **Elimination of the need for time derivatives**: Through the generalized momentum formulation, the need for time derivatives in the standard Newton - Lorentz treatment is avoided, thus retaining the full - time accuracy of the field solver.
2. **High - order spatial discretization**: A fifth - order spatial discretization method is adopted, which reduces the number of required grid points and further saves computational resources.
3. **Analytical calculation of spatial derivatives**: The integral solution method is used to analytically calculate the spatial derivatives of the field, ensuring that the convergence rate of these derivatives is the same as that of the field itself, improving numerical accuracy.
4. **Applicable to complex geometric structures**: This method can be applied to complex geometric structures without using cut - cells or staircase approximations, thereby improving geometric adaptability and computational efficiency.
### Application examples
The paper verifies the effectiveness of the new method through several classic benchmark test problems, including plasma sheath and relativistic particle beam. The results show that the new method performs well in these tests, especially when the grid spacing is close to the Debye length, it can still maintain stable numerical performance.
In summary, the new method proposed in this paper aims to improve the stability and accuracy of the existing PIC method in plasma simulations, especially in application scenarios dealing with complex geometric structures and high - precision requirements.