Authors: Roy C P Kerckhoffs Maxwell L Neal Quan Gu James B Bassingthwaighte Jeff H Omens Andrew D McCulloch
Publish Date: 2006/11/08
Volume: 35, Issue: 1, Pages: 1-18
Abstract
In this study we present a novel robust method to couple finite element FE models of cardiac mechanics to systems models of the circulation CIRC independent of cardiac phase For each time step through a cardiac cycle left and right ventricular pressures were calculated using ventricular compliances from the FE and CIRC models These pressures served as boundary conditions in the FE and CIRC models In succeeding steps pressures were updated to minimize cavity volume error FE minus CIRC volume using Newton iterations Coupling was achieved when a predefined criterion for the volume error was satisfied Initial conditions for the multiscale model were obtained by replacing the FE model with a varying elastance model which takes into account direct ventricular interactions Applying the coupling a novel multiscale model of the canine cardiovascular system was developed Global hemodynamics and regional mechanics were calculated for multiple beats in two separate simulations with a left ventricular ischemic region and pulmonary artery constriction respectively After the interventions global hemodynamics changed due to direct and indirect ventricular interactions in agreement with previously published experimental results The coupling method allows for simulations of multiple cardiac cycles for normal and pathophysiology encompassing levels from cell to systemThis work was supported by the National Biomedical Computation Resource NIH Grant P41 RR08605 to ADM National Science Foundation Grants BES0096492 and BES0506252 to ADM and BES0506477 to MLN NIH Grant HL32583 to JHO and NIH Grant EB001973 to JBB This investigation was conducted in a facility constructed with support from Research Facilities Improvement Program Grant Number C06 RR01758801 from the National Center for Research Resources National Institutes of Health ADM and JHO are cofounders of Insilicomed Inc a licensee of UCSDowned software used in this research Furthermore we are grateful to our programmers Sherief AbdelRahman Ryan Brown and Fred Lionetti for their excellent work on improving and extending ContinuityFrom the pressure and volume curves in Fig 2 it can be seen that ventricular cocompliances are pressuredependent ie the RV volume change for a LV pressure change is different at a constant low and high RV pressure Hence C LR and C RL in Eq A1 are written as a function of pressure see also Table A1
Keywords: