N Hooke | Search PubMed

PubMed1994

Linear algebraic transformations of the bidomain equations: implications for numerical methods.

A mathematical framework is presented for the treatment of the bidomain equations used to model propagation in cardiac tissue. This framework is independent of the model used to represent membrane ionic currents and incorporates boundary conditions and other constraints. By representing the bidomain equations in the operator notation L phi = F, various algebraic transformations can be expressed as PLQ-1 psi = PF, where P and Q are linear operators. The authors show how previous work fits into this framework and discuss the implications of various transformation for numerical methods of solution. Although such transformations allow many choices of independent variable, these results emphasize the fundamental importance of the transmembrane potential.

Animals↗

PubMed1992

Cardiac propagation simulation.

We have completed a range of membrane-based simulations of action potential propagation in two- and three-dimensional models of ventricular myocardium. The two-dimensional simulations included a bidomain representation of the myocardium which explicitly characterized the component volume conductors in the intracellular, interstitial, and extracellular spaces. With these simulations, we studied the contribution of the extracellular volume conductor to transmural myocardial propagation during depolarization. We also used two-dimensional bidomain simulations to study the effect of the interstitial volume conductor in the setting of planar myocardial depolarization with nominal and extreme tissue conductivities. Our three-dimensional simulations included a monodomain representation of the myocardium which characterized the three component volume conductors as a single lumped conductor. With these simulations, we examined the effects of the intramural rotation of the fiber axes on the timing and pattern of activation. To achieve practical solution times, we extended numerical techniques from previous reports and developed a range of new techniques applicable to this class of problems. Simulations of the depolarization wavefront used the nonlinear Ebihara and Johnson membrane equations for the fast sodium current as the membrane model. Simulations of the full action potential cycle combined the Ebihara and Johnson fast sodium current with the Beeler and Reuter membrane equations. Our results demonstrated that the individual volume conductors and the rotation of fiber axes have unique and identifiable consequences on the electrical activation in models of ventricular myocardium.

Action Potentials↗