Boundary optimized diagonal-norm SBP operators

Abstract By using non-equispaced grid points near boundaries, we derive boundary optimized first derivative finite difference operators, of orders up to twelve. The boundary closures are based on a diagonal-norm summation-by-parts (SBP) framework, thereby guaranteeing linear stability on piecewise curvilinear multi-block grids. The new operators lead to significantly more efficient numerical approximations, compared with traditional SBP operators on equidistant grids. We also show that the non-uniform grids make it possible to derive operators with fewer one-sided boundary stencils than their traditional counterparts. Numerical experiments with the 2D compressible Euler equations on a curvilinear multi-block grid demonstrate the accuracy and stability properties of the new operators.

[1]  Jan Nordström,et al.  On the order of Accuracy of Finite Difference Operators on Diagonal Norm Based Summation-by-Parts Form , 2018, SIAM J. Numer. Anal..

[2]  N. N. Kalitkin,et al.  Optimal first- to sixth-order accurate Runge-Kutta schemes , 2008 .

[3]  Steven H. Frankel,et al.  Entropy Stable Spectral Collocation Schemes for the Navier-Stokes Equations: Discontinuous Interfaces , 2014, SIAM J. Sci. Comput..

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

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

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

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

[8]  Michael A. Saunders,et al.  SNOPT: An SQP Algorithm for Large-Scale Constrained Optimization , 2005, SIAM Rev..

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

[10]  Magnus Svärd,et al.  On the convergence rates of energy-stable finite-difference schemes , 2019, J. Comput. Phys..

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

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

[13]  Nathan Albin,et al.  Existence of SBP operators with diagonal norm , 2014 .

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

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