Order Preserving Interpolation for Summation-by-Parts Operators at Non-Conforming Grid Interfaces

Martin Almquist,Siyang Wang,Jonatan Werpers

DOI: https://doi.org/10.1137/18M1191609

2018-06-06

Abstract:We study non-conforming grid interfaces for summation-by-parts finite difference methods applied to partial differential equations with second derivatives in space. To maintain energy stability, previous efforts have been forced to accept a reduction of the global convergence rate by one order, due to large truncation errors at the non-conforming interface. We avoid the order reduction by generalizing the interface treatment and introducing order preserving interpolation operators. We prove that, given two diagonal-norm summation-by-parts schemes, order preserving interpolation operators with the necessary properties are guaranteed to exist, regardless of the grid-point distributions along the interface. The new methods retain the stability and global accuracy properties of the underlying schemes for conforming interfaces.

Numerical Analysis

What problem does this paper attempt to address?

This paper aims to solve the problems that occur when using non - conforming grid interfaces for the numerical solution of partial differential equations. Specifically, when using partial differential equations with second - order spatial derivatives, the traditional Summation - by - Parts (SBP) finite - difference method will introduce large truncation errors at non - conforming grid interfaces, resulting in a one - order reduction in the global convergence rate. To solve this problem, the paper proposes Order Preserving Interpolation Operators to avoid the order reduction at non - conforming grid interfaces.

### Background and Problem Description of the Paper

In numerical simulation, Adaptive Mesh Refinement (AMR) is crucial for improving the efficiency of locally high - resolution areas. For wave - dominated phenomena, high - order Finite Difference (FD) methods usually have high computational efficiency but are not always robust. By combining Summation - by - Parts (SBP) operators and Simultaneous Approximation Terms (SATs), the SBP - SAT method can achieve energy - stable and conservative high - order finite - difference methods, which are suitable for multi - block and curvilinear grids.

However, when dealing with non - conforming grid interfaces, the SBP - SAT method faces a problem: at non - conforming grid interfaces, SATs will worsen the maximum local truncation error, thus resulting in a one - order reduction in the global convergence rate. This is because the order of the interpolation operator is limited by the upper limit of the order of the quadrature rule associated with each SBP operator.

### Solution

To overcome this problem, the paper proposes a new interpolation operator - Order Preserving Interpolation Operators (OP operators). The key feature of these operators is that they exist in pairs, and each pair of operators satisfies certain accuracy limitations, but through specific SATs, the accuracy reduction at the interface can be avoided.

Specifically, the paper proves that given two SBP operators based on diagonal norms, it is always possible to find matching OP interpolation operators. These new operators show the same global convergence rate as the conforming grid interface in numerical experiments and are applicable to the heat equation, the Schrödinger equation, and the second - order wave equation.

### Main Contributions

1. **Theoretical Results**: It is proved that given two SBP operators based on diagonal norms, it is always possible to find matching OP interpolation operators.
2. **Numerical Experiments**: It is shown that after using OP interpolation operators, the global convergence rates of the heat equation, the Schrödinger equation, and the second - order wave equation are the same as those of the conforming grid interface.
3. **Stability Analysis**: By introducing the concept of the Hilbert adjoint operator, the stability of the interpolation operator is ensured.

### Conclusion

The method proposed in the paper successfully solves the order - reduction problem at non - conforming grid interfaces and improves the overall stability and accuracy of the numerical method. These results are of great significance for the efficient numerical simulation of multi - scale problems.