摘要

In this article, we present and analyse an unfitted mesh method for the Poisson interface problem. By constructing a novel ansatz function in the vicinity of the interface, we are able to derive an extended Poisson problem whose interface fits a given quasi-uniform triangular mesh. Then we adopt a hybridizable discontinuous Galerkin method to solve the extended problem with an appropriate choice of flux for treating the jump conditions. In contrast with existing approaches, the ansatz function is designed through a delicate piecewise quadratic Hermite polynomial interpolation with a post-processing via a standard Lagrange polynomial interpolation. Such an explicit function offers a third-order approximation to the singular part of the underlying solution for interfaces of any shape. It is also essential for both stability and convergence of the proposed method. Moreover, we provide rigorous error analysis to showthat the scheme can achieve a second-order convergence rate for the approximation of the solution and its gradient. Ample numerical examples with complex interfaces demonstrate the expected convergence order and robustness of the method.

全文