ISBI'09 Proceedings of the Sixth IEEE international conference on Symposium on Biomedical Imaging: From Nano to Macro
STACOM'10/CESC'10 Proceedings of the First international conference on Statistical atlases and computational models of the heart, and international conference on Cardiac electrophysiological simulation challenge
Preconditioning the bidomain model with almost linear complexity
Journal of Computational Physics
SIAM Journal on Scientific Computing
Towards real-time computation of cardiac electrophysiology for training simulator
STACOM'12 Proceedings of the third international conference on Statistical Atlases and Computational Models of the Heart: imaging and modelling challenges
| Hi-index | 0.01 |
The bidomain model is a system of partial differential equations used to model the propagation of electrical potential waves in the myocardium. It is composed of coupled parabolic and elliptic partial differential equations, as well as at least one ordinary differential equation to model the ion activity through the cardiac cell membranes. The purpose of this paper is to propose and analyze several implicit, semi-implicit, and explicit time-stepping methods to solve that model, in particular to avoid the expensive resolution of a nonlinear system through the Newton-Raphson method. We identify necessary stability conditions on the time step $\Delta t$ for the proposed methods through a theoretical analysis based on energy estimates. We next compare the methods for one- and two-dimensional test cases, in terms of both stability and accuracy of the numerical solutions. The theoretical stability conditions are seen to be consistent with those observed in practice. Our analysis allows us to recommend using either the Crank-Nicolson/Adams-Bashforth method or the second order semi-implicit backward differention method. These semi-implicit methods produce a good numerical solution; unlike the explicit methods, their stability does not depend on the spatial grid size; and unlike the implicit methods, they do not require the resolution of a system of nonlinear equations.