Document Detail

Three-wall segment (TriSeg) model describing mechanics and hemodynamics of ventricular interaction.
Jump to Full Text
MedLine Citation:
PMID:  19718527     Owner:  NLM     Status:  MEDLINE    
Abstract/OtherAbstract:
A mathematical model (TriSeg model) of ventricular mechanics incorporating mechanical interaction of the left and right ventricular free walls and the interventricular septum is presented. Global left and right ventricular pump mechanics were related to representative myofiber mechanics in the three ventricular walls, satisfying the principle of conservation of energy. The walls were mechanically coupled satisfying tensile force equilibrium in the junction. Wall sizes and masses were rendered by adaptation to normalize mechanical myofiber load to physiological standard levels. The TriSeg model was implemented in the previously published lumped closed-loop CircAdapt model of heart and circulation. Simulation results of cardiac mechanics and hemodynamics during normal ventricular loading, acute pulmonary hypertension, and chronic pulmonary hypertension (including load adaptation) agreed with clinical data as obtained in healthy volunteers and pulmonary hypertension patients. In chronic pulmonary hypertension, the model predicted right ventricular free wall hypertrophy, increased systolic pulmonary flow acceleration, and increased right ventricular isovolumic contraction and relaxation times. Furthermore, septal curvature decreased linearly with its transmural pressure difference. In conclusion, the TriSeg model enables realistic simulation of ventricular mechanics including interaction between left and right ventricular pump mechanics, dynamics of septal geometry, and myofiber mechanics in the three ventricular walls.
Authors:
Joost Lumens; Tammo Delhaas; Borut Kirn; Theo Arts
Related Documents :
10496447 - Relation between negative u waves in precordial leads on the admission electrocardiogra...
3875437 - Non-q wave perioperative myocardial infarction: assessment of the incidence and severit...
22141747 - Right-left digit ratio (2d:4d) and maximal oxygen uptake.
223427 - Two dimensional echocardiographic quantification of infarct size alteration by pharmaco...
24490007 - Clinical manifestations of right ventricle involvement in inferior myocardial infarction.
3262487 - Interleukin-2 administration causes reversible hemodynamic changes and left ventricular...
Publication Detail:
Type:  Journal Article; Research Support, Non-U.S. Gov't     Date:  2009-08-29
Journal Detail:
Title:  Annals of biomedical engineering     Volume:  37     ISSN:  1573-9686     ISO Abbreviation:  Ann Biomed Eng     Publication Date:  2009 Nov 
Date Detail:
Created Date:  2010-01-29     Completed Date:  2010-03-09     Revised Date:  2013-05-31    
Medline Journal Info:
Nlm Unique ID:  0361512     Medline TA:  Ann Biomed Eng     Country:  United States    
Other Details:
Languages:  eng     Pagination:  2234-55     Citation Subset:  IM    
Affiliation:
Department of Physiology, Cardiovascular Research Institute Maastricht, Maastricht University, PO Box 616, 6200 MD, Maastricht, The Netherlands. j.lumens@fys.unimaas.nl
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Descriptor/Qualifier:
Heart Septum / physiology*
Heart Ventricles / ultrasonography
Hemodynamics / physiology*
Humans
Hypertension / physiopathology
Models, Cardiovascular
Models, Theoretical
Ventricular Function, Left / physiology*
Ventricular Function, Right / physiology*
Comments/Corrections

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): Ann Biomed Eng
ISSN: 0090-6964
ISSN: 1521-6047
Publisher: Springer US, Boston
Article Information
Download PDF
? The Author(s) 2009
Received Day: 26 Month: 11 Year: 2008
Accepted Day: 30 Month: 7 Year: 2009
Electronic publication date: Day: 29 Month: 8 Year: 2009
Print publication date: Month: 11 Year: 2009
Volume: 37 Issue: 11
First Page: 2234 Last Page: 2255
ID: 2758607
PubMed Id: 19718527
Publisher Id: 9774
DOI: 10.1007/s10439-009-9774-2
issue-copyright-statement: ? Biomedical Engineering Society 2009

Three-Wall Segment (TriSeg) Model Describing Mechanics and Hemodynamics of Ventricular Interaction
Joost Lumens12 Address: +31-43-3881658 +31-43-3884166 j.lumens@fys.unimaas.nl
Tammo Delhaas1
Borut Kirn23
Theo Arts2
1Department of Physiology, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, PO Box 616, 6200 MD Maastricht, The Netherlands
2Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, The Netherlands
3Institute of Physiology, School of Medicine, University of Ljubljana, Ljubljana, Slovenia

Introduction

The left (LV) and right (RV) ventricular cavities of the normal human heart are separated from the intrathoracic space by the left and right ventricular free walls, respectively (Fig.?1a). The cavities are mutually separated by the septal wall. In order to simulate LV and RV mechanics and hemodynamics, the close anatomic coupling between these two ventricles must be considered. Experimental data demonstrate that a change in loading condition of either ventricle directly influences pump function of the other ventricle.5,21,50,55,68 It is also well recognized that septal geometry and motion depend on transseptal pressure difference.36,37,51 This dependency is best illustrated by flattening of the septum with increase of RV pressure relative to LV pressure (Fig.?1b). Septal geometry and motion appeared to be of diagnostic and predictive significance in patients with pulmonary hypertension.32,35

Several mathematical models of ventricular mechanics have been developed to quantify the effect of ventricular interaction on cardiac function. In several lumped models of ventricular hemodynamics,7,13,41,43,48,53,56,60,61 ventricular interaction is described by empirically determined coupling coefficients, quantifying interventricular cross-talk of pressures and volumes. The model designed by Beyar et al.7 describes ventricular interaction in a more mechanistic way. In this model, LV and RV cavities are enclosed by three ventricular walls. For given ventricular pressures, the mechanical equilibrium of tensile forces in the junction of the walls is used as constraint to predict ventricular geometry. The model is restricted to description of passive mechanics of walls lacking contractile myofiber properties. More recently, 3D finite element models of the cardiac ventricles were used to simulate ventricular pump function and local tissue mechanics,34,47,66 inherently including ventricular interaction via the septum. In comparison with lumped models of global ventricular mechanics, finite element models allow description of regional wall mechanics and geometry. Consequently, these models are computationally demanding.

For study of beat-to-beat hemodynamics and mechanics of heart and blood vessels, Arts et al. previously developed the closed-loop CircAdapt model of heart and circulation.3 In this model, mechanical interaction of the LV and RV has been simulated by a common outer wall, having a transmural pressure equal to RV pressure, encapsulating an inner wall, which represents the left ventricle.1 The inner wall encapsulates the LV cavity only and has a transmural pressure equal to the difference between LV and RV pressure. Under normal ventricular loading conditions, this model enables realistic simulation of global LV and RV pump mechanics.1,3 However, this model setup presumes RV pressure to be substantially lower than LV pressure during the whole cardiac cycle. This condition is not satisfied with pulmonary hypertension or with left-to-right asynchrony of electrical activation. Furthermore, septal geometry cannot be simulated. Thus, this model cannot be used to interpret this measurable signal that contains important information about the difference between LV and RV pressure.35,51 Therefore, we designed the TriSeg model of ventricular mechanics that realistically incorporates ventricular interaction via the interventricular septum.

In the TriSeg model, three wall segments, i.e., the LV free wall (LW), septal wall (SW), and RV free wall (RW), meet in a junction margin, thus encapsulating the LV and RV cavities (Fig.?2a). LW and RW geometries are calculated starting from an initial estimate of SW geometry. For each wall segment, myofiber strain is calculated from midwall surface area relative to a reference area. From myofiber strain, myofiber stress is determined using constitutive equations describing sarcomere mechanics. Substitution of myofiber stress in the known geometry of the walls results in representative midwall tension, acting on the junction line. Next, septal geometry is adjusted so that equilibrium of tension is achieved, i.e., summed tension at the junction of the three walls is equal to zero. As a result, geometries of the three wall segments and the two cavities are known, together with stresses, tensions, and pressures.

The TriSeg model of ventricular mechanics is designed to be incorporated as a module in the existing CircAdapt model simulating mechanics and hemodynamics of the whole circulation.3 According to the principles of CircAdapt, size and mass of the wall segments are determined by adaptation so that mechanical load of the myofibers is normalized to physiological standard levels.3,4 The TriSeg model, as integrated in the CircAdapt model, has been tested by simulation of time-dependent LV and RV mechanics and hemodynamics under normal ventricular loading conditions as well as with acute and chronic pulmonary hypertension (PH). PH has been simulated to test whether the TriSeg model realistically relates septal geometry to transseptal pressure. Simulation results have been compared with previously published experimental data on ventricular hemodynamics9,62 and on the relation between septal geometry and transmural pressure17,51 in healthy volunteers and in patients with chronic PH.


Methods
General Design of the TriSeg Model

The design of the new TriSeg model of ventricular mechanics should enable incorporation as a module in the existing CircAdapt model that simulates mechanics and hemodynamics of the whole circulation.3 Therefore, the TriSeg model calculates LV and RV pressures (pLV and pRV, respectively) as functions of LV and RV cavity volumes (VLV and VRV, respectively). Mechanics of ventricular interaction is incorporated assuming a simplified ventricular composite geometry (Fig.?2a). The three ventricular walls LW, SW, and RW are modeled to be thick-walled and spherical with a common junction circle with midwall radius ym (Fig.?2b). Midwall is defined as the spherical surface that divides wall volume Vw in an inner and an outer shell of equal volume. Midwall volume Vm is the volume enclosed by the midwall surface and the plane of the junction circle (Fig.?2c). The center of the junction circle is the origin of the applied cylindrical coordinate system (Fig.?2b). The x-direction is perpendicular to the plane of the junction circle and is defined positive toward the RV free wall. The y-coordinate represents the radial distance to the x-axis. The spherical midwall surface of a wall segment intersects the x-axis at value xm (Fig.?2c).

In presenting the TriSeg model, the following levels are distinguished: (1) ventricular hemodynamics, (2) ventricular composite mechanics, (3) ventricular wall segment mechanics, (4) curved wall patch mechanics, and (5) myofiber mechanics. Starting from LV and RV volumes, the sequence of calculations needed to obtain LV and RV pressures has been illustrated in Fig.?3. The ventricular hemodynamics section serves as interface with the CircAdapt model of the whole circulation.3 In the ventricular composite section, LV and RV volumes are used to calculate LV and RV pressures, applying a model of wall mechanics to each ventricular wall segment. In the ventricular wall segment section, volume Vm and boundary radius ym of a spherical wall segment are used to calculate the axial and radial vector components Tx and Ty, respectively, at the common junction circle from representative midwall tension Tm, applying a model of the mechanics of a curved wall patch. Such patch is defined as a fraction of a spherical wall segment with an arbitrarily shaped boundary and with midwall area Am and curvature Cm. Curvature is defined as the reciprocal of radius of curvature, having the advantage of being well defined for a flat wall. In the section of curved wall patch mechanics, Am and Cm are used to calculate representative midwall tension Tm, applying a model of myofiber mechanics. In the myofiber mechanics section, natural myofiber strain ?f is converted to Cauchy myofiber stress ?f, using a constitutive law of the myofiber. This model of myofiber mechanics has been heuristically obtained from reported physiological experiments.

Ventricular Composite Mechanics

In this section, it is shown how cavity volumes VLV and VRV are used to calculate pressures pLV and pRV. For that purpose, a model of ventricular wall segment mechanics is applied that renders axial and radial tension components (Tx and Ty, respectively) from representative midwall tension as a function of the midwall volume (Vm) and junction radius (ym) of the spherical wall segment. This latter model will be derived in the next section discussing mechanics of the ventricular wall segment.

The geometry of the ventricular composite (Figs.?2b and 2c) is defined by the common radius ym of the wall junction and the enclosed midwall cap volumes Vm,LW, Vm,SW, and Vm,RW of the wall segments LW, SW, and RW, respectively. Septal midwall volume Vm,SW and junction radius ym are initially estimated by the solutions as obtained in the preceding time point, i.e., Vm,SW,est and ym,est, respectively:

[Formula ID: Equ1]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ y_{\text{m}} = y_{\text{m,est}} $$\end{document}]
[Formula ID: Equ2]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ V_{\text{m,SW}} = V_{\text{m,SW,est}} $$\end{document}]
The latter estimates, together with ventricular cavity volumes VLV and VRV as given directly by the CircAdapt model, render an estimate of ventricular composite geometry. Using the fact that ventricular cavity volume added to half of wall volume Vw of the enclosing wall segments renders the sum of midwall volumes, it is found:
[Formula ID: Equ3]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ V_{\text{m,LW}} = - V_{\text{LV}} - {\frac{1}{2}}V_{{{\text{w}},{\text{LW}}}} - {\frac{1}{2}}V_{{{\text{w}},{\text{SW}}}} + V_{{{\text{m}},{\text{SW}}}} $$\end{document}]
[Formula ID: Equ4]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ V_{{{\text{m}},{\text{RW}}}} = + V_{\text{RV}} + {\frac{1}{2}}V_{{{\text{w}},{\text{RW}}}} + {\frac{1}{2}}V_{{{\text{w}},{\text{SW}}}} + V_{{{\text{m}},{\text{SW}}}} $$\end{document}]
Note that the sign of midwall volume Vm is positive if wall curvature is convex to the positive x-direction (Fig.?2b). Thus, for a normal heart, the sign is negative for LW and positive for RW and SW.

For each wall segment i, Vm,i and ym,i are used to calculate axial and radial tension components Tx,i and Ty,i, as well as transmural pressure pTrans,i. This is done by applying a model of ventricular wall segment mechanics as discussed in the next section. The axial tension components of the three wall segments Tx,LW, Tx,SW, and Tx,RW are summed to calculate the net axial tension Tx,Tot in the junction. Similarly, the net radial tension Ty,Tot in the junction is obtained. The initial estimates of ym and Vm,SW (Eqs.?1 and 2) are adjusted numerically so that for the summed tension components it holds:

[Formula ID: Equ5]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ T_{\text{x,Tot}} (y_{\text{m}} ,V_{{{\text{m}},{\text{SW}}}} ) = 0 $$\end{document}]
[Formula ID: Equ6]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ T_{{{\text{y}},{\text{Tot}}}} (y_{\text{m}} ,V_{{{\text{m}},{\text{SW}}}} ) = 0 $$\end{document}]
Having solved Eqs. (5) and (6), and assuming the pressure surrounding the ventricular composite to be zero, for LV and RV pressures it holds:
[Formula ID: Equ7]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ p_{\text{LV}} = - p_{{{\text{Trans}},{\text{LW}}}} $$\end{document}]
[Formula ID: Equ8]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ p_{\text{RV}} = + p_{{{\text{Trans}},{\text{RW}}}} $$\end{document}]
Note that the minus sign for pLV is caused by the definition of pTrans, which is positive for a negative x-gradient of pressure in the wall segment.

