Patient-specific computer simulations can be a powerful tool in clinical applications, helping in diagnostics and the development of new treatments. However, its practical use depends on the reliability of the models. The construction of cardiac simulations involves several steps with inherent uncertainties, including model parameters, the generation of personalized geometry and fibre orientation assignment, which are semi-manual processes subject to errors. Thus, it is important to quantify how these uncertainties impact model predictions. The present work performs uncertainty quantification and sensitivity analyses to assess the variability in important quantities of interest (QoI). Clinical quantities are analysed in terms of overall variability and to identify which parameters are the major contributors. The analyses are performed for simulations of the left ventricle function during the entire cardiac cycle. Uncertainties are incorporated in several model parameters, including regional wall thickness, fibre orientation, passive material parameters, active stress and the circulatory model. The results show that the QoI are very sensitive to active stress, wall thickness and fibre direction, where ejection fraction and ventricular torsion are the most impacted outputs. Thus, to improve the precision of models of cardiac mechanics, new methods should be considered to decrease uncertainties associated with geometrical reconstruction, estimation of active stress and of fibre orientation. This article is part of the theme issue ‘Uncertainty quantification in cardiac and cardiovascular modelling and simulation’.

No abstract

No abstract

Computational modeling of the heart is a subject of substantial medical and scientific interest, which may contribute to increase the understanding of several phenomena associated with cardiac physiological and pathological states. Modeling the mechanics of the heart have led to considerable insights, but it still represents a complex and a demanding computational problem, especially in a strongly coupled electromechanical setting. Passive cardiac tissue is commonly modeled as hyperelastic and is characterized by quasi-incompressible, orthotropic, and nonlinear material behavior. These factors are known to be very challenging for the numerical solution of the model. The near-incompressibility is known to cause numerical issues such as the well-known locking phenomenon and ill-conditioning of the stiffness matrix. In this work, the augmented Lagrangian method is used to handle the nearly incompressible condition. This approach can potentially improve computational performance by reducing the condition number of the stiffness matrix and thereby improving the convergence of iterative solvers. We also improve the performance of iterative solvers by the use of an algebraic multigrid preconditioner. Numerical results of the augmented Lagrangian method combined with a preconditioned iterative solver for a cardiac mechanics benchmark suite are presented to show its improved performance.