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.

Constructing stable, high-order finite-difference operators on point clouds over complex geometries

Conservative and Stable Degree Preserving SBP Operators for Non-Conforming Meshes

A New Spectral Difference Method Using Hierarchical Polynomial Bases for Hyperbolic Conservation Laws

A dual-pairing summation-by-parts finite difference framework for nonlinear conservation laws

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

A High Order Stable Conservative Method for Solving Hyperbolic Conservation Laws on Arbitrarily Distributed Point Clouds.

An Entropy Stable High-Order Discontinuous Galerkin Method on Cut Meshes

Efficient discretization of the Laplacian on complex geometries

An optimization-based construction procedure for function space based summation-by-parts operators on arbitrary grids

On the robustness of high-order upwind summation-by-parts methods for nonlinear conservation laws

Accurate High-Order Tensor-Product Generalized Summation-By-Parts Discretizations of Hyperbolic Conservation Laws: General Curved Domains and Functional Superconvergence

Efficient Tensor-Product Spectral-Element Operators with the Summation-by-Parts Property on Curved Triangles and Tetrahedra

Positivity-preserving high-order schemes for conservation laws on arbitrarily distributed point clouds with a simple weno limiter

Compact finite-difference scheme for some Sobolev type equations with Dirichlet boundary conditions

A robust first order meshfree method for time-dependent nonlinear conservation laws

An optimization-based approach for high-order accurate discretization of conservation laws with discontinuous solutions

Generalized upwind summation-by-parts operators and their application to nodal discontinuous Galerkin methods

Multi-dimensional summation-by-parts operators for general function spaces: Theory and construction

Mortar-based entropy-stable discontinuous Galerkin methods on non-conforming quadrilateral and hexahedral meshes

A discontinuous Galerkin discretization of elliptic problems with improved convergence properties using summation by parts operators

Relative Efficiency of Finite-Difference and Discontinuous Spectral-Element Summation-by-Parts Methods on Distorted Meshes