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.

Efficient Spectral and Spectral Element Methods for Eigenvalue Problems of Schrödinger Equations with an Inverse Square Potential

Efficient Chebyshev Spectral Method for Solving Linear Elliptic Pdes Using Quasi-Inverse Technique

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

Exponential Collocation Methods Based on Continuous Finite Element Approximations for Efficiently Solving the Cubic Schrödinger Equation

An efficient spectral method for the fractional Schrödinger equation on the real line

Spectral Collocation Technique for Absorbing Boundary Conditions with Increasingly High Order Approximation

Efficient Legendre-Laguerre spectral element methods for problems on unbounded domains with diagonalization technique

The ultraspherical spectral element method

A spectral Galerkin exponential Euler time-stepping scheme for parabolic SPDEs on two-dimensional domains with a C2-boundary

Efficient Spectral Methods for Quasi-Equilibrium Closure Approximations of Symmetric Problems on Unit Circle and Sphere

A Legendre spectral method based on a hybrid format and its error estimation for fourth-order eigenvalue problems

Isogeometric analysis of the Laplace eigenvalue problem on circular sectors: Regularity properties, graded meshes & variational crimes

The Spectral Collocation Method for Efficiently Solving PDEs with Fractional Laplacian.

Highly efficient Gauss's law-preserving spectral algorithms for Maxwell's double-curl source and eigenvalue problems based on eigen-decomposition

An Efficient and Accurate Numerical Method for the Spectral Fractional Laplacian Equation

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

The Spectral Element Method for the Steklov Eigenvalue Problem

Possibility of magnetic field induced Composite Fermi sea in neutral Graphene

A Highly Efficient Spectral-Galerkin Method Based on Tensor Product for Fourth-Order Steklov Equation with Boundary Eigenvalue

An Adaptive Inverse Iteration for Maxwell Eigenvalue Problem Based on Edge Elements.

A new spectral element method for numerical solution of partial differential equations on annular-type regions