An oscillation-free fully staggered algorithm for velocity-dependent active models of cardiac mechanics

In this paper we address an unresolved problem in the numerical modeling of cardiac electromechanics, that is the onset of numerical oscillations due to the dependence of force generation models on the fibers shortening velocity. A way to avoid numerical oscillations is to use monolithic schemes for the solution of the coupled problem of active-passive mechanics. However, staggered strategies, which foresee the sequential solution of the models of force generation and of tissue mechanics, are preferable, due to their reduced computational cost and low implementation effort. In this paper we propose a cure for this issue, by introducing, with respect to the standard staggered scheme, a numerically consistent stabilization term. This term is derived in virtue of the identification of the cause of instability in the mismatch between macroscopic and microscopic strains, inconsistently expressed in Lagrangian and Eulerian coordinates, respectively. By considering a model problem of active mechanics we prove that the proposed scheme is unconditionally absolutely stable (i.e. it is stable for any time step size), yet within a fully staggered framework. As such, the new scheme removes the non-physical oscillations, as we prove by applying it to three force generation models, namely the Niederer-Hunter-Smith model, the model by Land and coworkers, and the mean-field force generation model that we have recently proposed. (C) 2020 The Author(s). Published by Elsevier B.V.

Automated summary

This paper addresses an unresolved problem in the numerical modeling of cardiac electromechanics, namely the onset of numerical oscillations due to the dependence of force generation models on the fibers shortening velocity. The proposed scheme introduces a numerically consistent stabilization term to remove non-physical oscillations in a fully staggered framework, proving its stability for any time step size and successful application to three force generation models.