zbMATH — the first resource for mathematics

Computational modeling of cardiac electrophysiology: a novel finite element approach. (English) Zbl 1171.92310
Summary: The key objective of this work is the design of an unconditionally stable, robust, efficient, modular, and easily expandable finite element-based simulation tool for cardiac electrophysiology. In contrast to existing formulations, we propose a global-local split of the system of equations in which the global variable is the fast action potential that is introduced as a nodal degree of freedom, whereas the local variable is the slow recovery variable introduced as an internal variable on the integration point level. Cell-specific excitation characteristics are thus strictly local and only affect the constitutive level. We illustrate the modular character of the model in terms of the FitzHugh-Nagumo model for oscillatory pacemaker cells and the Aliev-Panfilov model for non-oscillatory ventricular muscle cells. We apply an implicit Euler backward finite difference scheme for the temporal discretization and a finite element scheme for the spatial discretization. The resulting non-linear system of equations is solved with an incremental iterative Newton-Raphson solution procedure. Since this framework only introduces one single scalar-valued variable on the node level, it is extremely efficient, remarkably stable, and highly robust. The features of the general framework will be demonstrated by selected benchmark problems for cardiac physiology and a two-dimensional patient-specific cardiac excitation problem.

92C30 Physiology (general)
92C50 Medical applications (general)
68U20 Simulation (MSC2010)
Full Text: DOI
[1] Ganong, Review of Medical Physiology (2003)
[2] Berne, Cardiovascular Physiology (2001)
[3] Bers, Excitation-contraction Coupling and Cardiac Contractile Force (2001) · doi:10.1007/978-94-010-0658-3
[4] Durrer, Total excitation of the isolated heart, Circulation 41 pp 899– (1970) · doi:10.1161/01.CIR.41.6.899
[5] Opie, Heart Physiology: From Cell to Circulation (2003)
[6] Klabunde, Cardiovascular Physiology Concepts (2005)
[7] Plonsey, Bioelectricity. A Quantitative Approach (2007)
[8] Hodgkin, A quantitative description of membrane current and its application to conductance and excitation in nerve, Journal of Physiology 117 pp 500– (1952) · doi:10.1113/jphysiol.1952.sp004764
[9] FitzHugh, Impulses and physiological states in theoretical models of nerve membranes, Biophysical Journal 1 pp 445– (1961)
[10] Noble, A modification of the Hodgkin-Huxley equations applicable to Purkinje fibre action and pacemaker potentials, Journal of Physiology 160 pp 317– (1962) · doi:10.1113/jphysiol.1962.sp006849
[11] Beeler, Reconstruction of the action potential of ventricular myocardial fibers, Journal of Physiology 268 pp 177– (1977) · doi:10.1113/jphysiol.1977.sp011853
[12] Luo, A model of the ventricular cardiac action potential. Depolarization, repolarization, and their changes, Circulation Research 68 pp 1501– (1991) · doi:10.1161/01.RES.68.6.1501
[13] Keener, Mathematical Physiology (2004)
[14] Sachse, Computational Cardiology (2004) · doi:10.1007/b96841
[15] Clayton, A guide to modelling cardiac electrical activity in anatomically detailed ventricles, Progress in Biophysics and Molecular Biology 96 pp 19– (2008)
[16] ten Tuscher, Modelling of the ventricular conduction system, Progress in Biophysics and Molecular Biology 96 pp 152– (2008)
[17] Fenton, Multiple mechanisms of spiral wave breakup in a model of cardiac electrical activity, Chaos 12 pp 852– (2002)
[18] Aliev, A simple two-variable model of cardiac excitation, Chaos, Solitons and Fractals 7 pp 293– (1996)
[19] Fenton, Vortex dynamics in three-dimensional continuous myocardium with fiber rotation: filament instability and fibrillation, Chaos 8 pp 20– (1998) · Zbl 1069.92503
[20] Rogers, A collocation-Galerkin finite element model of cardiac action potential propagation, IEEE Transactions on Biomedical Engineering 41 pp 743– (1994)
[21] Rogers, Nonuniform muscle fiber orientation causes spiral wave drift in a finite element model of cardiac action potential propagation, Journal of Cardiovascular Electrophysiology 5 pp 496– (1994)
[22] Rogers, Wave front fragmentation due to ventricular geometry in a model of the rabbit heart, Chaos 12 pp 779– (2002)
[23] Hunter, Modelling the mechanical properties of cardiac muscle, Progress in Biophysics and Molecular Biology 69 pp 289– (1998)
[24] Nash, Electromechanical model of excitable tissue to study reentrant cardiac arrhytmias, Progress in Biophysics and Molecular Biology 85 pp 501– (2004)
[25] Panfilov, Self-organized pacemakers in a coupled reaction-diffusion-mechanics system, Physical Review Letters 95 pp 258104:1– (2005)
[26] Nickerson, Computational multiscale modeling in the IUPS physiome project: modeling cardiac electromechanics, IBM Journal of Research and Development 50 pp 617– (2006)
[27] Keldermann, Pacemakers in a reaction-diffusion mechanics system, Journal of Statistical Physics 128 pp 375– (2007) · Zbl 1115.92004
[28] Panfilov, Drift and breakup of spiral waves in reaction-diffusion-mechanics systems, Proceedings of the National Academy of Sciences 104 pp 7922– (2007)
[29] Tomlinson, A finite element method for an Eikonal equation model of myocardial excitation wavefront propagation, SIAM Journal on Applied Mathematics 63 pp 324– (2002) · Zbl 1014.92010
[30] van der Pol, On relaxation oscillations, Philosophical Magazine 2 pp 978– (1926) · JFM 52.0450.05 · doi:10.1080/14786442608564127
[31] Nagumo, Active pulse transmission line simulating nerve axon, Proceedings of the Institute of Radio Engineers 50 pp 2061– (1962)
[32] Gould, Biventricular pacing in heart failure: a review, Expert Review of Cardiovascular Therapy 4 pp 97– (2006)
[33] Hawkins, Selecting patients for cardiac resynchronization therapy: electrical or mechanical dyssynchronization, European Heart Journal 27 pp 1270– (2006)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.