Condition number analysis and preconditioning of the finite cell method

de Prenter, Frits; Verhoosel, C. V. Verhoosel; van Zwieten, G. J. van Zwieten; van Brummelen, E. H. van Brummelen

Condition number analysis and preconditioning of the finite cell method

Journal article (2017)

Authors

Frits de Prenter Eindhoven University of Technology

C. V. Verhoosel Verhoosel Eindhoven University of Technology

G. J. van Zwieten van Zwieten Eindhoven University of Technology

E. H. van Brummelen van Brummelen Eindhoven University of Technology

Affiliation

External organisation

Preconditioning Isogeometric Analysis Finite Cell Method Iterative solvers Condition number Immersed/fictitious domain methods

To reference this document use:

http://resolver.tudelft.nl/uuid:ba66a8e4-c7e9-4804-be46-fab1ff521700

More Info

expand_more

Published Date

2017

Language

English

Reuse Rights

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.

Affiliation

External organisation

Abstract

The (Isogeometric) Finite Cell Method–in which a domain is immersed in a structured background mesh–suffers from conditioning problems when cells with small volume fractions occur. In this contribution, we establish a rigorous scaling relation between the condition number of (I)FCM system matrices and the smallest cell volume fraction. Ill-conditioning stems either from basis functions being small on cells with small volume fractions, or from basis functions being nearly linearly dependent on such cells. Based on these two sources of ill-conditioning, an algebraic preconditioning technique is developed, which is referred to as Symmetric Incomplete Permuted Inverse Cholesky (SIPIC). A detailed numerical investigation of the effectivity of the SIPIC preconditioner in improving (I)FCM condition numbers and in improving the convergence speed and accuracy of iterative solvers is presented for the Poisson problem and for two- and three-dimensional problems in linear elasticity, in which Nitche's method is applied in either the normal or tangential direction. The accuracy of the preconditioned iterative solver enables mesh convergence studies of the finite cell method.