The bidomain equations are the most complete description of cardiac electrical activity. Their numerical solution is, however, computationally demanding, especially in three dimensions, because of the fine temporal and spatial sampling required. This paper methodically examines computational performance when solving the bidomain equations. Several techniques to speed up this computation are examined in this paper. The first step was to recast the equations into a parabolic part and an elliptic part. The parabolic part was solved by either the finite-element method (FEM) or the interconnected cable model model (ICCM). The elliptic equation was solved by FEM on a coarser grid than the parabolic problem and at a reduced frequency. The performance of iterative and direct linear equation system solvers was analyzed as well as the scalability and parallelizability of each method. Results indicate that the ICCM was twice as fast as the FEM for solving the parabolic problem, but when the total problem was considered, this resulted in only a 20% decrease in computation time. The elliptic problem could be solved on a coarser grid at one-quarter of the frequency at which the parabolic problem was solved and still maintain reasonable accuracy. Direct methods were faster than iterative methods by at least 50% when a good estimate of the extracellular potential was required. Parallelization over four processors was efficient only when the model comprised at least 500,000 nodes. Thus, it was possible to speed up solution of the bidomain equations by an order of magnitude with a slight decrease in accuracy.

We recently suggested that failure of implantable defibrillation therapy may be explained by the virtual electrode-induced phase singularity mechanism. The goal of this study was to identify possible mechanisms of vulnerability and defibrillation by externally applied shocks in vitro. We used bidomain simulations of realistic rabbit heart fibrous geometry to predict the passive polarization throughout the heart induced by external shocks. We also used optical mapping to assess anterior epicardium electrical activity during shocks in Langendorff-perfused rabbit hearts (n = 7). Monophasic shocks of either polarity (10-260 V, 8 ms, 150 microF) were applied during the T wave from a pair of mesh electrodes. Postshock epicardial virtual electrode polarization was observed after all 162 applied shocks, with positive polarization facing the cathode and negative polarization facing the anode, as predicted by the bidomain simulations. During arrhythmogenesis, a new wave front was induced at the boundary between the two regions near the apex but not at the base. It spread across the negatively polarized area toward the base of the heart and reentered on the other side while simultaneously spreading into the depth of the wall. Thus a scroll wave with a ribbon-shaped filament was formed during external shock-induced arrhythmia. Fluorescent imaging and passive bidomain simulations demonstrated that virtual electrode polarization-induced scroll waves underlie mechanisms of shock-induced vulnerability and failure of external defibrillation.

The purpose of this paper is to investigate the relationship between cardiac tissue structure, applied electric field, and the transmembrane potential induced in the process of defibrillation. It outlines a general understanding of the structural mechanisms that contribute to the outcome of a defibrillation shock. Electric shocks defibrillate by changing the transmembrane potential throughout the myocardium. In this process first and foremost the shock current must access the bulk of myocardial mass. The exogenous current traverses the myocardium along convoluted intracellular and extracellular pathways channeled by the tissue structure. Since individual fibers follow curved pathways in the heart, and the fiber direction rotates across the ventricular wall, the applied current perpetually engages in redistribution between the intra- and extracellular domains. This redistribution results in changes in transmembrane potential (membrane polarization): regions of membrane hyper- and depolarization of extent larger than a single cell are induced in the myocardium by the defibrillation shock. Tissue inhomogeneities also contribute to local membrane polarization in the myocardium which is superimposed over the large-scale polarization associated with the fibrous organization of the myocardium. The paper presents simulation results that illustrate various mechanisms by which cardiac tissue structure assists the changes in transmembrane potential throughout the myocardium. (c) 1998 American Institute of Physics.

No abstract

Background-Functional reentry in the heart takes the form of spiral waves. Drifting spiral waves can become pinned to anatomic obstacles and thus attain stability and persistence. Lidocaine is an antiarrhythmic agent commonly used to treat ventricular tachycardia clinically. We examined the ability of small obstacles to anchor spiral waves and the effect of lidocaine on their attachment. Methods and Results-Spiral waves were electrically induced in confluent monolayers of cultured, neonatal rat cardiomyocytes. Small, circular anatomic obstacles (0.6 to 2.6 mm in diameter) were situated in the center of the monolayers to provide an anchoring site. Eighty reentry episodes consisting of at least 4 revolutions were studied. In 36 episodes, the spiral wave attached to the obstacle and became stationary and sustained, with a shorter reentry cycle length and higher rate. Spiral waves could attach to obstacles as small as 0.6 mm, with a likelihood for attachment that increased with obstacle size. After attachment, both conduction velocity of the wave-front tip and wavelength near the obstacle adapted from their pre-reentry values and increased linearly with obstacle size. In contrast, reentry cycle length did not correlate significantly with obstacle size. Addition of lidocaine 90 mol/L depressed conduction velocity, increased reentry cycle length, and caused attached spiral waves to become quasiattached to the obstacle or terminate. Conclusions-Anchored spiral waves exhibit properties of both unattached spiral waves and anatomic reentry. Their behavior may be representative of functional reentry dynamics in cardiac tissue, particularly in the setting of monomorphic tachyarrhythmias.