Mechanics of the Ventricular Wall Segment

In this section, it is shown for a wall segment (Fig.?2c) how Vm and ym are used to calculate axial and radial midwall tension components Tx and Ty, and transmural pressure pTrans. For that purpose, a model of the mechanics of a curved wall patch is applied. Such curved wall patch is defined as a fraction of a spherical wall segment with midwall surface area Am and curvature Cm. This latter model renders representative midwall tension from midwall surface area Am and curvature Cm and will be derived in the next section discussing mechanics of the curved wall patch.

Wall segment geometry is most conveniently described by radius ym of the junction circle and distance xm between plane of junction circle and center of midwall surface (Fig.?2c). Vm, Am, and Cm depend on xm and ym by

[Formula ID: Equ9]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ V_{\text{m}} = {\frac{\pi }{6}}x_{\text{m}} (x_{\text{m}}^{2} + 3y_{\text{m}}^{2} ) $$\end{document}]
[Formula ID: Equ10]
10  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ A_{\text{m}} = \pi (x_{\text{m}}^{2} + y_{\text{m}}^{2} ) $$\end{document}]
[Formula ID: Equ11]
11  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ C_{\text{m}} = {\frac{{2x_{\text{m}} }}{{(x_{\text{m}}^{2} + y_{\text{m}}^{2} )}}} $$\end{document}]
First, xm is calculated from ym and Vm by solving Eq. (9). Second, Am and Cm are determined with Eqs.?(10) and (11). Next, Am and Cm are used to calculate representative midwall tension Tm, applying a model of the mechanics of a curved wall patch as discussed in the next section.

At the circular boundary of the midwall surface (Fig.?2c), geometric parameters xm and ym, and midwall tension Tm are used to calculate axial (Tx) and radial (Ty) tension components:

[Formula ID: Equ12]
12  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ T_{\text{x}} = T_{\text{m}} \cdot \sin \alpha \quad {\text{with}}\;\sin \alpha = {\frac{{2x_{\text{m}} y_{\text{m}} }}{{x_{\text{m}}^{2} + y_{\text{m}}^{2} }}} $$\end{document}]
[Formula ID: Equ13]
13  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ T_{\text{y}} = T_{\text{m}} \cdot \cos \alpha \quad {\text{with}}\;\cos \alpha = {\frac{{ - x_{\text{m}}^{2} + y_{\text{m}}^{2} }}{{x_{\text{m}}^{2} + y_{\text{m}}^{2} }}} $$\end{document}]
The symbol ? represents half opening angle of the circular wall segment (Fig.?2c).

Transmural pressure pTrans equals total axial force. Consequently, pTrans is calculated by multiplying Tx by the ratio of length of the junction contour (2?ym) to area of the junction circle (?ym2):

[Formula ID: Equ14]
14  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ p_{\text{Trans}} = {\frac{{2T_{\text{x}} }}{{y_{\text{m}} }}} $$\end{document}]

Mechanics of the Curved Wall Patch

In this section, it is shown how midwall surface area Am and curvature Cm are used to calculate representative midwall tension Tm. For that purpose, a constitutive model of myofiber mechanics is applied. This latter model renders Cauchy myofiber stress from natural myofiber strain and is presented in Appendix B.

In the ventricular module, originally implemented in the CircAdapt model,3 the one-fiber model developed by Arts et al.2 was used to relate ventricular pump mechanics, as described by cavity pressure and volume, to myofiber mechanics, as described by myofiber stress and strain. When assembling three wall segments to a ventricular composite with two cavities, contractile function of a wall segment, as described by representative midwall tension and area change, should be put between mechanics of cavity and myofiber. The following conditions should be satisfied. First, when folding a wall segment to a completely closed spherical surface, the relation between pump mechanics and myofiber mechanics should be equivalent to the equations of the one-fiber model. Second, like in the one-fiber model, where contractile myofiber work equals ventricular pump work, summed pump work of both cavities should be equal to summed work as generated by the three walls. Furthermore, within each wall segment, work as delivered by the wall through wall tension and changes of geometry, should be equal to the work generated by the myofibers.

In Appendix A, it is shown that such a model of wall segment mechanics can be found, although the derivation is quite complex. Because the shape of the midwall area is not relevant for the balance of work, the relations are extended to the more general case of a spherical wall patch with arbitrary shape of the boundary. Thus, for a wall patch it is found that natural myofiber strain ?f depends on midwall surface area Am and curvature Cm by the following approximation (relative error <1%), as derived in Appendix A:

[Formula ID: Equ15]
15  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \varepsilon_{\text{f}} \approx {\frac{1}{2}}\ln \left( {{\frac{{A_{\text{m}} }}{{A_{\text{m,ref}} }}}} \right) - {\frac{1}{12}}z^{2} - 0.019z^{4} \quad {\text{ with}}\quad \, z = {\frac{{3C_{\text{m}} V_{\text{w}} }}{{2A_{\text{m}} }}} $$\end{document}]
where Am,ref represents reference midwall surface area. Note that the dimensionless curvature parameter z is closely related to the dimensionless ratio of wall thickness to radius of curvature. With a constitutive model of the myofiber (Appendix B), natural myofiber strain ?f is used to calculate Cauchy myofiber stress ?f. In Appendix A, it is derived how myofiber stress and wall segment geometry are used to calculate representative midwall tension Tm. It holds by approximation (relative error <2%):
[Formula ID: Equ16]
16  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ T_{\text{m}} \approx {\frac{{V_{\text{w}} \sigma_{\text{f}} }}{{2A_{\text{m}} }}}\left( {1 + {\frac{{z^{2} }}{3}} + {\frac{{z^{4} }}{5}}} \right)\quad {\text{with}}\quad \, \sigma_{\text{f}} = f(\varepsilon_{\text{f}} ) $$\end{document}]
Equations (15) and (16) represent the first few terms of a Taylor series approximating the analytically derived expressions. These approximations have the advantage of being well defined for zero curvature. The analytical derivations cannot be used in numerical calculations, because zero by zero division occurs near zero curvature, while the ratio is analytically well defined.

Implementation of TriSeg Model in CircAdapt Model

The TriSeg model of ventricular mechanics is incorporated in the existing CircAdapt model of the whole circulation.3 This latter model supplies the required hemodynamic boundary conditions, i.e., LV and RV cavity volumes. The CircAdapt model is designed as a network of modules representing cardiac chambers, valves, large blood vessels, and peripheral resistances. Extended with the TriSeg model, the CircAdapt model allows beat-to-beat simulation of time-dependent ventricular mechanics and hemodynamics, e.g., ventricular cavity volumes and pressures, geometries and representative myofiber mechanics of ventricular walls, and flows through valves. An important feature of the existing CircAdapt model is that the number of independent parameters is reduced by incorporating adaptation of cavity size and wall mass of cardiac chambers and blood vessels to mechanical load so that stresses and strains in the walls of heart and blood vessels are normalized to tissue-specific physiological standard levels.

The model was initialized using realistic estimates of wall volume Vw, reference midwall surface area Am,ref, and reference sarcomere length Ls,ref (Table?1: Parameter values for model initialization). After model initialization, these parameters will be obtained for each wall segment by load adaptation as described in the ?Simulations? section. The actual parameter values chosen for initialization (Table?1) are noncritical due to load adaptation. However, for fast convergence to steady-state adaptation, these estimates were chosen in physiological range.

The set of differential equations that describe the CircAdapt model, including the TriSeg module, was solved by numerical integration with time steps of 2?ms using the ODE113 function in MATLAB 7.1.0 (MathWorks, Natick, MA). Simulation time of a single cardiac cycle was less than 7?s on a Windows XP? platform (version 2002) with a 2.00?GHz Intel? Core?2 Duo T7250 processor and 1?GB of RAM.

Simulations

The TriSeg model, as integrated in the CircAdapt model, was tested by simulation of human ventricular mechanics and hemodynamics with (1) normal ventricular loading conditions (NORM), (2) acute pulmonary hypertension (PHAc), and (3) chronic pulmonary hypertension including adaptation (PHCh). All simulations are presented with similar values of mean systemic blood flow (cardiac output), cardiac cycle time (heart rate), and mean systemic arterial blood pressure, simulating hemodynamics at rest (Table?1: Hemodynamics at rest).

NORM Simulation

The diameter of large blood vessels was assumed to be determined by adaptation to mean chronic circumstances, which state is most closely described by the condition of rest (Table?1: Hemodynamics at rest). However, geometry of the heart and wall thickness of large blood vessels were considered to be a result of adaptation to a state of moderate exercise, thus simulating the effect of relatively short periods of training by exercise. The applied simulation protocol of adaptation has been described earlier.3 Briefly recapitulating, with resting hemodynamics (Table?1: Hemodynamics at rest), diameters of the large blood vessels were adapted until mean blood flow velocity reached the set point value (Table?1: Set point values for adaptation). Next, a moderate state of exercise was simulated by tripling cardiac output and doubling heart rate (Table?1: Interventions with exercise). Under these circumstances, wall volume (Vw) and reference midwall surface area (Am,ref) of all cardiac wall segments were adapted until maximum and minimum sarcomere length as well as maximum passive myofiber stress reached the set point values of adaptation with exercise (Table?1). Also, wall thickness of the blood vessels was adapted until maximum wall stress was equal to the set point value (Table?1). Next, hemodynamics were returned to rest conditions and the above-mentioned adaptation protocol was repeated two or three times until steady-state geometry was reached (<1% deviation from set point values). In total, finding the steady-state adapted NORM simulation required simulation of about 100?200 cardiac cycles.

PHAc and PHCh Simulations

The NORM simulation was used as point of departure for simulation of ventricular mechanics and hemodynamics with increased pulmonary resistance. Steady-state simulations (PHAc1, PHAc2, and PHAc3), representing increasing degrees of acute pulmonary hypertension, were obtained by acute increase of mean pulmonary arteriovenous pressure drop in three steps without adaptation (Table?1: Interventions with PH simulations). Finally, adaptation was applied to the three steady-state PHAc simulations following the same protocol used for the NORM simulation. This resulted in three steady-state simulations (PHCh1, PHCh2, and PHCh3, respectively), representing gradually increasing degrees of chronic pulmonary hypertension at rest.

Simulation Data Analysis
Ventricular Pump Mechanics and Hemodynamics

Simulated time courses of normal (NORM) LV and RV pressures, volumes, and flows were compared with physiological data obtained in normal subjects.26,29,54,62,63 Area of ventricular pressure?volume relation was calculated to quantify ventricular pump stroke work (Wstroke) for the LV and RV. Time courses of blood flow velocities through the mitral, aortic, tricuspid, and pulmonary valves were used to study the effects of chronic pulmonary hypertension on ventricular hemodynamics. For comparison with clinical data obtained in healthy subjects and in patients with chronic pulmonary hypertension,9,62 several timing parameters of RV hemodynamics were quantified for the NORM and PHCh simulations. RV ejection time (ET) was quantified as the time from pulmonary valve opening to closure, acceleration time of pulmonary flow (AT) as the time from pulmonary valve opening to moment of maximal pulmonary flow velocity, RV isovolumic contraction time (ICT) as the time from tricuspid valve closure to pulmonary valve opening, and RV isovolumic relaxation time (IRT) as the time from pulmonary valve closure to tricuspid valve opening.

Tissue Mechanics

Area of myofiber stress?strain relation was calculated to quantify stroke work density (wstroke) for the LW, SW, and RW. Stroke work density was defined as contractile myofiber stroke work per unit of tissue volume.

Ventricular Wall Geometry

End-diastolic wall thicknesses (Hed) were calculated as wall volume divided by end-diastolic midwall surface area. End-diastole was defined as the moment of mitral valve closure. Furthermore, time courses of LW, SW, and RW midwall curvatures were used to assess acute and chronic effects of increase of pulmonary resistance on wall curvatures. The relation between septal geometry and the transseptal pressure gradient in the TriSeg model was compared with the relation found in patients with and without chronic pulmonary hypertension.17 For that purpose, curvature ratio (CR) was calculated as the SW/LW midwall curvature ratio at the moment of aortic valve closure, whereas transmural pressure ratio (PR) was defined as the difference between maximum LV and RV pressures divided by maximum LV pressure.


Results
Simulation of Normal Physiology (NORM)

Figure?4 shows simulated time courses of ventricular cavity pressures and volumes, aortic and pulmonary artery pressures, and flows through arterial and atrioventricular valves under normal ventricular loading conditions (NORM). These time courses show agreement with physiological data on the following aspects: (1) LV and RV cavity pressures and volumes rise at the end of diastole as a result of atrial contraction, (2) time courses of mitral and tricuspid valve flows show an early passive filling wave (E) followed by a clearly separated late filling wave (A) resulting from atrial contraction,26 (3) LV and RV end-diastolic volumes are not significantly different,54 (4) RV ejection starts earlier and ends later than LV ejection,29 and (5) durations of LV isovolumic contraction and relaxation exceed the corresponding RV time intervals (Fig.?4, bottom panel).62,63 Steady-state adaptation rendered LW, SW, and RW wall volumes of 76, 38, and 31?mL, respectively. Wall thicknesses Hed,LW and Hed,SW were almost equal and three times as large as Hed,RW (Table?2).

Simulations of Acute Pulmonary Hypertension (PHAc)

The upper panel row of Fig.?5 illustrates simulated ventricular pressure?volume loops under normal loading conditions (NORM) and with acutely increased pulmonary resistance (PHAc1, PHAc2, and PHAc3). With increasing pulmonary resistance, RV cavity volume, systolic and diastolic RV pressures, mean pulmonary artery pressure, and RV pump work (Table?3) increased, while LV cavity volume and pump work were relatively unaffected. The middle panel row of Fig.?5 shows stress?strain loops of the myofibers in the ventricular wall segments. Acute rise of pulmonary resistance caused increase of RW myofiber stress, myofiber strain, and stroke work density, while it caused decrease of LW and SW stroke work densities (Table?3). The lower panel row of Fig.?5 illustrates time courses of LW, SW, and RW midwall curvatures. With increasing pulmonary resistance, SW midwall curvature decreased predominantly during diastole. In the PHAc3 simulation, the septum even flattened to zero curvature and bulged shortly toward the LV free wall (negative curvature) during early diastole. LW and RW curvatures appeared relatively insensitive to acute changes of pulmonary resistance. With increasing pulmonary resistance (Table?2), Hed,RW decreased maximally by 26% in PHAc3 relative to NORM, Hed,LW decreased maximally by 8%, whereas Hed,SW increased only 1%.

Simulations of Chronic Pulmonary Hypertension (PHCh)

