High-order Finite Element--Integral Equation Coupling on Embedded Meshes

TL;DR

This paper introduces a high-order coupled finite element and integral equation method for interface problems on embedded meshes, achieving high accuracy and efficiency without modifying for jump conditions.

Contribution

It develops new FE-IE formulations for exterior and interface problems, enabling high-order convergence on embedded meshes with jump conditions.

Findings

Achieves high-order convergence close to embedded interfaces.

Efficient iterative solver with algebraic multigrid preconditioning.

Numerical results confirm high-order accuracy and efficiency.

Abstract

This paper presents a high-order method for solving an interface problem for the Poisson equation on embedded meshes through a coupled finite element and integral equation approach. The method is capable of handling homogeneous or inhomogeneous jump conditions without modification and retains high-order convergence close to the embedded interface. We present finite element-integral equation (FE-IE) formulations for interior, exterior, and interface problems. The treatments of the exterior and interface problems are new. The resulting linear systems are solved through an iterative approach exploiting the second-kind nature of the IE operator combined with algebraic multigrid preconditioning for the FE part. Assuming smooth continuations of coefficients and right-hand-side data, we show error analysis supporting high-order accuracy. Numerical evidence further supports our claims of efficiency and high-order accuracy for smooth data.