Inverses of SBP-SAT finite difference operators approximating the first and second derivative

04/19/2020
by   Sofia Eriksson, et al.
0

The scalar, one-dimensional advection equation and heat equation are considered. These equations are discretized in space, using a finite difference method satisfying summation-by-parts (SBP) properties. To impose the boundary conditions, we use a penalty method called simultaneous approximation term (SAT). Together, this gives rise to two semi-discrete schemes where the discretization matrices approximate the first and the second derivative operators, respectively. The discretization matrices depend on free parameters from the SAT treatment. We derive the inverses of the discretization matrices, interpreting them as discrete Green's functions. In this direct way, we also find out precisely which choices of SAT parameters that make the discretization matrices singular. In the second derivative case, it is shown that if the penalty parameters are chosen such that the semi-discrete scheme is dual consistent, the discretization matrix can become singular even when the scheme is energy stable. The inverse formulas hold for SBP-SAT operators of arbitrary order of accuracy. For second and fourth order accurate operators, the inverses are provided explicitly.

READ FULL TEXT
POST COMMENT

Comments

There are no comments yet.

Authors

page 1

page 2

page 3

page 4

03/02/2021

An energy-based summation-by-parts finite difference method for the wave equation in second order form

We develop an energy-based finite difference method for the wave equatio...
10/18/2021

A Stable FDTD Subgridding Scheme with SBP-SAT for Transient Electromagnetic Analysis

We proposed a provably stable FDTD subgridding method for accurate and e...
02/09/2021

Stability and Functional Superconvergence of Narrow-Stencil Second-Derivative Generalized Summation-By-Parts Discretizations

We analyze the stability and functional superconvergence of discretizati...
04/06/2021

Inverse properties of a class of pentadiagonal matrices related to higher order difference operators

This paper analyzes the convergence of fixed-point iterations of the for...
07/20/2019

Non-stiff narrow-stencil finite difference approximations of the Laplacian on curvilinear multiblock grids

The Laplacian appears in several partial differential equations used to ...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.