Figure?6 illustrates graphs of the same variables as presented in Fig.?5, but now rendered after adaptation to increased pulmonary resistance (PHCh1, PHCh2, and PHCh3). Like in absence of adaptation, systolic and diastolic RV pressures, mean pulmonary artery pressure, and RV pump work were elevated with chronic pulmonary hypertension (Table?3), while LV pump function remained relatively unaffected. In contrast to the PHAc simulations, RV cavity volume decreased 8% in the PHCh3 simulation (Fig.?6). The myofiber stress?strain loops in Fig.?6 show that adaptation resulted in a more homogeneous distribution of stroke work density in the ventricular walls (Table?3). In the PHCh3 simulation, RW stroke work density was increased (20%), whereas LW and SW stroke work densities were decreased (7 and 5%, respectively) relative to the NORM simulation. Increase of pulmonary resistance with load adaptation resulted in the following changes of ventricular geometry. Mean SW curvature (Fig.?6) decreased linearly with increasing pulmonary resistance, while the shape of the time course was preserved. LW and RW curvatures were practically unaffected. Relative to the NORM simulation, Hed,RW increased up to 120% in PHCh3, whereas changes of Hed,LW and Hed,SW were small (<5%) (Table?2).

In Fig.?7, time courses of LV and RV inflow and outflow velocities are shown for the NORM and PHCh simulations. With increasing pulmonary resistance, pulmonary and tricuspid valve flow velocities changed significantly relative to the NORM simulation, whereas aortic and mitral valve flow velocities remained relatively unaffected. Furthermore, in Fig.?7, RV ejection time (ET), acceleration time of pulmonary flow (AT), RV isovolumic contraction time (ICT), and RV isovolumic relaxation time (IRT) are indicated for the NORM simulation. In Fig.?8, changes of these timing parameters are compared with the changes as observed in clinical studies among normal individuals and patients with chronic pulmonary hypertension.9,62 In the simulations, chronic increase of pulmonary resistance caused decrease of AT and AT/ET-ratio, indicating increase of pulmonary flow curve skewness. This is in accordance with measurements as obtained in patients with chronic pulmonary hypertension.9 Furthermore, simulated ICT and IRT lengthened (maximally 163% and 47% in the PHCh3 simulation, respectively). The ratio of total RV isovolumic time (ICT?+?IRT) to ET almost doubled. Again, the simulated changes were similar to the changes as observed in patients with pulmonary hypertension.62 Furthermore, the tricuspid early filling wave narrowed and the tricuspid E/A-ratio of peak atrioventricular flow during early filling (E) to that during atrial contraction (A) decreased from 1.49 in the NORM simulation to unity in the PHCh3 simulation.

Figure?9 shows simulated time courses of SW/LW transmural pressure ratio and SW/LW midwall curvature ratio. SW transmural pressure was quantified as LV cavity pressure minus RV cavity pressure, while LW transmural pressure was quantified as LV cavity pressure. At each time point in the cardiac cycle, both ratios decreased linearly with increase of pulmonary resistance.

Both CR and PR decreased with increase of pulmonary resistance (Table?4). Figure?10 shows the relation between CR and PR, as derived from the NORM and PHCh simulations, as well as the relation found among patients with and without pulmonary hypertension.17 In this latter patient population, pulmonary hypertension was defined as an RV systolic pressure >40?mmHg. Linear regression analysis of data points as derived from the NORM and PHCh simulations revealed a strong correlation (r2?=?1.00, SEE?=?0.008) between CR and PR with similar slope as the relation found in patients (r2?=?0.73, SEE?=?0.044). In the simulations, however, an offset was found, shifting the relation to higher CR values (Fig.?10). It is noted (not shown) that although acute change (?20%) of mean systemic flow or mean arterial blood pressure affected CR and PR, their interdependence followed the linear relation as derived from the chronic simulations within 2% deviation. The same holds for the data points derived from the PHAc1 and PHAc2 simulations, while for the PHAc3 simulation the deviation was more than 10% (Table?4, not shown in Fig.?10).


Discussion

The newly designed TriSeg model of ventricular mechanics incorporates mechanical interaction of the LV free wall, RV free wall, and septal wall resulting in a strong coupling of LV and RV pump mechanics and hemodynamics. The TriSeg model was successfully implemented as a module in the available CircAdapt model, simulating cardiac mechanics and hemodynamics of the closed-loop circulation. Effects of ventricular interaction on cardiac mechanics and hemodynamics were assessed by simulation of pulmonary hypertension in the acute phase as well as in the chronic phase, the latter implying adaptation of ventricular geometry to mechanical load. For chronic pulmonary hypertension, simulated ventricular geometry, hemodynamics, and septal mechanics agreed surprisingly well with corresponding measurements in patients.

Model Assumptions

In the TriSeg model, ventricular geometry was approximated by three thick-walled spherical segments encapsulating the LV and RV cavities (Fig.?2a). In reality, the ventricular cavities are enclosed by truncated ellipsoidal muscular walls and the noncontractile basal sheet with valves.58 The simplification to spherical segments without a noncontractile sheet resulted in a general underestimation of ventricular dimensions. Although the fact that there are indications that bending stiffness of the myocardium may be important for simulation of septal geometry,7,21 we neglected this effect for simplicity. Despite these inaccuracies, relative changes of dimensions during the cardiac cycle as a result of adaptation were simulated realistically.

In the TriSeg model, myofiber strain was estimated from midwall curvature, area, and wall volume by application of the one-fiber model2 to a spherical wall segment (Eq.?15 and Appendix A). The one-fiber model has been shown to be insensitive to actual wall geometry by assuming conservation of energy and homogeneity of fiber stress in the wall.2 So, we expected that the present relation for transmural pressure as a function of midwall surface area and curvature was also applicable to the real, more irregular cardiac geometry, although this fact has not been proven.

The analytically derived dependency of myofiber strain on wall segment geometry (Eq.?A7) appeared continuous and differentiable around z?=?0. To avoid numerical inaccuracy near zero curvature, because of zero-division, a fourth-order Taylor series approximation (Eq.?15) was used instead. For similar reasons, Eq. (A4) for midwall tension was also approximated by a fourth-order Taylor series (Eq.?16). Within the physiological range of ventricular geometry, the errors of the approximations as compared to the analytically derived relations were smaller than 1% and 2% for strain and tension, respectively. With respect to total ventricular pump work, total myofiber stroke work was overestimated less than 2.6%, whereas total tensile stroke work at the midwall surface was underestimated less than 0.5%.

In the NORM and PHCh simulations, a load-controlling adaptation mechanism was applied to render size and mass of each wall segment. The applied adaptation rules required prescribed values for maximum and minimum sarcomere lengths and for maximum passive myofiber stress in each wall segment (Table?1). Maximum and minimum sarcomere lengths were derived from experiments on isolated cardiac muscle of the rat.15,31,64 Maximum passive myofiber stress was chosen as adaptation stimulus because experimental data obtained in dogs with chronic volume overload suggested that end-diastolic myofiber stress and ejection strain were important mechanical stimuli for hypertrophy, while peak systolic myofiber stress appeared irrelevant.22 However, when assuming similar values of maximum passive myofiber stress in the three ventricular walls of the TriSeg model, ventricular geometry did not develop anatomically accurate. Therefore, in our model, we adjusted levels of maximum diastolic stress per wall segment (Table?1) so that after adaptation the weight ratios for the walls agreed with findings in healthy volunteers.19,28

In this study, the external pressure surrounding the LV and RV free walls was assumed to be zero. The real heart is surrounded by the pericardium, which constrains increase of total heart volume during volume overload. Under resting conditions, the pericardium is believed to play a minor role, setting pericardial pressure close to zero.42,46,67 However, with acute increase of total heart volume, the pericardium affects cardiac hemodynamics and interaction of the cardiac chambers significantly.6,11,25 In our simulations of acute pulmonary hypertension (PHAc), the effect of the pericardium should be considered because of severe RV dilatation (Fig.?5). In the PHAc3 simulation, RV cavity volume even increased by 120% relative to the NORM simulation. Together with the changes of right atrial and LV cavity volumes (45% increase and 10% decrease, respectively), total heart volume increased by 25%. In the chronic PH simulations, the effect of the pericardium on ventricular mechanics is likely to be of minor importance, because measurements demonstrated adaptive dilatation of the pericardium in patients with chronic pulmonary hypertension.8,24

Comparison of Model Simulations with Measurements

The relation between curvature ratio CR and transmural pressure ratio PR as extracted from the NORM and PHCh simulations (Fig.?10) agreed quite well with the relation as found in a patient group consisting of patients with and without chronic pulmonary hypertension.17 The relations were about linear with equal slope. In the simulations, however, an offset was found, shifting the relation to higher CR values, implying overestimation of septal curvature for a given right-to-left pressure ratio.

The overestimation of septal curvature may have many causes. In the model, wall geometry was considered spherical having a clear unique radius of curvature by definition. The real ventricular walls are not spherical,58 implying that radii of curvature along the circumferential and base-to-apex direction are different. Also, the junction of the real ventricular walls is smoothed over the wall boundaries, thus further hindering a clear definition of curvature. Besides an error in the curvature ratio as derived from the simulations, the offset might originate from a systematic error in the pressure or curvature measurements obtained in patients as discussed by Dellegrottaglie et al.17 Nonetheless, the linear relationship between CR and PR might be a useful tool for noninvasive estimation of systolic RV cavity pressure in patients with RV pressure overload.17,35,51 In the simulations, the linear relation appeared also valid (<2% deviation) for the PHAc1 and PHAc2 simulations (Table?4) as well as after acute changes (?20%) of mean arterial pressure and cardiac output in the chronic pulmonary hypertension simulations (not shown in Fig.?10). The latter findings suggested that the linear relation between CR and PR is insensitive to acute changes of hemodynamic status.

Experimental studies with acute manipulation of transseptal pressure difference as well as clinical studies among patients with chronic pulmonary hypertension showed that septal curvature and position of the septum between the LV and RV free walls depend instantaneously on transseptal pressure difference.10,20,21,35,36,52 This dependency was found during systole as well as diastole. Similar dependencies were found in our simulations of pulmonary hypertension. Septal curvature decreased with increase of pulmonary resistance (Figs.?5 and 6). Furthermore, leftward shift of the septum, increase of RV volume, and decrease of LV volume with acute increase of pulmonary resistance (Fig.?5), as predicted by the TriSeg model, were also predicted by Kerckhoffs et al.34 using a finite element model of the ventricles coupled to a lumped circulation model.

Simulations and clinical observations were in agreement concerning changes in timing of cardiac flow events due to pulmonary hypertension (Figs.?7 and 8). For example, acceleration time of pulmonary flow velocity was significantly decreased and varied linearly with mean pulmonary artery pressure.9,14 Furthermore, RV isovolumic contraction and relaxation times increased while pulmonary ejection time decreased.62,70 Moreover, tricuspid E/A-ratio was decreased in patients with chronic pulmonary hypertension indicating deterioration of RV diastolic function.71

Ventricular wall volumes and end-diastolic wall thicknesses (Table?2) in the NORM simulation were all about 25% smaller than values measured in healthy volunteers.19,28,57 This difference most likely resulted from an overestimation of contractility used in the sarcomere mechanics model (Appendix B, Eq.?B7), causing walls to be thinner after adaptation to mechanical load. Furthermore, pulmonary acceleration and ejection times as derived from the NORM simulation (Fig.?8) were underestimated as compared to data obtained in healthy individuals.9,27 A probable cause is inaccuracy in the model description of sarcomere mechanics, which is primarily derived from experiments on isolated cardiac muscle of the rat.15,31,64 Human myocardium under in vivo conditions is likely to behave differently.

In the PHCh simulations, LV pump function and hemodynamics were relatively unaffected by increase of pulmonary resistance (Figs.?6 and 7). Also in rats with chronic pulmonary hypertension, resting LV pump function remains unaffected as long as the myocardium is able to compensate for increased tissue load by structural adaptation.23 Furthermore, in the PHCh simulations, the RV free wall hypertrophied and RV cavity volume decreased (Table?2 and Fig.?6). These geometric changes were in agreement with experimental observations in rats with mild chronic pulmonary hypertension.49 In the latter study, chronic pulmonary artery banding resulted in an increase of RV systolic pressure from 33 to 71?mmHg. This increase of RV afterload resulted in 76% increase of thickness of the RV free wall and 14% decrease of RV free wall area, which suggests a small reduction in RV cavity volume. In our simulations, RV systolic pressure increased from 27?mmHg in the NORM simulation to 72?mmHg in the PHCh3 simulation (Fig.?6), end-diastolic wall thickness increased 120% (Table?2), and RV cavity volume decreased 8%. In another animal study in rats,44 it was shown that structural adaptation to chronic increase of RV afterload was associated with concentric hypertrophy up to a certain level. Beyond this level, however, the myocardium could not fully compensate for further load increase. Consequently, the RV cavity dilated, mainly due to dilatation of the RV free wall. In patients with severe chronic pulmonary hypertension, RV free wall hypertrophy is found to occur together with RV dilatation.12,30 Secondary to RV failure, systolic and diastolic LV function deteriorate by ventricular interaction.38,39 As a result, cardiac performance deteriorates and exercise capacity reduces.12,45 These phenomena indicating RV failure were also found in our PHAc simulations.

In the PHAc simulations, the myocardial tissue was unable to compensate for increased RV afterload by load adaptation. As a result, mechanical myofiber load was inhomogeneously distributed over the ventricular walls (Table?3 and Fig.?5). RW stroke work density increased with increase of RV afterload, whereas LW and SW stroke work densities decreased. In the PHAc3 simulation, RW stroke work density was increased by almost 200% with respect to the NORM simulation. When assuming stroke work density to be correlated to oxygen consumption,59 this implies an increase of oxygen demand by the RV free wall. Since oxygen supply is not included in our model, the potential effect of perfusion imposed limitations on myocardial performance is unknown in our simulations of pulmonary hypertension.

In patients with pulmonary hypertension, tricuspid regurgitation often occurs due to RV and tricuspid annular dilatation.18,65 In our simulations, tricuspid valve regurgitation was not included. In the PHAc3 simulation, RV end-diastolic volume was increased by almost 100% with respect to the NORM simulation. It is likely that the absence of tricuspid regurgitation resulted in underestimation of RV volume overload in the PHAc simulations.

The TriSeg model was successfully integrated as a module in the CircAdapt model of the closed-loop cardiovascular system. As shown previously, the CircAdapt environment is flexible by its modular setup and enables realistic simulation of cardiovascular mechanics and hemodynamics under normal as well as various pathological conditions.3,33,40 A set of physiological adaptation rules, expressing structural adaptation of the system to mechanical load, makes the model self-structuring and reduces the number of independent model parameters. For each ventricular wall segment in the TriSeg model, wall volume, midwall area, and reference sarcomere length were varied so that mechanical myofiber load was normalized to a known physiological level, which was assumed to be the same for all ventricular walls. The field of application of the CircAdapt model was substantially enlarged by implementation of the TriSeg model. In the future, the combined model may be used to study fundamental research questions concerning ventricular interaction and its role in cardiac pathologies. For example, specific material properties may be changed per wall segment in order to study effects of heterogeneity of wall properties on cardiac mechanics and hemodynamics, e.g., asynchronous mechanical activation (and pacing) or a localized myocardial infarct.33,40

