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

TL;DR

This paper presents a spectral element numerical method for simulating streamer discharges in low-temperature atmospheric-pressure plasmas, combining spectral discretization, hierarchical solvers, and adaptive mesh refinement for efficiency and flexibility.

Contribution

The paper introduces a novel spectral element approach for streamer modeling, integrating the Poisson and electron continuity equations with advanced stabilization and adaptive techniques.

Findings

Method is validated against test problems and previous studies.

Simulation scheme is sufficiently fast for practical applications.

Flexibility allows extension to unstructured meshes and broader problems.

Abstract

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 Poincare-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.