Numerical simulation of mathematical heart model in COMSOL

  • Abstract
  • Keywords
  • References
  • PDF
  • Abstract

    Electrical activity is essential for the cardiac cell to perform its function. Mathematical modeling of cardiac electrical activity is performed from the cell, tissue and organ levels through to the body surface level. The electrical activity of the cardiac as a whole is thus characterized by a complex multiscale structure. The most complete model of such a complex setting is the anisotropic bidomain model that consists of a system of two degenerate parabolic reaction diffusion equations describing the intra and extracellular potentials in the cardiac muscle, coupled with a system of ordinary deferential equations describing the ionic currents flowing through the cellular membrane. This study describes an anatomically realistic 3D Bidomain model of whole-heart electrical activity. The heart was embedded in a human torso, incorporating spontaneous activation with heterogeneous action potential (AP) morphologies throughout the heart. The aim of this study is the development of a geometrically simple and computationally efficient 3D model of heart. In this paper a finite element formulation, model and simulation of Bidomain equation has been conducted. The FitzHugh-Nagumo (FHN) equations were incorporated into Bidomain model of cardiac electrical activity, which was comprised of a simplified geometry of the whole heart with the torso as an extracellular volume conductor. Laplace equation for the torso also considered. Simplified 3D cardiac model was implemented using COMSOL Multiphysics 5.0 finite element software. Electrical potential at different point on torso is measured. Propagation of electrical excitation on heart surface is also observed. This study represents the first stage toward the development of an accurate computer model of heart activation.




  • Keywords

    Bidomain Model; Electrical activity; FHN equation; Finite Element Method; Mathematical Model.

  • References

      [1] A. L. Hodgkin and A. F. Huxley, A quantitative description of membrane current and its application to conduction and excitation in nerve. J. Physiol. (Lond), 117(4):500–544, 1952.

      [2] Jacquemet V, van Oosterom A, Vesin J-M, Kappenberger L, Analysis of electrocardiograms during atrial fibrillation. Engineering in Medicine and Biology, p: 79-88, Magazine, IEEE 2006.

      [3] Malmivuo J, Plonsey R, Bioelectromagnetism: principles and applications of bioelectric and biomagnetic fields, Oxford University Press; 1995.

      [4] C. S. Henriquez, Simulating the electrical behavior of cardiac tissue using the Bidomain model, Crit. Rev. Biomedical Engineering 21, 1-77 (1993).

      [5] J.M.Rogers and A.D.McCulloch. A collocation – Galerkin finite element model of cardiac action potential propagation, IEEE Trans. Biomedical Engineering 41, 743-57, 1994.

      [6] Pullan AJ, Buist ML, Cheng LK, mathematically modeling the electrical activity of the heart: From Cell to Body Surface and Back Again, World Scientific; 2005.

      [7] Daryl L. Logan, AFirst Course in the Finite Element Method Fourth Edition,

      [8] J. Sundnes, G. T. Lines, X. Cai, B. F. Nielsen, K.-A. Mardal, and A. Tveito. Computing the electrical activity in the heart. Springer-Verlag, Berlin, 2006.

      [9] Socrates Dokos, Modeling Organs, Tissues, cells and Devices using MATLAB and COMSOL Multiphysics, Springer Publication 2016.

      [10] Socrates Dokos, Shaun L. Cloherty, Computational Model of Atrial Electrical Activation and Propagation, 29th Annual International Conference of the IEEE EMBS, 2007.




Article ID: 30052
DOI: 10.14419/ijet.v9i1.30052

Copyright © 2012-2015 Science Publishing Corporation Inc. All rights reserved.