A two-level Newton–Krylov–Schwarz (NKS) solver is constructed and analyzed for implicit time discretizations of the Bidomain reaction-diffusion system. 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. 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 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 2,000 processors for scaled and standard speedup tests in three dimensions.

A Two-Level Newton-Krylov-Schwarz Method for the Bidomain Model of Electrocardiology / M. Munteanu, L.F. Pavarino, S. Scacchi - In: Numerical mathematics and advanced applications 2009 / [a cura di] G. Kreiss, P. Lötstedt, A. Målqvist, M. Neytcheva. - Berlin : Springer, 2010. - ISBN 978-3-642-11795-4. - pp. 683-691 (( Intervento presentato al 8. convegno ENUMATH 2009 tenutosi a Uppsala nel 2009 [10.1007/978-3-642-11795-4_73].

A Two-Level Newton-Krylov-Schwarz Method for the Bidomain Model of Electrocardiology

L.F. Pavarino
Secondo
;
S. Scacchi
Ultimo
2010

Abstract

A two-level Newton–Krylov–Schwarz (NKS) solver is constructed and analyzed for implicit time discretizations of the Bidomain reaction-diffusion system. 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. 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 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 2,000 processors for scaled and standard speedup tests in three dimensions.
Settore MAT/08 - Analisi Numerica
Book Part (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/150577
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus ND
  • ???jsp.display-item.citation.isi??? 0
social impact