Although the TriSeg model can be easily modified to include a description of inhomogeneous myocardial wall properties, a finite element model of cardiac mechanics is more accurate and better suited to describe local inhomogeneities in mechanical load. For patient-specific modeling of the circulation, the model should simulate many cycles in order to find a best match with the available set of measurements. For that purpose, the TriSeg model is to be preferred because calculation effort was about 1,000 times less than that of a finite element approach. The finite element model is absolutely needed to estimate and evaluate possible errors introduced by the applied simplifications of the TriSeg model. Currently, the CircAdapt model with the TriSeg module will be evaluated in its possibility to simulate hemodynamics of pulmonary hypertension patient specifically.


Conclusions

We presented the TriSeg model of ventricular mechanics and hemodynamics incorporating mechanical interaction of the LV free wall, RV free wall, and septal wall, resulting in a strong coupling of LV and RV hemodynamics. The model enables calculation of LV and RV pressures given the respective cavity volumes. LV and RV hemodynamics are related to myofiber mechanics in the three ventricular walls, satisfying the principle of conservation of energy. The three ventricular walls are mechanically coupled satisfying equilibrium of tensile forces in their junction. After implementation as a module in the lumped closed-loop CircAdapt model of heart and circulation, the TriSeg model enables simulation of ventricular hemodynamics and wall mechanics as functions of time.

Simulations of cardiac mechanics and hemodynamics during normal ventricular loading, acute pulmonary hypertension, and chronic pulmonary hypertension after load adaptation agreed with clinical data as obtained in normal subjects and in chronic pulmonary hypertension patients. With increasing levels of chronic pulmonary hypertension, the TriSeg model predicted increase of systolic pulmonary flow acceleration, increase of isovolumic contraction and relaxation times, and linear decrease of septal-to-LV free wall curvature ratio. Summarizing, the TriSeg model realistically describes ventricular mechanics including the interaction between left and right ventricular pump mechanics, dynamics of septal geometry, and contractile myofiber function in the three ventricular walls.


We gratefully acknowledge the financial support of Actelion Pharmaceuticals Nederland B.V. (Woerden, The Netherlands) and the Netherlands Heart Foundation Grant 2007B203.

Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.

Appendix A: Derivations of Mechanics of the Curved Wall Patch

This appendix contains the analytical derivations of Eqs. (15) and (16) relating mechanics of the curved wall patch to myofiber mechanics. For the derivation of Eq. (16), calculating representative midwall tension Tm from myofiber stress, the logarithmic one-fiber model relating global ventricular pump mechanics to myofiber mechanics is used.2 Applying the definition of midwall surface as used in Eqs. (3) and (4), for the ratio of transmural pressure pTrans to myofiber stress ?f in a closed thick-walled sphere with midwall volume Vm and wall volume Vw, it holds:

[Formula ID: Equ17]
A1  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\frac{{p_{\text{Trans}} }}{{\sigma_{\text{f}} }}} = \left. {{\frac{1}{3}}\ln \left( {{\frac{1 + z}{1 - z}}} \right)} \right|_{\text{sphere}} \;{\text{with}}\;z = \left. {{\frac{{V_{\text{w}} }}{{2V_{\text{m}} }}}} \right|_{\text{sphere}} $$\end{document}]
where the definition of z is equivalent to the more general definition used in Eq. (15) for a curved wall patch. For a closed sphere with midwall surface area Am and curvature Cm, the relation between pTrans and representative midwall tension Tm is determined by conservation of energy, i.e.,
[Formula ID: Equ18]
A2  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ p_{\text{Trans}} dV_{\text{m}} = \left. {T_{\text{m}} dA_{\text{m}} } \right|_{\text{sphere}} $$\end{document}]
where dVm and dAm represent small increments of midwall volume and area, respectively. Substitution of Vm?=?(4/3)?r3, Am?=?4?r2, and r?=?1/Cm in Eq. (A2) yields:
[Formula ID: Equ19]
A3  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ p_{\text{Trans}} = 2T_{\text{m}} C_{\text{m}} $$\end{document}]

Note that Eq. (A3) represents Laplace?s law expressed in such a way that the applied scalar definition of Tm satisfies conservation of energy (Eq.?A2). Although we have derived Eq. (A3) for a complete sphere, it also holds for a fraction of the same sphere, since both curvature Cm and tension Tm are local parameters. Thus, Eq. (A3) may also be applied to a spherical wall segment (Fig.?2c). Substituting Eq. (A3) in Eq. (A1) and applying Cm?=?(2Amz)/(3Vw) renders for representative midwall tension Tm in a curved wall patch:

[Formula ID: Equ20]
A4  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ T_{\text{m}} = {\frac{{\sigma_{\text{f}} V_{\text{w}} }}{{2A_{\text{m}} }}} \cdot {\frac{1}{2z}}\ln \left( {{\frac{1 + z}{1 - z}}} \right) $$\end{document}]
Note that this latter function is continuous around zero curvature. To avoid numerical inaccuracy near z?=?0, Eq. (A4) is replaced by its close fourth-order polynomial approximation as specified in Eq. (16) (relative error <2%).

For the derivation of Eq. (15), myofiber strain ?f is considered a function of two variables, i.e., Am and Cm. Applying the principle of conservation of energy, work of the wall by a small area change dAm with tension Tm equals work of the myofibers by a change of natural myofiber strain d?f with Cauchy myofiber stress ?f:

[Formula ID: Equ21]
A5  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ T_{\text{m}} dA_{\text{m}} = V_{\text{w}} \sigma_{\text{f}} d\varepsilon_{\text{f}} $$\end{document}]
Rearranging Eq. (A5), substituting Tm with Eq. (A4), and keeping Cm constant renders the following partial derivative:
[Formula ID: Equ22]
A6  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\frac{{\partial \varepsilon_{\text{f}} }}{{\partial A_{\text{m}} }}} = {\frac{1}{{4A_{\text{m}} z}}}\ln \left( {{\frac{1 + z}{1 - z}}} \right)\quad {\text{with}}\;\;z = {\frac{{3C_{\text{m}} V_{\text{w}} }}{{2A_{\text{m}} }}} $$\end{document}]
A standard mathematical solution technique is used to find the general solution of the first-order partial differential equation in Eq. (A6).16 For that purpose, ?f is written as the sum of a strain term ?f,AC and a strain term ?f,C:
[Formula ID: Equ23]
A7  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \varepsilon_{\text{f}} = \varepsilon_{\text{f,AC}} (A_{\text{m}} ,C_{\text{m}} ) + \varepsilon_{\text{f,C}} (C_{\text{m}} ) $$\end{document}]

First, integration of Eq. (A6) with respect to Am renders strain term ?f,AC under the condition that Cm remains constant during integration:

[Formula ID: Equ24]
A8  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \varepsilon_{\text{f,AC}} (A_{\text{m}} ,C_{\text{m}} ) = {\frac{1}{2}}\ln \left( {{\frac{{A_{\text{m}} }}{{A_{\text{m,ref}} }}}} \right) + {\frac{1}{4z}}\ln \left( {{\frac{{(1 + z)^{1 + z} }}{{(1 - z)^{1 - z} }}}} \right) $$\end{document}]
where Am,ref represents a reference value of midwall surface area Am. Below, the term ?f,C will be derived for a relatively simple special case of inflating a closed sphere. The thus found expression is valid for all combinations of Am and Cm. For a closed thick-walled sphere myofiber strain ?f is known as a function of Vm. Using the principle that work as generated by the myofibers equals pump work, it holds:
[Formula ID: Equ25]
A9  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ V_{\text{w}} \sigma_{\text{f}} d\varepsilon_{\text{f}} = \left. {p_{\text{Trans}} dV_{\text{m}} } \right|_{\text{sphere}} $$\end{document}]

Substitution of Eq. (A1) in Eq. (A9) renders for a closed sphere:

[Formula ID: Equ26]
A10  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\frac{{d\varepsilon_{\text{f}} }}{{dV_{\text{m}} }}} = \left. {{\frac{1}{{3V_{\text{w}} }}}\ln \left( {{\frac{1 + z}{1 - z}}} \right)} \right|_{\text{sphere}} $$\end{document}]

Then, midwall volume Vm and curvature ratio z are expressed as functions of curvature Cm:

[Formula ID: Equ27]
A11  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ V_{\text{m}} = \left. {{\frac{4\pi }{{3C_{\text{m}}^{3} }}}} \right|_{\text{sphere}} $$\end{document}]
[Formula ID: Equ28]
A12  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ z = \left. {{\frac{{3V_{\text{w}} C_{\text{m}}^{3} }}{8\pi }}} \right|_{\text{sphere}} $$\end{document}]
Substitution of Eqs. (A11) and (A12) in Eq. (A10) results in a differential equation in Cm only. Integration of this differential equation with respect to Cm yields an expression of myofiber strain ?f for a closed sphere as function of Cm (full equation given in Appendix C):
[Formula ID: Equ29]
A13  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \varepsilon_{\text{f}} = \left. {f(C_{\text{m}} )} \right|_{\text{sphere}} $$\end{document}]

Applying Eq. (A8) to a closed sphere with curvature Cm and substituting the resulting expression as well as Eq. (A13) in Eq. (A7), results in the unknown strain term ?f,C (full equation given in Appendix C):

[Formula ID: Equ30]
A14  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \varepsilon_{\text{f,C}} (C_{\text{m}} ) = \left. {\varepsilon_{\text{f}} } \right|_{\text{sphere}} - \left. {\varepsilon_{\text{f,AC}} (A_{\text{m}} ,C_{\text{m}} )} \right|_{\text{sphere}} $$\end{document}]

Finally, substitution of Eqs. (A8) and (A14) in Eq. (A7) yields a complicated analytical expression of natural myofiber strain ?f in a curved wall patch (Appendix C). Mathematica 5.2 (Wolfram Research, Inc) was used to obtain a simpler fourth-order polynomial fit as specified in Eq. (15). This latter expression closely approximates the analytical derivation (relative error <1%).

Appendix B: Model of Myofiber Mechanics

In this appendix, the applied heuristic relation is presented that converts natural myofiber strain ?f to Cauchy myofiber stress ?f. This empirical representation of sarcomere contraction simulates experiments on isolated rat cardiac muscle.15,64 Parameter values are given in Table?B1.

Natural myofiber strain is converted to sarcomere length Ls by

[Formula ID: Equ31]
B1  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ L_{\text{s}} = L_{\text{s,ref}} e^{{\varepsilon_{\text{f}} }} $$\end{document}]

The sarcomere has been modeled as a passive element in parallel with a series combination of a contractile element and series elastic element. The time-dependent behavior is described by two state variables, i.e., contractile element length Lsc and mechanical activation C. This latter activation parameter is physiologically related to intracellular calcium concentration. The time derivative of Lsc depends linearly on length of the series elastic element (Ls???Lsc) and equals zero for isometric contraction:

[Formula ID: Equ32]
B2  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\frac{{dL_{\text{sc}} }}{dt}} = \left( {{\frac{{L_{\text{s}} - L_{\text{sc}} }}{{L_{\text{se,iso}} }}} - 1} \right)v_{\max } $$\end{document}]
where vmax represents velocity of sarcomere shortening with zero load and Lse,iso length of the isometrically stressed series elastic element (Table?B1). Dependence on vmax represents the myofiber force?velocity relation so that shortening velocity increases with applied external force.

The time derivative of C is heuristically obtained, having separate terms to describe rise and decay of mechanical activation C:

[Formula ID: Equ33]
B3  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\frac{dC}{dt}} = {\frac{1}{{\tau_{\text{R}} }}} \cdot C_{\text{L}} (L_{\text{sc}} ) \cdot F_{\text{rise}} (t) + {\frac{1}{{\tau_{\text{D}} }}} \cdot {\frac{{C_{\text{rest}} - C}}{{1 + e^{{{{\left( {T(L_{\text{sc}} ) - t} \right)} \mathord{\left/ {\vphantom {{\left( {T(L_{\text{sc}} ) - t} \right)} {\tau_{\text{D}} }}} \right. \kern-\nulldelimiterspace} {\tau_{\text{D}} }}}} }}} $$\end{document}]
Parameters ?R and ?D are scaling rise and decay time, respectively (Table?B1). Symbols t and Crest represent time and diastolic resting level of activation, respectively. Functions CL, Frise, and T describe increase of activation with sarcomere length, rise of mechanical activation, and decrease of activation duration with decrease of sarcomere length, respectively:
[Formula ID: Equ34]
B4  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ C_{\text{L}} (L_{\text{sc}} ) = \tanh \left( {4.0(L_{\text{sc}} - L_{{{\text{sc}}0}} )^{2} } \right) $$\end{document}]
[Formula ID: Equ35]
B5  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F_{\text{rise}} (t) = 0.02 \cdot x^{3} (8 - x)^{2} {\text{e}}^{ - x} \quad {\text{with}}\;\;x = \min \left( {8,\max (0,t/\tau_{\text{R}} )} \right) $$\end{document}]
[Formula ID: Equ36]
B6  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ T(L_{\text{sc}} ) = \tau_{\text{sc}} (0.29 + 0.3L_{\text{sc}} ) $$\end{document}]
where Lsc0 and ?sc represent contractile element length with zero load and a time factor scaling duration of contraction, respectively (Table?B1).

Active myofiber stress ?f,act depends on length of series elastic element, Lsc, and C:

[Formula ID: Equ37]
B7  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \sigma_{\text{f,act}} = \sigma_{\text{act}} \cdot C \cdot (L_{\text{sc}} - L_{\text{sc0}} ) \cdot {\frac{{L_{\text{s}} - L_{\text{sc}} }}{{L_{\text{se,iso}} }}} $$\end{document}]

In order to simulate the contributions of titin as well as collagen to passive myocardial stiffness,69 passive stress ?f,pas represents a soft behavior for low myofiber strain and a stiff behavior for large positive strain:

[Formula ID: Equ38]
B8  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \sigma_{\text{f,pas}} = \sigma_{\text{pas}} \left( {36 \cdot \max (0,\varepsilon_{\text{f}} - 0.1)^{2} + 0.1(\varepsilon_{\text{f}} - 0.1) + 0.0025{\text{e}}^{{30\varepsilon_{\text{f}} }} } \right) $$\end{document}]

The value of ?pas has been adapted so that at maximum sarcomere length with exercise, passive myofiber stress assumes the maximum value (Table?1: Set point values for adaptation). Total Cauchy myofiber stress ?f is the sum of active and passive stress:

