Pith. sign in

REVIEW

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

Not yet reviewed by Pith; the record is open.

This paper has not been read by Pith yet. Machine review is queued; the pith claim, tier, and objections will appear here once it completes.

SPECIMEN: schema-true, not a live event

T0 review · schema-true

One-sentence machine reading of the paper's core claim.

pith:XXXXXXXX · record.json · timestamp

arxiv 2104.07866 v1 pith:7FLA37KT submitted 2021-04-16 math.NA cs.NA

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

classification math.NA cs.NA
keywords gammainterfacepoissoncompactdifferencefiniteproblemscheme
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved
0 comments
read the original abstract

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

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.