Multidimensional Summation-by-Parts Operators: General Theory and Application to Simplex Elements

Summation-by-parts (SBP) finite-difference discretizations share many attractive properties with Galerkin finite-element methods (FEMs), including time stability and superconvergent functionals; however, unlike FEMs, SBP operators are not completely determined by a basis, so the potential exists to tailor SBP operators to meet different objectives. To date, application of high-order SBP discretizations to multiple dimensions has been limited to tensor product domains. This paper presents a definition for multi-dimensional SBP finite-difference operators that is a natural extension of one-dimensional SBP operators. Theoretical implications of the definition are investigated for the special case of a diagonal norm (mass) matrix. In particular, a diagonal-norm SBP operator exists on a given domain if and only if there is a cubature rule with positive weights on that domain and the polynomial-basis matrix has full rank when evaluated at the cubature nodes. Appropriate simultaneous-approximation terms are developed to impose boundary conditions weakly, and the resulting discretizations are shown to be time stable. Concrete examples of multi-dimensional SBP operators are constructed for the triangle and tetrahedron; similarities and differences with spectral-element and spectral-difference methods are discussed. An assembly process is described that builds diagonal-norm SBP operators on a global domain from element-level operators. Numerical results of linear advection on a doubly periodic domain demonstrate the accuracy and time stability of the simplex operators.

[1]  T. Koornwinder.  Two-Variable Analogues of the Classical Orthogonal Polynomials , 1975 .

[2]  P. Moin,et al.  Fully Conservative Higher Order Finite Difference Schemes for Incompressible Flow , 1998 .

[3]  J. Nordström,et al.  Summation by Parts Operators for Finite Difference Approximations of Second-Derivatives with Variable Coefficients , 2004, Journal of Scientific Computing.

[4]  Long Chen.  FINITE VOLUME METHODS , 2011 .

[5]  Robert H. Halstead,et al.  Matrix Computations , 2011, Encyclopedia of Parallel Computing.

[6]  R. McLachlan,et al.  SKEW-ADJOINT FINITE DIFFERENCE METHODS ON NONUNIFORM GRIDS , 2022 .

[7]  Moshe Dubiner.  Spectral methods on triangles and other domains , 1991 .

[8]  H. Kreiss,et al.  Finite Element and Finite Difference Methods for Hyperbolic Partial Differential Equations , 1974 .

[9]  Jing Gong,et al.  A stable and conservative high order multi-block method for the compressible Navier-Stokes equations , 2009, J. Comput. Phys..

[10]  Jason E. Hicken,et al.  Summation-by-parts operators and high-order quadrature , 2011, J. Comput. Appl. Math..

[11]  H. C. Yee.  DESIGNING ADAPTIVE LOW-DISSIPATIVE HIGH ORDER SCHEMES FOR LONG-TIME INTEGRATIONS , .

[12]  Jason E. Hicken,et al.  Dual consistency and functional accuracy: a finite-difference perspective , 2014, J. Comput. Phys..

[13]  B. Strand.  Summation by parts for finite difference approximations for d/dx , 1994 .

[14]  Qiqi Wang,et al.  A Conservative Mesh-Free Scheme and Generalized Framework for Conservation Laws , 2012, SIAM J. Sci. Comput..

[15]  Len Bos,et al.  On certain configurations of points in R n which are unisolvent for polynomial interpolation , 1991 .

[16]  H. C. Yee,et al.  Designing Adaptive Low-Dissipative High Order Schemes for Long-Time Integrations. Chapter 1 , 2001 .

[17]  Hui,et al.  A SET OF SYMMETRIC QUADRATURE RULES ON TRIANGLES AND TETRAHEDRA , 2009 .

[18]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

[19]  Magnus Svärd,et al.  On Coordinate Transformations for Summation-by-Parts Operators , 2004, J. Sci. Comput..

[20]  Marcel Vinokur,et al.  Spectral difference method for unstructured grids I: Basic formulation , 2006, J. Comput. Phys..

[21]  J. Hesthaven,et al.  Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications , 2007 .

[22]  Jan Nordström,et al.  Finite volume methods, unstructured meshes and strict stability for hyperbolic problems , 2003 .

[23]  Jan Nordström,et al.  Conservative Finite Difference Formulations, Variable Coefficients, Energy Estimates and Artificial Dissipation , 2006, J. Sci. Comput..

[24]  Marcel Vinokur,et al.  Exact Integrations of Polynomials and Symmetric Quadrature Formulas over Arbitrary Polyhedral Grids , 1998 .

[25]  Ken Mattsson,et al.  Summation by Parts Operators for Finite Difference Approximations of Second-Derivatives with Variable Coefficients , 2012, J. Sci. Comput..

[26]  D. Gottlieb,et al.  Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: methodology and application to high-order compact schemes , 1994 .

[27]  Gregor Gassner,et al.  A Skew-Symmetric Discontinuous Galerkin Spectral Element Discretization and Its Relation to SBP-SAT Finite Difference Methods , 2013, SIAM J. Sci. Comput..

[28]  Antony Jameson,et al.  A Conservative Meshless Scheme: General Order Formulation and Application to Euler Equations , 2011 .

[29]  David C. Del Rey Fernández,et al.  A generalized framework for nodal first derivative summation-by-parts operators , 2014, J. Comput. Phys..

[30]  Wim A. Mulder,et al.  HIGHER-ORDER MASS-LUMPED FINITE ELEMENTS FOR THE WAVE EQUATION , 2001 .

[31]  Jason E. Hicken.  Output error estimation for summation-by-parts finite-difference schemes , 2012, J. Comput. Phys..

[32]  Freddie D. Witherden,et al.  On the identification of symmetric quadrature rules for finite element methods , 2014, Comput. Math. Appl..

[33]  Stephen P. Boyd,et al.  Convex Optimization , 2004, Algorithms and Theory of Computation Handbook.

[34]  H. C. Yee,et al.  Entropy Splitting and Numerical Dissipation , 2000 .

[35]  Francis X. Giraldo,et al.  A diagonal-mass-matrix triangular-spectral-element method based on cubature points , 2007 .

[36]  D. César,et al.  Generalized Summation-by-parts Operators for First and Second Derivatives , 2015 .

[37]  David Gottlieb,et al.  Spectral Methods on Arbitrary Grids , 1995 .

[38]  Magnus Svärd,et al.  Steady-State Computations Using Summation-by-Parts Operators , 2005, J. Sci. Comput..

[39]  Jean E. Roberts,et al.  Higher Order Triangular Finite Elements with Mass Lumping for the Wave Equation , 2000, SIAM J. Numer. Anal..

[40]  D. Gottlieb,et al.  A Stable and Conservative Interface Treatment of Arbitrary Spatial Accuracy , 1999 .

[41]  Jason E. Hicken,et al.  Parallel Newton-Krylov Solver for the Euler Equations Discretized Using Simultaneous-Approximation Terms , 2008 .

[42]  LiuYen,et al.  Spectral difference method for unstructured grids I , 2006 .

[43]  D. Gottlieb,et al.  A new method of imposing boundary conditions in pseudospectral approximations of hyperbolic equations , 1988 .