| This work considers the approximation of the cardiac bidomain equations, either isolated or coupled with the torso, via first order semi-implicit time-marching schemes involving a fully decoupled computation of the unknown fields (ionic state, transmembrane potential, extracellular and torso potentials). For the isolated bidomain system, we show that the Gauss-Seidel and Jacobi like splittings do not compromise energy stability; they simply alter the energy norm. Within the framework of the numerical simulation of electrocardiograms (ECG), these bidomain splittings are combined with an explicit Robin-Robin treatment of the heart-torso coupling conditions. We show that the resulting schemes allow a fully decoupled (energy) stable computation of the heart and torso fields, under an additional hyperbolic-CFL like condition. The accuracy and convergence rate of the considered schemes are investigated numerically with a series of numerical experiments. |
BibTeX,EndNote,...
1Nejib Zemzemi 1
