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.

On the Interface Identification of Free Boundary Problem by Method of Fundamental Solution.

Asymptotic Behavior of Solutions of a Free Boundary Problem Modelling the Growth of Tumors with Stokes Equations

A Volume of Fluid Method Algorithm for Simulation of Surface Tension Dominant Two-Phase Flows

Development of a Volume of Fluid Method for Computing Interfacial Incompressible Fluid Flows

On self-similar solutions of a multi-phase Stefan problem in the half-line

A Kind of Extended Korteweg-De Vries Equation and Solitary Wave Solutions for Interfacial Waves in A Two-Fluid System

Numerical studies of a class of reaction–diffusion equations with Stefan conditions

A stable and accurate scheme for solving the Stefan problem coupled with natural convection using the Immersed Boundary Smooth Extension method

A Cartesian grid-based boundary integral method for moving interface problems

Analytical solution of Stefan-type problems

One-phase Stefan problem with a latent heat depending on the position of the free boundary and its rate of change

An Extended Method of Fundamental Solutions to the Problems of Inverse Boundary Value of the Heat-conduction Equation Involving One Specific Type of Heat Source

Kernel-free boundary integral method for two-phase Stokes equations with discontinuous viscosity on staggered grids

A Deep Neural Network/Meshfree Method for Solving Dynamic Two-phase Interface Problems

Self-similar diffuse boundary method for phase boundary driven flow

Smoothed Boundary Method for Solving Partial Differential Equations with General Boundary Conditions on Complex Boundaries

Extended Smoothed Boundary Method for Solving Partial Differential Equations with General Boundary Conditions on Complex Boundaries

On the two-phase Navier-Stokes equations with surface tension

A Novel Supermesh Method for Computing Solutions to the Multi-material Stefan Problem with Complex Deforming Interfaces and Microstructure

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

A half boundary method for two dimensional unsteady convection–diffusion equations