A Meshfree Method for Simulating Myocardial Electrical Activity

An element-free Galerkin method (EFGM) is proposed to simulate the propagation of myocardial electrical activation without explicit mesh constraints using a monodomain model. In our framework the geometry of myocardium is first defined by a meshfree particle representation that is, a sufficient number of sample nodes without explicit connectivities are placed in and inside the surface of myocardium. Fiber orientations and other material properties of myocardium are then attached to sample nodes according to their geometrical locations, and over the meshfree particle representation spatial variation of these properties is approximated using the shape function of EFGM. After the monodomain equations are converted to their Galerkin weak form and solved using EFGM, the propagation of myocardial activation can be simulated over the meshfree particle representation. The derivation of this solution technique is presented along a series of numerical experiments and a solution of monodomain model using a FitzHugh-Nagumo (FHN) membrane model in a canine ventricular model and a human-heart model which is constructed from digitized virtual Chinese dataset.

[1]  Chris P. Bradley,et al.  A coupled cubic hermite finite element/boundary element procedure for electrocardiographic problems , 1996 .

[2]  Indra Vir Singh,et al.  Parallel implementation of the EFG Method for heat transfer and fluid flow problems , 2004 .

[3]  T. Belytschko,et al.  Element‐free Galerkin methods , 1994 .

[4]  P. Hunter,et al.  A generalized finite difference method for modeling cardiac electrical activation on arbitrary, irregular computational meshes. , 2005, Mathematical biosciences.

[5]  Huafeng Liu,et al.  Noninvasive Computational Imaging of Cardiac Electrophysiology for 3-D Infarct , 2011, IEEE Transactions on Biomedical Engineering.

[6]  Huafeng Liu,et al.  Maximum a Posteriori Strategy for the Simultaneous Motion and Material Property Estimation of the Heart , 2009, IEEE Transactions on Biomedical Engineering.

[7]  J. Keener.  An eikonal-curvature equation for action potential propagation in myocardium , 1991, Journal of mathematical biology.

[8]  Ted Belytschko,et al.  Numerical integration of the Galerkin weak form in meshfree methods , 1999 .

[9]  P J Hunter,et al.  Analytical models of propagation in excitable cells. , 1975, Progress in biophysics and molecular biology.

[10]  Andrew J. Pullan,et al.  Mathematically Modelling the Electrical Activity of the Heart: From Cell to Body Surface and Back Again , 2005 .

[11]  Tong Wang,et al.  [Research report of experimental database establishment of digitized virtual Chinese No.1 female]. , 2003, Di 1 jun yi da xue xue bao = Academic journal of the first medical college of PLA.

[12]  Bernice E. Rogowitz,et al.  Mapping, modeling, and visual exploration of structure-function relationships in the heart , 2001, IBM Syst. J..

[13]  D. Durrer,et al.  Total Excitation of the Isolated Human Heart , 1970, Circulation.

[14]  J. Westwater,et al.  The Mathematics of Diffusion. , 1957 .

[15]  A. M. Scher,et al.  Influence of Cardiac Fiber Orientation on Wavefront Voltage, Conduction Velocity, and Tissue Resistivity in the Dog , 1979, Circulation research.

[16]  C S Henriquez,et al.  A space-time adaptive method for simulating complex cardiac dynamics. , 2000, Physical review letters.

[17]  Pengcheng Shi,et al.  An H(∞) approach for elasticity properties reconstruction. , 2012, Medical physics.

[18]  P Le Beux,et al.  An interactive 3D anisotropic cellular automata model of the heart. , 1998, Computers and biomedical research, an international journal.

[19]  Huafeng Liu,et al.  Physiological Fusion of Functional and Structural Images for Cardiac Deformation Recovery , 2011, IEEE Transactions on Medical Imaging.

[20]  J L Lehr,et al.  Effects of paddle placement and size on defibrillation current distribution: a three-dimensional finite element model. , 1993, IEEE transactions on bio-medical engineering.

[21]  Huafeng Liu,et al.  Stochastic finite element framework for simultaneous estimation of cardiac kinematic functions and material parameters , 2003, Medical Image Anal..

[22]  A. McCulloch,et al.  A collocation-Galerkin finite element model of cardiac action potential propagation , 1994, IEEE Transactions on Biomedical Engineering.

[23]  Mark A Fleming,et al.  ENRICHED ELEMENT-FREE GALERKIN METHODS FOR CRACK TIP FIELDS , 1997 .

[24]  Andrew J. Pullan,et al.  A Finite Element Method for an Eikonal Equation Model of Myocardial Excitation Wavefront Propagation , 2002, SIAM J. Appl. Math..

[25]  Pengcheng Shi,et al.  Physiome-model-based state-space framework for cardiac deformation recovery. , 2007, Academic radiology.

[26]  Santiago Falcón,et al.  Penalty functions in constrained variational principles for element free Galerkin method , 2000 .

[27]  Guirong Liu.  Mesh Free Methods: Moving Beyond the Finite Element Method , 2002 .

[28]  Frank B. Sachse,et al.  Computational Cardiology , 2004, Lecture Notes in Computer Science.

[29]  J. Lehr,et al.  Effect of paddle placement and size on defibrillation current distribution: a three-dimensional finite element model , 1993, IEEE Transactions on Biomedical Engineering.

[30]  H. Raper,et al.  Principles of Human Physiology , 1946, Nature.

[31]  T. Belytschko,et al.  Adaptivity for structured meshfree particle methods in 2D and 3D , 2005 .

[32]  Nicolas P. Smith,et al.  Mathematical modelling of the ischaemic heart , 2001 .

[33]  C R Johnson,et al.  A computer model for the study of electrical current flow in the human thorax. , 1992, Computers in biology and medicine.

[34]  Pengcheng Shi,et al.  A Meshfree Method for Solving Cardiac Electrical Propagation , 2005, 2005 IEEE Engineering in Medicine and Biology 27th Annual Conference.

[35]  P. C. Franzone,et al.  Spreading of excitation in 3-D models of the anisotropic cardiac tissue. I. Validation of the eikonal model. , 1993, Mathematical biosciences.

[36]  Huafeng Liu,et al.  State-Space Analysis of Cardiac Motion With Biomechanical Constraints , 2007, IEEE Transactions on Image Processing.

[37]  Mark A Fleming,et al.  Meshless methods: An overview and recent developments , 1996 .