A Block ILUT Smoother for Multipatch Geometries in Isogeometric Analysis
Roel Tielen (TU Delft - Numerical Analysis)
Matthias Moller (TU Delft - Numerical Analysis)
K. Vuik (TU Delft - Delft Institute of Applied Mathematics)
More Info
expand_more
Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.
Abstract
Since its introduction in [20], Isogeometric Analysis (IgA) has established itself as a viable alternative to the Finite Element Method (FEM). Solving the resulting linear systems of equations efficiently remains, however, challenging when high-order B-spline basis functions of order p> 1 are adopted for approximation. The use of Incomplete LU (ILU) type factorizations, like ILU(k) or ILUT, as a preconditioner within a Krylov method or as a smoother within a multigrid method is very effective, but costly [37]. In this paper, we investigate the use of a block ILUT smoother within a p-multigrid method, where the coarse grid correction is obtained at p= 1, and compare it to a global ILUT smoother in case of multipatch geometries. A spectral analysis indicates that the use of the block ILUT smoother improves the overall convergence rate of the resulting p-multigrid method. Numerical results, obtained for a variety of two dimensional benchmark problems, illustrate the potential of this block ILUT smoother for multipatch geometries.