Simulating cardiac electrophysiology using anisotropic mesh adaptivity

Abstract The simulation of cardiac electrophysiology requires small time steps and a fine mesh in order to resolve very sharp, but highly localized, wavefronts. The use of very high resolution meshes containing large numbers of nodes results in a high computational cost, both in terms of CPU hours and memory footprint. In this paper an anisotropic mesh adaptivity technique is implemented in the Chaste physiological simulation library in order to reduce the mesh resolution away from the depolarization front. Adapting the mesh results in a reduction in the number of degrees of freedom of the system to be solved by an order of magnitude during propagation and 2–3 orders of magnitude in the subsequent plateau phase. As a result, a computational speedup by a factor of between 5 and 12 has been obtained with no loss of accuracy, both in a slab-like geometry and for a realistic heart mesh with a spatial resolution of 0.125 mm.

[1]  G Plank,et al.  Computational tools for modeling electrical activity in cardiac tissue. , 2003, Journal of electrocardiology.

[2]  Matthew D. Piggott,et al.  Conservative interpolation between unstructured meshes via supermesh construction , 2009 .

[3]  C. Henriquez,et al.  A finite volume model of cardiac propagation , 1997, Annals of Biomedical Engineering.

[4]  C. C. Pain,et al.  Adjoint goal-based error norms for adaptive mesh ocean modelling , 2006 .

[5]  David Gavaghan,et al.  Generation of histo-anatomically representative models of the individual heart: tools and application , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[6]  P. Hunter,et al.  Integration from proteins to organs: the Physiome Project , 2003, Nature Reviews Molecular Cell Biology.

[7]  Patrick E. Farrell,et al.  Galerkin projection of discrete fields via supermesh construction , 2009 .

[8]  Bruce Smaill,et al.  A Finite Volume Method for Modeling Discontinuous Electrical Activation in Cardiac Tissue , 2005, Annals of Biomedical Engineering.

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

[10]  Elizabeth M Cherry,et al.  Efficient simulation of three-dimensional anisotropic cardiac tissue using an adaptive mesh refinement method. , 2003, Chaos.

[11]  Alexander G. Fletcher,et al.  Chaste: A test-driven approach to software development for biological modelling , 2009, Comput. Phys. Commun..

[12]  James P. Keener,et al.  Mathematical physiology , 1998 .

[13]  J. Trangenstein,et al.  Operator splitting and adaptive mesh refinement for the Luo-Rudy I model , 2004 .

[14]  Natalia A. Trayanova,et al.  Computational techniques for solving the bidomain equations in three dimensions , 2002, IEEE Transactions on Biomedical Engineering.

[15]  Gustavo C. Buscaglia,et al.  Anisotropic mesh optimization and its application in adaptivity , 1997 .

[16]  Pras Pathmanathan,et al.  Chaste: using agile programming techniques to develop computational biology software , 2008, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[17]  Jonathan P. Whiteley,et al.  Physiology Driven Adaptivity for the Numerical Solution of the Bidomain Equations , 2007, Annals of Biomedical Engineering.

[18]  P. Deuflhard,et al.  Adaptive finite element simulation of ventricular fibrillation dynamics , 2009 .

[19]  C. Luo,et al.  A model of the ventricular cardiac action potential. Depolarization, repolarization, and their interaction. , 1991, Circulation research.

[20]  Jens Lang,et al.  Konrad-zuse-zentrum F ¨ Ur Informationstechnik Berlin Adaptivity in Space and Time for Reaction-diffusion Systems in Electrocardiology Adaptivity in Space and Time for Reaction-diffusion Systems in Electrocardiology , 2022 .

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

[22]  J. Keener,et al.  A numerical method for the solution of the bidomain equations in cardiac tissue. , 1998, Chaos.

[23]  R. B. Simpson.  Anisotropic mesh transformations and optimal error control , 1994 .

[24]  G Plank,et al.  Solvers for the cardiac bidomain equations. , 2008, Progress in biophysics and molecular biology.

[25]  Frédéric Hecht,et al.  Anisotropic unstructured mesh adaption for flow simulations , 1997 .

[26]  O. C. Zienkiewicz,et al.  Adaptive remeshing for compressible flow computations , 1987 .

[27]  C. C. Pain,et al.  h, r, and hr adaptivity with applications in numerical ocean modelling , 2005 .

[28]  Paul-Louis George,et al.  Delaunay triangulation and meshing : application to finite elements , 1998 .

[29]  C.R.E. de Oliveira,et al.  Tetrahedral mesh optimisation and adaptivity for steady-state and transient finite element calculations , 2001 .

[30]  Youssef Belhamadia,et al.  A Time-Dependent Adaptive Remeshing for Electrical Waves of the Heart , 2008, IEEE Transactions on Biomedical Engineering.

[31]  Youssef Belhamadia,et al.  Towards accurate numerical method for monodomain models using a realistic heart geometry. , 2009, Mathematical biosciences.

[32]  Pascal Frey,et al.  Anisotropic mesh adaptation for CFD computations , 2005 .