Uncertainty in cardiac myofiber orientation and stiffnesses dominate the variability of left ventricle deformation response

Computational cardiac modelling is a mature area of biomedical computing and is currently evolving from a pure research tool to aiding in clinical decision making. Assessing the reliability of computational model predictions is a key factor for clinical use, and uncertainty quantification (UQ) and sensitivity analysis are important parts of such an assessment. In this study, we apply UQ in computational heart mechanics to study uncertainty both in material parameters characterizing global myocardial stiffness and in the local muscle fiber orientation that governs tissue anisotropy. The uncertainty analysis is performed using the polynomial chaos expansion (PCE) method, which is a nonintrusive meta‐modeling technique that surrogates the original computational model with a series of orthonormal polynomials over the random input parameter space. In addition, in order to study variability in the muscle fiber architecture, we model the uncertainty in orientation of the fiber field as an approximated random field using a truncated Karhunen‐Loéve expansion. The results from the UQ and sensitivity analysis identify clear differences in the impact of various material parameters on global output quantities. Furthermore, our analysis of random field variations in the fiber architecture demonstrate a substantial impact of fiber angle variations on the selected outputs, highlighting the need for accurate assignment of fiber orientation in computational heart mechanics models. Sample of random fiber orientation field generated by a Gaussian random field with a standard deviation σKL = 0.5 radians and correlation length equals 3 cm (left panel). This loss of the helical arrangements of the myofiber orientation has a substantial impact on global mechanical properties of the ventricle, as is shown for the inner cavity volume (right panel).