A Curved Boundary Treatment for Discontinuous Galerkin Method Applied to Euler Equations on Triangular and Tetrahedral Grids.
Junhui Yin,Li Xu,Peng Xie,Lan Zhu,Shucheng Huang,Hangxin Liu,Zhonghai Yang,Bin Li
DOI: https://doi.org/10.1016/j.cpc.2020.107549
IF: 4.717
2021-01-01
Computer Physics Communications
Abstract:For Euler equations defined in a domain Ω with boundary conditions specified on a curved boundary Γ, we consider a discontinuous Galerkin (DG) method with high order polynomial basis functions on a geometry fitting triangular and tetrahedral grids. We derive a modification to the DG scheme defined on a boundary triangle and tetrahedron based on the general curvilinear element DG method. In the modified DG scheme, even though integration along the curved boundary is still necessary, integrals over any curvilinear element are avoided. To further simplify the calculation of integration along the curved boundary, the modified DG scheme is constructed by the relationship between normal vector of geometric boundary and surface Jacobian. In the new modified DG scheme, not only integrals over any curvilinear element are avoided, but also the line and face integrations along the curved boundary are not required. Numerical tests suggest that such a new modified DG scheme is as accurate as the basic modified DG scheme and full curvilinear DG scheme and more efficient than latter, and is stable on reasonably coarse, as well as more refined, grids. It is shown that this new modified DG scheme can be extremely useful to solve the problems defined in a domain with boundary conditions specified on a curved boundary.