Generalized Summation-by-Parts Operators for the Second Derivative

The generalization of summation-by-parts operators for the first derivative of Del Rey Fernandez, Boom, and Zingg [J. Comput. Phys., 266 (2014), pp. 214--239] is extended to approximations of second derivatives with a constant or variable coefficient. This enables the construction of second-derivative operators with one or more of the following characteristics: (i) nonrepeating interior point operators, (ii) nonuniform nodal distributions, and (iii) exclusion of one or both boundary nodes. Definitions are proposed that give rise to generalized summation-by-parts operators that result in consistent, conservative, and stable discretizations of partial differential equations with or without mixed derivatives. It is proven that approximations to the second derivative with a variable coefficient can be constructed using the constituent matrices of the constant-coefficient operator. Moreover, for operators with a repeating interior point operator, a decomposition is proposed that makes the application of such o...

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

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

[3]  Jan Nordström,et al.  High-order finite difference methods, multidimensional linear problems, and curvilinear coordinates , 2001 .

[4]  Ramji Kamakoti,et al.  High-Order Narrow Stencil Finite-Difference Approximations of Second-Order Derivatives Involving Variable Coefficients , 2009, SIAM J. Sci. Comput..

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

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

[7]  A. Chertock,et al.  Strict Stability of High-Order Compact Implicit Finite-Difference Schemes: The Role of Boundary Conditions for Hyperbolic PDEs, II , 2000 .

[8]  I. Babuska,et al.  A DiscontinuoushpFinite Element Method for Diffusion Problems , 1998 .

[9]  Bradley K. Alpert,et al.  Hybrid Gauss-Trapezoidal Quadrature Rules , 1999, SIAM J. Sci. Comput..

[10]  David C. Del Rey Fernández,et al.  High-Order Compact-Stencil Summation-By-Parts Operators for the Compressible Navier-Stokes Equations , 2013 .

[11]  Jie Shen,et al.  Spectral Methods: Algorithms, Analysis and Applications , 2011 .

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

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

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

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

[16]  Ken Mattsson,et al.  Optimal diagonal-norm SBP operators , 2014, J. Comput. Phys..

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

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

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

[20]  Magnus Svärd,et al.  Review of summation-by-parts schemes for initial-boundary-value problems , 2013, J. Comput. Phys..

[21]  Jan Nordström,et al.  Revisiting and Extending Interface Penalties for Multi-domain Summation-by-Parts Operators , 2010, J. Sci. Comput..

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

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

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

[25]  Magnus Svärd,et al.  On the order of accuracy for difference approximations of initial-boundary value problems , 2006, J. Comput. Phys..

[26]  David C. Del Rey Fernández,et al.  Review of summation-by-parts operators with simultaneous approximation terms for the numerical solution of partial differential equations , 2014 .

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

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

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

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

[31]  J. Oden,et al.  A discontinuous hp finite element method for convection—diffusion problems , 1999 .

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

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