Score: 0

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

Published: October 27, 2025 | arXiv ID: 2510.23231v1

By: Mirco Ciallella

Potential Business Impact:

Fixes computer simulations of moving shapes.

Business Areas:
Embedded Systems Hardware, Science and Engineering, Software

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.

Page Count
24 pages

Category
Mathematics:
Numerical Analysis (Math)