Numerical simulations for a hybrid model of kinetic ions and mass-less fluid electrons in canonical formulations
We study the structure-preserving discretizations of a hybrid model with kinetic ions and mass-less electrons. Different from most existing works in the literature, we conduct the discretizations based on two equivalent formulations with vector potentials in different gauges, and the distribution functions depend on canonical momentum (not velocity). Particle-in-cell methods are used for the distribution functions, and vector potentials are discretized by finite element methods in the framework of finite element exterior calculus. Splitting methods are used for time discretizations. For the first formulation, filters are used to reduce the noises from particles and are shown to improve the numerical results significantly. The schemes of the second formulation show good stability and accuracy because of the use of symplectic methods for canonical Hamiltonian systems. Magnetic fields obtained from the vector potentials are divergence-free naturally. Some numerical experiments are conducted to validate and compare the two discretizations.
READ FULL TEXT