Year: 2026
Type: article
Abstract: In the present study, we consider the Extra-Membrane-Intra (EMI) model for the simulation of excitable tissues at the cellular level. We provide the (possibly large) system of partial differential equations (PDEs), equipped with ad hoc boundary conditions, relevant for modelling portions of excitable tissues, composed of several cells. In particular, we study two geometrical settings in the context of computational cardiology and neuroscience. The Galerkin approximations of the considered system of PDEs lead to large linear systems of algebraic equations, where the coefficient matrices depend on the number N of cells and the fineness parameters. We present a structural and spectral analysis of the related matrix-sequences as the fineness parameters tends to zero. Based on the theoretical results, we propose preconditioners and specific multilevel solvers. Numerical experiments are presented and critically discussed, showing that a monolithic multilevel solver is efficient and robust with respect to all the problem and discretization parameters. In particular, we include numerical results for an increasing number of cells N , both for idealized geometries (with N exceeding $$10^5$$ 10 5 ) and for realistic, densely populated 3D tissue reconstructions.
Source: BIT Numerical Mathematics
Language: English
Topic: Matrix Theory and Algorithms
Subfield: Computational Theory and Mathematics
Field: Computer Science
Domain: Physical Sciences
Open Access status: hybrid
Estimated APC paid: $3,190
OpenAlex ID: W4403758534