PIs: Christian Wieners, Axel Loewe
Aim:
This project aims to extend the limits of multi-physics simulations of the human heart by developing and evaluating robust coupling schemes for models of the cardiac electro-mechanical system with high biophysical accuracy across multiple scales and dimensions. By combining our profound interdisciplinary experience in modeling each of the cardiac physical systems separately with our preliminary work in coupling them, we are well positioned to construct and analyze an efficient and accurate parallel finiteelement realization of the fully coupled system.
As specific goals, we will develop and provide
(1) a formulation for a robustly coupled model for the whole heart comprising active electrophysiology, cardiac mechanics and circulatory function;
2) an efficient and scalable numerical implementation of this model and a comprehensive convergence analysis of the approximations in space and time;
(3) a user-friendly integration of the fully coupled model into the established openCARP simulation framework to foster community adoption of this research tool and pave the way for clinical translation.
Description:
Our multi-physics model consists of several model components that describe the electromechanical contraction of the heart at different spatial scales. Ion concentrations within the cell are modeled via ion channels of the cardiomyocyte cell membrane. Based on the resulting calcium concentration, a contraction force is obtained using a sarcomere model. The activation of the myocardium is based on the electrical propagation across the organ, which – like the deformation of the heart – is calculated using the finite element method. A circulatory model is coupled to the deformation of the heart via the volumes and pressures of the atria and ventricles.
Building on our current work, we will further develop our multi-physics model by extending it with stable and efficient implicit-explicit time integration and finite element formulations and making the fully coupled system scalable. Subsequently, we will use this advanced model to study the effects of local electrical dyssynchrony and the application of cardiac resynchronization therapy in collaboration with our clinical partners. In addition, we will identify application scenarios where the highly complex full model can be reduced and quantify the resulting uncertainty.
Involved Institutions:
Karlsruhe Institute of Technology
Institute of Biomedical Engineering
Computational Cardiac Modeling
Karlsruhe Institute of Technology
Institute for Applied and Numerical Mathematics
Scientific Computing
Applicants:
Publications
Gerach, Tobias; Loewe, Axel
In: The Journal of Physiology, 2024.
@article{Gerach2024,
title = {Differential effects of mechano-electric feedback mechanisms on whole-heart activation, repolarization, and tension},
author = {Tobias Gerach and Axel Loewe},
editor = {John Wiley The Physiological Society},
url = {https://physoc.onlinelibrary.wiley.com/doi/10.1113/JP285022},
doi = {https://doi.org/10.1113/JP285022},
year = {2024},
date = {2024-01-07},
urldate = {2024-01-07},
journal = {The Journal of Physiology},
abstract = {The human heart is subject to highly variable amounts of strain during day-to-day activities and needs to adapt to a wide range of physiological demands. This adaptation is driven by an autoregulatory loop that includes both electrical and the mechanical components. In particular, mechanical forces are known to feed back into the cardiac electrophysiology system, which can result in pro- and anti-arrhythmic effects. Despite the widespread use of computational modelling and simulation for cardiac electrophysiology research, the majority of in silico experiments ignore this mechano-electric feedback entirely due to the high computational cost associated with solving cardiac mechanics. In this study, we therefore use an electromechanically coupled whole-heart model to investigate the differential and combined effects of electromechanical feedback mechanisms with a focus on their physiological relevance during sinus rhythm. In particular, we consider troponin-bound calcium, the effect of deformation on the tissue diffusion tensor, and stretch-activated channels. We found that activation of the myocardium was only significantly affected when including deformation into the diffusion term of the monodomain equation. Repolarization, on the other hand, was influenced by both troponin-bound calcium and stretch-activated channels and resulted in steeper repolarization gradients in the atria. The latter also caused afterdepolarizations in the atria. Due to its central role for tension development, calcium bound to troponin affected stroke volume and pressure. In conclusion, we found that mechano-electric feedback changes activation and repolarization patterns throughout the heart during sinus rhythm and lead to a markedly more heterogeneous electrophysiological substrate.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Fröhlich, Jonathan; Gerach, Tobias; Krauß, Jonathan; Loewe, Axel; Stengel, Laura; Wieners, Christian
In: GAMM Mitteilungen, Bd. 46, Ausg. 2, 2024, ISSN: 1522-2608.
@article{Fröhlich2024,
title = {Numerical evaluation of elasto-mechanical and visco-elastic electro-mechanical models of the human heart},
author = {Jonathan Fröhlich and Tobias Gerach and Jonathan Krauß and Axel Loewe and Laura Stengel and Christian Wieners},
editor = {Gesellschaft Mechanik},
url = {https://onlinelibrary.wiley.com/doi/10.1002/gamm.202370010},
doi = {https://doi.org/10.1002/gamm.202370010},
issn = {1522-2608},
year = {2024},
date = {2024-01-02},
urldate = {2024-01-02},
journal = {GAMM Mitteilungen},
volume = {46},
issue = {2},
abstract = {We investigate the properties of static mechanical and dynamic electro-mechanical models for the deformation of the human heart. Numerically this is realized by a staggered scheme for the coupled partial/ordinary differential equation (PDE-ODE) system. First, we consider a static and purely mechanical benchmark configuration on a realistic geometry of the human ventricles. Using a penalty term for quasi-incompressibility, we test different parameters and mesh sizes and observe that this approach is not sufficient for lowest order conforming finite elements. Then, we compare the approaches of active stress and active strain for cardiac muscle contraction. Finally, we compare in a coupled anatomically realistic electro-mechanical model numerical Newmark damping with a visco-elastic model using Rayleigh damping. Nonphysiological oscillations can be better mitigated using viscosity.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}