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

by   Qiwei Feng, et al.

In this paper, we develop sixth-order hybrid finite difference methods (FDMs) for the elliptic interface problem -∇·( a∇ u)=f in Ω\Γ, where Γ is a smooth interface inside Ω. The variable scalar coefficient a>0 and source f are possibly discontinuous across Γ. The hybrid FDMs utilize a 9-point compact stencil at any interior regular point of the grid and a 13-point stencil at irregular points near Γ. For interior regular points away from Γ, we obtain a sixth-order 9-point compact FDM satisfying the M-matrix property. Consequently, for the elliptic problem without interface (i.e., Γ is empty), our compact FDM satisfies the discrete maximum principle, which guarantees the theoretical sixth-order convergence. We also derive sixth-order compact (4-point for corners and 6-point for edges) FDMs having the M-matrix property at any boundary point subject to (mixed) Dirichlet/Neumann/Robin boundary conditions. For irregular points near Γ, we propose fifth-order 13-point FDMs, whose stencil coefficients can be effectively calculated by recursively solving several small linear systems. Theoretically, the proposed high order FDMs use high order (partial) derivatives of the coefficient a, the source term f, the interface curve Γ, the two jump functions along Γ, and the functions on ∂Ω. Numerically, we always use function values to approximate all required high order (partial) derivatives in our hybrid FDMs without losing accuracy. Our proposed FDMs are independent of the choice representing Γ and are also applicable if the jump conditions on Γ only depend on the geometry (e.g., curvature) of the curve Γ. Our numerical experiments confirm the sixth-order convergence in the l_∞ norm of the proposed hybrid FDMs for the elliptic interface problem.


page 14

page 15

page 16


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

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

A High Order Compact Finite Difference Scheme for Elliptic Interface Problems with Discontinuous and High-Contrast Coefficients

The elliptic interface problems with discontinuous and high-contrast coe...

Can 4th-order compact schemes exist for flux type BCs?

In this paper new innovative fourth order compact schemes for Robin and ...

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

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

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

In this paper we develop finite difference schemes for elliptic problems...

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...

Please sign up or login with your details

Forgot password? Click here to reset