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

TL;DR

This paper compares finite element and finite difference methods for solving elliptic cross-interface problems, revealing their similarities under certain conditions and significant differences in high-contrast, high-frequency scenarios, with implications for multiscale modeling.

Contribution

It provides a detailed comparison of FEM and FDM for elliptic cross-interface problems, highlighting conditions where solutions diverge, especially in high-contrast, high-frequency cases, and offers practical implementation details.

Findings

FEM and FDM solutions are similar for low-frequency oscillations.

Solutions differ significantly in high-contrast, high-frequency cases.

Results are relevant for multiscale methods and benchmark problems.

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 https://www.spe.org/web/csp/datasets/set02.htm) 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.