Comments to the manuscript:
This manuscript describes methods for simulating mechanics of the cardiac ventricles. The authors use finite deformation hyperelasticity with an orthotropic constitutive relation, and linearize and solve the governing nonlinear equations using total-Lagrangian system dynamics. They present a limited validation study in order to verify their implementation, followed by a simulation of ventricular systole that is only qualitatively compared to a limited set of MR images. Intro/Methods/Results sections are presented rather well, however Discussion of the findings and limitations are overly brief or non-existent.
The authors claim to be the first to implement an orthotropic exponential constitutive relation (“Costa law”) at the ventricular level [Intro: last parag], however the UCSD group have published several studies on the use of orthotropic exponential constitutive relations for whole ventricle mechanics (e.g. Usyk et al. Journal of Elasticity 61:143–164, 2000; which also uses the penalty method for incompressibility). Throughout the manuscript, the authors repeatedly highlight some issues regarding one previous study (ref [9]), such as the choice of constitutive relation (“pole-zero law”) and coordinate system (“prolate spheroidal coordinates”). Instead, they would do better to concentrate on the/novel aspects of their own work – noting that models using Cartesian coordinates (e.g. Vetter et al., Annals of Biomedical Engineering, 28:781–792, 2000), and orthotropic material relations (e.g. Usyk et al., as above) have been widely studied.
The authors also claim that their representation of the constraint that the pericardial sac places on wall motion is novel. However, the use of this choice of boundary conditions is remains to be quantitatively validated against experimental recordings or medical images (more on this below). Furthermore, a similar constraint was previously described by the INRIA group (Sermesant et al.).
Perhaps the main novelty of this study is the use of total-Lagrangian system dynamics numerical methods for cardiac electro-mechanics, but this requires further validation (more on this below). It would be interesting to quantitatively compare over multiple heart cycles the accuracy of this technique against that of the more traditional implicit finite element methods.
The validation studies presented here are inadequate. A deformation study involving simple shear is presented. Analytical homogenous data are compared against numerical (finite element) non-homogeneous model predictions (Fig 1a shows non-homogeneous deformation). Thus two different types of deformation are being compared. Also of concern is: “yet the errors increase when the deformation becomes large” (parag. below Eq.14) – this requires further investigation (e.g. did errors reduce with smaller step size?). Furthermore, this test is based on just one mode of deformation, which involves just 2 of the 6 independent components of the strain tensor. A full validation should include multiple modes of deformation involving the shear and normal components of strain. Perhaps a more appropriate and comprehensive closed-form solution is that of the pressure inflation, axial extension and torsion of a homogeneous, isotropic cylinder (eg. Rivlin 1950, Phil. Trans. A242: 173–195).
Finally, another major shortcoming of this manuscript is the distinct lack of quantitative validation against experimental/clinical measures of ventricular wall motion. Just one set of medical images are presented and the comparison against model predictions is limited to a few qualitative statements (see parag. below Fig 4). Such gross observations have been previously reproduced using axi-symmetric truncated ellipsoids (e.g. the Dutch groups of Arts, Bovendeerd, et al.). The use of a geometrically detailed biventricular model should be justified and validated using regional measures of wall motion (e.g. from tagged MRI).
Specific minor comments:
- In several places, “In consequence” should be changed to -> “As a consequence”.
- Page 2, Intro, line 3: “At the macroscopic”
- Page 2, Intro, line 10: “apart from” -> “as well as”
- Page 2, bottom: last line “Furthermore, …” needs rewording.
- Page 3, bottom: last line “The f-n-s, …” needs rewording.
- Page 4, “Elasticity Tensor”, line 3: “stain” -> “strain”
- Page 5, eq (8): there are other nearly incompressible formulations. Cite a ref, and explain why this one was selected.
- Page 5, Sec 2.3: it is unclear why you would want to “impose transmural shear resistance…”. This sentence is very unclear. Note that another (simpler?) way to approximate the effects would be to add an external traction constraint on the epicardium. Discuss pros/cons of this.
- Page 9, line 1: “20 frames” of “5 frames” as pictured?
- Page 9, bottom parag: the first 5-6 lines describes methods (and should thus be moved to Methods section) and takes material constants from a pig study. But at the bottom, there is a suggestion that these parameters are too stiff (“stroke volume is very small”). Thus, it would perhaps be better to use parameters for canine myocardium (e.g. from Usyk et al. Journal of Elasticity 61: 143–164, 2000).
- Page 10, bottom: “only the blood filling phase is included”: this point is confusing and probably needs rewording. Suggests that the authors only studied diastole (however results for systole were also presented).