A Second-Order Algorithm for Solving Dynamic Cell Membrane Equations

This paper describes an extension of the so-called Rush-Larsen scheme, which is a widely used numerical method for solving dynamic models of cardiac cell electrophysiology. The proposed method applies a local linearization of nonlinear terms in combination with the analytical solution of linear ordinary differential equations to obtain a second-order accurate numerical scheme. We compare the error and computational load of the second-order scheme to the original Rush-Larsen method and a second-order Runge-Kutta (RK) method. The numerical results indicate that the new method outperforms the original Rush-Larsen scheme for all the test cases. The comparison with the RK solver reveals that the new method is more efficient for stiff problems.

[1]  Leslie Tung,et al.  A bi-domain model for describing ischemic myocardial d-c potentials , 1978 .

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

[3]  A. Tveito,et al.  An operator splitting method for solving the bidomain equations coupled to a volume conductor model for the torso. , 2005, Mathematical biosciences.

[4]  L. Shampine,et al.  A 3(2) pair of Runge - Kutta formulas , 1989 .

[5]  Uri M. Ascher,et al.  Computer methods for ordinary differential equations and differential-algebraic equations , 1998 .

[6]  A V Panfilov,et al.  A guide to modelling cardiac electrical activity in anatomically detailed ventricles. , 2008, Progress in biophysics and molecular biology.

[7]  E. Hairer,et al.  Stiff differential equations solved by Radau methods , 1999 .

[8]  Aslak Tveito,et al.  An unconditionally stable numerical method for the Luo-Rudy 1 model used in simulations of defibrillation. , 2007, Mathematical biosciences.

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

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

[11]  Xing Cai,et al.  On the Computational Complexity of the Bidomain and the Monodomain Models of Electrophysiology , 2006, Annals of Biomedical Engineering.

[12]  M. Courtemanche,et al.  Ionic mechanisms underlying human atrial action potential properties: insights from a mathematical model. , 1998, The American journal of physiology.

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

[14]  R. Winslow,et al.  Mechanisms of altered excitation-contraction coupling in canine tachycardia-induced heart failure, II: model studies. , 1999, Circulation research.