References & Citations
Computer Science > Computational Engineering, Finance, and Science
Title: A splitting, discontinuous Galerkin solver for the cell-by-cell electroneutral Nernst-Planck framework
(Submitted on 12 Apr 2024 (v1), last revised 25 Apr 2024 (this version, v2))
Abstract: Mathematical models for excitable tissue with explicit representation of individual cells are highly detailed and can, unlike classical homogenized models, represent complex cellular geometries and local membrane variations. However, these cell-based models are challenging to approximate numerically, partly due to their mixed-dimensional nature with unknowns both in the bulk and at the lower-dimensional cellular membranes. We here develop and evaluate a novel solution strategy for the cell-based KNP-EMI model describing ionic electrodiffusion in and between intra- and extracellular compartments with explicit representation of individual cells. The strategy is based on operator splitting, a multiplier-free formulation of the coupled dynamics across sub-regions, and a discontinuous Galerkin discretization. In addition to desirable theoretical properties, such as local mass conservation, the scheme is practical as it requires no specialized functionality in the finite element assembly and order optimal solvers for the resulting linear systems can be realized with black-box algebraic multigrid preconditioners. Numerical investigations show that the proposed solution strategy is accurate, robust with respect to discretization parameters, and that the parallel scalability of the solver is close to optimal - both for idealized and realistic two and three dimensional geometries.
Submission history
From: Ada Johanne Ellingsrud [view email][v1] Fri, 12 Apr 2024 08:27:06 GMT (5190kb,D)
[v2] Thu, 25 Apr 2024 11:55:33 GMT (5190kb,D)
Link back to: arXiv, form interface, contact.