A novel strategy for deriving high-order stable boundary closures based on global conservation, I: Basic formulas

Abstract It is well known that for initial boundary value problems one may need to decrease the order of boundary closures to ensure strict stability (also called time stability) of finite difference schemes. For high-order schemes, this often reduces the global convergence rate significantly. To cure this issue, we derive high-order accurate boundary closures that satisfy a global conservation property. We show that the boundary closures cannot be both high-order accurate and conservative in the l 2 -inner product, if the grid is equidistant. Therefore, we derive particular sets of nonuniform solution points and quadrature weights near the domain boundaries. The new schemes are high-order accurate on the modified grids and conservative in the derived quadrature rules. Although we have no theoretical stability proof, numerical eigenvalue analysis suggests that the new schemes are strictly stable up to eleventh-order global accuracy for scalar linear advection equations. Numerical experiments with one- and two-dimensional linear equations corroborate the design orders of the schemes.

[1]  Xiaolin Zhong,et al.  Derivation of high-order compact finite difference schemes for non-uniform grid using polynomial interpolation , 2005 .

[2]  Matteo Parsani,et al.  Entropy Stable Staggered Grid Discontinuous Spectral Collocation Methods of any Order for the Compressible Navier-Stokes Equations , 2016, SIAM J. Sci. Comput..

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

[4]  Jan Nordström,et al.  The SBP-SAT technique for initial value problems , 2014, J. Comput. Phys..

[5]  Ken Mattsson,et al.  Compatible diagonal-norm staggered and upwind SBP operators , 2018, J. Comput. Phys..

[6]  Dan Kosloff,et al.  A modified Chebyshev pseudospectral method with an O(N –1 ) time step restriction , 1993 .

[7]  Changhoon Lee,et al.  A new compact spectral scheme for turbulence simulations , 2002 .

[8]  Miguel R. Visbal,et al.  On the use of higher-order finite-difference schemes on curvilinear and deforming meshes , 2002 .

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

[10]  Taku Nonomura,et al.  Geometric interpretations and spatial symmetry property of metrics in the conservative form for high-order finite-difference schemes on moving and deforming grids , 2014, J. Comput. Phys..

[11]  Miguel Hermanns,et al.  Stable high‐order finite‐difference methods based on non‐uniform grid point distributions , 2008 .

[12]  Xiaogang Deng,et al.  Developing high-order weighted compact nonlinear schemes , 2000 .

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

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

[15]  Thomas Hagstrom,et al.  Grid stabilization of high-order one-sided differencing I: First-order hyperbolic systems , 2007, J. Comput. Phys..

[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]  Mahidhar Tatineni,et al.  High-order non-uniform grid schemes for numerical simulation of hypersonic boundary-layer stability and transition , 2003 .

[19]  P. Thomas,et al.  Geometric Conservation Law and Its Application to Flow Computations on Moving Grids , 1979 .

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

[21]  Alina Chertock,et al.  Strict Stability of High-Order Compact Implicit Finite-Difference Schemes , 2000 .

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

[23]  Huayong Liu,et al.  Further studies on Geometric Conservation Law and applications to high-order finite difference schemes with stationary grids , 2013, J. Comput. Phys..

[24]  Bertil Gustafsson,et al.  The convergence rate for difference approximations to general mixed initial boundary value problems , 1981 .

[25]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[26]  Thomas Hagstrom,et al.  Grid stabilization of high-order one-sided differencing II: Second-order wave equations , 2012, J. Comput. Phys..

[27]  Hiroshi Maekawa,et al.  Compact High-Order Accurate Nonlinear Schemes , 1997 .

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

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

[30]  Rémi Abgrall,et al.  High‐order CFD methods: current status and perspective , 2013 .

[31]  Taku Nonomura,et al.  Freestream and vortex preservation properties of high-order WENO and WCNS on curvilinear grids , 2010 .

[32]  Huayong Liu,et al.  Geometric conservation law and applications to high-order finite difference schemes with stationary grids , 2011, J. Comput. Phys..

[33]  Soogab Lee,et al.  Grid-optimized dispersion-relation-preserving schemes on general geometries for computational aeroacoustics , 2001 .

[34]  Chi-Wang Shu,et al.  High order finite difference and finite volume WENO schemes and discontinuous Galerkin methods for CFD , 2001 .

[35]  Xiaogang Deng,et al.  A family of hybrid cell-edge and cell-node dissipative compact schemes satisfying geometric conservation law , 2015 .

[36]  John H. Kolias,et al.  A CONSERVATIVE STAGGERED-GRID CHEBYSHEV MULTIDOMAIN METHOD FOR COMPRESSIBLE FLOWS , 1995 .