Next-generation In Silico Cardiac Electrophysiology Through Immersed Grid Meshfree Modeling: Application To Simulation Of Myocardial Infarction

State-of-art simulators of cardiac tissue electrophysiology are commonly based on the Finite Element Method (FEM). FEM is known to be a robust and accurate numerical method, but its accuracy highly depends on the quality of the mesh. Generating a good-quality mesh may be cumbersome and time consuming for models with complex geometries, such as those representing the anatomy of human organs. This limitation restricts the clinical application of FEM. To overcome this challenge, we propose the use of a meshfree method, the Moving Kriging Mixed Collocation (MKMC) method for in silico cardiac electrophysiology. MKMC is a purely meshfree method requiring the definition of a point cloud rather than a mesh. We propose the construction of the point cloud as an immersed grid of points generated automatically from image data. In simulations on a swine biventricular model, both under baseline and myocardial infarction conditions, we demonstrate the capability of the MKMC method to generate results in very good agreement with FEM while alleviating the mesh requirement. Differences in local activation time and action potential duration between MKMC and FEM are in mean below 3%. The proposed MKMC method represents a promising alternative to FEM for cardiac in silico investigations with the potential to be integrated in the clinic.

[1]  Zhilin Qu,et al.  Cardiac fibrosis and arrhythmogenesis: the road to repair is paved with perils. , 2014, Journal of molecular and cellular cardiology.

[2]  A. Garfinkel,et al.  An advanced algorithm for solving partial differential equation in cardiac conduction , 1999, IEEE Transactions on Biomedical Engineering.

[3]  Konstantinos A. Mountris,et al.  The Radial Point Interpolation Mixed Collocation (RPIMC) Method for the Solution of Transient Diffusion Problems , 2020, Engineering Analysis with Boundary Elements.

[4]  Ellen Kuhl,et al.  Generating Purkinje networks in the human heart. , 2016, Journal of biomechanics.

[5]  Nicolas Duchateau,et al.  Development of a Swine Model of Left Bundle Branch Block for Experimental Studies of Cardiac Resynchronization Therapy , 2013, Journal of Cardiovascular Translational Research.

[6]  Bernadette A. Thomas,et al.  Global, regional, and national age–sex specific all-cause and cause-specific mortality for 240 causes of death, 1990–2013: a systematic analysis for the Global Burden of Disease Study 2013 , 2015, The Lancet.

[7]  Oscar Camara,et al.  Smoothed Particle Hydrodynamics for Electrophysiological Modeling: An Alternative to Finite Element Methods , 2017, FIMH.

[8]  Gabriel Bernardino,et al.  A rule‐based method to model myocardial fiber orientation in cardiac biventricular geometries with outflow tracts , 2018, International journal for numerical methods in biomedical engineering.

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

[10]  L. Gu,et al.  Moving kriging interpolation and element‐free Galerkin method , 2003 .

[11]  Martin Bishop,et al.  Personalized Cardiac Computational Models: From Clinical Data to Simulation of Infarct-Related Ventricular Tachycardia , 2019, Front. Physiol..

[12]  Esther Pueyo,et al.  A Novel Paradigm for In Silico Simulation of Cardiac Electrophysiology Through the Mixed Collocation Meshless Petrov-Galerkin Method , 2019, 2019 Computing in Cardiology (CinC).