Implicit Eulerian Fluid-Structure Interaction Method for the Modeling of Highly Deformable Elastic Membranes

This paper is concerned with the development of a fully implicit and purely Eulerian fluid-structure interaction method tailored for the modeling of the large deformations of elastic membranes in a surrounding Newtonian fluid. We consider a simplified model for the mechanical properties of the membrane, in which the surface strain energy depends on the membrane stretching. The fully Eulerian description is based on the advection of a modified surface tension tensor, and the deformations of the membrane are tracked using a level set strategy. The resulting nonlinear problem is solved by a Newton-Raphson method, featuring a quadratic convergence behavior. A monolithic solver is implemented, and we report several numerical experiments aimed at model validation and illustrating the accuracy of the presented method. We show that stability is maintained for significantly larger time steps.

[1]  David Salac,et al.  A level set projection model of lipid vesicles in general flows , 2011, J. Comput. Phys..

[2]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[3]  Ming-Chih Lai,et al.  Simulating the dynamics of inextensible vesicles by the penalty immersed boundary method , 2010, J. Comput. Phys..

[4]  Harald Garcke,et al.  Finite element approximation for the dynamics of asymmetric fluidic biomembranes , 2016, Math. Comput..

[5]  Saroja Ramanujan,et al.  Deformation of liquid capsules enclosed by elastic membranes in simple shear flow: large deformations and the effect of fluid viscosities , 1998, Journal of Fluid Mechanics.

[6]  I. C. Howard,et al.  On the opening mechanism of the aortic valve: some observations from simulations , 2003, Journal of medical engineering & technology.

[7]  Jean-Yves L'Excellent,et al.  MUMPS: A Multifrontal Massively Parallel Solver , 2002 .

[8]  Alfio Quarteroni,et al.  Numerical modeling of heart valves using resistive Eulerian surfaces , 2016, International journal for numerical methods in biomedical engineering.

[9]  Zhilin Li,et al.  An immersed interface method for solving incompressible viscous flows with piecewise constant viscosity across a moving elastic membrane , 2008, J. Comput. Phys..

[10]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[11]  S. Schmid,et al.  GTPase Cycle of Dynamin Is Coupled to Membrane Squeeze and Release, Leading to Spontaneous Fission , 2008, Cell.

[12]  Gábor Székely,et al.  Eulerian finite element method for the numerical modeling of fluid dynamics of natural and pathological aortic valves , 2017, J. Comput. Appl. Math..

[13]  Boyce E. Griffith,et al.  Immersed boundary model of aortic heart valve dynamics with physiological driving and loading conditions , 2012, International journal for numerical methods in biomedical engineering.

[14]  S. Osher,et al.  Algorithms Based on Hamilton-Jacobi Formulations , 1988 .

[15]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[16]  Message Passing Interface Forum.  MPI: A message - passing interface standard , 1994 .

[17]  George Biros,et al.  Author ' s personal copy Dynamic simulation of locally inextensible vesicles suspended in an arbitrary two-dimensional domain , a boundary integral method , 2010 .

[18]  Pierre Saramito,et al.  Computing the dynamics of biomembranes by combining conservative level set and adaptive finite element methods , 2014, J. Comput. Phys..

[19]  Thomas Williams,et al.  Gnuplot 4.4: an interactive plotting program , 2010 .

[20]  Dominique Barthès-Biesel,et al.  Large deformations and burst of a capsule freely suspended in an elongational flow , 1988, Journal of Fluid Mechanics.

[21]  Ming-Chih Lai,et al.  An immersed boundary method for simulating vesicle dynamics in three dimensions , 2016, J. Comput. Phys..

[22]  Udo Seifert,et al.  Configurations of fluid membranes and vesicles , 1997 .

[23]  Morton E. Gurtin,et al.  A continuum theory of elastic material surfaces , 1975 .

[24]  Boyce E. Griffith,et al.  A finite strain nonlinear human mitral valve model with fluid-structure interaction , 2014, International journal for numerical methods in biomedical engineering.

[25]  M. Thubrikar,et al.  Stress analysis of the aortic valve with and without the sinuses of valsalva. , 2001, The Journal of heart valve disease.

[26]  Pierre Saramito,et al.  An adaptive finite element method for the modeling of the equilibrium of red blood cells , 2016 .

[27]  Alfio Quarteroni,et al.  Fully Eulerian finite element approximation of a fluid‐structure interaction problem in cardiac cells , 2013 .

[28]  S. Hysing,et al.  A new implicit surface tension implementation for interfacial flows , 2006 .

[29]  高等学校計算数学学報編輯委員会編.  高等学校計算数学学報 = Numerical mathematics , 1979 .

[30]  P. Saramito.  Efficient C++ finite element computing with Rheolef , 2012 .

[31]  Dominique Barthès Biesel.  Motion and Deformation of Elastic Capsules and Vesicles in Flow , 2016 .

[32]  F P T Baaijens,et al.  A three-dimensional computational analysis of fluid-structure interaction in the aortic valve. , 2003, Journal of biomechanics.

[33]  G. Cottet,et al.  A level-set formulation of immersed boundary methods for fluid–structure interaction problems , 2004 .

[34]  M. Lai,et al.  An Immersed Boundary Method with Formal Second-Order Accuracy and Reduced Numerical Viscosity , 2000 .

[35]  Patrick Amestoy,et al.  A Fully Asynchronous Multifrontal Solver Using Distributed Dynamic Scheduling , 2001, SIAM J. Matrix Anal. Appl..