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.