[Formula ID: Equ39]
B9  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \sigma_{\text{f}} (\varepsilon_{\text{f}} ) = \sigma_{\text{f,pas}} + \sigma_{\text{f,act}} $$\end{document}]
The latter result has been substituted in Eq. (16).

Appendix C: Polynomial Series Approximations

The polynomial series approximation given by Eq. (15) is obtained from the analytical solution expressing natural myofiber strain in terms of Am and Cm (Eq.?A7) by applying the following steps:

  1. Expression of myofiber strain ?f for a closed sphere (use Am?=?4?/Cm2) as function of Cm (Eq.?A13) yields:

[Formula ID: Equ40]
C1a  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left. {\varepsilon_{\text{f}} } \right|_{\text{sphere}} = - \ln (C_{\text{m}} ) + {\frac{\ln (8\pi )}{3}} + {\frac{G(x)}{6}} $$\end{document}]
with
[Formula ID: Equ41]
C1b  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ G(x) = {\frac{1}{x}} \cdot \ln \left( {{\frac{{(1 + x)^{1 + x} }}{{(1 - x)^{1 - x} }}}} \right) $$\end{document}]
[Formula ID: Equ42]
C1c  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x = \left. z \right|_{\text{sphere}} = {\frac{{3C_{m}^{3} V_{w} }}{8\pi }} $$\end{document}]
  1. Applying Eq. (A8) to a closed sphere (use Am?=?4?/Cm2) yields:
[Formula ID: Equ43]
C2  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left. {\varepsilon_{\text{f,AC}} } \right|_{\text{sphere}} = - \ln (C_{\text{m}} ) + {\frac{1}{2}} \cdot \ln \left( {{\frac{4\pi }{{A_{\text{m,ref}} }}}} \right) + {\frac{G(x)}{4}} $$\end{document}]
Substitution of Eqs. (C1a) and (C2) in Eq. (A14) yields for strain term ?f,C:
[Formula ID: Equ44]
C3  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \varepsilon_{\text{f,C}} = {\frac{\ln (8\pi )}{3}} - {\frac{1}{2}} \cdot \ln \left( {{\frac{4\pi }{{A_{\text{m,ref}} }}}} \right) - {\frac{G(x)}{12}} $$\end{document}]

The thus found expression for ?f,C, which holds for the special case of inflating a sphere, is valid for all combinations of Am and Cm.

  1. Substitution of Eqs. (C3) and (A8) in Eq. (A7) yields the general analytical solution of natural myofiber strain ?f in a curved wall patch:
[Formula ID: Equ45]
C4  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \varepsilon_{\text{f}} = {\frac{\ln (8\pi )}{3}} - {\frac{1}{2}} \cdot \ln \left( {{\frac{4\pi }{{A_{\text{m,ref}} }}}} \right) + {\frac{1}{2}} \cdot \ln \left( {{\frac{{A_{\text{m}} }}{{A_{\text{m,ref}} }}}} \right) - {\frac{G(x)}{12}} + {\frac{G(z)}{4}} $$\end{document}]
  1. In order to facilitate series approximation, z?=?K/A and x?=?K3 are substituted in Eq. (C4):
[Formula ID: Equ46]
C5a  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \begin{aligned} \varepsilon_{\text{f}} & = {\frac{\ln (8\pi )}{3}} - {\frac{1}{2}} \cdot \ln \left( {{\frac{4\pi }{{A_{\text{m,ref}} }}}} \right) + {\frac{1}{2}} \cdot \ln \left( {{\frac{{A_{\text{m}} }}{{A_{\text{m,ref}} }}}} \right) \\ & \quad - {\frac{1}{{12K^{3} }}} \cdot \ln \left( {{\frac{{(1 + K^{3} )^{{1 + K^{3} }} }}{{(1 - K^{3} )^{{1 - K^{3} }} }}}} \right) + {\frac{A}{4K}} \cdot \ln \left( {{\frac{{\left( {1 + {\frac{K}{A}}} \right)^{{1 + {\frac{K}{A}}}} }}{{\left( {1 - {\frac{K}{A}}} \right)^{{1 - {\frac{K}{A}}}} }}}} \right) \\ \end{aligned} $$\end{document}]
with A and K representing normalized area and curvature, respectively:
[Formula ID: Equ47]
C5b  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ A = {\frac{{A_{\text{m}} }}{{\root{3}\of{{{9\pi V_{\text{w}}^{2} }}}}}} $$\end{document}]
[Formula ID: Equ48]
C5c  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ K = {\frac{{C_{\text{m}} }}{2}}\root{3}\of{{{{\frac{{3V_{\text{w}} }}{\pi }}}}} $$\end{document}]
  1. The first two terms of Eq. (C5a) are constant strain offsets and are therefore eliminated from the strain relation, which describes the transition from a reference to a new state with state variables Am and Cm. Taylor series approximation of this simplified version of Eq. (C5a) yields:
[Formula ID: Equ49]
C6  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \varepsilon_{\text{f}} = {\frac{1}{2}} \cdot \ln \left( {{\frac{{A_{\text{m}} }}{{A_{\text{m,ref}} }}}} \right) - {\frac{{K^{2} }}{{12A^{2} }}} - {\frac{{K^{4} }}{{40A^{4} }}} + {\rm O}(K^{6} ) $$\end{document}]
where ?(K6) represents the higher order terms.
  1. It is used that z?=?K/A. Furthermore, the fourth-order polynomial approximation of natural myofiber strain in a curved wall patch (Eq.?15) is obtained after slight adjustment of the fourth-order term of Eq. (C6) for best accuracy within the whole range of operation (|z|?<?0.8).

The polynomial series approximation given by Eq. (16) is obtained from the analytical solution expressing representative midwall tension in terms of Am, Cm, and ?f (Eq.?A4) by applying the following steps:

  1. Substitution of [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ z = 3C_{\text{m}} V_{\text{w}} /2A_{\text{m}} $$\end{document}] in Eq. (A4) yields:
[Formula ID: Equ50]
C7  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\frac{{T_{\text{m}} }}{{\sigma_{\text{f}} }}} = {\frac{1}{{6C_{\text{m}} }}}\ln \left( {{\frac{{2A_{\text{m}} + 3C_{\text{m}} V_{\text{w}} }}{{2A_{\text{m}} - 3C_{\text{m}} V_{\text{w}} }}}} \right) $$\end{document}]
  1. Similar to what has been done to obtain Eq. (C5a), A and K are substituted in Eq. (C7) to facilitate series approximation:
[Formula ID: Equ51]
C8a  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\frac{{T_{\text{m}} }}{{\sigma_{\text{f}} }}} = {\frac{L}{6K}}\ln \left( {{\frac{{1 + {\frac{K}{A}}}}{{1 - {\frac{K}{A}}}}}} \right) $$\end{document}]
with
[Formula ID: Equ52]
C8b  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ L = {\frac{1}{2}}\root{3}\of{{{{\frac{{3V_{\text{w}} }}{\pi }}}}} $$\end{document}]
  1. Taylor series approximation of Eq. (C8a) yields:
[Formula ID: Equ53]
C9  [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\frac{{T_{\text{m}} }}{{\sigma_{f} }}} = {\frac{L}{3A}} + {\frac{{LK^{2} }}{{9A^{3} }}} + {\frac{{LK^{4} }}{{15A^{5} }}} + {\rm O}(K^{6} ) $$\end{document}]
where ?(K6) represents the higher order terms, which are neglected in the TriSeg model.
  1. Finally, the fourth-order polynomial approximation of representative midwall tension (Eq.?16), is obtained by rearranging Eq. (C9) and using z?=?K/A.


