Exact and Inexact IETI-DP Solvers for Continuous and Discontinuous Galerkin IgA Equations

MSc Christoph Hofer

Jan. 31, 2017, 3:30 p.m. S2 059

In this talk, we construct and investigate fast solvers for large-scale linear systems of algebraic equations arising from isogeometric analysis (IgA) of diffusion problems with heterogeneous diffusion coefficient on multipatch domains. In particular, we investigate the adaption of the Dual-Primal Finite Element Tearing and Interconnecting (FETI-DP) method to IgA, called Dual-Primal IsogEometric Tearing and Interconnecting (IETI-DP) method, [HL:HoferLanger:2016b]. We consider the cases where we have matching and non-matching meshes on the interfaces. In the latter case we use a discontinuous Galerkin (dG) method to couple the different patches. This requires a special extension of the IETI-DP method to the dG-IgA formulation, see [HL:HoferLanger:2016a]. We use ideas from the finite element case in order to formulate the corresponding IETI-DP method, called dG-IETI-DP. We design the dG-IETI-DP method in such a way that it can be seen as a IETI-DP method on an extended discrete interface space. Furthermore, we extend the dG-IETI-DP method to the case of non-matching interfaces due to incorrect segmentation producing gaps and overlaps in the domain decomposition, [HL:HoferLangerToulopoulos:2016b]. We investigate inexact versions utilizing multigrid methods, cf. [HL:HofreitherTakacs:2016a], for the formulation with an energy minimizing primal subspace.

These methods are highly suited for parallelization. We investigate the scaling behaviour up to 1024 cores. We present numerical results for complicated two and three dimensional domains. We observe that the condition number $\kappa$ behaves like $O((1 + \log(H/h))^2)$, and is robust with respect to jumping diffusion coefficients and changing mesh-sizes across patch interfaces. We also study the dependence of $\kappa$ on the underlying polynomial degree $p$ of the NURBS used.

[HL:HoferLanger:2016a]
C. Hofer and U. Langer.
Dual-Primal Isogeometric Tearing and Interconnecting solvers for multipatch {dG-IgA} equations.
Computer Methods in Applied Mechanics and Engineering, 2016, DOI 10.1016/j.cma.2016.03.031,
in Press.

[HL:HoferLanger:2016b]
C. Hofer, U. Langer.
Dual-Primal Isogeometric Tearing and Interconnecting Methods.
in "Contributions to PDE for Applications", ed. by P. Neittanmakki et al, Springer-ECCOMAS series “Computational Methods in Applied Sciences”, to appear, 2016,
also available as NFN Report 2015-39.

[HL:HofreitherTakacs:2016a]
C. Hofreither, S. Takacs.
Robust Multigrid for Isogeometric Analysis Based on Stable Splittings of Spline Spaces.
available as NFN Report 2016-45.

[HL:HoferLangerToulopoulos:2016b]
C. Hofer and U. Langer and I. Toulopoulos,
Discontinuous Galerkin Isogeometric Analysis on non-matching Segmentation: error estimates and efficient solvers.
available as RICAM-Report 2015-44.

Abstract