6.1 Synopsis

The CardiacEPSolver is used to model the electrophysiology of cardiac tissue, specifically using the monodomain or bidomain model. These models are continuum models and represent an average of the electrical activity over many cells. The system is a reaction-diffusion system, with the reaction term modeling the flow of current in and out of the cells using a separate set of ODEs.

6.1.1 Bidomain Model

The Bidomain model is given by the following PDEs,

gix∂2V
---2i
 ∂x + giy∂2V
---i2
∂y = χ[   ∂(V -  V )             ]
 Cm ---i----e-+ Gm (Vi - Ve)
        ∂t
gex∂2Ve-
∂x2 + gey∂2Ve-
∂y2 = -χ[                           ]
     ∂(Vi --Ve)
 Cm     ∂t     + Gm (Vi - Ve).
However, when solving numerically, one often rewrites these equations in terms of the transmembrane potential and extracellular potential,
χ[   ∂Vm        ]
 Cm -∂t- + Jion = gex∂2Ve
-∂x2- + gey∂2Ve
-∂y2-
(gix + gex) 2
∂-Ve-
 ∂x2 + (giy + gey) 2
∂-Ve-
∂y2 = -gix 2
∂-Vm-
 ∂x2 - giy 2
∂-Vm-
 ∂y2

6.1.2 Monodomain Model

In the case where the intracellular and extracellular conductivities are proportional, that is gix = kgex for some k, then the above two PDEs can be reduced to a single PDE:

χ[              ]
     ∂Vm-
 Cm   ∂t + Jion = ∇⋅ (σ∇V m)

6.1.3 Cell Models

The action potential of a cardiac cell can be modelled at either a biophysical level of detail, including a number of transmembrane currents, or as a phenomenological model, to reproduce the features of the action potential, with fewer variables. Each cell model will include a unique system of ODEs to represent the gating variables of that model.

A number of ionic cell models are currently supported by the solver including:

Phenomological cell models are also supported:

It is important to ensure that the units of the voltage and currents from the cell model are consistent with the units expected by the tissue level solver (monodomain/bidomain). We will show as an example the Courtemanche, Ramirez, Nattel, 1998 human atrial model.

The monodomain equation:

χ[              ]
     ∂Vm-
 Cm   ∂t + Jion = ∇⋅ (σ∇V m)