Diagonal-norm upwind SBP operators

High-order accurate first derivative finite difference operators are derived that naturally introduce artificial dissipation. The boundary closures are based on the diagonal-norm summation-by-parts (SBP) framework and the boundary conditions are imposed using a penalty (SAT) technique, to guarantee linear stability for a large class of initial boundary value problems. These novel first derivative SBP operators have a non-central difference stencil in the interior, and come in pairs (for each order of accuracy). The resulting SBP-SAT approximations lead to fully explicit ODE systems. The accuracy and stability properties are demonstrated for linear first- and second-order hyperbolic problems in 1D, and for the compressible Euler equations in 2D. The newly derived first derivative SBP operators lead to significantly more robust and accurate numerical approximations, compared with the exclusive usage of (previously derived central) non-dissipative first derivative SBP operators.

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

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

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

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

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

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

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

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

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

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

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

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

[13]  Jason E. Hicken.  Output error estimation for summation-by-parts , 2014 .

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

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

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

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

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

[19]  Magnus Svärd,et al.  Well-Posed Boundary Conditions for the Navier-Stokes Equations , 2005, SIAM J. Numer. Anal..

[20]  N. Anders Petersson,et al.  Discretizing singular point sources in hyperbolic wave propagation problems , 2016, J. Comput. Phys..

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

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

[23]  Olsson,et al.  SUMMATION BY PARTS, PROJECTIONS, AND STABILITY. I , 2010 .

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

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

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

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

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

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

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

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

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

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

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

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

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

[37]  Leonid Dovgilovich,et al.  High-accuracy finite-difference schemes for solving elastodynamic problems in curvilinear coordinates within multiblock approach , 2014 .

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

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