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

Mostafa Bendahmane, Raimund Bürger, Ricardo Ruiz-Baier

Research output: Contribution to journalArticleResearchpeer-review

24 Citations (Scopus)

Abstract

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 subcase 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 nonsignificant 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.

Original languageEnglish
Pages (from-to)1377-1404
Number of pages28
JournalNumerical Methods for Partial Differential Equations
Volume26
Issue number6
DOIs
Publication statusPublished - 1 Nov 2010
Externally publishedYes

Keywords

  • Bidomain model
  • Electrocardiology
  • Fully adaptive multiresolution schemes
  • Local time stepping
  • Parabolic-elliptic system
  • Runge-Kutta-Fehlberg

Cite this