Comparison of coupling techniques in a high-order discontinuous Galerkin-based particle-in-cell solver

Highly rarefied plasma flows in technical devices are physically modelled by the Maxwell–Lorentz equations. They combine the solution of the Maxwell equations, where the electric field E and magnetic induction B are determined, with the Lorentz system, accounting for the movement of charged particles due to the electromagnetic forces. To solve these equations for complex-shaped domains, a fully electromagnetic particle-in-cell (PIC) code has been developed using high-order discontinuous Galerkin methods for the Maxwell equations on a computational mesh, coupled with a Lorentz solver on the basis of a second-order leapfrog scheme, acting on the particles at their current positions. Since the particles move freely in space, the mesh-based and the mesh-free values have to be coupled. This coupling includes the deposition of the charge and current densities from the current particle positions onto the mesh as well as the interpolation of the electromagnetic fields from the mesh to the actual particle positions. Both steps have to be computed with appropriate accuracy. Different approaches to particle-grid coupling within the PIC solver have been investigated. In this paper, these concepts are described and corresponding simulation results with respect to accuracy and computational demand are presented.