Adaptive Mesh Refinement and Adaptive Time Integration for Electrical Wave Propagation on the Purkinje System

A both space and time adaptive algorithm is presented for simulating electrical wave propagation in the Purkinje system of the heart. The equations governing the distribution of electric potential over the system are solved in time with the method of lines. At each timestep, by an operator splitting technique, the space-dependent but linear diffusion part and the nonlinear but space-independent reactions part in the partial differential equations are integrated separately with implicit schemes, which have better stability and allow larger timesteps than explicit ones. The linear diffusion equation on each edge of the system is spatially discretized with the continuous piecewise linear finite element method. The adaptive algorithm can automatically recognize when and where the electrical wave starts to leave or enter the computational domain due to external current/voltage stimulation, self-excitation, or local change of membrane properties. Numerical examples demonstrating efficiency and accuracy of the adaptive algorithm are presented.

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

[2]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[3]  Craig S. Henriquez,et al.  A Brief History of Tissue Models For Cardiac Electrophysiology , 2014, IEEE Transactions on Biomedical Engineering.

[4]  Kenneth B. Huber.  Department of Mathematics , 1894 .

[5]  C. Brezinski,et al.  Extrapolation methods , 1992 .

[6]  D. Fayuk,et al.  The Journal of Physiology , 1978, Medical History.

[7]  Wenjun Ying,et al.  A multilevel adaptive approach for computational cardiology , 2005 .

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

[9]  J. Gillis,et al.  Matrix Iterative Analysis , 1961 .

[10]  Morten Hjorth Jensen,et al.  COMPUTATIONAL Physics , 2015 .

[11]  Wenjun Ying,et al.  Efficient Fully Implicit Time Integration Methods for Modeling Cardiac Dynamics , 2008, IEEE Transactions on Biomedical Engineering.

[12]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

[13]  A. Tveito,et al.  Numerical solution of the bidomain equations , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[14]  J. Rogers.  Chaos , 1876 .

[15]  Alan Garny,et al.  A numerical guide to the solution of the bi-domain equations of cardiac electrophysiology. , 2010, Progress in biophysics and molecular biology.

[16]  G. W. Beeler,et al.  Reconstruction of the action potential of ventricular myocardial fibres , 1977, The Journal of physiology.

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

[18]  A. Hodgkin,et al.  A quantitative description of membrane current and its application to conduction and excitation in nerve , 1952, The Journal of physiology.

[19]  October I.  Physical Review Letters , 2022 .

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

[21]  Edward J. Vigmond,et al.  Construction of a Computer Model to Investigate Sawtooth Effects in the Purkinje System , 2007, IEEE Transactions on Biomedical Engineering.

[22]  M. Berger,et al.  Adaptive mesh refinement for hyperbolic partial differential equations , 1982 .

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