References
1.. Arts T,Bovendeerd P,Delhaas T,Prinzen F. Modeling the relation between cardiac pump function and myofiber mechanicsJ. Biomech. 2003;36:731–736. [doi: 10.1016/S0021-9290(02)00451-7]
1.. Arts, T., P. Bovendeerd, T. Delhaas, and F. Prinzen. Modeling the relation between cardiac pump function and myofiber mechanics. J. Biomech. 36:731?736, 2003. [pmid: 12695003]
2.. Arts T,Bovendeerd PH,Prinzen FW,Reneman RS. Relation between left ventricular cavity pressure and volume and systolic fiber stress and strain in the wallBiophys. J. 1991;59:93–102. [doi: 10.1016/S0006-3495(91)82201-9]
2.. Arts, T., P. H. Bovendeerd, F. W. Prinzen, and R. S. Reneman. Relation between left ventricular cavity pressure and volume and systolic fiber stress and strain in the wall. Biophys. J. 59:93?102, 1991. [pmid: 2015392]
3.. Arts T,Delhaas T,Bovendeerd P,Verbeek X,Prinzen FW. Adaptation to mechanical load determines shape and properties of heart and circulation: the CircAdapt modelAm. J. Physiol. Heart Circ. Physiol. 2005;288:H1943–H1954. [doi: 10.1152/ajpheart.00444.2004]
3.. Arts, T., T. Delhaas, P. Bovendeerd, X. Verbeek, and F. W. Prinzen. Adaptation to mechanical load determines shape and properties of heart and circulation: the CircAdapt model. Am. J. Physiol. Heart Circ. Physiol. 288:H1943?H1954, 2005. [pmid: 15550528]
4.. Arts T,Prinzen FW,Snoeckx LH,Rijcken JM,Reneman RS. Adaptation of cardiac structure by mechanical feedback in the environment of the cell: a model studyBiophys. J. 1994;66:953–961. [doi: 10.1016/S0006-3495(94)80876-8]
4.. Arts, T., F. W. Prinzen, L. H. Snoeckx, J. M. Rijcken, and R. S. Reneman. Adaptation of cardiac structure by mechanical feedback in the environment of the cell: a model study. Biophys. J. 66:953?961, 1994. [pmid: 8038399]
5.. Baker AE,Dani R,Smith ER,Tyberg JV,Belenkie I. Quantitative assessment of independent contributions of pericardium and septum to direct ventricular interactionAm. J. Physiol. 1998;275:H476–H483.
5.. Baker, A. E., R. Dani, E. R. Smith, J. V. Tyberg, and I. Belenkie. Quantitative assessment of independent contributions of pericardium and septum to direct ventricular interaction. Am. J. Physiol. 275:H476?H483, 1998. [pmid: 9683435]
6.. Belenkie I,Dani R,Smith ER,Tyberg JV. The importance of pericardial constraint in experimental pulmonary embolism and volume loadingAm. Heart J. 1992;123:733–742. [doi: 10.1016/0002-8703(92)90514-V]
6.. Belenkie, I., R. Dani, E. R. Smith, and J. V. Tyberg. The importance of pericardial constraint in experimental pulmonary embolism and volume loading. Am. Heart J. 123:733?742, 1992. [pmid: 1539525]
7.. Beyar R,Dong SJ,Smith ER,Belenkie I,Tyberg JV. Ventricular interaction and septal deformation: a model compared with experimental dataAm. J. Physiol. 1993;265:H2044–H2056.
7.. Beyar, R., S. J. Dong, E. R. Smith, I. Belenkie, and J. V. Tyberg. Ventricular interaction and septal deformation: a model compared with experimental data. Am. J. Physiol. 265:H2044?H2056, 1993. [pmid: 8285243]
8.. Blanchard DG,Dittrich HC. Pericardial adaptation in severe chronic pulmonary hypertension. An intraoperative transesophageal echocardiographic studyCirculation 1992;85:1414–1422.
8.. Blanchard, D. G., and H. C. Dittrich. Pericardial adaptation in severe chronic pulmonary hypertension. An intraoperative transesophageal echocardiographic study. Circulation 85:1414?1422, 1992. [pmid: 1532541]
9.. Bossone E,Avelar E,Bach DS,Gillespie B,Rubenfire M,Armstrong WF. Diagnostic value of resting tricuspid regurgitation velocity and right ventricular ejection flow parameters for the detection of exercise induced pulmonary arterial hypertensionInt. J. Card. Imaging 2000;16:429–436. [doi: 10.1023/A:1010604913656]
9.. Bossone, E., E. Avelar, D. S. Bach, B. Gillespie, M. Rubenfire, and W. F. Armstrong. Diagnostic value of resting tricuspid regurgitation velocity and right ventricular ejection flow parameters for the detection of exercise induced pulmonary arterial hypertension. Int. J. Card. Imaging 16:429?436, 2000. [pmid: 11482708]
10.. Brinker JA,Weiss JL,Lappe DL,Rabson JL,Summer WR,Permutt S,Weisfeldt ML. Leftward septal displacement during right ventricular loading in manCirculation 1980;61:626–633.
10.. Brinker, J. A., J. L. Weiss, D. L. Lappe, J. L. Rabson, W. R. Summer, S. Permutt, and M. L. Weisfeldt. Leftward septal displacement during right ventricular loading in man. Circulation 61:626?633, 1980. [pmid: 7353253]
11.. Calvin JE. Optimal right ventricular filling pressures and the role of pericardial constraint in right ventricular infarction in dogsCirculation 1991;84:852–861.
11.. Calvin, J. E. Optimal right ventricular filling pressures and the role of pericardial constraint in right ventricular infarction in dogs. Circulation 84:852?861, 1991. [pmid: 1713537]
12.. Chemla D,Castelain V,Herve P,Lecarpentier Y,Brimioulle S. Haemodynamic evaluation of pulmonary hypertensionEur. Respir. J. 2002;20:1314–1331. [doi: 10.1183/09031936.02.00068002]
12.. Chemla, D., V. Castelain, P. Herve, Y. Lecarpentier, and S. Brimioulle. Haemodynamic evaluation of pulmonary hypertension. Eur. Respir. J. 20:1314?1331, 2002. [pmid: 12449189]
13.. Chung DC,Niranjan SC,Clark JW Jr,Bidani A,Johnston WE,Zwischenberger JB,Traber DL. A dynamic model of ventricular interaction and pericardial influenceAm. J. Physiol. 1997;272:H2942–H2962.
13.. Chung, D. C., S. C. Niranjan, J. W. Clark, Jr., A. Bidani, W. E. Johnston, J. B. Zwischenberger, and D. L. Traber. A dynamic model of ventricular interaction and pericardial influence. Am. J. Physiol. 272:H2942?H2962, 1997. [pmid: 9227574]
14.. Dabestani A,Mahan G,Gardin JM,Takenaka K,Burn C,Allfie A,Henry WL. Evaluation of pulmonary artery pressure and resistance by pulsed Doppler echocardiographyAm. J. Cardiol. 1987;59:662–668. [doi: 10.1016/0002-9149(87)91189-1]
14.. Dabestani, A., G. Mahan, J. M. Gardin, K. Takenaka, C. Burn, A. Allfie, and W. L. Henry. Evaluation of pulmonary artery pressure and resistance by pulsed Doppler echocardiography. Am. J. Cardiol. 59:662?668, 1987. [pmid: 3825910]
15.. Tombe PP,ter Keurs HE. Force and velocity of sarcomere shortening in trabeculae from rat heart. Effects of temperatureCirc. Res. 1990;66:1239–1254.
15.. de Tombe, P. P., and H. E. ter Keurs. Force and velocity of sarcomere shortening in trabeculae from rat heart. Effects of temperature. Circ. Res. 66:1239?1254, 1990. [pmid: 2335024]
16.. Debnath L. Nonlinear Partial Differential Equations for Scientists and Engineers (2) 20052. Boston: Birkh?user; :737.
16.. Debnath, L. Nonlinear Partial Differential Equations for Scientists and Engineers (2nd ed.). Boston: Birkh?user, p. 737, 2005.
17.. Dellegrottaglie S,Sanz J,Poon M,Viles-Gonzalez JF,Sulica R,Goyenechea M,Macaluso F,Fuster V,Rajagopalan S. Pulmonary hypertension: accuracy of detection with left ventricular septal-to-free wall curvature ratio measured at cardiac MRRadiology 2007;243:63–69. [doi: 10.1148/radiol.2431060067]
17.. Dellegrottaglie, S., J. Sanz, M. Poon, J. F. Viles-Gonzalez, R. Sulica, M. Goyenechea, F. Macaluso, V. Fuster, and S. Rajagopalan. Pulmonary hypertension: accuracy of detection with left ventricular septal-to-free wall curvature ratio measured at cardiac MR. Radiology 243:63?69, 2007. [pmid: 17392248]
18.. Dittrich HC,McCann HA,Blanchard DG. Cardiac structure and function in chronic thromboembolic pulmonary hypertensionAm. J. Card. Imaging 1994;8:18–27.
18.. Dittrich, H. C., H. A. McCann, and D. G. Blanchard. Cardiac structure and function in chronic thromboembolic pulmonary hypertension. Am. J. Card. Imaging 8:18?27, 1994. [pmid: 8130612]
19.. Doherty NE 3rd,Fujita N,Caputo GR,Higgins CB. Measurement of right ventricular mass in normal and dilated cardiomyopathic ventricles using cine magnetic resonance imagingAm. J. Cardiol. 1992;69:1223–1228. [doi: 10.1016/0002-9149(92)90940-Z]
19.. Doherty, 3rd, N. E., N. Fujita, G. R. Caputo, and C. B. Higgins. Measurement of right ventricular mass in normal and dilated cardiomyopathic ventricles using cine magnetic resonance imaging. Am. J. Cardiol. 69:1223?1228, 1992. [pmid: 1575194]
20.. Dong SJ,Crawley AP,MacGregor JH,Petrank YF,Bergman DW,Belenkie I,Smith ER,Tyberg JV,Beyar R. Regional left ventricular systolic function in relation to the cavity geometry in patients with chronic right ventricular pressure overload. A three-dimensional tagged magnetic resonance imaging studyCirculation 1995;91:2359–2370.
20.. Dong, S. J., A. P. Crawley, J. H. MacGregor, Y. F. Petrank, D. W. Bergman, I. Belenkie, E. R. Smith, J. V. Tyberg, and R. Beyar. Regional left ventricular systolic function in relation to the cavity geometry in patients with chronic right ventricular pressure overload. A three-dimensional tagged magnetic resonance imaging study. Circulation 91:2359?2370, 1995. [pmid: 7729022]
21.. Dong SJ,Smith ER,Tyberg JV. Changes in the radius of curvature of the ventricular septum at end diastole during pulmonary arterial and aortic constrictions in the dogCirculation 1992;86:1280–1290.
21.. Dong, S. J., E. R. Smith, and J. V. Tyberg. Changes in the radius of curvature of the ventricular septum at end diastole during pulmonary arterial and aortic constrictions in the dog. Circulation 86:1280?1290, 1992. [pmid: 1394934]
22.. Donker DW,Volders PG,Arts T,Bekkers BC,Hofstra L,Spatjens RL,Beekman JD,Borgers M,Crijns HJ,Vos MA. End-diastolic myofiber stress and ejection strain increase with ventricular volume overload?serial in vivo analyses in dogs with complete atrioventricular blockBasic Res. Cardiol. 2005;100:372–382. [doi: 10.1007/s00395-005-0525-8]
22.. Donker, D. W., P. G. Volders, T. Arts, B. C. Bekkers, L. Hofstra, R. L. Spatjens, J. D. Beekman, M. Borgers, H. J. Crijns, and M. A. Vos. End-diastolic myofiber stress and ejection strain increase with ventricular volume overload?serial in vivo analyses in dogs with complete atrioventricular block. Basic Res. Cardiol. 100:372?382, 2005. [pmid: 15795796]
23.. Faber MJ,Dalinghaus M,Lankhuizen IM,Steendijk P,Hop WC,Schoemaker RG,Duncker DJ,Lamers JM,Helbing WA. Right and left ventricular function after chronic pulmonary artery banding in rats assessed with biventricular pressure?volume loopsAm. J. Physiol. Heart Circ. Physiol. 2006;291:H1580–H1586. [doi: 10.1152/ajpheart.00286.2006]
23.. Faber, M. J., M. Dalinghaus, I. M. Lankhuizen, P. Steendijk, W. C. Hop, R. G. Schoemaker, D. J. Duncker, J. M. Lamers, and W. A. Helbing. Right and left ventricular function after chronic pulmonary artery banding in rats assessed with biventricular pressure?volume loops. Am. J. Physiol. Heart Circ. Physiol. 291:H1580?H1586, 2006. [pmid: 16679397]
24.. Freeman GL,LeWinter MM. Pericardial adaptations during chronic cardiac dilation in dogsCirc. Res. 1984;54:294–300.
24.. Freeman, G. L., and M. M. LeWinter. Pericardial adaptations during chronic cardiac dilation in dogs. Circ. Res. 54:294?300, 1984. [pmid: 6230171]
25.. Freeman GL,LeWinter MM. Role of parietal pericardium in acute, severe mitral regurgitation in dogsAm. J. Cardiol. 1984;54:217–219. [doi: 10.1016/0002-9149(84)90332-1]
25.. Freeman, G. L., and M. M. LeWinter. Role of parietal pericardium in acute, severe mitral regurgitation in dogs. Am. J. Cardiol. 54:217?219, 1984. [pmid: 6741816]
26.. Fujii J,Yazaki Y,Sawada H,Aizawa T,Watanabe H,Kato K. Noninvasive assessment of left and right ventricular filling in myocardial infarction with a two-dimensional Doppler echocardiographic methodJ. Am. Coll. Cardiol. 1985;5:1155–1160.
26.. Fujii, J., Y. Yazaki, H. Sawada, T. Aizawa, H. Watanabe, and K. Kato. Noninvasive assessment of left and right ventricular filling in myocardial infarction with a two-dimensional Doppler echocardiographic method. J. Am. Coll. Cardiol. 5:1155?1160, 1985. [pmid: 3989126]
27.. Gardin JM,Burn CS,Childs WJ,Henry WL. Evaluation of blood flow velocity in the ascending aorta and main pulmonary artery of normal subjects by Doppler echocardiographyAm. Heart J. 1984;107:310–319. [doi: 10.1016/0002-8703(84)90380-6]
27.. Gardin, J. M., C. S. Burn, W. J. Childs, and W. L. Henry. Evaluation of blood flow velocity in the ascending aorta and main pulmonary artery of normal subjects by Doppler echocardiography. Am. Heart J. 107:310?319, 1984. [pmid: 6695664]
28.. Hajduczok ZD,Weiss RM,Stanford W,Marcus ML. Determination of right ventricular mass in humans and dogs with ultrafast cardiac computed tomographyCirculation 1990;82:202–212.
28.. Hajduczok, Z. D., R. M. Weiss, W. Stanford, and M. L. Marcus. Determination of right ventricular mass in humans and dogs with ultrafast cardiac computed tomography. Circulation 82:202?212, 1990. [pmid: 2364512]
29.. Hirschfeld S,Meyer R,Schwartz DC,Korfhagen J,Kaplan S. Measurement of right and left ventricular systolic time intervals by echocardiographyCirculation 1975;51:304–309.
29.. Hirschfeld, S., R. Meyer, D. C. Schwartz, J. Korfhagen, and S. Kaplan. Measurement of right and left ventricular systolic time intervals by echocardiography. Circulation 51:304?309, 1975. [pmid: 1112010]
30.. Horan LG,Flowers NC,Havelda CJ. Relation between right ventricular mass and cavity size: an analysis of 1500 human heartsCirculation 1981;64:135–138.
30.. Horan, L. G., N. C. Flowers, and C. J. Havelda. Relation between right ventricular mass and cavity size: an analysis of 1500 human hearts. Circulation 64:135?138, 1981. [pmid: 7237711]
31.. Hunter PJ,McCulloch AD,ter Keurs HE. Modelling the mechanical properties of cardiac muscleProg. Biophys. Mol. Biol. 1998;69:289–331. [doi: 10.1016/S0079-6107(98)00013-3]
31.. Hunter, P. J., A. D. McCulloch, and H. E. ter Keurs. Modelling the mechanical properties of cardiac muscle. Prog. Biophys. Mol. Biol. 69:289?331, 1998. [pmid: 9785944]
32.. Kaul S. The interventricular septum in health and diseaseAm. Heart J. 1986;112:568–581. [doi: 10.1016/0002-8703(86)90523-5]
32.. Kaul, S. The interventricular septum in health and disease. Am. Heart J. 112:568?581, 1986. [pmid: 3529903]
33.. Kerckhoffs RC,Lumens J,Vernooy K,Omens JH,Mulligan LJ,Delhaas T,Arts T,McCulloch AD,Prinzen FW. Cardiac resynchronization: insight from experimental and computational modelsProg. Biophys. Mol. Biol. 2008;97:543–561. [doi: 10.1016/j.pbiomolbio.2008.02.024]
33.. Kerckhoffs, R. C., J. Lumens, K. Vernooy, J. H. Omens, L. J. Mulligan, T. Delhaas, T. Arts, A. D. McCulloch, and F. W. Prinzen. Cardiac resynchronization: insight from experimental and computational models. Prog. Biophys. Mol. Biol. 97:543?561, 2008. [pmid: 18417196]
34.. Kerckhoffs RC,Neal ML,Gu Q,Bassingthwaighte JB,Omens JH,McCulloch AD. Coupling of a 3D finite element model of cardiac ventricular mechanics to lumped systems models of the systemic and pulmonic circulationAnn. Biomed. Eng. 2007;35:1–18. [doi: 10.1007/s10439-006-9212-7]
34.. Kerckhoffs, R. C., M. L. Neal, Q. Gu, J. B. Bassingthwaighte, J. H. Omens, and A. D. McCulloch. Coupling of a 3D finite element model of cardiac ventricular mechanics to lumped systems models of the systemic and pulmonic circulation. Ann. Biomed. Eng. 35:1?18, 2007. [pmid: 17111210]
35.. King ME,Braun H,Goldblatt A,Liberthson R,Weyman AE. Interventricular septal configuration as a predictor of right ventricular systolic hypertension in children: a cross-sectional echocardiographic studyCirculation 1983;68:68–75.
35.. King, M. E., H. Braun, A. Goldblatt, R. Liberthson, and A. E. Weyman. Interventricular septal configuration as a predictor of right ventricular systolic hypertension in children: a cross-sectional echocardiographic study. Circulation 68:68?75, 1983. [pmid: 6851056]
36.. Kingma I,Tyberg JV,Smith ER. Effects of diastolic transseptal pressure gradient on ventricular septal position and motionCirculation 1983;68:1304–1314.
36.. Kingma, I., J. V. Tyberg, and E. R. Smith. Effects of diastolic transseptal pressure gradient on ventricular septal position and motion. Circulation 68:1304?1314, 1983. [pmid: 6640880]
37.. Little WC,Reeves RC,Arciniegas J,Katholi RE,Rogers EW. Mechanism of abnormal interventricular septal motion during delayed left ventricular activationCirculation 1982;65:1486–1491.
37.. Little, W. C., R. C. Reeves, J. Arciniegas, R. E. Katholi, and E. W. Rogers. Mechanism of abnormal interventricular septal motion during delayed left ventricular activation. Circulation 65:1486?1491, 1982. [pmid: 7074805]
38.. Louie EK,Lin SS,Reynertson SI,Brundage BH,Levitsky S,Rich S. Pressure and volume loading of the right ventricle have opposite effects on left ventricular ejection fractionCirculation 1995;92:819–824.
38.. Louie, E. K., S. S. Lin, S. I. Reynertson, B. H. Brundage, S. Levitsky, and S. Rich. Pressure and volume loading of the right ventricle have opposite effects on left ventricular ejection fraction. Circulation 92:819?824, 1995. [pmid: 7641362]
39.. Louie EK,Rich S,Levitsky S,Brundage BH. Doppler echocardiographic demonstration of the differential effects of right ventricular pressure and volume overload on left ventricular geometry and fillingJ. Am. Coll. Cardiol. 1992;19:84–90.
39.. Louie, E. K., S. Rich, S. Levitsky, and B. H. Brundage. Doppler echocardiographic demonstration of the differential effects of right ventricular pressure and volume overload on left ventricular geometry and filling. J. Am. Coll. Cardiol. 19:84?90, 1992. [pmid: 1729350]
40.. Lumens, J., T. Delhaas, B. Kirn, and T. Arts. Modeling ventricular interaction: a multiscale approach from sarcomere mechanics to cardiovascular system hemodynamics. Pac. Symp. Biocomput. 378?389, 2008.
41.. Luo C,Ware DL,Zwischenberger JB,Clark JW Jr. Using a human cardiopulmonary model to study and predict normal and diseased ventricular mechanics, septal interaction, and atrio-ventricular blood flow patternsCardiovasc. Eng. 2007;7:17–31. [doi: 10.1007/s10558-007-9025-9]
41.. Luo, C., D. L. Ware, J. B. Zwischenberger, and J. W. Clark, Jr. Using a human cardiopulmonary model to study and predict normal and diseased ventricular mechanics, septal interaction, and atrio-ventricular blood flow patterns. Cardiovasc. Eng. 7:17?31, 2007. [pmid: 17334942]
42.. Mangano DT,Dyke DC,Hickey RF,Ellis RJ. Significance of the pericardium in human subjects: effects on left ventricular volume, pressure and ejectionJ. Am. Coll. Cardiol. 1985;6:290–295.
42.. Mangano, D. T., D. C. Van Dyke, R. F. Hickey, and R. J. Ellis. Significance of the pericardium in human subjects: effects on left ventricular volume, pressure and ejection. J. Am. Coll. Cardiol. 6:290?295, 1985. [pmid: 4019916]
43.. Maughan WL,Sunagawa K,Sagawa K. Ventricular systolic interdependence: volume elastance model in isolated canine heartsAm. J. Physiol. 1987;253:H1381–H1390.
43.. Maughan, W. L., K. Sunagawa, and K. Sagawa. Ventricular systolic interdependence: volume elastance model in isolated canine hearts. Am. J. Physiol. 253:H1381?H1390, 1987. [pmid: 3425740]
44.. Minami S,Onodera T,Okazaki F,Miyazaki H,Ohsawa S,Mochizuki S. Myocyte morphological characteristics differ between the phases of pulmonary hypertension-induced ventricular hypertrophy and failureInt. Heart J. 2006;47:629–637. [doi: 10.1536/ihj.47.629]
44.. Minami, S., T. Onodera, F. Okazaki, H. Miyazaki, S. Ohsawa, and S. Mochizuki. Myocyte morphological characteristics differ between the phases of pulmonary hypertension-induced ventricular hypertrophy and failure. Int. Heart J. 47:629?637, 2006. [pmid: 16960417]
45.. Miyamoto S,Nagaya N,Satoh T,Kyotani S,Sakamaki F,Fujita M,Nakanishi N,Miyatake K. Clinical correlates and prognostic significance of six-minute walk test in patients with primary pulmonary hypertension. Comparison with cardiopulmonary exercise testingAm. J. Respir. Crit. Care Med. 2000;161:487–492.
45.. Miyamoto, S., N. Nagaya, T. Satoh, S. Kyotani, F. Sakamaki, M. Fujita, N. Nakanishi, and K. Miyatake. Clinical correlates and prognostic significance of six-minute walk test in patients with primary pulmonary hypertension. Comparison with cardiopulmonary exercise testing. Am. J. Respir. Crit. Care Med. 161:487?492, 2000. [pmid: 10673190]
46.. Morris-Thurgood JA,Frenneaux MP. Diastolic ventricular interaction and ventricular diastolic fillingHeart Fail. Rev. 2000;5:307–323. [doi: 10.1023/A:1026555012135]
46.. Morris-Thurgood, J. A., and M. P. Frenneaux. Diastolic ventricular interaction and ventricular diastolic filling. Heart Fail. Rev. 5:307?323, 2000. [pmid: 11498644]
47.. Nash MP,Hunter PJ. Computational mechanics of the heart: from tissue structure to ventricular functionJ. Elast. 2000;61:113–141. [doi: 10.1023/A:1011084330767]
47.. Nash, M. P., and P. J. Hunter. Computational mechanics of the heart: from tissue structure to ventricular function. J. Elast. 61:113?141, 2000.
48.. Olansen JB,Clark JW,Khoury D,Ghorbel F,Bidani A. A closed-loop model of the canine cardiovascular system that includes ventricular interactionComput. Biomed. Res. 2000;33:260–295. [doi: 10.1006/cbmr.2000.1543]
48.. Olansen, J. B., J. W. Clark, D. Khoury, F. Ghorbel, and A. Bidani. A closed-loop model of the canine cardiovascular system that includes ventricular interaction. Comput. Biomed. Res. 33:260?295, 2000. [pmid: 10944405]
49.. Olivetti G,Ricci R,Lagrasta C,Maniga E,Sonnenblick EH,Anversa P. Cellular basis of wall remodeling in long-term pressure overload-induced right ventricular hypertrophy in ratsCirc. Res. 1988;63:648–657.
49.. Olivetti, G., R. Ricci, C. Lagrasta, E. Maniga, E. H. Sonnenblick, and P. Anversa. Cellular basis of wall remodeling in long-term pressure overload-induced right ventricular hypertrophy in rats. Circ. Res. 63:648?657, 1988. [pmid: 2970334]
50.. Olsen CO,Tyson GS,Maier GW,Spratt JA,Davis JW,Rankin JS. Dynamic ventricular interaction in the conscious dogCirc. Res. 1983;52:85–104.
50.. Olsen, C. O., G. S. Tyson, G. W. Maier, J. A. Spratt, J. W. Davis, and J. S. Rankin. Dynamic ventricular interaction in the conscious dog. Circ. Res. 52:85?104, 1983. [pmid: 6848213]
51.. Reisner SA,Azzam Z,Halmann M,Rinkevich D,Sideman S,Markiewicz W,Beyar R. Septal/free wall curvature ratio: a noninvasive index of pulmonary arterial pressureJ. Am. Soc. Echocardiogr. 1994;7:27–35.
51.. Reisner, S. A., Z. Azzam, M. Halmann, D. Rinkevich, S. Sideman, W. Markiewicz, and R. Beyar. Septal/free wall curvature ratio: a noninvasive index of pulmonary arterial pressure. J. Am. Soc. Echocardiogr. 7:27?35, 1994. [pmid: 8155331]
52.. Roeleveld RJ,Marcus JT,Faes TJ,Gan TJ,Boonstra A,Postmus PE,Vonk-Noordegraaf A. Interventricular septal configuration at mr imaging and pulmonary arterial pressure in pulmonary hypertensionRadiology 2005;234:710–717. [doi: 10.1148/radiol.2343040151]
52.. Roeleveld, R. J., J. T. Marcus, T. J. Faes, T. J. Gan, A. Boonstra, P. E. Postmus, and A. Vonk-Noordegraaf. Interventricular septal configuration at mr imaging and pulmonary arterial pressure in pulmonary hypertension. Radiology 234:710?717, 2005. [pmid: 15634939]
53.. Santamore WP,Burkhoff D. Hemodynamic consequences of ventricular interaction as assessed by model analysisAm. J. Physiol. 1991;260:H146–H157.
53.. Santamore, W. P., and D. Burkhoff. Hemodynamic consequences of ventricular interaction as assessed by model analysis. Am. J. Physiol. 260:H146?H157, 1991. [pmid: 1992793]
54.. Sechtem U,Pflugfelder PW,Gould RG,Cassidy MM,Higgins CB. Measurement of right and left ventricular volumes in healthy individuals with cine MR imagingRadiology 1987;163:697–702.
54.. Sechtem, U., P. W. Pflugfelder, R. G. Gould, M. M. Cassidy, and C. B. Higgins. Measurement of right and left ventricular volumes in healthy individuals with cine MR imaging. Radiology 163:697?702, 1987. [pmid: 3575717]
55.. Slinker BK,Glantz SA. End-systolic and end-diastolic ventricular interactionAm. J. Physiol. 1986;251:H1062–H1075.
55.. Slinker, B. K., and S. A. Glantz. End-systolic and end-diastolic ventricular interaction. Am. J. Physiol. 251:H1062?H1075, 1986. [pmid: 3777194]
56.. Smith BW,Chase JG,Shaw GM,Nokes RI. Simulating transient ventricular interaction using a minimal cardiovascular system modelPhysiol. Meas. 2006;27:165–179. [doi: 10.1088/0967-3334/27/2/007]
56.. Smith, B. W., J. G. Chase, G. M. Shaw, and R. I. Nokes. Simulating transient ventricular interaction using a minimal cardiovascular system model. Physiol. Meas. 27:165?179, 2006. [pmid: 16400203]
57.. Stolzmann P,Scheffel H,Leschka S,Schertler T,Frauenfelder T,Kaufmann PA,Marincek B,Alkadhi H. Reference values for quantitative left ventricular and left atrial measurements in cardiac computed tomographyEur. Radiol. 2008;18:1625–1634. [doi: 10.1007/s00330-008-0939-4]
57.. Stolzmann, P., H. Scheffel, S. Leschka, T. Schertler, T. Frauenfelder, P. A. Kaufmann, B. Marincek, and H. Alkadhi. Reference values for quantitative left ventricular and left atrial measurements in cardiac computed tomography. Eur. Radiol. 18:1625?1634, 2008. [pmid: 18446346]
58.. Streeter DD Jr. Berne RM,Sperelakis N,Geiger SRGross morphology and fiber geometry of the heartHandbook of Physiology 1979Bethesda, MD: American Physiological Society; :61–112.
58.. Streeter Jr., D. D. Gross morphology and fiber geometry of the heart. In: Handbook of Physiology, edited by R. M. Berne, N. Sperelakis, and S. R. Geiger. Bethesda, MD: American Physiological Society, 1979, pp. 61?112.
59.. Suga H,Hayashi T,Shirahata M. Ventricular systolic pressure?volume area as predictor of cardiac oxygen consumptionAm. J. Physiol. 1981;240:H39–H44.
59.. Suga, H., T. Hayashi, and M. Shirahata. Ventricular systolic pressure?volume area as predictor of cardiac oxygen consumption. Am. J. Physiol. 240:H39?H44, 1981. [pmid: 7457620]
60.. Sun Y,Beshara M,Lucariello RJ,Chiaramida SA. A comprehensive model for right-left heart interaction under the influence of pericardium and baroreflexAm. J. Physiol. 1997;272:H1499–H1515.
60.. Sun, Y., M. Beshara, R. J. Lucariello, and S. A. Chiaramida. A comprehensive model for right-left heart interaction under the influence of pericardium and baroreflex. Am. J. Physiol. 272:H1499?H1515, 1997. [pmid: 9087629]
61.. Taher MF,Santamore WP,Bogen DK. Ventricular interaction is described by three coupling coefficientsAm. J. Physiol. 1994;266:H228–H234.
61.. Taher, M. F., W. P. Santamore, and D. K. Bogen. Ventricular interaction is described by three coupling coefficients. Am. J. Physiol. 266:H228?H234, 1994. [pmid: 8304503]
62.. Tei C,Dujardin KS,Hodge DO,Bailey KR,McGoon MD,Tajik AJ,Seward SB. Doppler echocardiographic index for assessment of global right ventricular functionJ. Am. Soc. Echocardiogr. 1996;9:838–847. [doi: 10.1016/S0894-7317(96)90476-9]
62.. Tei, C., K. S. Dujardin, D. O. Hodge, K. R. Bailey, M. D. McGoon, A. J. Tajik, and S. B. Seward. Doppler echocardiographic index for assessment of global right ventricular function. J. Am. Soc. Echocardiogr. 9:838?847, 1996. [pmid: 8943444]
63.. Tei C,Dujardin KS,Hodge DO,Kyle RA,Tajik AJ,Seward JB. Doppler index combining systolic and diastolic myocardial performance: clinical value in cardiac amyloidosisJ. Am. Coll. Cardiol. 1996;28:658–664. [doi: 10.1016/0735-1097(96)00202-1]
63.. Tei, C., K. S. Dujardin, D. O. Hodge, R. A. Kyle, A. J. Tajik, and J. B. Seward. Doppler index combining systolic and diastolic myocardial performance: clinical value in cardiac amyloidosis. J. Am. Coll. Cardiol. 28:658?664, 1996. [pmid: 8772753]
64.. ter Keurs HE,Rijnsburger WH,Heuningen R,Nagelsmit MJ. Tension development and sarcomere length in rat cardiac trabeculae. Evidence of length-dependent activationCirc. Res. 1980;46:703–714.
64.. ter Keurs, H. E., W. H. Rijnsburger, R. van Heuningen, and M. J. Nagelsmit. Tension development and sarcomere length in rat cardiac trabeculae. Evidence of length-dependent activation. Circ. Res. 46:703?714, 1980. [pmid: 7363419]
65.. Thistlethwaite PA,Jamieson SW. Tricuspid valvular disease in the patient with chronic pulmonary thromboembolic diseaseCurr. Opin. Cardiol. 2003;18:111–116. [doi: 10.1097/00001573-200303000-00007]
65.. Thistlethwaite, P. A., and S. W. Jamieson. Tricuspid valvular disease in the patient with chronic pulmonary thromboembolic disease. Curr. Opin. Cardiol. 18:111?116, 2003. [pmid: 12652215]
66.. Usyk TP,Grice IJ,McCulloch AD. Computational model of three-dimensional cardiac electromechanicsComput. Vis. Sci. 2002;4:249–257. [doi: 10.1007/s00791-002-0081-9]
66.. Usyk, T. P., I. J. Le Grice, and A. D. McCulloch. Computational model of three-dimensional cardiac electromechanics. Comput. Vis. Sci. 4:249?257, 2002.
67.. Watkins MW,LeWinter MM. Physiologic role of the normal pericardiumAnnu. Rev. Med. 1993;44:171–180. [doi: 10.1146/annurev.me.44.020193.001131]
67.. Watkins, M. W., and M. M. LeWinter. Physiologic role of the normal pericardium. Annu. Rev. Med. 44:171?180, 1993. [pmid: 8476238]
68.. Weber KT,Janicki JS,Shroff S,Fishman AP. Contractile mechanics and interaction of the right and left ventriclesAm. J. Cardiol. 1981;47:686–695. [doi: 10.1016/0002-9149(81)90556-7]
68.. Weber, K. T., J. S. Janicki, S. Shroff, and A. P. Fishman. Contractile mechanics and interaction of the right and left ventricles. Am. J. Cardiol. 47:686?695, 1981. [pmid: 7008570]
69.. Wu Y,Cazorla O,Labeit D,Labeit S,Granzier H. Changes in titin and collagen underlie diastolic stiffness diversity of cardiac muscleJ. Mol. Cell. Cardiol. 2000;32:2151–2162. [doi: 10.1006/jmcc.2000.1281]
69.. Wu, Y., O. Cazorla, D. Labeit, S. Labeit, and H. Granzier. Changes in titin and collagen underlie diastolic stiffness diversity of cardiac muscle. J. Mol. Cell. Cardiol. 32:2151?2162, 2000. [pmid: 11112991]
70.. Yeo TC,Dujardin KS,Tei C,Mahoney DW,McGoon MD,Seward JB. Value of a Doppler-derived index combining systolic and diastolic time intervals in predicting outcome in primary pulmonary hypertensionAm. J. Cardiol. 1998;81:1157–1161. [doi: 10.1016/S0002-9149(98)00140-4]
70.. Yeo, T. C., K. S. Dujardin, C. Tei, D. W. Mahoney, M. D. McGoon, and J. B. Seward. Value of a Doppler-derived index combining systolic and diastolic time intervals in predicting outcome in primary pulmonary hypertension. Am. J. Cardiol. 81:1157?1161, 1998. [pmid: 9605059]
71.. Yu CM,Sanderson JE,Chan S,Yeung L,Hung YT,Woo KS. Right ventricular diastolic dysfunction in heart failureCirculation 1996;93:1509–1514.
71.. Yu, C. M., J. E. Sanderson, S. Chan, L. Yeung, Y. T. Hung, and K. S. Woo. Right ventricular diastolic dysfunction in heart failure. Circulation 93:1509?1514, 1996. [pmid: 8608618]

