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.

翻译：暂无翻译