Authors: Wobker, Hilmar
Title: Efficient multilevel solvers and high performance computing techniques for the finite element simulation of large-scale elasticity problems
Language (ISO): en
Abstract: 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.
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.
Subject Headings: Iterativer Löser
Multilevel
Mehrgitter
Gebietszerlegung
Mehrgitter-Krylov
Nicht-konformes Mehrgitter
ScaRC
Adaptive Grobgitterkorrektur
Minimale Überlappung
Sattelpunkt-Problem
Schurkomplement-Vorkonditionierer
Vanka
Gedämpftes Newton-Raphson
Globales Newton-Raphson
Inexaktes Newton-Raphson
Liniensuche
FEAST
Hardware-orientiert
Großskalig
Paralleles Rechnen
Parallele Effizienz
Finite-Elemente-Methode
Gemischte Formulierung
LBB Stabilisierung
Irreguläres Gitter
Festkörpermechanik
Strukturmechanik
Elastizität
Elastostatisch
Elastodynamisch
Zeitabhängig
Inkompressibles Material
Finite Deformation
Große Deformation
Volumenversteifung
Schubversteifung
Iterative solver
Multilevel
Multigrid
Domain decomposition
Multigrid-Krylov
Nonconforming multigrid
ScaRC
Adaptive coarse grid correction
Minimal overlap
Saddle point problem
Schur complement preconditioning
Vanka
Newton-Raphson
Damped Newton-Raphson
Global Newton-Raphson
Inexact Newton-Raphson
Line-search
FEAST
High performance computing
Hardware-oriented
Large-scale
Parallel computing
Parallel efficiency
Finite element method
Mixed formulation
LBB stabilisation
Equal-order finite elements
Irregular grids
Solid mechanics
Structural mechanics
Elasticity
Elastostatic
Elastodynamic
Transient
Incompressible material
Finite deformation
Large deformation
Volume locking
Shear locking
URI: http://hdl.handle.net/2003/26998
http://dx.doi.org/10.17877/DE290R-497
Issue Date: 2010-03-24T10:08:12Z
Appears in Collections:Lehrstuhl III Angewandte Mathematik und Numerik

Files in This Item:
File Description SizeFormat 
Wobker_Dissertation_03-2010.pdfDNB7.43 MBAdobe PDFView/Open


This item is protected by original copyright



This item is protected by original copyright rightsstatements.org