Simulation of a pulsatile total artificial heart: Development of a partitioned Fluid Structure Interaction model

Abstract Heart disease is one of the leading causes of death in the world. Due to a shortage in donor organs artificial hearts can be a bridge to transplantation or even serve as a destination therapy for patients with terminal heart insufficiency. A pusher plate driven pulsatile membrane pump, the Total Artificial Heart (TAH) ReinHeart, is currently under development at the Institute of Applied Medical Engineering of RWTH Aachen University. This paper presents the methodology of a fully coupled three-dimensional time-dependent Fluid Structure Interaction (FSI) simulation of the TAH using a commercial partitioned block-Gauss–Seidel coupling package. Partitioned coupling of the incompressible fluid with the slender flexible membrane as well as a high fluid/structure density ratio of about unity led inherently to a deterioration of the stability (‘artificial added mass instability’). The objective was to conduct a stable simulation with high accuracy of the pumping process. In order to achieve stability, a combined resistance and pressure outlet boundary condition as well as the interface artificial compressibility method was applied. An analysis of the contact algorithm and turbulence condition is presented. Independence tests are performed for the structural and the fluid mesh, the time step size and the number of pulse cycles. Because of the large deformation of the fluid domain, a variable mesh stiffness depending on certain mesh properties was specified for the fluid elements. Adaptive remeshing was avoided. Different approaches for the mesh stiffness function are compared with respect to convergence, preservation of mesh topology and mesh quality. The resulting mesh aspect ratios, mesh expansion factors and mesh orthogonalities are evaluated in detail. The membrane motion and flow distribution of the coupled simulations are compared with a top-view recording and stereo Particle Image Velocimetry (PIV) measurements, respectively, of the actual pump.

[1]  Katsutoshi Sato,et al.  Analysis of flow patterns in a ventricular assist device: a comparative study of particle image velocimetry and computational fluid dynamics. , 2009, Artificial organs.

[2]  Carsten Lund,et al.  Estimating flow distributions from sampled flow statistics , 2003, SIGCOMM '03.

[3]  Toshiaki Hisada,et al.  Analysis of fluid-structure interaction problems with structural buckling and large domain changes by ALE finite element method , 2001 .

[4]  Yang Liu,et al.  Unsteady flow and diaphragm motion in total artificial heart , 2007 .

[5]  D. Steinman,et al.  Simulation of non-Newtonian blood flow in an end-to-side anastomosis. , 1994, Biorheology.

[6]  K. Hameyer,et al.  New Linear Motor Concepts for Artificial Hearts , 2008, IEEE Transactions on Magnetics.

[7]  Mikko Lyly,et al.  FLUID-STRUCTURE INTERACTION BOUNDARY CONDITIONS BY ARTIFICIAL COMPRESSIBILITY , 2001 .

[8]  D Kwak,et al.  Computational approach for probing the flow through artificial heart devices. , 1997, Journal of biomechanical engineering.

[9]  S Najarian,et al.  SIMULATION OF BLOOD FLOW IN THE SAC-TYPE VENTRICULAR ASSIST DEVICE USING COMPUTATIONAL FLUID DYNAMIC , 2005 .

[10]  Thomas J. R. Hughes,et al.  Fast projection algorithm for unstructured meshes , 1992 .

[11]  T. Laursen Computational Contact and Impact Mechanics , 2003 .

[12]  Thomas J. R. Hughes,et al.  Encyclopedia of computational mechanics , 2004 .

[13]  Shmuel Einav,et al.  Numerical model of flow in a sac-type ventricular assist device. , 2006, Artificial organs.

[14]  K. Bathe Finite Element Procedures , 1995 .

[15]  G Rosenberg,et al.  Fluid dynamics of a pediatric ventricular assist device. , 2000, Artificial organs.

[16]  Byoung Goo Min,et al.  Numerical analysis of the three-dimensional blood flow in the korean artificial heart. , 2003, Artificial organs.

