A multiresolution space-time adaptive scheme for the bidomain model in electrocardiology

The bidomain model of electrical activity of myocardial tissue consists of a possibly degenerate parabolic PDE coupled with an elliptic PDE for the transmembrane and extracellular potentials, respectively. This system of two scalar PDEs is supplemented by a time-dependent ODE modeling the evolution of the gating variable. In the simpler sub-case of the monodomain model, the elliptic PDE reduces to an algebraic equation. Since typical solutions of the bidomain and monodomain models exhibit wavefronts with steep gradients, we propose a finite volume scheme enriched by a fully adaptive multiresolution method, whose basic purpose is to concentrate computational effort on zones of strong variation of the solution. Time adaptivity is achieved by two alternative devices, namely locally varying time stepping and a Runge-Kutta-Fehlberg-type adaptive time integration. A series of numerical examples demonstrates that these methods are efficient and sufficiently accurate to simulate the electrical activity in myocardial tissue with affordable effort. In addition, the optimal choice of the threshold for discarding non-significant information in the multiresolution representation of the solution is addressed, and the numerical efficiency and accuracy of the method is measured in terms of CPU time speed-up, memory compression, and errors in different norms.


Published in:
Numerical Methods for Partial Differential Equations, 26, 1377-1404
Year:
2010
Publisher:
Wiley-Blackwell
ISSN:
0749-159X
Keywords:
Laboratories:




 Record created 2009-02-16, last modified 2018-03-17

n/a:
Download fulltext
PDF

Rate this document:

Rate this document:
1
2
3
 
(Not yet reviewed)