A novel two-level Newton-Krylov-Schwarz (NKS) solver is constructed and analyzed for implicit time discretizations of the bidomain reaction-diffusion system in three dimensions. This multiscale system describes the bioelectrical activity of the heart by coupling two degenerate parabolic equations with several ordinary differential equations at each point in space. Together with a finite element discretization in space, the proposed NKS Bidomain solver employs an outer inexact Newton iteration to solve the nonlinear finite element system originating at each time step of the decoupled implicit discretization. The Jacobian update during the Newton iteration is solved by a Krylov method employing a two-level overlapping Schwarz preconditioner. A convergence rate estimate is proved for the resulting preconditioned operator, showing that its condition number is independent of the number of subdomains (scalability) and bounded by the ratio of the subdomains characteristic size and the overlap size. This theoretical result is confirmed by several parallel simulations employing up to more than 2000 processors for scaled and standard speedup tests in three dimensions. The results show the scalability of the proposed NKS Bidomain solver in terms of both nonlinear and linear iterations, in both Cartesian slabs and ellipsoidal cardiac domains.

A scalable Newton-Krylov-Schwarz method for the Bidomain reaction-diffusion system / M. Munteanu, L. F. Pavarino, S. Scacchi. - In: SIAM JOURNAL ON SCIENTIFIC COMPUTING. - ISSN 1064-8275. - 31:5(2009), pp. 3861-3883. [10.1137/08074355X]

A scalable Newton-Krylov-Schwarz method for the Bidomain reaction-diffusion system

M. Munteanu
Primo
;
L.F. Pavarino
Secondo
;
S. Scacchi
Ultimo
2009

Abstract

A novel two-level Newton-Krylov-Schwarz (NKS) solver is constructed and analyzed for implicit time discretizations of the bidomain reaction-diffusion system in three dimensions. This multiscale system describes the bioelectrical activity of the heart by coupling two degenerate parabolic equations with several ordinary differential equations at each point in space. Together with a finite element discretization in space, the proposed NKS Bidomain solver employs an outer inexact Newton iteration to solve the nonlinear finite element system originating at each time step of the decoupled implicit discretization. The Jacobian update during the Newton iteration is solved by a Krylov method employing a two-level overlapping Schwarz preconditioner. A convergence rate estimate is proved for the resulting preconditioned operator, showing that its condition number is independent of the number of subdomains (scalability) and bounded by the ratio of the subdomains characteristic size and the overlap size. This theoretical result is confirmed by several parallel simulations employing up to more than 2000 processors for scaled and standard speedup tests in three dimensions. The results show the scalability of the proposed NKS Bidomain solver in terms of both nonlinear and linear iterations, in both Cartesian slabs and ellipsoidal cardiac domains.
Domain decomposition overlapping Schwarz preconditioners; Monodomain and bidomain models; Newton-Krylov-Schwarz methods; Nonlinear reaction-diffusion systems; Parallel numerical simulations
Settore MAT/08 - Analisi Numerica
Article (author)
File in questo prodotto:
Non ci sono file associati a questo prodotto.
Pubblicazioni consigliate

Caricamento 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/143706
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 41
  • ???jsp.display-item.citation.isi??? 40
social impact