A multilevel hybrid Newton-Krylov-Schwarz (NKS) method is constructed and studied numerically for implicit time discretizations of the Bidomain reaction-diffusion system in three dimensions. This model describes the bioelectrical activity of the heart by coupling two degenerate parabolic equations with a stiff system of ordinary differential equations. The 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 multilevel hybrid overlapping Schwarz preconditioner, additive within the levels and multiplicative among the levels. Several parallel tests on Linux clusters are performed, showing that the convergence of the method is independent of the number of subdomains (scalability), the discretization parameters and the number of levels (optimality).
A multilevel hybrid Newton-Krylov-Schwarz method for the Bidomain model of electrocardiology / S. Scacchi. - In: COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING. - ISSN 0045-7825. - 200:5-8(2011), pp. 717-725. [10.1016/j.cma.2010.09.016]
A multilevel hybrid Newton-Krylov-Schwarz method for the Bidomain model of electrocardiology
S. Scacchi
2011
Abstract
A multilevel hybrid Newton-Krylov-Schwarz (NKS) method is constructed and studied numerically for implicit time discretizations of the Bidomain reaction-diffusion system in three dimensions. This model describes the bioelectrical activity of the heart by coupling two degenerate parabolic equations with a stiff system of ordinary differential equations. The 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 multilevel hybrid overlapping Schwarz preconditioner, additive within the levels and multiplicative among the levels. Several parallel tests on Linux clusters are performed, showing that the convergence of the method is independent of the number of subdomains (scalability), the discretization parameters and the number of levels (optimality).File | Dimensione | Formato | |
---|---|---|---|
1-s2.0-S0045782510002768-main.pdf
accesso riservato
Tipologia:
Publisher's version/PDF
Dimensione
910.87 kB
Formato
Adobe PDF
|
910.87 kB | Adobe PDF | Visualizza/Apri Richiedi una copia |
Pubblicazioni consigliate
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.