Nomenclature General
LV Left ventricle/ventricular
RV Right ventricle/ventricular
LW Left ventricular free wall
SW Septal wall
RW Right ventricular free wall
Geometry-Related Parameters
VLV Left ventricular cavity volume (m3)
VRV Right ventricular cavity volume (m3)
Vw Wall volume of ventricular wall segment (m3)
Vm Volume of spherical cap, formed by midwall surface of wall segment (m3)
Am Midwall surface area of curved wall segment (m2)
Cm Curvature of midwall surface (reciprocal of radius) (m?1)
xm Maximal axial distance from midwall surface to origin (m)
ym Radius of midwall junction circle (m)
z Ratio of wall thickness to midwall radius of curvature of curved wall segment
? Half the opening angle of spherical midwall surface
?f Natural myofiber strain
Force-Related Parameters
pLV Left ventricular cavity pressure (Pa)
pRV Right ventricular cavity pressure (Pa)
pTrans Transmural pressure difference across curved wall segment (Pa)
?f Cauchy myofiber stress (Pa)
Tm Representative midwall tension (N?m?1)
Tx Axial midwall tension component (N?m?1)
Ty Radial midwall tension component (N?m?1)

Figures

[Figure ID: Fig1]
Figure?1 

Two typical 2D echocardiographic images showing end-diastolic parasternal short-axis views of the hearts of a normal individual (a) and a pulmonary hypertension patient (b). Note that the septum is flattened with pulmonary hypertension. Cross-sections of left (LV) and right (RV) ventricular cavities, LV free wall (LW), septal wall (SW), and RV free wall (RW) are indicated. The bright dots indicate the cross-section of the junction margin of the walls. In the right upper corner, image plane orientation (gray plane) is indicated in a schematic representation of the cardiac ventricles



