Multilevel preconditioners for embedded enriched partition of unity approximations

• Published in: Advanced modeling and simulation in engineering sciences Vol. 5; no. 1; pp. 1 - 16
• Main Authors: Schweitzer, Marc Alexander, Ziegenhagel, Albert
• Format: Journal Article
• Language: English
• Published: Cham Springer International Publishing 18.05.2018; Springer; SpringerOpen
• Series: Non-intrusive computational techniques and related methods Non-intrusive computational techniques and related methods
• Subjects: Classical and Continuum Physics; Computational Science and Engineering; Domain decompositon; Engineering; Engineering Sciences; Generalized finite element method; Multilevel solver; Non-intrusive computational techniques and related methods; Partition of unity method; Research Article; Subspace correction; Theoretical and Applied Mechanics; Multilevel solver; Subspace correction; Primary 65N55; Domain decompositon; Generalized finite element method; 65N30; Partition of unity method
• ISSN: 2213-7467; 2213-7467
• DOI: 10.1186/s40323-018-0107-6

In this paper we are concerned with the non-invasive embedding of enriched partition of unity approximations in classical finite element simulations and the efficient solution of the resulting linear systems. The employed embedding is based on the partition of unity approach introduced in Schweitzer and Ziegenhagel (Embedding enriched partition of unity approximations in finite element simulations. In: Griebel M, Schweitzer MA, editors. Meshfree methods for partial differential equations VIII. Lecture notes in science and engineering, Cham, Springer International Publishing; 195–204, 2017 ) which is applicable to any finite element implementation and thus allows for a stable enrichment of e.g. commercial finite element software to improve the quality of its approximation properties in a non-invasive fashion. The major remaining challenge is the efficient solution of the arising linear systems. To this end, we apply classical subspace correction techniques to design non-invasive efficient multilevel solvers by blending a non-invasive algebraic multigrid method (applied to the finite element components) with a (geometric) multilevel solver (Griebel and Schweitzer in SIAM J Sci Comput 24:377–409, 2002 ; Schweitzer in Numer Math 118:307–28, 2011 ) (applied to the enriched embedded components). We present first numerical results in two and three space dimensions which clearly show the (close to) optimal performance of the proposed solver.