Spherical topology in cardiac simulations

Computational simulations of the electrodynamics of cardiac fibrillation yield a great deal of useful data and provide a framework for theoretical explanations of heart behavior. Extending the application of these mathematical models to defibrillation studies requires that a simulation should sustain fibrillation without defibrillation intervention. In accordance with the critical mass hypothesis, the simulated tissue should be of a large enough size. The choice of biperiodic boundary conditions sustains fibrillation for a longer duration than no-flux boundary conditions for a given area, and so is commonly invoked. Here, we show how this leads to a boundary condition artifact that may complicate the analysis of defibrillation efficacy; we implement an alternative coordinate scheme that utilizes spherical shell topology and mitigates singularities in the Laplacian found with the usual spherical curvilinear coordinate system

A Model for Multi-site Pacing of Fibrillation Using Nonlinear Dynamics Feedback

Traditionally, cardiac defibrillation requires a strong electric shock. Many unwanted side effects of this shock could be eliminated if defibrillation were performed using weak stimuli applied to several locations throughout the heart. Such multi-site pacing algorithms have been shown to defibrillate both experimentally (Pak et al., Am J Physiol 285:H2704–H2711, 2003) and theoretically (Puwal and Roth, J Biol Systems 14:101–112, 2006). Gauthier et al. (Chaos, 12:952–961, 2002) proposed a method to pace the heart using an algorithm based on nonlinear dynamics feedback applied through a single electrode. Our study applies a related but simpler algorithm, which essentially configures each electrode as a demand pacemaker, to simulate the multi-site pacing of fibrillating cardiac tissue. We use the numerical model developed by Fenton et al. (Chaos, 12:852–892, 2002) as the reaction term in a reaction–diffusion equation that we solve over a two-dimensional sheet of tissue. The defibrillation rate after pacing for 3 s is about 30%, which is significantly higher than the spontaneous defibrillation rate and is higher than observed in previous experimental and theoretical studies. Tuning the algorithm period can increase this rate to 45%

Fourier-based magnetic induction tomography for mapping resistivity

Magnetic induction tomography is used as an experimental tool for mapping the passive electromagnetic properties of conductors, with the potential for imaging biological tissues. Our numerical approach to solving the inverse problem is to obtain a Fourier expansion of the resistivity and the stream functions of the magnetic fields and eddy current density. Thus, we are able to solve the inverse problem of determining the resistivity from the applied and measured magnetic fields for a two-dimensional conducting plane. When we add noise to the measured magnetic field, we find the fidelity of the measured to the true resistivity is quite robust for increasing levels of noise and increasing distances of the applied and measured field coils from the conducting plane, when properly filtered. We conclude that Fourier methods provide a reliable alternative for solving the inverse problem