Efficient multilevel solvers and high performance computing techniques for the finite element simulation of large-scale elasticity problems
Creators: name: Wobker, Hilmar
Publisher: name: Technische Universität Dortmund
Year: 2010
Subjects: subject: Iterativer Löser
subject: Multilevel
subject: Mehrgitter
subject: Gebietszerlegung
subject: Mehrgitter-Krylov
subject: Nicht-konformes Mehrgitter
subject: ScaRC
subject: Adaptive Grobgitterkorrektur
subject: Minimale Überlappung
subject: Sattelpunkt-Problem
subject: Schurkomplement-Vorkonditionierer
subject: Vanka
subject: Gedämpftes Newton-Raphson
subject: Globales Newton-Raphson
subject: Inexaktes Newton-Raphson
subject: Liniensuche
subject: FEAST
subject: Hardware-orientiert
subject: Großskalig
subject: Paralleles Rechnen
subject: Parallele Effizienz
subject: Finite-Elemente-Methode
subject: Gemischte Formulierung
subject: LBB Stabilisierung
subject: Irreguläres Gitter
subject: Festkörpermechanik
subject: Strukturmechanik
subject: Elastizität
subject: Elastostatisch
subject: Elastodynamisch
subject: Zeitabhängig
subject: Inkompressibles Material
subject: Finite Deformation
subject: Große Deformation
subject: Volumenversteifung
subject: Schubversteifung
subject: Iterative solver
subject: Multigrid
subject: Domain decomposition
subject: Multigrid-Krylov
subject: Nonconforming multigrid
subject: Adaptive coarse grid correction
subject: Minimal overlap
subject: Saddle point problem
subject: Schur complement preconditioning
subject: Newton-Raphson
subject: Damped Newton-Raphson
subject: Global Newton-Raphson
subject: Inexact Newton-Raphson
subject: Line-search
subject: High performance computing
subject: Hardware-oriented
subject: Large-scale
subject: Parallel computing
subject: Parallel efficiency
subject: Finite element method
subject: Mixed formulation
subject: LBB stabilisation
subject: Equal-order finite elements
subject: Irregular grids
subject: Solid mechanics
subject: Structural mechanics
subject: Elasticity
subject: Elastostatic
subject: Elastodynamic
subject: Transient
subject: Incompressible material
subject: Finite deformation
subject: Large deformation
subject: Volume locking
subject: Shear locking
subject: 510

