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

by   Qiwei Feng, et al.

Let Γ be a smooth curve inside a two-dimensional rectangular region Ω. In this paper, we consider the Poisson interface problem -∇^2 u=f in Ω∖Γ with Dirichlet boundary condition such that f is smooth in Ω∖Γ and the jump functions [u] and [∇ u·n⃗] across Γ are smooth along Γ. This Poisson interface problem includes the weak solution of -∇^2 u=f+gδ_Γ in Ω as a special case. Because the source term f is possibly discontinuous across the interface curve Γ and contains a delta function singularity along the curve Γ, both the solution u of the Poisson interface problem and its flux ∇ u·n⃗ are often discontinuous across the interface. To solve the Poisson interface problem with singular sources, in this paper we propose a sixth order compact finite difference scheme on uniform Cartesian grids. Our proposed compact finite difference scheme with explicitly given stencils extends the immersed interface method (IIM) to the highest possible accuracy order six for compact finite difference schemes on uniform Cartesian grids, but without the need to change coordinates into the local coordinates as in most papers on IIM in the literature. Also in contrast with most published papers on IIM, we explicitly provide the formulas for all involved stencils. The coefficient matrix A in the resulting linear system Ax=b, following from the proposed scheme, is independent of any source term f, jump condition gδ_Γ, interface curve Γ and Dirichlet boundary conditions. Our numerical experiments confirm the sixth accuracy order of the proposed compact finite difference scheme on uniform meshes for the Poisson interface problems with various singular sources.


page 1

page 2

page 3

page 4


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

Sixth Order Compact Finite Difference Method for 2D Helmholtz Equations with Singular Sources and Reduced Pollution Effect

Due to its highly oscillating solution, the Helmholtz equation is numeri...

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

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

A convergent hybrid neural network and finite difference scheme for Stokes interface problems

In this paper, we present a novel hybrid method for solving a Stokes int...

Decoupling Numerical Method Based on Deep Neural Network for Nonlinear Degenerate Interface Problems

Interface problems depict many fundamental physical phenomena and widely...

Please sign up or login with your details

Forgot password? Click here to reset