RECENT DEVELOPMENT OF IMMERSED FEM FOR ELLIPTIC AND ELASTIC INTERFACE PROBLEMS

Abstract

We survey a recently developed immersed finite element method (IFEM) for the interface problems. The IFEM uses structured grids such as uniform grids, even if the interface is a smooth curve. Instead of fitting the curved interface, the bases are modified so that they satisfy the jump conditions along the interface. The early versions of IFEM [1, 2] were suboptimal in convergence order [3]. Later, the consistency terms were added to the bilinear forms [4, 5], thus the scheme became optimal and the error estimates were proven. For elasticity problems with interfaces, we modify the Crouzeix-Raviart based element to satisfy the traction conditions along the interface [6], but the consistency terms are not needed. To satisfy the Korn's inequality, we add the stabilizing terms to the bilinear form. The optimal error estimate was shown for a triangular grid. Lastly, we describe the multigrid algorithms for the discretized system arising from IFEM. The prolongation operators are designed so that the prolongated function satisfy the flux continuity condition along the interface. The W-cycle convergence was proved, and the number of V-cycle is independent of the mesh size.