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 t...

Full description

Saved in:
Bibliographic Details
Published in:Computer methods in biomechanics and biomedical engineering 2014, Vol.17 (6), p.604-615
Main Authors: Sundnes, J., Wall, S., Osnes, H., Thorvaldsen, T., McCulloch, A.D.
Format: Article
Language:English
Subjects:
Algorithms
cardiac electro-mechanics
Computer Simulation
Elasticity
Electrophysiological Phenomena
Heart - physiology
Models, Cardiovascular
Myocardial Contraction
operator splitting
strongly coupled simulation
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary: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.
ISSN:1025-5842
1476-8259
DOI:10.1080/10255842.2012.704368