Abstract:In this paper, we develop sixth-order hybrid finite difference methods (FDMs) for the elliptic interface problem $-\nabla \cdot( a\nabla u)=f$ in $\Omega\backslash \Gamma$, where $\Gamma$ is a smooth interface inside $\Omega$. The variable scalar coefficient $a>0$ and source $f$ are possibly discontinuous across $\Gamma$. The hybrid FDMs utilize a $9$-point compact stencil at any interior regular points of the grid and a $13$-point stencil at irregular points near $\Gamma$. For interior regular points away from $\Gamma$, we obtain a sixth-order $9$-point compact FDM satisfying the sign and sum conditions for ensuring the M-matrix property. We also derive sixth-order compact ($4$-point for corners and $6$-point for edges) FDMs satisfying the sign and sum conditions for the M-matrix property at any boundary point subject to (mixed) Dirichlet/Neumann/Robin boundary conditions. Thus, for the elliptic problem without interface (i.e., $\Gamma$ is empty), our compact FDM has the M-matrix property for any mesh size $h>0$ and consequently, satisfies the discrete maximum principle, which guarantees the theoretical sixth-order convergence. For irregular points near $\Gamma$, we propose fifth-order $13$-point FDMs, whose stencil coefficients can be effectively calculated by recursively solving several small linear systems. Theoretically, the proposed high order FDMs use high order (partial) derivatives of the coefficient $a$, the source term $f$, the interface curve $\Gamma$, the two jump functions along $\Gamma$, and the functions on $\partial \Omega$. Numerically, we always use function values to approximate all required high order (partial) derivatives in our hybrid FDMs without losing accuracy. Our numerical experiments confirm the sixth-order convergence in the $l_{\infty}$ norm of the proposed hybrid FDMs for the elliptic interface problem.

A Combined Finite Element Method for Elliptic Problems Posted in Domains with Rough Boundaries

An adaptive non-symmetric finite volume and boundary element coupling method for a fluid mechanics interface problem

Stable non-symmetric coupling of the finite volume and the boundary element method for convection-dominated parabolic-elliptic interface problems

A CutFE-LOD Method for the Multiscale Elliptic Problems on Complex Domains

An Adaptive Scaled Boundary Finite Element Method by Subdividing Subdomains for Elastodynamic Problems

Distinct Numerical Solutions for Elliptic Cross-Interface Problems Using Finite Element and Finite Difference Methods

A Nonconforming Extended Virtual Element Method for Elliptic Interface Problems

An extended mixed finite element method for elliptic interface problems

Sixth-Order Hybrid Finite Difference Methods for Elliptic Interface Problems with Mixed Boundary Conditions

An extended virtual element method for elliptic interface problems

Solution of the Elliptic Interface Problem by a Hybrid Mixed Finite Element Method

A Combined Finite Element and Multiscale Finite Element Method for the Multiscale Elliptic Problems

A direct finite element method for elliptic interface problems

A Combined Multiscale Finite Element Method Based on the LOD Technique for the Multiscale Elliptic Problems with Singularities

Boundary integro-differential equations of elliptic boundary value problems and their numerical solutions

Immersed Virtual Element Methods for Elliptic Interface Problems in Two Dimensions

Biquadratic element discrete duality finite volume method for solving elliptic equations on quadrilateral mesh

Superconvergence of Modified Nonconforming Cut Finite Element Method for Elliptic Problems

Adaptive Finite Element Method for a Nonlinear Helmholtz Equation with High Wave Number

Adaptive FEM for Helmholtz Equation with Large Wavenumber

An arbitrarily high order unfitted finite element method for elliptic interface problems with automatic mesh generation, Part II. Piecewise-smooth interfaces