Semi-implicit formulation of the immersed finite element method

Abstract

The immersed finite element method (IFEM) is a novel numerical approach to solve fluid---structure interaction types of problems that utilizes non-conforming meshing concept. The fluid and the solid domains are represented independently. The original algorithm of the IFEM follows the interpolation process as illustrated in the original immersed boundary method where the fluid velocity and the interaction force are explicitly coupled. However, the original approach presents many numerical difficulties when the fluid and solid physical properties have large mismatches, such as when the density difference is large and when the solid is a very stiff material. Both situations will lead to divergent or unstable solutions if not handled properly. In this paper, we develop a semi-implicit formulation of the IFEM algorithm so that several terms of the interfacial forces are implicitly evaluated without going through the force distribution process. Based on the 2-D and 3-D examples that we study in this paper, we show that the semi-implicit approach is robust and is capable of handling these highly discontinuous physical properties quite well without any numerical difficulties.