Simulating the fluid dynamics of natural and prosthetic heart valves using the immersed boundary method

The immersed boundary method is both a general mathematical framework and a particular numerical approach to problems of fluid-structure interaction. In the present work, we describe the application of the immersed boundary method to the simulation of the fluid dynamics of heart valves, including a model of a natural aortic valve and a model of a chorded prosthetic mitral valve. Each valve is mounted in a semi-rigid flow chamber. In the case of the mitral valve, the flow chamber is a circular pipe, and in the case of the aortic valve, the flow chamber is a model of the aortic root. The model valves and flow chambers are immersed in a viscous incompressible fluid, and realistic fluid boundary conditions are prescribed at the upstream and downstream ends of the chambers. To connect the immersed boundary models to the boundaries of the fluid domain, we introduce a novel modification of the standard immersed boundary scheme. In particular, near the outer boundaries of the fluid domain, we modify the construction of the regularized delta function which mediates fluid-structure coupling in the immersed boundary method, whereas in the interior of the fluid domain, we employ a standard four-point delta function which is frequently used with the immersed boundary method. The standard delta function is used wherever possible, and the modified delta function continuously transitions to the standard delta function away from the outer boundaries of the fluid domain. Three-dimensional computational results are presented to demonstrate the capabilities of our immersed boundary approach to simulating the fluid dynamics of heart valves.

[1]  Michael S Sacks,et al.  On the biaxial mechanical properties of the layers of the aortic valve leaflet. , 2005, Journal of biomechanical engineering.

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

[3]  C. Peskin,et al.  Implicit second-order immersed boundary methods with boundary mass , 2008 .

[4]  Phillip Colella,et al.  A cell-centered adaptive projection method for the incompressible Navier-Stokes equations in three dimensions , 2007, J. Comput. Phys..

[5]  Charles S. Peskin,et al.  2-D Parachute Simulation by the Immersed Boundary Method , 2006, SIAM J. Sci. Comput..

[6]  The penalty immersed boundary method and its application to aerodynamics , 2003 .

[7]  Boyce E. Griffith,et al.  Parallel and Adaptive Simulation of Cardiac Fluid Dynamics , 2009 .

[8]  N. Stergiopulos,et al.  Total arterial inertance as the fourth element of the windkessel model. , 1999, American journal of physiology. Heart and circulatory physiology.

[9]  C. Peskin,et al.  A three-dimensional computer model of the human heart for studying cardiac fluid dynamics , 2000, SIGGRAPH 2000.

[10]  Xiaoyu Luo,et al.  Effect of ventricle motion on the dynamic behaviour of chorded mitral valves , 2008 .

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

[12]  A. Chorin.  Numerical Solution of the Navier-Stokes Equations* , 1989 .

[13]  Joao A. C. Lima,et al.  Transesophageal magnetic resonance imaging of the aortic arch and descending thoracic aorta in patients with aortic atherosclerosis. , 2001, Journal of the American College of Cardiology.

[14]  E. Lansac,et al.  Aortic and pulmonary root: are their dynamics similar? , 2002, European journal of cardio-thoracic surgery : official journal of the European Association for Cardio-thoracic Surgery.

[15]  Aaron L. Fogelson,et al.  Stability of approximate projection methods on cell-centered grids , 2005 .

[16]  C. Peskin,et al.  Fluid Dynamics of the Heart and its Valves , 1996 .

[17]  Boyce E. Griffith,et al.  Simulating the blood-muscle-valve mechanics of the heart by an adaptive and parallel version of the immersed boundary method , 2005 .

[18]  Yongsam Kim,et al.  Penalty immersed boundary method for an elastic boundary with mass , 2007 .

[19]  F. Baaijens,et al.  Collagen fibers reduce stresses and stabilize motion of aortic valve leaflets during systole. , 2004, Journal of biomechanics.

[20]  P. Colella,et al.  A Conservative Adaptive Projection Method for the Variable Density Incompressible Navier-Stokes Equations , 1998 .

[21]  M. Minion,et al.  Accurate projection methods for the incompressible Navier—Stokes equations , 2001 .

[22]  Charles S. Peskin,et al.  Flow patterns around heart valves: a digital computer method for solving the equations of motion , 1973 .

[23]  Daniel F. Martin,et al.  A Cell-Centered Adaptive Projection Method for the Incompressible Euler Equations , 2000 .

[24]  P. Colella,et al.  A second-order projection method for the incompressible navier-stokes equations , 1989 .

[25]  C. Peskin,et al.  Heart Simulation by an Immersed Boundary Method with Formal Second-order Accuracy and Reduced Numerical Viscosity , 2001 .

[26]  M. Sacks,et al.  Time-dependent biaxial mechanical behavior of the aortic heart valve leaflet. , 2007, Journal of biomechanics.

[27]  M. Berger,et al.  An Adaptive Version of the Immersed Boundary Method , 1999 .

