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 $C^0$ finite element algorithm for the sixth order problem with simply supported boundary conditions

A $C^0$ finite element method for the biharmonic problem with Dirichlet boundary conditions in a polygonal domain

A mixed finite element method for a sixth order elliptic problem

A New Stabilized Nonconforming-Mixed Finite Element Method for the Second Order Elliptic Boundary Value Problem

Mixed Finite Element Method and Higher-Order Local Artificial Boundary Conditions for Exterior 3-D Poisson Equation

Two families of $n$-rectangle nonconforming finite elements for sixth-order elliptic equations

The Numerical Solution of Linear Sixth Order Boundary Value Problems with Quartic B-Splines.

A robust lower order mixed finite element method for a strain gradient elasticity model

The polygonal scaled boundary thin plate element based on the discrete Kirchhoff theory

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

A New Least-Squares Nonconforming Mixed Finite Element Analysis for Second Order Elliptic Problems

C~0-Continuous Time Stepping Linear Finite Element Method for Second Order Linear Hyperbolic Equations

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

A sixth-order wavelet integral collocation method for solving nonlinear boundary value problems in three dimensions

A new mixed finite-element method for $H^2$ elliptic problems

A recovery-based linear C0 finite element method for a fourth-order singularly perturbed Monge-Ampère equation

An Arbitrary-Order Discontinuous Galerkin Method with One Unknown Per Element

A new recovery based C 0 element method for fourth-order equations

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

Geometrical mixed finite element methods for fourth order obstacle problems in linearised elasticity

Adaptive discontinuous Galerkin finite element methods for the Allen-Cahn equation on polygonal meshes