The focus of this work is to provide an extensive numerical study of the parallel efficiency and robustness of a staggered dual-primal Newton–Krylov deluxe solver for implicit time discretizations of the Bidomain model. This model describes the propagation of the electrical impulse in the cardiac tissue, by means of a system of parabolic reaction-diffusion partial differential equations. This system is coupled to a system of ordinary differential equations, modeling the ionic currents dynamics. A staggered approach is employed for the solution of a fully implicit time discretization of the problem, where the two systems are solved successively. The arising nonlinear algebraic system is solved with a Newton–Krylov approach, preconditioned by a dual-primal Domain Decomposition algorithm in order to improve convergence. The theoretical analysis and numerical validation of this strategy has been carried out in Huynh et al. (SIAM J. Sci. Comput. 44, B224–B249, 2022) considering only simple ionic models. This paper extends this study to include more complex biophysical ionic models, as well as the presence of ischemic regions, described mathematically by heterogeneous diffusion coefficients with possible discontinuities between subregions. The results of several numerical experiments show robustness and scalability of the proposed parallel solver.

Scalable and Robust Dual-Primal Newton–Krylov Deluxe Solvers for Cardiac Electrophysiology with Biophysical Ionic Models / N.M.M. Huynh, L.F. Pavarino, S. Scacchi. - In: VIETNAM JOURNAL OF MATHEMATICS. - ISSN 2305-221X. - 50:4(2022), pp. 1029-1052. [10.1007/s10013-022-00576-1]

Scalable and Robust Dual-Primal Newton–Krylov Deluxe Solvers for Cardiac Electrophysiology with Biophysical Ionic Models

S. Scacchi
Ultimo
2022

Abstract

The focus of this work is to provide an extensive numerical study of the parallel efficiency and robustness of a staggered dual-primal Newton–Krylov deluxe solver for implicit time discretizations of the Bidomain model. This model describes the propagation of the electrical impulse in the cardiac tissue, by means of a system of parabolic reaction-diffusion partial differential equations. This system is coupled to a system of ordinary differential equations, modeling the ionic currents dynamics. A staggered approach is employed for the solution of a fully implicit time discretization of the problem, where the two systems are solved successively. The arising nonlinear algebraic system is solved with a Newton–Krylov approach, preconditioned by a dual-primal Domain Decomposition algorithm in order to improve convergence. The theoretical analysis and numerical validation of this strategy has been carried out in Huynh et al. (SIAM J. Sci. Comput. 44, B224–B249, 2022) considering only simple ionic models. This paper extends this study to include more complex biophysical ionic models, as well as the presence of ischemic regions, described mathematically by heterogeneous diffusion coefficients with possible discontinuities between subregions. The results of several numerical experiments show robustness and scalability of the proposed parallel solver.
Bidomain system; Deluxe scaling; Domain decomposition; FETI-DP and BDDC preconditioners; Implicit time discretizations
Settore MAT/08 - Analisi Numerica
Article (author)
File in questo prodotto:
File Dimensione Formato  
huynhPS_2022_viet.pdf

accesso aperto

Tipologia: Publisher's version/PDF
Dimensione 1.55 MB
Formato Adobe PDF
1.55 MB Adobe PDF Visualizza/Apri
Pubblicazioni consigliate

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/2434/944030
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 0
  • ???jsp.display-item.citation.isi??? 0
social impact