English

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

Numerical Analysis 2021-04-19 v1 Numerical Analysis

Abstract

Let Γ\Gamma be a smooth curve inside a two-dimensional rectangular region Ω\Omega. In this paper, we consider the Poisson interface problem 2u=f-\nabla^2 u=f in ΩΓ\Omega\setminus \Gamma with Dirichlet boundary condition such that ff is smooth in ΩΓ\Omega\setminus \Gamma and the jump functions [u][u] and [un][\nabla u\cdot \vec{n}] across Γ\Gamma are smooth along Γ\Gamma. This Poisson interface problem includes the weak solution of 2u=f+gδΓ-\nabla^2 u=f+g\delta_\Gamma in Ω\Omega as a special case. Because the source term ff is possibly discontinuous across the interface curve Γ\Gamma and contains a delta function singularity along the curve Γ\Gamma, both the solution uu of the Poisson interface problem and its flux un\nabla u\cdot \vec{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 AA in the resulting linear system Ax=bAx=b, following from the proposed scheme, is independent of any source term ff, jump condition gδΓg\delta_\Gamma, interface curve Γ\Gamma 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.

Keywords

Cite

@article{arxiv.2104.07866,
  title  = {Sixth Order Compact Finite Difference Scheme for Poisson Interface Problem with Singular Sources},
  author = {Qiwei Feng and Bin Han and Peter Minev},
  journal= {arXiv preprint arXiv:2104.07866},
  year   = {2021}
}

Comments

27 pages, 12 figures, 8 tables

R2 v1 2026-06-24T01:13:42.359Z