Application of Efficient Numerical Methods in Solution of Ordinary Differential Equations for Modeling Electrical Activity in Cardiac Cells

Yu Zhang,Ling Xia,Yinglan Gong
DOI: https://doi.org/10.1007/978-3-540-74282-1_49
2007-01-01
Abstract:There is a large number of ordinary differential equations (ODEs) characterize the electrical behavior generated by ionic movements in human myocardial cell. In this paper, several approaches were investigated in order to improve the efficiency of solving the ODE systems for ten Tusscher et al.'s ionic model of human ventricular tissue. By using non-standard finite difference (NSFD) scheme, the stiffness of the ODEs system will be successfully reduced, so a larger step-size can be used. A popular multi-step method called backward differentiation formulation (BDF) was also incorporated into the computational model for testing the largest possible time steps. The results show that NSFD can be as much as 10 times more efficient than standard forward Euler in single cell model simulation while maintaining an acceptable level of accuracy. The investigation of BDF method shows that a large step size is not recommended for the single cell simulation application. All solvers were coupled to the partial differential equations for the complete simulation of heart tissue, and such computation scheme may be a good calculation technique in heart modeling and simulation.
What problem does this paper attempt to address?