Semi-implicit formulation of the immersed finite element method

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.

[1]  K. Bube,et al.  The Immersed Interface Method for Nonlinear Differential Equations with Discontinuous Coefficients and Singular Sources , 1998 .

[2]  Luca Heltai,et al.  On the CFL condition for the finite element immersed boundary method , 2007 .

[3]  Wing Kam Liu,et al.  Modelling and simulation of fluid structure interaction by meshfree and FEM , 2003 .

[4]  T. Tezduyar.  Computation of moving boundaries and interfaces and stabilization parameters , 2003 .

[5]  Wing Kam Liu,et al.  Extended immersed boundary method using FEM and RKPM , 2004 .

[6]  C. Peskin,et al.  Mechanical equilibrium determines the fractal fiber architecture of aortic heart valve leaflets. , 1994, The American journal of physiology.

[7]  Yaling Liu,et al.  Rheology of red blood cell aggregation by computer simulation , 2006, J. Comput. Phys..

[8]  C. Peskin.  Acta Numerica 2002: The immersed boundary method , 2002 .

[9]  L. Heltai,et al.  A finite element approach to the immersed boundary method , 2003 .

[10]  H. Othmer,et al.  Case Studies in Mathematical Modeling: Ecology, Physiology, and Cell Biology , 1997 .

[11]  Lucy T. Zhang,et al.  Stent modeling using immersed finite element method , 2006 .

[12]  Randall J. LeVeque,et al.  An Immersed Interface Method for Incompressible Navier-Stokes Equations , 2003, SIAM J. Sci. Comput..

[13]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .

[14]  Howard H. Hu,et al.  Direct numerical simulations of fluid-solid systems using the arbitrary Langrangian-Eulerian technique , 2001 .

[15]  C S Peskin,et al.  Computer-assisted design of pivoting disc prosthetic mitral valves. , 1983, The Journal of thoracic and cardiovascular surgery.

[16]  C S Peskin,et al.  Cardiac fluid dynamics. , 1992, Critical reviews in biomedical engineering.

[17]  C. Peskin.  The immersed boundary method , 2002, Acta Numerica.

[18]  Randall J. LeVeque,et al.  Immersed Interface Methods for Stokes Flow with Elastic Boundaries or Surface Tension , 1997, SIAM J. Sci. Comput..

[19]  James P. Keener,et al.  Immersed Interface Methods for Neumann and Related Problems in Two and Three Dimensions , 2000, SIAM J. Sci. Comput..

[20]  David Farrell,et al.  Immersed finite element method and its applications to biological systems. , 2006, Computer methods in applied mechanics and engineering.

[21]  Lucy T. Zhang,et al.  Interpolation functions in the immersed boundary and finite element methods , 2010 .

[22]  Hiroshi Terashima,et al.  A front-tracking/ghost-fluid method for fluid interfaces in compressible flows , 2009, J. Comput. Phys..

[23]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[24]  R. LeVeque,et al.  A comparison of the extended finite element method with the immersed interface method for elliptic equations with discontinuous coefficients and singular sources , 2006 .

[25]  Antonio Huerta,et al.  Viscous flow with large free surface motion , 1988 .

[26]  C. Peskin,et al.  A three-dimensional computational method for blood flow in the heart. 1. Immersed elastic fibers in a viscous incompressible fluid , 1989 .

[27]  C. Peskin.  Numerical analysis of blood flow in the heart , 1977 .

[28]  Lucy T. Zhang.  Immersed finite element method for fluid-structure interactions , 2007 .

[29]  Charles S. Peskin,et al.  A vortex-grid method for blood flow through heart valves , 1980 .

[30]  Wing Kam Liu,et al.  Lagrangian-Eulerian finite element formulation for incompressible viscous flows☆ , 1981 .

[31]  Zhilin Li,et al.  The immersed interface method for the Navier-Stokes equations with singular forces , 2001 .

[32]  Chi-Wang Shu,et al.  Total variation diminishing Runge-Kutta schemes , 1998, Math. Comput..

[33]  Wing Kam Liu,et al.  Computer implementation aspects for fluid-structure interaction problems , 1982 .

[34]  Lucy T. Zhang,et al.  Immersed finite element method , 2004 .

[35]  Ted Belytschko,et al.  Arbitrary Lagrangian-Eulerian Petrov-Galerkin finite elements for nonlinear continua , 1988 .

[36]  Andreas Wiegmann,et al.  The Explicit-Jump Immersed Interface Method: Finite Difference Methods for PDEs with Piecewise Smooth Solutions , 2000, SIAM J. Numer. Anal..