A spectral element method for modelling streamer discharges in low-temperature atmospheric-pressure plasmas

Streamers are ionization fronts that occur in gases at atmospheric and sub-atmospheric pressures. Numerical studies of streamers are important for practical applications but are challenging due to the multiscale nature of this discharge type. This paper introduces a spectral element method for modelling streamer discharges. The method is developed for Cartesian grids but can be extended to be used on unstructured meshes. The streamer model is based on the Poisson equation for the electric potential and the electron continuity equation. The Poisson equation is discretized via a spectral method based on the integral representation of the solution. The hierarchical Poincaré-Steklov (HPS) scheme is used to solve the resulting set of equations. The electron continuity equation is solved by means of the discontinuous Galerkin spectral element method (DGSEM). The DGSEM is extended by an alternative definition of the diffusion flux. A subcell finite volume method is used to stabilize the DGSEM scheme, if required. The entire simulation scheme is validated by solving a number of test problems and reproducing the results of previous studies. Adaptive mesh refinement is used to reduce the number of unknowns. The proposed method is found to be sufficiently fast for being used in practical applications. The flexibility of the method provides an interesting opportunity to broaden the range of problems that can be addressed in numerical studies of low-temperature plasma discharges.