PARALLEL MULTILEVEL SCHWARZ AND BLOCK PRECONDITIONERS FOR THE BIDOMAIN PARABOLIC-PARABOLIC AND PARABOLIC-ELLIPTIC FORMULATIONS

• Published in: SIAM journal on scientific computing Vol. 33; no. 3-4; pp. 1897 - 1919
• Main Authors: PAVARINO, L. F, SCACCHI, S
• Format: Journal Article
• Language: English
• Published: Philadelphia, PA Society for Industrial and Applied Mathematics 01.01.2011
• Subjects: Exact sciences and technology; Linux; Mathematical analysis; Mathematics; Methods of scientific computing (including symbolic computation, algebraic computation); Numerical analysis; Numerical analysis. Scientific computation; Partial differential equations; Partial differential equations, boundary value problems; Partial differential equations, initial value problems and time-dependant initial-boundary value problems; Sciences and techniques of general use; Simulation; Differential equation; Decomposition method; Optimal estimation; Partial differential equation; Non linear equation; Parabolic equation; 65F10; Multigrid; Convergence rate; 65M55; 65N55; 92C30; Transcendental equation; block preconditioners; Convergence acceleration; Initial value problem; domain decomposition methods; Optimal convergence; Reaction diffusion equation; Non linear system; Bidomain system; Equation system; Numerical analysis; Boundary value problem; Scientific computation; Algebraic equation; Optimal approximation; Domain decomposition
• ISSN: 1064-8275; 1095-7197
• DOI: 10.1137/100808721

The aim of this work is to develop parallel multilevel and block preconditioners for the Bidomain model of electrocardiology. The Bidomain model describes the electrical activity of the heart tissue and consists of a system of two parabolic nonlinear partial differential equations (PDEs) of reaction-diffusion type (PP formulation) or alternatively of a system of a parabolic nonlinear PDE and an elliptic linear PDE (PE formulation). In both formulations, the PDEs are coupled with a system of ordinary differential equations, modeling the cellular membrane ionic currents. The first goal of the present study is to construct, analyze, and numerically test a multilevel additive Schwarz preconditioner for the PE formulation of the Bidomain model, extending previous results obtained for the PP formulation. Optimal convergence rate estimates are established and confirmed by 3D numerical test on Linux clusters. The second goal of the present study is to analyze the scalability of multilevel Schwarz block-diagonal and block-factorized preconditioners for both PP and PE formulations of the Bidomain model and to compare them with multilevel Schwarz coupled preconditioners. The 3D parallel numerical tests show that block preconditioners for the PP formulation are not scalable, while they are scalable for the PE formulation, but less efficient than the coupled preconditioners.