Application of Parallel Processing to the Simulation of Heart Mechanics

Simulations of the mechanics of the left ventricle of the heart with fluid-structure interaction benefit greatly from the parallel processing power of a high performance computing cluster, such as HPCVL. The objective of this paper is to describe the computational requirements for our simulations. Results of parallelization studies show that, as expected, increasing the number of threads per job reduces the total wall clock time for the simulations. Further, the speed-up factor increases with increasing problem size. Comparative simulations with different computational meshes and time steps show that our numerical solutions are nearly independent of the mesh density in the solid wall (myocardium) and the time step duration. The results of these tests allow our simulations to continue with the confidence that we are optimizing our computational resources while minimizing errors due to choices in spatial or temporal resolution.

[1]  A D McCulloch,et al.  Transmural distribution of three-dimensional strain in the isolated arrested canine left ventricle. , 1991, The American journal of physiology.

[2]  D. D. Streeter,et al.  Engineering Mechanics for Successive States in Canine Left Ventricular Myocardium: II. Fiber Angle and Sarcomere Length , 1973, Circulation research.

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

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

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

[6]  J. Holmes Determinants of left ventricular shape change during filling. , 2004, Journal of biomechanical engineering.

[7]  P. Hunter,et al.  Computational Mechanics of the Heart , 2000 .

[8]  H N Sabbah,et al.  Pressure‐Diameter Relations during Early Diastole in Dogs: Incompatibility with the Concept of Passive Left Ventricular Filling , 1981, Circulation research.

[9]  Miguel A. Fernández,et al.  Numerical Simulation of the Electromechanical Activity of the Heart , 2009, FIMH.

[10]  P. Hunter,et al.  Mathematical model of geometry and fibrous structure of the heart. , 1991, The American journal of physiology.

[11]  Toshiaki Hisada,et al.  Multiphysics simulation of left ventricular filling dynamics using fluid-structure interaction finite element method. , 2004, Biophysical journal.

[12]  D. D. Streeter Engineering Mechanics for Successive States in Canine Left Ventricular Myocardium: I. CAVITY AND WALL GEOMETRY , 1973, Circulation research.

[13]  Hervé Delingette,et al.  Functional Imaging and Modeling of the Heart, 5th International Conference, FIMH 2009, Nice, France, June 3-5, 2009. Proceedings , 2009, Functional Imaging and Modeling of the Heart.

[14]  J. Ross,et al.  Fiber Orientation in the Canine Left Ventricle during Diastole and Systole , 1969, Circulation research.

[15]  Hou Zhang,et al.  Direct and iterative computing of fluid flows fully coupled with structures , 2001 .

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

[17]  D N Firmin,et al.  Subject-specific computational simulation of left ventricular flow based on magnetic resonance imaging , 2008, Proceedings of the Institution of Mechanical Engineers. Part H, Journal of engineering in medicine.

[18]  R. Judd,et al.  Compressibility of perfused passive myocardium. , 1996, The American journal of physiology.

[19]  Y. Bourgault,et al.  Adaptation of a rabbit myocardium material model for use in a canine left ventricle simulation study. , 2010, Journal of biomechanical engineering.

[20]  Liang Ge,et al.  The use of computational flow modeling (CFD) to determine the effect of left ventricular shape on blood flow in the left ventricle. , 2009, The Annals of thoracic surgery.

[21]  F. Yin,et al.  A multiaxial constitutive law for mammalian left ventricular myocardium in steady-state barium contracture or tetanus. , 1998, Journal of biomechanical engineering.

[22]  D. H. van Campen,et al.  A two-phase finite element model of the diastolic left ventricle. , 1991, Journal of biomechanics.

[23]  P. Hunter,et al.  Computational mechanics of the heart : From tissue structure to ventricular function , 2000 .