Efficient Multigrid based solvers for isogeometric analysis

Research output: Chapter in Book/Conference proceedings/Edited volumeConference contributionScientificpeer-review

13 Downloads (Pure)

Abstract

Introduced in [1], 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 [2]. 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 [3], 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.
Original languageEnglish
Title of host publicationProceedings of the 6th European Conference on Computational Mechanics (ECCM 6) and 7th European Conference on Computational Fluid Dynamics (ECFD 7)
EditorsR. Owen , R. de Borst, J. Reese
Place of PublicationBarcelona
PublisherECCOMAS
Pages1-12
Number of pages12
Publication statusPublished - 2018
Event6th European Conference on Computational Mechanics ( Solids, Structures and Coupled Problems) and 7th European Conference on Computational Fluid Dynamics - Glasgow, United Kingdom
Duration: 11 Jun 201815 Jun 2018
http://congress.cimne.com/eccm_ecfd2018/frontal/introduction.asp

Conference

Conference6th European Conference on Computational Mechanics ( Solids, Structures and Coupled Problems) and 7th European Conference on Computational Fluid Dynamics
Abbreviated titleECCM 6 ECFD 7
Country/TerritoryUnited Kingdom
CityGlasgow
Period11/06/1815/06/18
Internet address

Bibliographical note

Green 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.

Keywords

  • Isogeometric Analysis
  • p-Multigrid
  • Iterative Solvers

Fingerprint

Dive into the research topics of 'Efficient Multigrid based solvers for isogeometric analysis'. Together they form a unique fingerprint.

Cite this