Sixth order compact finite difference schemes for Poisson interface problems with singular sources

• Published in: Computers & mathematics with applications (1987) Vol. 99; pp. 2 - 25
• Main Authors: Feng, Qiwei, Han, Bin, Minev, Peter
• Format: Journal Article
• Language: English
• Published: Oxford Elsevier Ltd 01-10-2021; Elsevier BV
• Subjects: Algorithms; Boundary conditions; Cartesian coordinates; Delta function; Delta source functions along curves; Dirichlet problem; Discontinuous and singular source terms; Finite difference method; Flow control; High order compact finite difference schemes; Mathematical analysis; Matrices (mathematics); Maximum principle; Piecewise smooth solutions; Poisson interface equations; The convergence proof; The convergence proof; Piecewise smooth solutions; Discontinuous and singular source terms; High order compact finite difference schemes; Delta source functions along curves; Poisson interface equations
• ISSN: 0898-1221; 1873-7668
• DOI: 10.1016/j.camwa.2021.07.020

Let Γ be a smooth curve inside a two-dimensional rectangular region Ω. In this paper, we consider the Poisson interface problem −∇2u=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 −∇2u=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 and therefore, our proposed scheme can be easily implemented and is of interest to practitioners dealing with Poisson interface problems. Note that the curve Γ splits Ω into two disjoint subregions Ω+ and Ω−. 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, while only b depends on these factors and is explicitly given, according to the configuration of the nine stencil points in Ω+ or Ω−. The constant coefficient matrix A facilitates the parallel implementation of the algorithm in case of a large size matrix and only requires the update of the right hand side vector b for different Poisson interface problems. Due to the flexibility and explicitness of the proposed scheme, it can be generalized to obtain the highest order compact finite difference scheme for non-uniform grids as well. We prove the order 6 convergence for the proposed scheme using the discrete maximum principle. 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.