Efficient integration of a realistic two-dimensional cardiac tissue model by domain decomposition

The size of realistic cardiac tissue models has been limited by their high computational demands. In particular, the Luo-Rudy phase II membrane model, used to simulate a thin sheet of ventricular tissue with arrays of coupled ventricular myocytes, is usually limited to 100/spl times/100 arrays. The authors introduce a new numerical method based on domain decomposition and a priority queue integration scheme which reduces the computational cost by a factor of 3-17. In the standard algorithm all the nodes advance with the same time step /spl Delta/t, whose size Is limited by the time scale of activation. However, at any given time, many regions may he inactive and do not require the same small /spl Delta/t and consequent extensive computations. Hence, adjusting /spl Delta/t locally is a key factor in improving computational efficiency, since most of the computing time is spent calculating ionic currents. This paper proposes an efficient adaptive numerical scheme for integrating a two-dimensional (2-D) propagation model, by incorporating local adjustments of /spl Delta/t. In this method, alternating direction Cooley-Dodge and Rush-Larsen methods were used for numerical integration. Between consecutive integrations over the whole domain using an implicit method, the model was spatially decomposed into many subdomains, and /spl Delta/t adjusted locally. The Euler method was used for numerical integration in the subdomains. Local boundary values were determined from the boundary mesh elements of the neighboring subdomains using linear interpolation. Because /spl Delta/t was defined locally, a priority queue was used to store and order nest update times for each subdomain. The subdomain with the earliest update time was given the highest priority and advanced first. This new method yielded stable solutions with relative errors less than 1% and reduced computation time by a factor of 3-17 and will allow much larger (e.g. 500/spl times/500) models based on realistic membrane kinetics and realistic dimensions to simulate reentry, triggered activity, and their interactions.

[1]  Bruce M. Steinhaus,et al.  Action Potential Collision in Heart Tissue-Computer Simulations and Tissue Expenrments , 1985, IEEE Transactions on Biomedical Engineering.

[2]  Ronald W. Joyner,et al.  Simulation of Action Potential Propagation in an Inhomogeneous Sheet of Coupled Excitable Cells , 1975, Circulation research.

[3]  W. Baxter,et al.  Spiral waves of excitation underlie reentrant activity in isolated cardiac muscle. , 1993, Circulation research.

[4]  C. Henriquez,et al.  Cardiac propagation simulation. , 1992, Critical reviews in biomedical engineering.

[5]  M. Fishbein,et al.  Anisotropic repolarization in ventricular tissue. , 1997, The American journal of physiology.

[6]  Michael D. Lesh,et al.  Cellular Uncoupling Can Unmask Dispersion of Action Potential Duration in Ventricular Myocardium A Computer Modeling Study , 1989, Circulation research.

[7]  R. Barr,et al.  Propagation of excitation in idealized anisotropic two-dimensional tissue. , 1984, Biophysical journal.

[8]  F A Roberge,et al.  Directional characteristics of action potential propagation in cardiac muscle. A model study. , 1991, Circulation research.

[9]  S. Rush,et al.  A Practical Algorithm for Solving Dynamic Membrane Equations , 1978, IEEE Transactions on Biomedical Engineering.

[10]  D Durrer,et al.  Computer Simulation of Arrhythmias in a Network of Coupled Excitable Elements , 1980, Circulation research.

[11]  D. Noble,et al.  A model of sino-atrial node electrical activity based on a modification of the DiFrancesco-Noble (1984) equations , 1984, Proceedings of the Royal Society of London. Series B. Biological Sciences.

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

[13]  Walmor C. De Mello,et al.  Chapter 12 – The Healing-Over Process in Cardiac and Other Muscle Fibers , 1972 .

[14]  J. Cooley,et al.  Digital computer solutions for excitation and propagation of the nerve impulse. , 1966, Biophysical journal.

[15]  Richard L. Burden,et al.  Numerical analysis: 4th ed , 1988 .

[16]  H. H. Rachford,et al.  The Numerical Solution of Parabolic and Elliptic Differential Equations , 1955 .

[17]  F A Roberge,et al.  Reconstruction of Propagated Electrical Activity with a Two‐Dimensional Model of Anisotropic Heart Muscle , 1986, Circulation research.

[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]  R. W. Joyner,et al.  Effects of the Discrete Pattern of Electrical Coupling on Propagation through an Electrical Syncytium , 1982, Circulation research.

[20]  B. Victorri,et al.  Numerical integration in the reconstruction of cardiac action potentials using Hodgkin-Huxley-type models. , 1985, Computers and biomedical research, an international journal.

[21]  D DiFrancesco,et al.  A model of cardiac electrical activity incorporating ionic pumps and concentration changes. , 1985, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[22]  Yoram Rudy,et al.  A Model Study of the Effects of the Discrete Cellular Structure on Electrical Propagation in Cardiac Tissue , 1987, Circulation research.

[23]  Gul'ko Fb,et al.  Mechanism of formation of closed propagation pathways in excitable media , 1972 .

[24]  C. S. Henriquez.  An examination of a computationally efficient algorithm for modeling propagation in cardiac tissue , 1989 .

[25]  F. Roberge,et al.  Structural complexity effects on transverse propagation in a two-dimensional model of myocardium , 1991, IEEE Transactions on Biomedical Engineering.

[26]  C. Luo,et al.  A dynamic model of the cardiac ventricular action potential. I. Simulations of ionic currents and concentration changes. , 1994, Circulation research.

[27]  C. Canuto.  Spectral methods in fluid dynamics , 1991 .

[28]  J. Clark,et al.  A mathematical model of electrophysiological activity in a bullfrog atrial cell. , 1990, The American journal of physiology.

[29]  C. Luo,et al.  A dynamic model of the cardiac ventricular action potential. II. Afterdepolarizations, triggered activity, and potentiation. , 1994, Circulation research.