AbstractComputational models of the cardiovascular system and heart function are currently being investigated as analytic tools to assist medical practice and clinical trials. Recent technological advances allow for finite element models of heart ventricles and atria to be customized to medical images and to assimilate electrical and hemodynamic measurements. Optimizing model parameters to physiological data is, however, challenging due to the computational complexity of finite element models. Metaheuristic algorithms and other optimization strategies typically require sampling hundreds of points in the model parameter space before converging to optimal solutions. Similarly, resolving uncertainty of model outputs to input assumptions is difficult for finite element models due to their computational cost. In this paper, we present a novel, multifidelity strategy for model order reduction of 3-D finite element models of ventricular mechanics. Our approach is centered around well established findings on the similarity between contraction of an isolated muscle and the whole ventricle. Specifically, we demonstrate that simple linear transformations between sarcomere strain (tension) and ventricular volume (pressure) are sufficient to reproduce global pressure-volume outputs of 3-D finite element models even by a reduced model with just a single myocyte unit. We further develop a procedure for congruency training of a surrogate low-order model from multi-scale finite elements, and we construct an example of parameter optimization based on medical images. We discuss how the presented approach might be employed to process large datasets of medical images as well as databases of echocardiographic reports, paving the way towards application of heart mechanics models in the clinical practice.
Stable model order reduction for time-domain exterior vibro-acoustic finite element simulations
