A numerical guide to the solution of the bidomain equations of cardiac electrophysiology

Simulation of cardiac electrical activity using the bidomain equations can be a massively computationally demanding problem. This study provides a comprehensive guide to numerical bidomain modelling. Each component of bidomain simulations-discretisation, ODE-solution, linear system solution, and parallelisation-is discussed, and previously-used methods are reviewed, new methods are proposed, and issues which cause particular difficulty are highlighted. Particular attention is paid to the choice of stimulus currents, compatibility conditions for the equations, the solution of singular linear systems, and convergence of the numerical scheme. © 2010 Elsevier Ltd.

Pathmanathan P, Bernabeu MO, Bordas R, Cooper J, Garny A, et al. (2010) A numerical guide to the solution of the bidomain equations of cardiac electrophysiology. Progress in Biophysics and Molecular Biology 102: 136–155. Available: http://dx.doi.org/10.1016/j.pbiomolbio.2010.05.006.

PP is supported by the EPSRC-funded OXMOS project New frontiers in the mathematics of solids, (grant reference EP/D048400/1). JC is supported by the European Community's Seventh Framework Programme [FP7/2007-2013] under grant agreement 223920 (VPH NoE). AG is funded through the preDiCT and euHeart projects (numbers 224381 and 224495, respectively) which are both supported by the European Commission, DG Information Society, through the Seventh Framework Programme of Information and Communication Technologies. JPW is supported by Award No. KUK-C1-013-04, made by King Abdullah University of Science and Technology (KAUST).

Elsevier BV

Progress in Biophysics and Molecular Biology


PubMed ID

Permanent link to this record