Tissue Electrical Response

The electrical response of the tissue is characterized by an action potential, ϕ, and recovery variable, r. Hurtado and Kuhl describes the action potential and recovery variable.

The global electrical analysis assumes a monodomain response:

ϕ˙+div(q(ϕ))=fϕ(ϕ,r),
where the flux term, q, characterizes the propagating nature of the electrical waves:
q=Dϕ.

The term, D, is a second-order diffusion tensor, which can account for anisotropic diffusion.

The source term, fϕ, characterizes the local action potential profile:

fϕ(ϕ,r)=cϕ[ϕα][1ϕ]rϕ.

The local biochemical portion of the analysis is modeled through a temporal evolution of the recovery variable, r:

r˙=fr(ϕ,r).

The source term, fr, characterizes the slow features of the action potential:

fr(ϕ,r)=[γ+rγ¯(ϕ)][rcϕ[ϕb1]],
where
γ¯(ϕ)=μ1μ2+ϕ
with parameters as defined in the table below.
Table 1. Constitutive Parameters for the Electrical Response
Parameters Description
c Scaling parameter for the source term, fϕ
α Oscillation threshold (positive values characterize stable nonpacemaker cells, negative values characterize oscillatory pacemaker cells)
γ Refractoriness that governs the time it takes the tissue to repolarize (that is, the time to get back to the resting potential)
b A phenomenological scaling parameter
μ1,μ2 Scaling factors that control the shape of the restitution curve

The electrical properties are calibrated to provide physiologically observed activation times; see Conduction System Tutorial. The bundle of His and Purkinje fibers are assigned conduction parameters that generate the physically observed wave propagation pattern within the heart in which the electrical signal first travels down the ventricular septum to the apex and then up the ventricular side walls.