Abstract: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 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 A x = 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.

A locally second order symmetric method for discontinuous solution of Poisson's equation on uniform cartesian grids

An Approach to Numerically Solving the Poisson Equation

High order solution of Poisson problems with piecewise constant coefficients and interface jumps

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

An efficient neural-network and finite-difference hybrid method for elliptic interface problems with applications

A computational approach to solving a second-order singularly perturbed Fredholm integro-differential equation with discontinuous source term

An optimization-based approach for high-order accurate discretization of conservation laws with discontinuous solutions

Sixth-order compact difference scheme and multigrid method for solving the Poisson equation

Local Discontinuous Galerkin Methods for Solving Convection-Diffusion and Cahn-Hilliard Equations on Surfaces

A direct method for the solution of poisson's equation with neumann boundary conditions on a staggered grid of arbitrary size

A finite difference - discontinuous Galerkin method for the wave equation in second order form

A hybrid approach for singularly perturbed parabolic problem with discontinuous data

A finite difference scheme for two-dimensional singularly perturbed convection-diffusion problem with discontinuous source term

A Graded Mesh Refinement for 2D Poisson's Equation on Non Convex Polygonal Domains

Efficient finite element method for 2D singularly perturbed parabolic convection diffusion problems with discontinuous source term

$C^0$ discontinuous Galerkin finite element methods for second order linear elliptic partial differential equations in non-divergence form

Solution of 1D Poisson Equation with Neumann-Dirichlet and Dirichlet-Neumann Boundary Conditions, Using the Finite Difference Method

A selective immersed discontinuous Galerkin method for elliptic interface problems

A Fast Algorithm for Solving the Poisson Equation on a Nested Grid

Numerical treatment for the solution of singularly perturbed pseudo-parabolic problem on an equidistributed grid

Stability and Error Analysis of an Efficient Numerical Method for Convection Dominated Parabolic PDEs with Jump Discontinuity in Source Function on Modified Layer-Adapted Mesh