Abstract:In this paper, we discuss the second-order finite element method (FEM) and finite difference method (FDM) for numerically solving elliptic cross-interface problems characterized by vertical and horizontal straight lines, piecewise constant coefficients, two homogeneous jump conditions, continuous source terms, and Dirichlet boundary conditions. For brevity, we consider a 2D simplified version where the intersection points of the interface lines coincide with grid points in uniform Cartesian grids. Our findings reveal interesting and important results: (1) When the coefficient functions exhibit either high jumps with low-frequency oscillations or low jumps with high-frequency oscillations, the finite element method and finite difference method yield similar numerical solutions. (2) However, when the interface problems involve high-contrast and high-frequency coefficient functions, the numerical solutions obtained from the finite element and finite difference methods differ significantly. Given that the widely studied SPE10 benchmark problem (see <a class="link-external link-https" href="https://www.spe.org/web/csp/datasets/set02.htm" rel="external noopener nofollow">this https URL</a>) typically involves high-contrast and high-frequency permeability due to varying geological layers in porous media, this phenomenon warrants attention. Furthermore, this observation is particularly important for developing multiscale methods, as reference solutions for these methods are usually obtained using the standard second-order finite element method with a fine mesh, and analytical solutions are not available. We provide sufficient details to enable replication of our numerical results, and the implementation is straightforward. This simplicity ensures that readers can easily confirm the validity of our findings.

Tailored Finite Point Method for the Interface Problem

Tailored Finite Point Method for a Singular Perturbation Problem with Variable Coefficients in Two Dimensions

An accurate finite element method for elliptic interface problems

Tailored Finite Point Method for First Order Wave Equation

A high order compact finite difference scheme for elliptic interface problems with discontinuous and high-contrast coefficients

A Tailored Finite Point Method for the Helmholtz Equation with High Wave Numbers in Heterogeneous Medium

A class of HOC finite difference method for elliptic interface problems with imperfect contact

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

The Tailored Finite Point Method

A high order finite difference method for the elastic wave equation in bounded domains with nonconforming interfaces

Higher-order finite element methods for elliptic problems with interfaces

A Nonconforming Extended Virtual Element Method for Elliptic Interface Problems

A Tailored Finite Point Method for a Singular Perturbation Problem on an Unbounded Domain

Immersed Virtual Element Methods for Elliptic Interface Problems in Two Dimensions

A selective immersed discontinuous Galerkin method for elliptic interface problems

High order finite difference methods for the wave equation with non-conforming grid interfaces

The Immersed Weak Galerkin and Continuous Galerkin Finite Element Method for Elliptic Interface Problem

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

An Interface Penalty Finite Element Method for Elliptic Interface Problems on Piecewise Meshes

Multiscale Tailored Finite Point Method for Second Order Elliptic Equations with Rough or Highly Oscillatory Coefficients

A direct finite element method for elliptic interface problems