Improved discretisation and linearisation of active tension in strongly coupled cardiac electro-mechanics simulations.

Mathematical models of cardiac electro-mechanics typically consist of three tightly coupled parts: systems of ordinary differential equations describing electro-chemical reactions and cross-bridge dynamics in the muscle cells, a system of partial differential equations modelling the propagation of the electrical activation through the tissue and a nonlinear elasticity problem describing the mechanical deformations of the heart muscle. The complexity of the mathematical model motivates numerical methods based on operator splitting, but simple explicit splitting schemes have been shown to give severe stability problems for realistic models of cardiac electro-mechanical coupling. The stability may be improved by adopting semi-implicit schemes, but these give rise to challenges in updating and linearising the active tension. In this paper we present an operator splitting framework for strongly coupled electro-mechanical simulations and discuss alternative strategies for updating and linearising the active stress component. Numerical experiments demonstrate considerable performance increases from an update method based on a generalised Rush-Larsen scheme and a consistent linearisation of active stress based on the first elasticity tensor.

Uncertainty analysis of ventricular mechanics using the probabilistic collocation method.

Uncertainty and variability in material parameters are fundamental challenges in computational biomechanics. Analyzing and quantifying the resulting uncertainty in computed results with parameter sweeps or Monte Carlo methods has become very computationally demanding. In this paper, we consider a stochastic method named the probabilistic collocation method, and investigate its applicability for uncertainty analysis in computing the passive mechanical behavior of the left ventricle. Specifically, we study the effect of uncertainties in material input parameters upon response properties such as the increase in cavity volume, the elongation of the ventricle, the increase in inner radius, the decrease in wall thickness, and the rotation at apex. The numerical simulations conducted herein indicate that the method is well suited for the problem of consideration, and is far more efficient than the Monte Carlo simulation method for obtaining a detailed uncertainty quantification. The numerical experiments also give interesting indications on which material parameters are most critical for accurately determining various global responses.