[17]  Steven Deutsch,et al.  Diaphragm motion affects flow patterns in an artificial heart. , 2003, Artificial organs.

[18]  F. Menter Two-equation eddy-viscosity turbulence models for engineering applications , 1994 .

[19]  Charbel Farhat,et al.  Partitioned analysis of coupled mechanical systems , 2001 .

[20]  Emmanuel Guilmineau,et al.  Numerical simulation with a DES approach for automotive flows , 2011 .

[21]  Zhi Tao,et al.  Turbulence modeling of deep dynamic stall at relatively low Reynolds number , 2012 .

[22]  Jan Vierendeels,et al.  Simulation of fluid–structure interaction with the interface artificial compressibility method , 2010 .

[23]  Yves Bourgault,et al.  Numerical simulations of blood flow in artificial and natural hearts with fluid-structure interaction. , 2008, Artificial organs.

[24]  Helge I. Andersson,et al.  Experimental study on the wake behind tapered circular cylinders , 2011 .

[25]  H. Saunders Book Reviews : NUMERICAL METHODS IN FINITE ELEMENT ANALYSIS K.-J. Bathe and E.L. Wilson Prentice-Hall, Inc, Englewood Cliffs, NJ , 1978 .

[26]  E. Ramm,et al.  Artificial added mass instabilities in sequential staggered coupling of nonlinear structures and incompressible viscous flows , 2007 .

[27]  Ulrich Steinseifer,et al.  Implementation of cerebral autoregulation into computational fluid dynamics studies of cardiopulmonary bypass. , 2012, Artificial organs.

[28]  Ali Cemal Benim,et al.  Simulation of blood flow in human aorta with emphasis on outlet boundary conditions , 2011 .

[29]  J. C. Simo,et al.  An augmented lagrangian treatment of contact problems involving friction , 1992 .

[30]  Roland Wüchner,et al.  A framework for stabilized partitioned analysis of thin membrane–wind interaction , 2007 .

[31]  A. Huerta,et al.  Arbitrary Lagrangian–Eulerian Methods , 2004 .

[32]  Timothy J. Barth,et al.  The design and application of upwind schemes on unstructured meshes , 1989 .

[33]  H. Udaykumar,et al.  Application of large-eddy simulation to the study of pulsatile flow in a modeled arterial stenosis. , 2001, Journal of biomechanical engineering.

[34]  Marvin J Slepian,et al.  Cardiac replacement with a total artificial heart as a bridge to transplantation. , 2004, The New England journal of medicine.

[35]  M.H. Moosavi,et al.  Flow Simulation of a Diaphragm-type Ventricular Assist Device with Structural Interactions , 2007, 2007 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[36]  G. Hou,et al.  Numerical Methods for Fluid-Structure Interaction — A Review , 2012 .

[37]  Fabio Nobile,et al.  Added-mass effect in the design of partitioned algorithms for fluid-structure problems , 2005 .

[38]  R D Kamm,et al.  Performance analysis of a cardiac assist device in counterpulsation. , 1998, Journal of biomechanical engineering.

[39]  Eiji Okamoto,et al.  Design of a miniature implantable left ventricular assist device using CAD/CAM technology , 2003, Journal of Artificial Organs.

[40]  C S König,et al.  Investigation of unsteady flow in a model of a ventricular assist device by numerical modelling and comparison with experiment. , 1999, Medical engineering & physics.

[41]  J. C. Simo,et al.  Algorithmic symmetrization of coulomb frictional problems using augmented lagrangians , 1993 .

[42]  B. Norrving,et al.  Global atlas on cardiovascular disease prevention and control. , 2011 .

[43]  Jan Vierendeels,et al.  Stability of a coupling technique for partitioned solvers in FSI applications , 2008 .

[44]  R. Körfer,et al.  The MiniACcor: Constructive Redesign of an Implantable Total Artificial Heart, Initial Laboratory Testing and Further Steps , 2007 .