[28]  Mei-Lin Lai,et al.  A Projection Method for Reacting Flow in the Zero Mach Number Limit , 1994 .

[29]  Alexandre Joel Chorin,et al.  On the Convergence of Discrete Approximations to the Navier-Stokes Equations* , 1989 .

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

[31]  D J Wheatley,et al.  Dynamic modelling of prosthetic chorded mitral valves using the immersed boundary method. , 2007, Journal of biomechanics.

[32]  John B. Bell,et al.  A Numerical Method for the Incompressible Navier-Stokes Equations Based on an Approximate Projection , 1996, SIAM J. Sci. Comput..

[33]  William J. Rider,et al.  Accurate monotonicity- and extrema-preserving methods through adaptive nonlinear hybridizations , 2007, J. Comput. Phys..

[34]  Karim Azer,et al.  A One-dimensional Model of Blood Flow in Arteries with Friction and Convection Based on the Womersley Velocity Profile , 2007, Cardiovascular engineering.

[35]  Robert Michael Kirby,et al.  Unconditionally stable discretizations of the immersed boundary equations , 2007, J. Comput. Phys..

[36]  Andrea Prosperetti,et al.  A second-order boundary-fitted projection method for free-surface flow computations , 2006, J. Comput. Phys..

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

[38]  A. Guyton,et al.  Textbook of Medical Physiology , 1961 .

[39]  Charles S. Peskin,et al.  Shared-Memory Parallel Vector Implementation of the Immersed Boundary Method for the Computation of Blood Flow in the Beating Mammalian Heart , 2004, The Journal of Supercomputing.

[40]  Daniel J. Bodony,et al.  Analysis of sponge zones for computational fluid mechanics , 2006, J. Comput. Phys..

[41]  D. J. Hart.  Fluid-structure interaction in the aortic heart valve : a three-dimensional computational analysis , 2002 .

[42]  Alexandre J. Chorin,et al.  On the Convergence of Discrete Approximations to the Navier-Stokes Equations , 1969 .

[43]  Aahj Fons Sauren.  The mechanical behaviour of the aortic valve , 1981 .

[44]  C. Peskin,et al.  Simulation of a Flapping Flexible Filament in a Flowing Soap Film by the Immersed Boundary Method , 2002 .

[45]  Boyce E. Griffith,et al.  An adaptive, formally second order accurate version of the immersed boundary method , 2007, J. Comput. Phys..

[46]  John B. Bell,et al.  Approximate Projection Methods: Part I. Inviscid Analysis , 2000, SIAM J. Sci. Comput..

[47]  B. Wetton,et al.  Analysis and computation of immersed boundaries, with application to pulp fibres , 1997 .

[48]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

[49]  W. M. Swanson,et al.  Dimensions and Geometric Relationships of the Human Aortic Value as a Function of Pressure , 1974, Circulation research.

[50]  Robert M. Kirby,et al.  A Comparison of Implicit Solvers for the Immersed Boundary Equations , 2008 .

[51]  C. Peskin,et al.  Interaction of two flapping filaments in a flowing soap film , 2003 .

[52]  J. Szmelter.  Incompressible flow and the finite element method , 2001 .

[53]  F P T Baaijens,et al.  A computational fluid-structure interaction analysis of a fiber-reinforced stentless aortic valve. , 2003, Journal of biomechanics.

[54]  Boyce E. Griffith,et al.  An accurate and efficient method for the incompressible Navier-Stokes equations using the projection method as a preconditioner , 2009, J. Comput. Phys..

[55]  Boyce E. Griffith,et al.  On the order of accuracy of the immersed boundary method: Higher order convergence rates for sufficiently smooth problems , 2005 .

[56]  R. Sani,et al.  Incompressible Flow and the Finite Element Method, Volume 1, Advection-Diffusion and Isothermal Laminar Flow , 1998 .

[57]  Yongsam Kim,et al.  On various techniques for computer simulation of boundaries with mass , 2003 .

[58]  N. Westerhof,et al.  Aortic Input Impedance in Normal Man: Relationship to Pressure Wave Forms , 1980, Circulation.

[59]  T. Schmitz-Rode,et al.  The geometry of the aortic root in health, at valve disease and after valve replacement. , 1990, Journal of biomechanics.

[60]  Charles S. Peskin,et al.  Dynamics of a Closed Rod with Twist and Bend in Fluid , 2008, SIAM J. Sci. Comput..

[61]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[62]  A. Chorin.  Numerical solution of the Navier-Stokes equations , 1968 .

[63]  Michael S Sacks,et al.  The flexural rigidity of the aortic valve leaflet in the commissural region. , 2006, Journal of biomechanics.

[64]  David Saloner,et al.  Asymmetric mechanical properties of porcine aortic sinuses. , 2008, The Annals of thoracic surgery.