Abstract We formulate and analyse interior penalty discontinuous Galerkin methods for coupled elliptic partial differential equations (PDEs) modelling excitable tissue, represented by intracellular and extracellular domains sharing a common interface. The PDEs are coupled through a dynamic boundary condition, posed on the interface, that relates the normal gradients of the solutions to the time derivative of their jump. This system is referred to as the extracellular membrane intracellular model or the cell-by-cell model. Due to the dynamic nature of the interface condition and to the presence of corner singularities the analysis of discontinuous Galerkin methods is nonstandard. We prove the existence and uniqueness of solutions by a reformulation of the problem to one posed on the membrane. Convergence is shown by utilizing face-to-element lifting operators and notions of weak consistency suitable for solutions with low spatial regularity. Further, we present parameter-robust preconditioned iterative solvers. Numerical examples in idealized geometries demonstrate our theoretical findings, and simulations in multiple cells portray the robustness of the method.
Masri et al. (Thu,) studied this question.