A comparison of coupled and uncoupled solvers for the cardiac Bidomain model

Abstract

Abstract. The aim of this work is to compare a new uncoupled solver for the cardiac Bidomain model with a usual coupled solver. The Bidomain model describes the bioelectric activity of the cardiac tissue and consists of a system of a non-linear parabolic reaction-diffusion partial differential equation (PDE) and an elliptic linear PDE. This system models at macroscopic level the evolution of the transmembrane and extracellular electric potentials of the anisotropic cardiac tissue. The evolution equation is coupled through the non-linear reaction term with a stiff system of ordinary differential equations (ODEs), the so-called membrane model, describing the ionic currents through the cellular membrane. A novel uncoupled solver for the Bidomain system is here introduced, based on solving twice the parabolic PDE and once the elliptic PDE at each time step, and it is compared with a usual coupled solver. Three-dimensional numerical tests have been performed in order to show that the proposed uncoupled method has the same accuracy of the coupled strategy. Parallel numerical tests on structured meshes have also shown that the uncoupled technique is as scalable as the coupled one. Moreover, the conjugate gradient method preconditioned by Multilevel Hybrid Schwarz preconditioners converges faster for the linear systems deriving from the uncoupled method than from the coupled one. Finally, in all parallel numerical tests considered, the uncoupled technique proposed is always about two or three times faster than the coupled approach.Mathematics Subject Classification. 65N55, 65M55, 65F10, 92C30.

Role of infarct scar dimensions, border zone repolarization properties and anisotropy in the origin and maintenance of cardiac reentry

“…The space discretization of system (2) is performed by employing hexahedral isoparametric Q 1 finite elements, while the time discretization is based on the following double operator splitting procedure: a) split the ODEs from the PDEs and b) split the elliptic PDE from the parabolic one; for further details see reference [25]. We refer to [80,82,57] for other numerical strategies adopted in the literature.…”

Isogeometric Analysis of the electrophysiology in the human heart: Numerical simulation of the bidomain equations on the atria

“…Under the assumptions that v in ∈ W, w in ∈ [W] nw , and that I s i,e ∈ L 2 (Ω × (0, T )) satisfy Eq. ( 4), the weak formulation of the bidomain equations [25] reads…”

BPX preconditioners for the Bidomain model of electrocardiology

“…Alternatively, most previous works have considered semi-implicit (IMEX) time discretizations and/or operator splitting schemes, where the reaction and diffusion terms are treated separately, see e.g. [5,6,7,8,9,10,11,12,13,14,15]. The advantage of IMEX and operator splitting schemes is that they only require the solution of linear systems at each time step.…”