Efficient Entropy-Stable Discontinuous Spectral-Element Methods Using Tensor-Product Summation-by-Parts Operators on Triangles and Tetrahedra

Abstract

We present a new class of efficient and robust discontinuous spectral-element methods of arbitrary order for nonlinear hyperbolic systems of conservation laws on curved triangular and tetrahedral unstructured grids. Such discretizations employ a recently introduced family of sparse tensor-product summation-by-parts (SBP) operators in collapsed coordinates within an entropy-stable modal formulation. The proposed algorithms exploit the structure of such SBP operators alongside that of the Proriol-Koornwinder-Dubiner polynomial basis, and a weight-adjusted approximation is used to efficiently invert the local mass matrix for curvilinear elements. Using such techniques, the number of required entropy-conservative two-point flux evaluations between pairs of quadrature nodes is significantly reduced relative to existing entropy-stable formulations using (non-tensor-product) multidimensional SBP operators, particularly for high polynomial degrees, with an improvement in time complexity from $\mathcal{O}(p^{2d})$ to $\mathcal{O}(p^{d+1})$, where $p$ is the polynomial degree of the approximation and $d$ is the number of spatial dimensions. In numerical experiments involving smooth solutions to the compressible Euler equations, the proposed tensor-product schemes demonstrate similar levels of accuracy for a given mesh and polynomial degree to those using multidimensional SBP operators based on symmetric quadrature rules. Furthermore, both operator families are shown to give rise to entropy-stable methods which exhibit excellent robustness for under-resolved problems. Such results suggest that the algorithmic advantages resulting from the use of tensor-product operators are obtained without compromising accuracy or robustness, enabling the efficient extension of the benefits of entropy stability to higher polynomial degrees than previously considered for triangular and tetrahedral elements.

38 pages, 9 figures

References in corpus (13)

• Review of Summation-by-parts schemes for initial-boundary-value problems
• Split Form Nodal Discontinuous Galerkin Schemes with Summation-By-Parts Property for the Compressible Euler Equations
• Nodal Discontinuous Galerkin Methods on Graphics Processors
• Kinetic energy preserving and entropy stable finite volume schemes for compressible Euler and Navier-Stokes equations
• On discretely entropy conservative and entropy stable discontinuous Galerkin methods
• Comparison of some Entropy Conservative Numerical Fluxes for the Euler Equations
• Subcell limiting strategies for discontinuous Galerkin spectral element methods
• Preventing pressure oscillations does not fix local linear stability issues of entropy-based split-form high-order schemes
• Analysis and entropy stability of the line-based discontinuous Galerkin method
• A positivity preserving strategy for entropy stable discontinuous Galerkin discretizations of the compressible Euler and Navier-Stokes equations
• Efficient implementation of modern entropy stable and kinetic energy preserving discontinuous Galerkin methods for conservation laws
• A unifying algebraic framework for discontinuous Galerkin and flux reconstruction methods based on the summation-by-parts property
• Efficient Tensor-Product Spectral-Element Operators with the Summation-by-Parts Property on Curved Triangles and Tetrahedra

Cited by in corpus (2)

• Generalized upwind summation-by-parts operators and their application to nodal discontinuous Galerkin methods
• Finite-difference compatible entropy-conserving schemes for the compressible Euler equations