Explicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems
- School of Nuclear Science and Technology and Department of Modern Physics, University of Science and Technology of China, Hefei, Anhui 230026 (China)
- LSEC, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, P.O. Box 2719, Beijing 100190 (China)
Explicit high-order non-canonical symplectic particle-in-cell algorithms for classical particle-field systems governed by the Vlasov-Maxwell equations are developed. The algorithms conserve a discrete non-canonical symplectic structure derived from the Lagrangian of the particle-field system, which is naturally discrete in particles. The electromagnetic field is spatially discretized using the method of discrete exterior calculus with high-order interpolating differential forms for a cubic grid. The resulting time-domain Lagrangian assumes a non-canonical symplectic structure. It is also gauge invariant and conserves charge. The system is then solved using a structure-preserving splitting method discovered by He et al. [preprint http://arxiv.org/abs/arXiv:1505.06076 (2015)], which produces five exactly soluble sub-systems, and high-order structure-preserving algorithms follow by combinations. The explicit, high-order, and conservative nature of the algorithms is especially suitable for long-term simulations of particle-field systems with extremely large number of degrees of freedom on massively parallel supercomputers. The algorithms have been tested and verified by the two physics problems, i.e., the nonlinear Landau damping and the electron Bernstein wave.
- OSTI ID:
- 22489840
- Journal Information:
- Physics of Plasmas, Journal Name: Physics of Plasmas Journal Issue: 11 Vol. 22; ISSN PHPAEN; ISSN 1070-664X
- Country of Publication:
- United States
- Language:
- English
Similar Records
Gauge properties of the guiding center variational symplectic integrator
Field theory and a structure-preserving geometric particle-in-cell algorithm for drift wave instability and turbulence