A numerically efficient model for simulation of defibrillation in an active bidomain sheet of myocardium

Presented here is an efficient algorithm for solving the bidomain equations describing myocardial tissue with active membrane kinetics. An analysis of the accuracy shows advantages of this numerical technique over other simple and therefore popular approaches. The modular structure of the algorithm provides the critical flexibility needed in simulation studies: fiber orientation and membrane kinetics can be easily modified. The computational tool described here is designed specifically to simulate cardiac defibrillation, i.e., to allow modeling of strong electric shocks applied to the myocardium extracellularly. Accordingly, the algorithm presented also incorporates modifications of the membrane model to handle the high transmembrane voltages created in the immediate vicinity of the defibrillation electrodes.