Summation-by-parts approximations of the second derivative: Pseudoinverses of singular operators and revisiting the sixth order accurate narrow-stencil operator

by   Sofia Eriksson, et al.

We consider finite difference approximations of the second derivative, exemplified in Poisson's equation, the heat equation and the wave equation. The finite difference operators satisfy a summation-by-parts property, which mimics the integration-by-parts. Since the operators approximate the second derivative, they are singular by construction. To impose boundary conditions, these operators are modified using Simultaneous Approximation Terms. This makes the modified matrices non-singular, for most choices of boundary conditions. Recently, inverses of such matrices were derived. However, when considering Neumann boundary conditions on both boundaries, the modified matrix is still singular. For such matrices, we have derived an explicit expression for the Moore-Penrose pseudoinverse, which can be used for solving elliptic problems and some time-dependent problems. The condition for this new pseudoinverse to be valid, is that the modified matrix does not have more than one zero eigenvalue. We have reconstructed the sixth order accurate narrow-stencil operator with a free parameter and show that more than one zero eigenvalue can occur. We have performed a detailed analysis on the free parameter to improve the properties of the second derivative operator. We complement the derivations by numerical experiments to demonstrate the improvements of the new second derivative operator.



There are no comments yet.


page 1

page 2

page 3

page 4


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

The scalar, one-dimensional advection equation and heat equation are con...

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

We analyze the stability and functional superconvergence of discretizati...

Numerical analysis of a self-similar turbulent flow in Bose–Einstein condensates

We study a self-similar solution of the kinetic equation describing weak...

Local Compatibility Boundary Conditions for High-Order Accurate Finite-Difference Approximations of PDEs

We describe a new approach to derive numerical approximations of boundar...

Finite difference approach to fourth-order linear boundary-value problems

Discrete approximations to the equation L_contu = u^(4) + D(x) u^(3) ...

Bi-Parametric Operator Preconditioning

We extend the general operator preconditioning framework [R. Hiptmair, C...

On computing derivatives of transfer operators and linear responses in higher dimensions

We show that the derivative of the transfer operator with respect to per...
This week in AI

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