@article{GoebelHuynhCheginietal.2025, author = {G{\"o}bel, Fritz and Huynh, Ngoc Mai Monica and Chegini, Fatemeh and Pavarino, Luca and Weiser, Martin and Scacchi, Simone and Anzt, Hartwig}, title = {A BDDC Preconditioner for the Cardiac EMI Model in three Dimensions}, journal = {SIAM J. Sci. Comput.}, arxiv = {http://arxiv.org/abs/2502.07722}, year = {2025}, abstract = {We analyze a Balancing Domain Decomposition by Constraints (BDDC) preconditioner for the solution of three dimensional composite Discontinuous Galerkin discretizations of reaction-diffusion systems of ordinary and partial differential equations arising in cardiac cell-by-cell models like the Extracellular space, Membrane and Intracellular space (EMI) Model. These microscopic models are essential for the understanding of events in aging and structurally diseased hearts which macroscopic models relying on homogenized descriptions of the cardiac tissue, like Monodomain and Bidomain models, fail to adequately represent. The modeling of each individual cardiac cell results in discontinuous global solutions across cell boundaries, requiring the careful construction of dual and primal spaces for the BDDC preconditioner. We provide a scalable condition number bound for the precondition operator and validate the theoretical results with extensive numerical experiments.}, language = {en} } @article{HuynhCheginiPavarinoetal.2023, author = {Huynh, Ngoc and Chegini, Fatemeh and Pavarino, Luca and Weiser, Martin and Scacchi, Simone}, title = {Convergence analysis of BDDC preconditioners for hybrid DG discretizations of the cardiac cell-by-cell model}, volume = {45}, journal = {SIAM Journal on Scientific Computing}, number = {6}, arxiv = {http://arxiv.org/abs/2212.12295}, pages = {A2836 -- A2857}, year = {2023}, abstract = {A Balancing Domain Decomposition by Constraints (BDDC) preconditioner is constructed and analyzed for the solution of hybrid Discontinuous Galerkin discretizations of reaction-diffusion systems of ordinary and partial differential equations arising in cardiac cell-by-cell models. The latter are different from the classical Bidomain and Monodomain cardiac models based on homogenized descriptions of the cardiac tissue at the macroscopic level, and therefore they allow the representation of individual cardiac cells, cell aggregates, damaged tissues and nonuniform distributions of ion channels on the cell membrane. The resulting discrete cell-by-cell models have discontinuous global solutions across the cell boundaries, hence the proposed BDDC preconditioner is based on appropriate dual and primal spaces with additional constraints which transfer information between cells (subdomains) without influencing the overall discontinuity of the global solution. A scalable convergence rate bound is proved for the resulting BDDC cell-by-cell preconditioned operator, while numerical tests validate this bound and investigate its dependence on the discretization parameters.}, language = {en} } @inproceedings{CheginiFroehlyHuynhetal.2023, author = {Chegini, Fatemeh and Froehly, Algiane and Huynh, Ngoc Mai Monica and Pavarino, Luca and Potse, Mark and Scacchi, Simone and Weiser, Martin}, title = {Efficient numerical methods for simulating cardiac electrophysiology with cellular resolution}, booktitle = {10th Int. Conf. Computational Methods for Coupled Problems in Science and Engineering 2023}, doi = {10.23967/c.coupled.2023.004}, year = {2023}, abstract = {The cardiac extracellular-membrane-intracellular (EMI) model enables the precise geometrical representation and resolution of aggregates of individual myocytes. As a result, it not only yields more accurate simulations of cardiac excitation compared to homogenized models but also presents the challenge of solving much larger problems. In this paper, we introduce recent advancements in three key areas: (i) the creation of artificial, yet realistic grids, (ii) efficient higher-order time stepping achieved by combining low-overhead spatial adaptivity on the algebraic level with progressive spectral deferred correction methods, and (iii) substructuring domain decomposition preconditioners tailored to address the complexities of heterogeneous problem structures. The efficiency gains of these proposed methods are demonstrated through numerical results on cardiac meshes of different sizes.}, language = {en} }