[Figure ID: Fig2]
Figure?2 

TriSeg model of ventricular mechanics. (a) The TriSeg model (gray shading) incorporated in the modular CircAdapt model of the systemic (syst) and pulmonary (pulm) circulations.3 Three thick-walled spherical segments representing LV free wall (LW), RV free wall (RW), and septal wall (SW) form a ventricular composite consisting of a left ventricular (LV) cavity, with volume VLV and pressure pLV, and a right ventricular (RV) cavity, with volume VRV and pressure pRV. (b) Cross-section of the ventricular composite through the axis of rotational symmetry (x). The ventricular walls are coupled mechanically at the common junction circle where their midwall surfaces (thick lines) meet. This junction circle perpendicularly intersects the plane of drawing (open dots). The midwall surface divides the wall with wall volume Vw in two shells of equal volume. The center of the junction circle is the origin (0) of the cylindrical coordinate system applied. Note that x-direction is defined positive toward the RV free wall. (c) Cross-section of a single wall segment through the axis of rotational symmetry. Note that only midwall geometry is shown. The shaded area indicates midwall volume Vm enclosed between the midwall surface and the junction plane. Midwall volume Vm, area Am, and curvature Cm (reciprocal of radius of curvature) depend on distance xm and radius ym of the boundary circle. Each wall is loaded by transmural pressure difference resulting in a representative midwall tension Tm. With half opening angle ?, Tm is resolved in an axial (Tx) and a radial (Ty) component



[Figure ID: Fig3]
Figure?3 

Flowchart of the TriSeg model of left (LV) and right (RV) ventricular mechanics. Applied equations are indicated by their numbers. The ventricular hemodynamics section serves as interface with the CircAdapt model of the whole circulation. In the ventricular composite section, for each wall i, midwall boundary radius (ym,i) and midwall volume (Vm,i) are calculated from LV and RV volumes (VLV and VRV, respectively) and septal geometry. Septal geometry is initially defined by estimates of boundary radius (ym,est) and midwall volume (Vm,SW,est). LV free wall, septum, and RV free wall are indicated by i?=?LW, SW, and RW, respectively. The first dashed arrow indicates calculation of axial (Tx,i) and radial (Ty,i) wall tension components and transmural pressure (pTrans,i). LV and RV pressure (pLV and pRV, respectively) and a new estimate of septal geometry is obtained from the summed axial (Tx,Tot) and radial (Ty,Tot) tension components. If septal geometry is sufficiently accurate by iteration, the solution is accepted. Because direct calculation along the dashed arrow is not possible, the sequence of calculations has to follow the solid lines. For each wall segment, midwall area (Am,i) and curvature (Cm,i) are calculated. Wall tension (Tm,i) is calculated following the second dashed arrow and renders wall tension components and transmural pressure. Again, direct calculation is not possible. Therefore, in the myofiber section, midwall area and curvature are used to calculate myofiber strain (?f,i). Myofiber strain is converted to myofiber stress (?f,i) using a model of sarcomere mechanics. Myofiber stress and wall segment geometry are used to calculate wall tension. By cascading the different sections, calculations can be performed sequentially along the pathways marked by solid lines



[Figure ID: Fig4]
Figure?4 

Simulated time courses of normal (NORM) LV and RV hemodynamics at rest. The bar plot (bottom panel) displays RV and LV ejection, filling, and isovolumic contraction and relaxation periods. Note that simulation duration amounts to one and a half cardiac cycle. Symbols: A, late active ventricular filling wave; E, early passive ventricular filling wave; IC, isovolumic contraction; IR, isovolumic relaxation; LV, left ventricular; pao, aortic pressure; pLV, left ventricular pressure; ppa, pulmonary artery pressure; pRV, right ventricular pressure; qao, aortic valve flow; qmi, mitral valve flow; qpv, pulmonary valve flow; qtr, tricuspid valve flow; RV, right ventricular; VLV, left ventricular cavity volume; and VRV, right ventricular cavity volume



[Figure ID: Fig5]
Figure?5 

Simulated ventricular pressure?volume loops, myofiber stress?strain loops, and time courses of midwall curvatures under normal (NORM) and acute pulmonary hypertensive (PHAc1, PHAc2, and PHAc3) ventricular loading conditions. Acute increase of pulmonary resistance resulted in RV dilatation, increased RV pump work, increased contractile work density in the RV free wall, and decreased septal midwall curvature. Simulation duration amounts to one cardiac cycle. LV ejection is indicated in the curvature plots. Abbreviations: LV, left ventricular; LW, left ventricular free wall; RV, right ventricular; RW, right ventricular free wall; and SW, septal wall. *Note that the absolute value of LW midwall curvature is used for graphical representation, while this curvature is defined to be negative in the TriSeg model



[Figure ID: Fig6]
Figure?6 

Simulated ventricular pressure?volume loops, myofiber stress?strain loops, and time courses of midwall curvatures under normal (NORM) and chronic pulmonary hypertensive (PHCh1, PHCh2, and PHCh3) ventricular loading conditions. Chronic increase of pulmonary resistance including adaptation resulted in increase of RV pump work, homogeneity of contractile work densities in the walls, and linear decrease of septal midwall curvature. Simulation duration amounts to one cardiac cycle. LV ejection is indicated in the curvature plots. Abbreviations: LV, left ventricular; LW, left ventricular free wall; RV, right ventricular; RW, right ventricular free wall; and SW, septal wall. *Note that the absolute value of LW midwall curvature is used for graphical representation, while this curvature is defined to be negative in the TriSeg model



[Figure ID: Fig7]
Figure?7 

Simulated time courses of pulmonary and tricuspid valve flow velocities (upper two panels) show lengthening of RV isovolumic contraction, shortening of RV ejection, and increase of pulmonary flow acceleration after adaptation to increased pulmonary resistance (PHCh1, PHCh2, and PHCh3). Aortic and mitral valve flow velocities (lower two panels) appear relatively insensitive to increase of pulmonary resistance. Simulation duration amounts to one cardiac cycle. The following time intervals are indicated for the NORM simulation: AT, acceleration time of pulmonary valve flow; ET, RV ejection time; ICT, RV isovolumic contraction time; and IRT, RV isovolumic relaxation time. Symbols: A, late active ventricular filling wave; and E, early passive ventricular filling wave



[Figure ID: Fig8]
Figure?8 

Comparison of simulated and patient data on timing of RV hemodynamics. Simulated changes of AT, ET, ICT, and IRT as chronic response to increased pulmonary resistance agree with measurements obtained in healthy normals and patients with pulmonary hypertension [adapted from Bossone et al.9 (upper row) and Tei et al.62 (lower row)]. Patient data are presented by mean values with standard deviations indicated by error bars. Symbols: AT, pulmonary flow acceleration time; ET, RV ejection time; ICT, RV isovolumic contraction time; and IRT, RV isovolumic relaxation time



[Figure ID: Fig9]
Figure?9 

Simulated time courses of SW/LW transmural pressure ratio (upper panel) and SW/LW midwall curvature ratio (lower panel) appear similar revealing a linear relation between both ratios at each time point in the cardiac cycle. Simulation duration amounts to one cardiac cycle. Duration of LV ejection is indicated by vertical dashed lines



[Figure ID: Fig10]
Figure?10 

Simulation results on the correlation between curvature ratio CR and transmural pressure ratio PR (black data points), see text for definitions, as compared with the relation extracted from data (gray data points) obtained in patients with and without pulmonary hypertension (PH?=?RV systolic pressure?>?40?mmHg) [adapted from Dellegrottaglie et al.17]. Linear regression analysis of data points as derived from the NORM and PHCh simulations revealed a strong correlation (bold black line, r2?=?1.00, SEE?=?0.008) between CR and PR with similar slope as the relation found in patients (thick gray dashed line, r2?=?0.73, SEE?=?0.044). The thin dashed lines indicate the 95% confidence interval of linear regression fit to the patient data. Linear regression formulas are shown for simulated as well as measured data points. Symbols: r, regression residual and SEE, standard error of the estimate



Tables
[TableWrap ID: Tab1] Table?1 

Input parameter values for NORM and PH simulations


Parameter Unit Value
Hemodynamics at rest
Mean systemic arterial blood pressure kPa 12.2
Mean systemic blood flow mL?s?1 85
Cardiac cycle time s 0.850
Mean pulmonary arteriovenous pressure drop kPa 1.5 (NORM)
Interventions with exercise
Mean systemic blood flow mL?s?1 255
Cardiac cycle time s 0.425
Set point values for adaptation
Mean flow velocity in large blood vessels (rest) m?s?1 0.17
Maximum vascular wall stress (exercise) kPa 500
Maximum sarcomere length (exercise) ?m 2.2
Minimum sarcomere length (exercise) ?m 1.75
Maximum passive myofiber stress (exercise) kPa 10 (LW), 8 (SW), 20 (RW)
Interventions with PH simulations
Mean pulmonary arteriovenous pressure drop kPa 3.0 (PHAc1 and PHCh1), 4.5 (PHAc2 and PHCh2), 6.0 (PHAc3 and PHCh3)
Parameter values for model initialization
Wall volume (Vw) mL 75 (LW), 40 (SW), 30 (RW)
Reference midwall surface area (Am,ref) cm2 80 (LW), 45 (SW), 100 (RW)
Reference sarcomere length (Ls,ref) ?m 2.0 (LW, SW, and RW)
Septal midwall volume (Vm,SW) mL 42
Radius of midwall junction circle (ym) cm 3.3

LW, left ventricular free wall; NORM, normal simulation; PHAc, acute pulmonary hypertension simulation; PHCh, chronic pulmonary hypertension simulation; RW, right ventricular free wall; SW, septal wall


[TableWrap ID: Tab2] Table?2 

Simulation results: ventricular wall thickness


Parameter Unit NORM PHAc PHCh
1 2 3 1 2 3
End-diastolic wall thickness
Hed,LW mm 7.9 7.8 7.6 7.3 8.3 8.2 8.3
Hed,SW mm 6.8 6.9 7.0 6.9 6.7 6.9 6.5
Hed,RW mm 2.3 2.1 2.0 1.7 3.3 4.2 5.1

Hed, end-diastolic wall thickness; LV, left ventricular; LW, left ventricular free wall; NORM, normal simulation; PHAc, acute pulmonary hypertension simulation; PHCh, chronic pulmonary hypertension simulation; RV, right ventricular; RW, right ventricular free wall; SW, septal wall


[TableWrap ID: Tab3] Table?3 

Simulation results: pulmonary artery pressures and ventricular work


Parameter Unit NORM PHAc PHCh
1 2 3 1 2 3
Ventricular pump stroke work
Wstroke,LV Pa?m3 1.04 1.04 1.03 1.02 1.05 1.04 1.04
Wstroke,RV Pa?m3 0.22 0.32 0.42 0.51 0.35 0.48 0.61
Stroke work density
wstroke,LW kPa 9.41 9.54 9.32 8.66 8.88 8.90 8.78
wstroke,SW kPa 8.94 8.28 7.54 6.15 9.19 8.66 8.49
wstroke,RW kPa 7.91 11.66 15.94 21.59 8.59 9.13 9.48
Mean pulmonary artery pressure
ppa kPa 2.04 3.56 5.11 6.77 3.53 5.00 6.47

LV, left ventricular; LW, left ventricular free wall; NORM, normal simulation; PHAc, acute pulmonary hypertension simulation; PHCh, chronic pulmonary hypertension simulation; ppa, mean pulmonary artery pressure; RV, right ventricular; RW, right ventricular free wall; SW, septal wall; wstroke, stroke work density; Wstroke, ventricular pump stroke work


[TableWrap ID: Tab4] Table?4 

Simulation results: curvature ratio and transmural pressure ratio


Parameter Unit NORM PHAc PHCh
1 2 3 1 2 3
CR ? 0.94 0.85 0.72 0.54 0.81 0.67 0.55
PR ? 0.78 0.68 0.57 0.46 0.65 0.52 0.40

CR, septal-to-LV free wall midwall curvature ratio at the moment of aortic valve closure; NORM, normal simulation; PHAc, acute pulmonary hypertension simulation; PHCh, chronic pulmonary hypertension simulation; PR, systolic septal-to-LV free wall transmural pressure ratio


[TableWrap ID: Tab5] Table B1 

Parameter values of the sarcomere mechanics model


Symbol Unit Value Parameter description
Crest ? 0.02 Diastolic resting level of activation
Lsc0 ?m 1.51 Contractile element length with zero active stress
Ls,ref ?m 2.0 Reference sarcomere length at zero strain
Lse,iso ?m 0.04 Length of isometrically stressed series elastic element
vmax ?m?s?1 7 Sarcomere shortening velocity with zero load
?D ms 32 Factor scaling contraction decay time
?R ms 48 Factor scaling contraction rise time
?sc ms 425 Factor scaling duration of contraction
?act kPa 120 Factor scaling active myofiber stress
?pas kPa 7 Factor scaling passive myofiber stress


Article Categories:
  • Article

Keywords: Keywords Pulmonary hypertension, Septal motion, Adaptation, Stress, Strain, Myofiber, Cardiac mechanics.

Previous Document:  Ethnic health care advisors: a good strategy to improve the access to health care and social welfare...
Next Document:  Pediatric surgical oncology in China: present and future.