Multilevel spectral domain decomposition

Highly heterogeneous, anisotropic coefficients, e.g. in the simulation of carbon-fibre composite components, can lead to extremely challenging finite element systems. Direct solvers for the resulting large and sparse linear systems suffer from severe memory requirements and limited parallel scalabil...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Bastian, Peter (VerfasserIn) , Scheichl, Robert (VerfasserIn) , Seelinger, Linus (VerfasserIn) , Strehlow, Arne (VerfasserIn)
Dokumenttyp: Article (Journal) Kapitel/Artikel
Sprache:Englisch
Veröffentlicht: 15 Jun 2021
Ausgabe:Version 2
In: Arxiv
Year: 2021, Pages: 1-25
DOI:10.48550/arXiv.2106.06404
Online-Zugang:Resolving-System, kostenfrei, Volltext: https://doi.org/10.48550/arXiv.2106.06404
Verlag, kostenfrei, Volltext: http://arxiv.org/abs/2106.06404
Volltext
Verfasserangaben:Peter Bastian, Robert Scheichl, Linus Seelinger, and Arne Strehlow
Beschreibung
Zusammenfassung:Highly heterogeneous, anisotropic coefficients, e.g. in the simulation of carbon-fibre composite components, can lead to extremely challenging finite element systems. Direct solvers for the resulting large and sparse linear systems suffer from severe memory requirements and limited parallel scalability, while iterative solvers in general lack robustness. Two-level spectral domain decomposition methods can provide such robustness for symmetric positive definite linear systems, by using coarse spaces based on independent generalized eigenproblems in the subdomains. Rigorous condition number bounds are independent of mesh size, number of subdomains, as well as coefficient contrast. However, their parallel scalability is still limited by the fact that (in order to guarantee robustness) the coarse problem is solved via a direct method. In this paper, we introduce a multilevel variant in the context of subspace correction methods and provide a general convergence theory for its robust convergence for abstract, elliptic variational problems. Assumptions of the theory are verified for conforming, as well as for discontinuous Galerkin methods applied to a scalar diffusion problem. Numerical results illustrate the performance of the method for two- and three-dimensional problems and for various discretization schemes, in the context of scalar diffusion and linear elasticity.
Beschreibung:Identifizierung der Ressource nach: 15 Jun 2021
Gesehen am 23.09.2022
Beschreibung:Online Resource
DOI:10.48550/arXiv.2106.06404