A hybrid multilevel Schwarz method is studied numerically for the anisotropic Bidomain model in both two and three dimensions. This multiscale system models the electrical activity of the heart and it consists of two degenerate parabolic non-linear reaction–diffusion equations, coupled with a stiff system of ordinary differential equations. The numerical discretization of the whole system by finite elements in space and semi-implicit methods in time generates ill-conditioned linear systems that must be solved at each time step. The multilevel algorithm studied employs a hierarchy of nested meshes with overlapping Schwarz preconditioners on each level and is additive within the levels and multiplicative among the levels. We perform several parallel tests on two Linux clusters, showing that the convergence of the method is independent of the number of subdomains (scalability), the discretization parameters and the number of levels (optimality). Moreover the comparison with the traditional Block Jacobi ILU parallel preconditioner and the V-cycle Multigrid parallel preconditioner shows that, on a whole heart cycle simulation, the proposed method attains the best performances in terms of CPU times.

A hybrid multilevel Schwarz method for the bidomain model / S. Scacchi. - In: COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING. - ISSN 0045-7825. - 197:45-48(2008), pp. 4051-4061.

A hybrid multilevel Schwarz method for the bidomain model

S. Scacchi
Primo
2008

Abstract

A hybrid multilevel Schwarz method is studied numerically for the anisotropic Bidomain model in both two and three dimensions. This multiscale system models the electrical activity of the heart and it consists of two degenerate parabolic non-linear reaction–diffusion equations, coupled with a stiff system of ordinary differential equations. The numerical discretization of the whole system by finite elements in space and semi-implicit methods in time generates ill-conditioned linear systems that must be solved at each time step. The multilevel algorithm studied employs a hierarchy of nested meshes with overlapping Schwarz preconditioners on each level and is additive within the levels and multiplicative among the levels. We perform several parallel tests on two Linux clusters, showing that the convergence of the method is independent of the number of subdomains (scalability), the discretization parameters and the number of levels (optimality). Moreover the comparison with the traditional Block Jacobi ILU parallel preconditioner and the V-cycle Multigrid parallel preconditioner shows that, on a whole heart cycle simulation, the proposed method attains the best performances in terms of CPU times.
Multilevel Schwarz preconditioner ; Domain decomposition ; Bidomain model ; Computational electrocardiology
Settore MAT/08 - Analisi Numerica
2008
Article (author)
File in questo prodotto:
Non ci sono file associati a questo prodotto.
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/55748
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 47
  • ???jsp.display-item.citation.isi??? 46
social impact