A P-Adaptive Hybridizable Discontinuous Galerkin Spectral Element Method for Electrostatic Particle-in-Cell Simulations

TL;DR

This paper introduces a p-adaptive hybridizable discontinuous Galerkin spectral element method for electrostatic plasma simulations, improving efficiency by local polynomial degree refinement in particle-in-cell schemes.

Contribution

It develops a novel p-adaptive HDG-SEM approach integrated into PICLas, enabling efficient, high-order electrostatic plasma simulations with local refinement capabilities.

Findings

Reduces degrees of freedom compared to uniform high-order methods

Successfully validated with benchmark test cases

Demonstrates efficiency in complex plasma modeling

Abstract

This paper presents a p-adaptive high-order hybridizable discontinuous Galerkin spectral element method (HDG-SEM) for solving the Poisson equation in electrostatic plasma simulations using particle-in-cell (PIC) schemes. This approach enables element-local refinement of the polynomial degree, concentrating computational effort specifically in regions with strong gradients. Thus, the method significantly reduces the global number of degrees of freedom compared to uniform high-order methods. The proposed method is implemented in the open-source framework PICLas and validated through a series of benchmark test cases, including a dielectric sphere and a one-dimensional plasma sheath. Finally, a two-dimensional axisymmetric simulation of an ion optic demonstrates the method's capability to efficiently model complex plasma phenomena but also highlights current limitations.