We present a multi-level massively parallel additive Schwarz preconditioner for Isogeometric Analysis, a FEM-like numerical analysis for PDEs that permits exact geometry representation and high regularity basis functions. Two model problems are considered: the scalar elliptic equation and the advection-diffusion equation. Theoretical analysis proves that the adoption of a coarse correction grid is crucial in order to have the condition number of the preconditioned stiffness matrix independent from the number of subdomains, whenever the ratio between the coarse mesh size and the fine mesh size is kept fixed. Numerical tests for the scalar elliptic equation (in 2D and 3D on trivial and non-trivial domains) confirm the theory. The preconditioner is then applied to the advection-diffusion equation in 2D and 3D. Again, the numerical results shows that the condition number of the preconditioned linear system scales with the number of subdomains up to 8100 processors, eventually with SUPG stabilization. The tests are implemented in C programming language on the top of PETSc library.
PARALLEL ADDITIVE SCHWARZ PRECONDITIONING FOR ISOGEOMETRIC ANALYSIS
MARINI, FEDERICO
2015
Abstract
We present a multi-level massively parallel additive Schwarz preconditioner for Isogeometric Analysis, a FEM-like numerical analysis for PDEs that permits exact geometry representation and high regularity basis functions. Two model problems are considered: the scalar elliptic equation and the advection-diffusion equation. Theoretical analysis proves that the adoption of a coarse correction grid is crucial in order to have the condition number of the preconditioned stiffness matrix independent from the number of subdomains, whenever the ratio between the coarse mesh size and the fine mesh size is kept fixed. Numerical tests for the scalar elliptic equation (in 2D and 3D on trivial and non-trivial domains) confirm the theory. The preconditioner is then applied to the advection-diffusion equation in 2D and 3D. Again, the numerical results shows that the condition number of the preconditioned linear system scales with the number of subdomains up to 8100 processors, eventually with SUPG stabilization. The tests are implemented in C programming language on the top of PETSc library.File | Dimensione | Formato | |
---|---|---|---|
phd_unimi_R09930.pdf
accesso aperto
Dimensione
1.57 MB
Formato
Adobe PDF
|
1.57 MB | Adobe PDF | Visualizza/Apri |
I documenti in UNITESI sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.
https://hdl.handle.net/20.500.14242/74025
URN:NBN:IT:UNIMI-74025