Stability analysis of discontinuous Galerkin with a high order embedded boundary treatment for linear hyperbolic equations

TL;DR

This paper investigates the stability of high order discontinuous Galerkin methods with embedded boundary corrections for linear hyperbolic equations, revealing stability limitations with explicit time integration and solutions with implicit schemes.

Contribution

It provides a detailed stability analysis of high order embedded boundary treatments using shifted boundary polynomial correction in DG methods, highlighting stability constraints and implicit solution benefits.

Findings

Explicit time integration limits stability for high order methods.

Implicit time integration achieves unconditional stability.

Stability constraints increase with polynomial degree.

Abstract

Embedded, or immersed, approaches have the goal of reducing to the minimum the computational costs associated with the generation of body-fitted meshes by only employing fixed, possibly Cartesian, meshes over which complex boundaries can move freely. However, this boundary treatment introduces a geometrical error of the order of the mesh size that, if not treated properly, can spoil the global accuracy of a high order discretization, herein based on discontinuous Galerkin. The shifted boundary polynomial correction was proposed as a simplified version of the shifted boundary method, which is an embedded boundary treatment based on Taylor expansions to deal with unfitted boundaries. It is used to accordingly correct the boundary conditions imposed on a non-meshed boundary to compensate the aforementioned geometrical error, and reach high order accuracy. In this paper, the stability analysis of discontinuous Galerkin methods coupled with the shifted boundary polynomial correction is conducted in depth for the linear advection equation, by visualizing the eigenvalue spectrum of the high order discretized operators. The analysis considers a simplified one-dimensional setting by varying the degree of the polynomials and the distance between the real boundary and the closest mesh interface. The main result of the analysis shows that the considered high order embedded boundary treatment introduces a limitation to the stability region of high order discontinuous Galerkin methods with explicit time integration, which becomes more and more important when using higher order methods. The implicit time integration is also studied, showing that the implicit treatment of the boundary condition allows one to overcome such limitation and achieve an unconditionally stable high order embedded boundary treatment.