A Particle-in-Cell solver based on a high-order hybridizable discontinuous Galerkin spectral element method on unstructured curved meshes

Abstract A high-order hybridizable discontinuous Galerkin spectral element method (HDGSEM) for Particle-In-Cell (PIC) schemes is presented for the simulation of electrostatic applications on three-dimensional unstructured curved meshes. The electrostatic Poisson equation is solved and optionally a Boltzmann relation for the electron species can be used which leads to non-linear source terms. The hybridizable formulation reduces the total number of unknowns of the field solver, allowing the simulation of large problems. The implementation of the HDGSEM solver in a PIC code is described and validated using several test cases with successively increasing complexity. It is shown that the high-order convergence properties are retained on curvilinear meshes, likewise when material jumps are introduced. The simulation of an ion optic illustrates the applicability of the presented method for complex geometries and large problem sizes.