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.

Spectral methods on a triangle and W-systems

Stable Spectral Methods for Time-Dependent Problems and the Preservation of Structure

A Triangular Spectral Method for the Stokes Eigenvalue Problem by the Stream Function Formulation

Spectral quadrangulation with boundary conformation

A framework for stable spectral methods in $d$-dimensional unit balls

A triangular spectral element method; applications to the incompressible Navier-Stokes equations

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

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

Mathematical foundations of spectral methods for time-dependent PDEs

Spectral Quadrangulation with Orientation and Alignment Control.

Spectral Methods for Solution of Differential and Functional Equations

Anti-symmetric and Positivity Preserving Formulation of a Spectral Method for Vlasov-Poisson Equations

Spectral Galerkin Approximation and Error Analysis Based on a Mixed Scheme for Fourth‐Order Problems in Complex Regions

The numerical solution of fractional integral equations via orthogonal polynomials in fractional powers

Spectral Analysis of High Order Continuous FEM for Hyperbolic PDEs on Triangular Meshes: Influence of Approximation, Stabilization, and Time-Stepping

Numerical evolutions of fields on the 2-sphere using a spectral method based on spin-weighted spherical harmonics

Multiwavelets over a Triangular Domain

Spectral stability and perturbation results for kernel differentiation matrices on the sphere

A sparse spectral method for fractional differential equations in one-spatial dimension

Efficient and Accurate Spectral Method for Solving Fractional Differential Equations on the Half Line Using Orthogonal Generalized Rational Jacobi Functions

Weyl-Type Theorems of Upper Triangular Relation Matrices