An Asymptotic-Preserving and Energy-Conserving Particle-In-Cell Method for Vlasov-Maxwell Equations
arXiv:2209.08227 · doi:10.1063/5.0126389
Abstract
In this paper, we develop an asymptotic-preserving and energy-conserving (APEC) Particle-In-Cell (PIC) algorithm for the Vlasov-Maxwell system. This algorithm not only guarantees that the asymptotic limiting of the discrete scheme is a consistent and stable discretization of the quasi-neutral limit of the continuous model, but also preserves Gauss's law and energy conservation at the same time, thus it is promising to provide stable simulations of complex plasma systems even in the quasi-neutral regime. The key ingredients for achieving these properties include the generalized Ohm's law for electric field such that the asymptotic-preserving discretization can be achieved, and a proper decomposition of the effects of the electromagnetic fields such that a Lagrange multiplier method can be appropriately employed for correcting the kinetic energy. We investigate the performance of the APEC method with three benchmark tests in one dimension, including the linear Landau damping, the bump-on-tail problem and the two-stream instability. Detailed comparisons are conducted by including the results from the classical explicit leapfrog and the previously developed asymptotic-preserving PIC schemes. Our numerical experiments show that the proposed APEC scheme can give accurate and stable simulations both kinetic and quasi-neutral regimes, demonstrating the attractive properties of the method crossing scales.
arXiv comment21 pages, 30 figures
References in corpus (14)
- SMILEI: a collaborative, open-source, multi-purpose particle-in-cell code for plasma simulation
- An Energy- and Charge-conserving, Implicit, Electrostatic Particle-in-Cell Algorithm
- An Asymptotic-Preserving all-speed scheme for the Euler and Navier-Stokes equations
- Variational formulation of particle algorithms for kinetic plasma simulations
- A Hamiltonian splitting for the Vlasov-Maxwell system
- Energy-conserving discontinuous Galerkin methods for the Vlasov-Ampère system
- Asymptotic-Preserving methods and multiscale models for plasma physics
- Numerical approximation of the Euler-Maxwell model in the quasineutral limit
- A semi-implicit, energy- and charge-conserving particle-in-cell algorithm for the relativistic Vlasov-Maxwell equations
- Numerical approximation of the Euler-Poisson-Boltzmann model in the quasineutral limit
- The multi-dimensional Hermite-discontinuous Galerkin method for the Vlasov-Maxwell equations
- Quasi-Helmholtz Decomposition, Gauss' Laws and Charge Conservation for Finite Element Particle-in-Cell
- Numerical Study of the Two-Species Vlasov-Ampère System: Energy-Conserving Schemes and the Current-Driven Ion-Acoustic Instability
- Low-frequency whistler waves excited by relativistic laser pulses