Introduced in , Isogeometric Analysis (IgA) has become widely accepted in academia and industry. However, solving the resulting linear systems remains a challenging task. For instance, the condition number of the Poisson operator scales quadratically with the mesh width h, but, in contrast to standard Finite Elements, exponentially with the order of the approximation p . The performance of (standard) iterative solvers thus decreases fast for higher values of p. In this talk we propose an efficient solution strategy for IgA discretizations that is based on p- multigrid techniques used both as a solver and as a preconditioner in a Krylov subspace iteration method. The approach makes use of a hierarchy of B-spline based discretizations of different approximation orders, which is in contrast to (geometric) h-multigrid methods, where a hierarchy of coarser and finer meshes is constructed. The `coarse grid' correction is determined at level p = 1, which enables us to use established solution techniques developed for low-order Lagrange finite elements. Prolongation and restriction operators are defined as mappings between arbitrary spline spaces, solely determined by the generating knot vectors, allowing us to combine coarsening in both h and p, leading to a flexible hp-multigrid. Prelimenary numerical results are presented for different two-dimensional benchmark problems on non-trivial geometries. It follows from a Local Fourier Analysis , that the coarse grid correction and the smoothing procedure complement each other quite well. Moreover, the obtained convergence rates indicate that p-multigrid methods have the potential to efficiently solve IgA discretizations.
|Title of host publication||Proceedings of the 6th European Conference on Computational Mechanics (ECCM 6) and 7th European Conference on Computational Fluid Dynamics (ECFD 7)|
|Editors||R. Owen , R. de Borst, J. Reese|
|Place of Publication||Barcelona|
|Number of pages||12|
|Publication status||Published - 2018|
|Event||6th European Conference on Computational Mechanics ( Solids, Structures and Coupled Problems) and 7th European Conference on Computational Fluid Dynamics - Glasgow, United Kingdom|
Duration: 11 Jun 2018 → 15 Jun 2018
|Conference||6th European Conference on Computational Mechanics ( Solids, Structures and Coupled Problems) and 7th European Conference on Computational Fluid Dynamics|
|Abbreviated title||ECCM 6 ECFD 7|
|Period||11/06/18 → 15/06/18|
Bibliographical noteGreen Open Access added to TU Delft Institutional Repository ‘You share, we take care!’ – Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public.
- Isogeometric Analysis
- Iterative Solvers