A study on the simulation of HL-1 unicellular model based on numerical methods

摘要

With rapid development in imaging modalities and computer power, computational model of cardiac electrophysiology field has been integrated with more and more powerful and sophisticated numerical solvers, which aims to give people a more accurate understanding of the heart in pathological conditions. However, due to the limitation of the accuracy that the numerical methods could reach, the computational models may not yet be so accurate. In this paper, an implementation of Takeuchi-HL-1 unicellular model based on the first-order Euler method and the fourth-order Runge-Katta RK4 method is introduced. Besides accuracy, our work is also mainly focused on the optimization of data utilization based on the Runge-Kutta RK4 method. Through our work, it is shown that the higher order of accuracy method does give us a more accurate simulation result, while the time consumption is generally acceptable.