description: Bei der Simulation realistischer strukturmechanischer Probleme können Gleichungssysteme mit mehreren hundert Millionen Unbekannten entstehen. Für die effiziente Lösung solcher Systeme sind parallele Multilevel-Methoden unerlässlich, die in der Lage sind, die Leistung moderner Hardware-Technologien auszuschöpfen. Die Finite-Elemente- und Löser-Toolbox FEAST, die auf die Behandlung skalarer Gleichungen ausgelegt ist, verfolgt genau dieses Ziel. FEAST kombiniert Hardware-orientierte Implementierungstechniken mit einer Multilevel-Gebietszerlegungsmethode namens ScaRC. In der vorliegenden Arbeit wird ein Konzept entwickelt, multivariate Elastizitätsprobleme basierend auf der FEAST-Bibliothek zu lösen. Die generelle Herangehensweise besteht darin, die Lösung multivariater Probleme auf die Lösung einer Reihe von skalaren Problemen zurückzuführen. Dieser Ansatz ermöglicht eine strikte Trennung von skalaren "low level" Kernfunktionalitäten (in Form der FEAST-Bibliothek) und multivariatem "high level" Anwendungscode (in Form des Elastizitätsproblems), was aus Sicht der Softwareentwicklungstechnik sehr vorteilhaft ist: Alle Bemühungen zur Verbesserung der Hardware-Effizienz, sowie Anpassungen an zukünftige technologische Entwicklungen können auf skalare Operationen beschränkt werden, während die multivariate Anwendung automatisch von diesen Erweiterungen profitiert. Im ersten Teil der Arbeit werden substantielle Verbesserungen der skalaren ScaRC-Löser entwickelt, die dann als essentielle Bausteine zur Lösung multivariater Elastizitätsprobleme eingesetzt werden. Ausführliche numerische Untersuchungen zeigen, wie sich die Effizienz der skalaren FEAST-Bibliothek auf den multivariaten Lösungsprozess überträgt. Die Löserstrategie wird dann auf nichtlineare Probleme der Elastizität mit finiter Deformation angewandt. Durch Einsatz einer Liniensuche-Methode wird die Robustheit des Newton-Raphson-Verfahrens signifikant erhöht. Es werden verschiedene Strategien miteinander verglichen, wie genau die linearen Probleme innerhalb der nichtlinearen Iteration zu lösen sind. Zur Behandlung der wichtigen Klasse von (fast) inkompressiblen Materialien wird eine gemischte Verschiebung/Druck-Formulierung gewählt, die mit Hilfe von stabilisierten bilinearen finiten Elementen (Q1/Q1) diskretisiert wird. Eine erweiterte Version der klassischen "Druck-Poisson"-Stabilisierung wird präsentiert, die auch für hochgradig irreguläre Gitter geeignet ist. Es werden Vor- und Nachteile der Q1/Q1-Diskretisierung erörtert, insbesondere in Bezug auf zeitabhängige Rechnungen. Zwei Löser-Klassen zur Behandlung der entstehenden Sattelpunkt-Probleme werden beschrieben und miteinander verglichen: einerseits verschiedene Arten von (beschleunigten) entkoppelten Lösern (Uzawa, Druck-Schurkomplement-Methoden, Block-Vorkonditionierer), andererseits gekoppelte Mehrgitter-Verfahren mit Vanka-Glättern. Effiziente Schurkomplement-Vorkonditionierer, die für die erste Löser-Klasse notwendig sind, werden im Rahmen statischer und zeitabhängiger Probleme besprochen. Die zentrale Strategie, die Lösung multivariater Systeme auf die Lösung skalarer Systeme zu reduzieren, ist nur im Falle der entkoppelten Lösungsmethoden anwendbar. Es wird gezeigt, dass für die Klasse der Elastizitätsprobleme, die in dieser Arbeit behandelt werden, die entkoppelten Löser den gekoppelten hinsichtlich numerischer und paralleler Effizienz deutlich überlegen sind.
description: In the simulation of realistic solid mechanical problems, linear equation systems with hundreds of million unknowns can arise. For the efficient solution of such systems, parallel multilevel methods are mandatory that are able to exploit the capabilities of modern hardware technologies. The finite element and solution toolbox FEAST, which is designed to solve scalar equations, pursues exactly this goal. It combines hardware-oriented implementation techniques with a multilevel domain decomposition method called ScaRC that achieves high numerical and parallel efficiency. In this thesis a concept is developed to solve multivariate elasticity problems based on the FEAST library. The general strategy is to reduce the solution of multivariate problems to the solution of a series of scalar problems. This approach facilitates a strict separation of 'low level' scalar kernel functionalities (in the form of the FEAST library) and 'high level' multivariate application code (in the form of the elasticity problem), which is very attractive from a software-engineering point of view: All efforts to improve hardware-efficiency and adaptations to future technology trends can be restricted to scalar operations, and the multivariate application automatically benefits from these enhancements. In the first part of the thesis, substantial improvements of the scalar ScaRC solvers are developed, which are then used as essential building blocks for the efficient solution of multivariate elasticity problems. Extensive numerical studies demonstrate how the efficiency of the scalar FEAST library transfers to the multivariate solution process. The solver strategy is then applied to treat nonlinear problems of finite deformation elasticity. A line-search method is used to significantly increase the robustness of the Newton-Raphson method, and different strategies are compared how to choose the accuracy of the linear system solves within the nonlinear iteration. In order to treat the important class of (nearly) incompressible material, a mixed displacement/pressure formulation is used which is discretised with stabilised bilinear finite elements (Q1/Q1). An enhanced version of the classical 'pressure Poisson' stabilisation is presented which is suitable for highly irregular meshes. Advantages and disadvantages of the Q1/Q1 discretisation are discussed, especially in the context of transient computations. Two solver classes for the resulting saddle point systems are described and compared: on the one hand various kinds of (accelerated) segregated solvers (Uzawa, pressure Schur complement methods, block preconditioners), and on the other hand coupled multigrid solvers with Vanka-smoothers. Efficient Schur complement preconditioners, which are required for the former class, are discussed for the static and the transient case. The main strategy to reduce the solution of multivariate systems to the solution of scalar systems is only applicable in the case of segregated methods. It is shown that for the class of elasticity problems considered in this thesis, segregated solvers are clearly superior to Vanka-type solvers in terms of numerical and parallel efficiency.