An improved projection method

Abstract Strictly stable high-order accurate finite difference approximations are derived, for linear initial boundary value problems. The boundary closures are based on the diagonal-norm summation-by-parts framework and the boundary conditions are imposed using an improved projection method. The improved projection method removes some of the numerical issues with the previously derived (here referred to as traditional ) projection method. The finite difference approximations lead to fully explicit ODE systems. The accuracy and stability properties are demonstrated for linear hyperbolic and parabolic problems in 1D and 2D. In particular we study the problem where initial and boundary data are inconsistent, that causes the traditional projection method to fail. Comparisons are made also with weak (penalty) enforcement of boundary conditions.

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

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

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

[4]  Ken Mattsson,et al.  Diagonal-norm upwind SBP operators , 2017, J. Comput. Phys..

[5]  Tomas Edvinsson,et al.  High-fidelity numerical solution of the time-dependent Dirac equation , 2014, J. Comput. Phys..

[6]  Gianluca Iaccarino,et al.  Stable and accurate wave-propagation in discontinuous media , 2008, J. Comput. Phys..

[7]  Matteo Parsani,et al.  Entropy stable discontinuous interfaces coupling for the three-dimensional compressible Navier-Stokes equations , 2015, J. Comput. Phys..

[8]  Magnus Svärd,et al.  High-order accurate computations for unsteady aerodynamics , 2007 .

[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.  Entropy-Stable Schemes for the Euler Equations with Far-Field and Wall Boundary Conditions , 2014, J. Sci. Comput..

[11]  Jan Nordström,et al.  Superconvergent functional output for time-dependent problems using finite differences on summation-by-parts form , 2012, J. Comput. Phys..

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

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

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

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

[16]  David W. Zingg,et al.  Higher-order spatial discretization for turbulent aerodynamic computations , 2001 .

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

[18]  S. Abarbanel,et al.  Asymptotically Stable Fourth-Order Accurate Schemes for the Diffusion Equation on Complex Shapes , 1997 .

[19]  P. Olsson.  Summation by parts, projections, and stability. II , 1995 .

[20]  Ken Mattsson,et al.  Diagonal-norm summation by parts operators for finite difference approximations of third and fourth derivatives , 2014, J. Comput. Phys..

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

[22]  Eli Turkel,et al.  A fourth-order accurate finite-difference scheme for the computation of elastic waves , 1986 .

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

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

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

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

[27]  Anna Nissen,et al.  Stable Difference Methods for Block-Oriented Adaptive Grids , 2014, J. Sci. Comput..

[28]  Jan Nordström,et al.  Stable, high order accurate adaptive schemes for long time, highly intermittent geophysics problems , 2014, J. Comput. Appl. Math..

[29]  Ken Mattsson,et al.  A high-order accurate embedded boundary method for first order hyperbolic equations , 2017, J. Comput. Phys..

[30]  John C. Strikwerda.  HIGH‐ORDER‐ACCURATE SCHEMES FOR INCOMPRESSIBLE VISCOUS FLOW , 1997 .

[31]  Jan Nordström,et al.  High order finite difference methods for wave propagation in discontinuous media , 2006, J. Comput. Phys..

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

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

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

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

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

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

[38]  Ken Mattsson,et al.  Acoustic Wave Propagation in Complicated Geometries and Heterogeneous Media , 2014, J. Sci. Comput..

[39]  Kenneth Duru,et al.  Stable and high order accurate difference methods for the elastic wave equation in discontinuous media , 2014, J. Comput. Phys..

[40]  Matteo Parsani,et al.  Entropy stable wall boundary conditions for the three-dimensional compressible Navier-Stokes equations , 2014, J. Comput. Phys..

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

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

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

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

[45]  Ken Mattsson,et al.  High-fidelity numerical simulation of solitons in the nerve axon , 2016, J. Comput. Phys..

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