A fully implicit parallel algorithm for simulating the non‐linear electrical activity of the heart

In this paper, we study a fully implicit parallel Newton–Krylov–Schwarz method (NKS) for solving the bidomain equations describing the electrical excitation process of the heart. NKS has been used successfully for many non-linear problems, but this is the first attempt to use this method for the bidomain model which consists of a system of time dependent partial differential equations of mixed type. Our experiments on parallel computers show that the method is scalable and robust with respect to many of the parameters in the bidomain model. In the outer layer of the algorithm, we use a non-linearly implicit backward Euler method to discretize the time derivative, and the resulting systems of large sparse non-linear equations are solved using an inexact Newton method. The Jacobian system required to solve in each Newton iteration is solved with a GMRES method preconditioned by a new component-wise restricted additive Schwarz preconditioner. The efficiency and robustness of the overall method depend heavily on what preconditioner we use. By comparing several preconditioners, we found our new restricted additive Schwarz method offers the best performance. Our parallel software is developed using the PETSc package of Argonne National Laboratory. Numerical results obtained on an IBM SP will be reported. Copyright © 2004 John Wiley & Sons, Ltd.

[1]  John A. Board,et al.  A modular simulation system for the bidomain equations , 1999 .

[2]  William Gropp,et al.  Parallel Newton-Krylov-Schwarz Algorithms for the Transonic Full Potential Equation , 1996, SIAM J. Sci. Comput..

[3]  D. Rose,et al.  Linear algebraic transformations of the bidomain equations: implications for numerical methods. , 1994, Mathematical biosciences.

[4]  N. Hooke.  Efficient simulation of action potential propagation in a bidomain , 1992 .

[5]  Xiao-Chuan Cai,et al.  A Restricted Additive Schwarz Preconditioner for General Sparse Linear Systems , 1999, SIAM J. Sci. Comput..

[6]  Babskii Eb,et al.  ELECTRONIC DEVICE FOR SIMULATING THE ELECTRICAL ACTIVITY OF THE HEART. , 1964 .

[7]  A. Hodgkin,et al.  A quantitative description of membrane current and its application to conduction and excitation in nerve , 1990 .

[8]  N.V. Thakor,et al.  Initiation and termination of spiral waves in a two-dimensional bidomain model of cardiac tissue , 1996, Computers in Cardiology 1996.

[9]  C. Pappas,et al.  Application of a 3-D ischemic heart model derived from MRI data to the simulation of the electrical activity of the heart , 1996, Computers in Cardiology 1996.

[10]  Hans Othmer,et al.  A computational study of wave propagation in a model for anisotropic cardiac ventricular tissue , 1995, HPCN Europe.

[11]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[12]  V. Simoncini,et al.  Efficient algebraic solution of reaction-diffusion systems for the cardiac excitation process , 2002 .

[13]  Olof B. Widlund,et al.  Domain Decomposition Algorithms with Small Overlap , 1992, SIAM J. Sci. Comput..

[14]  Homer F. Walker,et al.  Globally Convergent Inexact Newton Methods , 1994, SIAM J. Optim..

[15]  I. Ruttkay-Nedecký,et al.  Physiology of the Cardiac Electric Field , 1993 .

[16]  B. Roth,et al.  Electrical stimulation of cardiac tissue by a bipolar electrode in a conductive bath , 1998, IEEE Transactions on Biomedical Engineering.