An Entropy Stable Discontinuous Galerkin Spectral Element Method on Heterogeneous Grids
Please login to view abstract download link
The discontinuous Galerkin spectral element method (DGSEM) has proven to be a highly efficient numerical scheme for the approximation of smooth multiscale problems, e.g., turbulence in the present case. Key aspects for the high efficiency are its excellent parallelization properties and the use of a tensor-product ansatz, where the multidimensional space operator is decomposed into a sequence of one-dimensional operations. Efficient and robust application-oriented large eddy simulations (LES) of turbulence were enabled by the use of entropy-stable DGSEM in a flux differencing formulation on quadrilateral and hexahedral elements. However, for the discretization of complex geometries, more flexible element types, e.g., simplexes, are in general more favorable. In a recent work, Montoya and Zingg presented an efficient entropy stable flux differencing DGSEM formulation on triangular and tetrahedral elements. To benefit from the computational efficiency of the tensor-product ansatz, the authors used a collapsed coordinate transformation on Legendre–Gauss-Lobatto and Legendre-Gauss- Radau nodes. In practical applications such as the flow around an airfoil in the majority of the computational domain, the use of, e.g., hexahedral elements can be favorable and only in a small subset, e.g., simplexes are required. For this, the approach in is extended to triangular prism, pyramid and tetrahedral elements on Legendre–Gauss nodes to enable an efficient entropy stable DGSEM on heterogeneous grids. In the present talk, an entropy stable DGSE scheme on heterogeneous meshes is presented. A particular focus is placed on the collapsed coordinates transformation, the interelement flux computation between different elements types, the modal time-update of the non-hexahedral elements, and the load balancing strategy. The free-stream preservation and the convergence properties are demonstrated based on numerical experiments and the excellent scaling properties of the DGSEM are shown. This talk concludes with several application-oriented numerical problems.
