Compact 9-Point Finite Difference Methods with High Accuracy Order and/or M-Matrix Property for Elliptic Cross-Interface Problems

by   Qiwei Feng, et al.

In this paper we develop finite difference schemes for elliptic problems with piecewise continuous coefficients that have (possibly huge) jumps across fixed internal interfaces. In contrast with such problems involving one smooth non-intersecting interface, that have been extensively studied, there are very few papers addressing elliptic interface problems with intersecting interfaces of coefficient jumps. It is well known that if the values of the permeability in the four subregions around a point of intersection of two such internal interfaces are all different, the solution has a point singularity that significantly affects the accuracy of the approximation in the vicinity of the intersection point. In the present paper we propose a fourth-order 9-point finite difference scheme on uniform Cartesian meshes for an elliptic problem whose coefficient is piecewise constant in four rectangular subdomains of the overall two-dimensional rectangular domain. Moreover, for the special case when the intersecting point of the two lines of coefficient jumps is a grid point, such a compact scheme, involving relatively simple formulas for computation of the stencil coefficients, can even reach sixth order of accuracy. Furthermore, we show that the resulting linear system for the special case has an M-matrix, and prove the theoretical sixth order convergence rate using the discrete maximum principle. Our numerical experiments demonstrate the fourth (for the general case) and sixth (for the special case) accuracy orders of the proposed schemes. In the general case, we derive a compact third-order finite difference scheme, also yielding a linear system with an M-matrix. In addition, using the discrete maximum principle, we prove the third order convergence rate of the scheme for the general elliptic cross-interface problem.


page 1

page 2

page 3

page 4


Hybrid Finite Difference Schemes for Elliptic Interface Problems with Discontinuous and High-Contrast Variable Coefficients

For elliptic interface problems with discontinuous coefficients, the max...

High order two-grid finite difference methods for interface and internal layer problems

Second order accurate Cartesian grid methods have been well developed fo...

Sixth-Order Hybrid FDMs and/or the M-Matrix Property for Elliptic Interface Problems with Mixed Boundary Conditions

In this paper, we develop sixth-order hybrid finite difference methods (...

Sixth Order Compact Finite Difference Scheme for Poisson Interface Problem with Singular Sources

Let Γ be a smooth curve inside a two-dimensional rectangular region Ω. I...

A Novel HOC-Immersed Interface Approach For Elliptic Problems

We present a new higher-order accurate finite difference explicit jump I...

A unified framework to generate optimized compact finite difference schemes

A unified framework to derive optimized compact schemes for a uniform gr...

A Compact Coupling Interface Method with Accurate Gradient Approximation for Elliptic Interface Problems

Elliptic interface boundary value problems play a major role in numerous...

Please sign up or login with your details

Forgot password? Click here to reset