A generalized framework for nodal first derivative summation-by-parts operators

A generalized framework is presented that extends the classical theory of finite-difference summation-by-parts (SBP) operators to include a wide range of operators, where the main extensions are (i) non-repeating interior point operators, (ii) nonuniform nodal distribution in the computational domain, (iii) operators that do not include one or both boundary nodes. Necessary and sufficient conditions are proven for the existence of nodal approximations to the first derivative with the SBP property. It is proven that the positive-definite norm matrix of each SBP operator must be associated with a quadrature rule; moreover, given a quadrature rule there exists a corresponding SBP operator, where for diagonal-norm SBP operators the weights of the quadrature rule must be positive. The generalized framework gives a straightforward means of posing many known approximations to the first derivative as SBP operators; several are surveyed, such as discontinuous Galerkin discretizations based on the Legendre-Gauss quadrature points, and shown to be SBP operators. Moreover, the new framework provides a method for constructing SBP operators by starting from quadrature rules; this is illustrated by constructing novel SBP operators from known quadrature rules. To demonstrate the utility of the generalization, the Legendre-Gauss and Legendre-Gauss-Radau quadrature points are used to construct SBP operators that do not include one or both boundary nodes.

[1]  Nico M. Temme,et al.  Numerical methods for special functions , 2007 .

[2]  Edward N. Tinoco,et al.  Summary of Data from the Fifth AIAA CFD Drag Prediction Workshop , 2013 .

[3]  Erik Schnetter,et al.  Optimized High-Order Derivative and Dissipation Operators Satisfying Summation by Parts, and Applications in Three-dimensional Multi-block Evolutions , 2005, J. Sci. Comput..

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

[5]  Alfio Quarteroni,et al.  Numerical Mathematics (Texts in Applied Mathematics) , 2006 .

[6]  K. Mattsson.  Accurate and stable finite volume operators for unstructured flow solvers , 2022 .

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

[8]  Kai Hormann,et al.  On the Lebesgue constant of barycentric rational interpolation at equidistant nodes , 2012, Numerische Mathematik.

[9]  Francis Jack Smith,et al.  Error Estimation in the Clenshaw-Curtis Quadrature Formula , 1968, Comput. J..

[10]  Jean-Paul Berrut,et al.  Convergence rates of derivatives of a family of barycentric rational interpolants , 2011 .

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

[12]  Ken Mattsson,et al.  A solution to the stability issues with block norm summation by parts operators , 2013, J. Comput. Phys..

[13]  Jean-Paul Berrut,et al.  Linear barycentric rational quadrature , 2012 .

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

[15]  Nail K. Yamaleev,et al.  A systematic methodology for constructing high-order energy stable WENO schemes , 2009, J. Comput. Phys..

[16]  Michael T. Heath,et al.  Energy stable numerical methods for hyperbolic partial differential equations using overlapping domain decomposition , 2012, J. Comput. Phys..

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

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

[19]  Jan S. Hesthaven,et al.  A Stable Penalty Method for the Compressible Navier-Stokes Equations: I. Open Boundary Conditions , 1996, SIAM J. Sci. Comput..

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

[21]  Ken Mattsson,et al.  Boundary Procedures for Summation-by-Parts Operators , 2003, J. Sci. Comput..

[22]  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..

[23]  Lloyd N. Trefethen,et al.  Is Gauss Quadrature Better than Clenshaw-Curtis? , 2008, SIAM Rev..

[24]  Jan Nordström,et al.  Boundary and Interface Conditions for High-Order Finite-Difference Methods Applied to the Euler and Navier-Stokes Equations , 1999 .

[25]  Jan S. Hesthaven,et al.  A Stable Penalty Method for the Compressible Navier-Stokes Equations: III. Multidimensional Domain Decomposition Schemes , 1998, SIAM J. Sci. Comput..

[26]  Jan S. Hesthaven,et al.  A Stable Penalty Method for the Compressible Navier-Stokes Equations: II. One-Dimensional Domain Decomposition Schemes , 1997, SIAM J. Sci. Comput..

[27]  H. Kreiss,et al.  Comparison of accurate methods for the integration of hyperbolic equations , 1972 .

[28]  Jing Gong,et al.  Interface procedures for finite difference approximations of the advection-diffusion equation , 2011, J. Comput. Appl. Math..

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

[30]  Magnus Svärd,et al.  Stable and Accurate Artificial Dissipation , 2004, J. Sci. Comput..

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

[32]  Nail K. Yamaleev,et al.  Boundary closures for fourth-order energy stable weighted essentially non-oscillatory finite-difference schemes , 2011, J. Comput. Phys..

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

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

[35]  Kai Hormann,et al.  Barycentric rational interpolation with no poles and high rates of approximation , 2007, Numerische Mathematik.

[36]  Magnus Svärd,et al.  Stability of finite volume approximations for the Laplacian operator on quadrilateral and triangular grids , 2004 .

[37]  Magnus Svärd,et al.  A stable high-order finite difference scheme for the compressible Navier-Stokes equations, far-field boundary conditions , 2007, J. Comput. Phys..

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

[39]  H. Kreiss,et al.  Initial-Boundary Value Problems and the Navier-Stokes Equations , 2004 .

[40]  Jason E. Hicken,et al.  The Role of Dual Consistency in Functional Accuracy: Error Estimation and Superconvergence , 2011 .

[41]  Magnus Svärd,et al.  Stable and accurate schemes for the compressible Navier-Stokes equations , 2008, J. Comput. Phys..

[42]  B. Gustafsson.  The convergence rate for difference approximations to mixed initial boundary value problems , 1975 .

[43]  Georges Klein,et al.  Applications of linear barycentric rational interpolation , 2012 .

[44]  H. Kreiss,et al.  Time-Dependent Problems and Difference Methods , 1996 .

[45]  Burton Wendroff,et al.  The Relative Efficiency of Finite Difference and Finite Element Methods. I: Hyperbolic Problems and Splines , 1974 .

[46]  B. Gustafsson.  High Order Difference Methods for Time Dependent PDE , 2008 .

[47]  Jason E. Hicken,et al.  Superconvergent Functional Estimates from Summation-By-Parts Finite-Difference Discretizations , 2011, SIAM J. Sci. Comput..

[48]  Magnus Svärd,et al.  A stable high-order finite difference scheme for the compressible Navier-Stokes equations: No-slip wall boundary conditions , 2008, J. Comput. Phys..

[49]  Nail K. Yamaleev,et al.  Third-Order Energy Stable WENO Scheme , 2008 .

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

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

[52]  Jan Nordström,et al.  Finite volume approximations and strict stability for hyperbolic problems , 2001 .

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