A Cell-Based Framework for Numerical Modeling of Electrical Conduction in Cardiac Tissue

In this paper, we study a mathematical model of cardiac tissue based on explicit representation of individual cells. In this EMI model, the extracellular (E) space, the cell membrane (M), and the intracellular (I) space are represented as separate geometrical domains. This representation introduces modeling flexibility needed for detailed representation of the properties of cardiac cells including their membrane. In particular, we will show that the model allows ion channels to be non-uniformly distributed along the membrane of the cell. Such features are difficult to include in classical homogenized models like the monodomain and bidomain models frequently used in computational analyses of cardiac electrophysiology. The EMI model is solved using a finite difference method (FDM) and two variants of the finite element method (FEM). We compare the three schemes numerically, reporting on CPU-efforts and convergence rates. Finally, we illustrate the distinctive capabilities of the EMI model compared to classical models by simulating monolayers of cardiac cells with heterogeneous distributions of ionic channels along the cell membrane. Because of the detailed representation of every cell, the computational problems that result from using the EMI model are much larger than for the classical homogenized models, and thus represent a computational challenge. However, our numerical simulations indicate that the FDM scheme is optimal in the sense that the computational complexity increases proportionally to the number of cardiac cells in the model. Moreover, we present simulations, based on systems of equations involving 117 million unknowns, representing up to 16,000 cells. We conclude that collections of cardiac cells can be simulated using the EMI model, and that the EMI model enable greater modeling flexibility than the classical monodomain and bidomain models.