Full metadata record
DC FieldValueLanguage
dc.contributor.advisorTurek, S.-
dc.contributor.authorWobker, Hilmar-
dc.date.accessioned2010-03-24T10:08:12Z-
dc.date.available2010-03-24T10:08:12Z-
dc.date.issued2010-03-24T10:08:12Z-
dc.identifier.urihttp://hdl.handle.net/2003/26998-
dc.identifier.urihttp://dx.doi.org/10.17877/DE290R-497-
dc.description.abstractBei 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.de
dc.description.abstractIn 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.en
dc.language.isoen-
dc.subjectIterativer Löserde
dc.subjectMultilevelde
dc.subjectMehrgitterde
dc.subjectGebietszerlegungde
dc.subjectMehrgitter-Krylovde
dc.subjectNicht-konformes Mehrgitterde
dc.subjectScaRCde
dc.subjectAdaptive Grobgitterkorrekturde
dc.subjectMinimale Überlappungde
dc.subjectSattelpunkt-Problemde
dc.subjectSchurkomplement-Vorkonditioniererde
dc.subjectVankade
dc.subjectGedämpftes Newton-Raphsonde
dc.subjectGlobales Newton-Raphsonde
dc.subjectInexaktes Newton-Raphsonde
dc.subjectLiniensuchede
dc.subjectFEASTde
dc.subjectHardware-orientiertde
dc.subjectGroßskaligde
dc.subjectParalleles Rechnende
dc.subjectParallele Effizienzde
dc.subjectFinite-Elemente-Methodede
dc.subjectGemischte Formulierungde
dc.subjectLBB Stabilisierungde
dc.subjectIrreguläres Gitterde
dc.subjectFestkörpermechanikde
dc.subjectStrukturmechanikde
dc.subjectElastizitätde
dc.subjectElastostatischde
dc.subjectElastodynamischde
dc.subjectZeitabhängigde
dc.subjectInkompressibles Materialde
dc.subjectFinite Deformationde
dc.subjectGroße Deformationde
dc.subjectVolumenversteifungde
dc.subjectSchubversteifungde
dc.subjectIterative solveren
dc.subjectMultilevelen
dc.subjectMultigriden
dc.subjectDomain decompositionen
dc.subjectMultigrid-Kryloven
dc.subjectNonconforming multigriden
dc.subjectScaRCen
dc.subjectAdaptive coarse grid correctionen
dc.subjectMinimal overlapen
dc.subjectSaddle point problemen
dc.subjectSchur complement preconditioningen
dc.subjectVankaen
dc.subjectNewton-Raphsonen
dc.subjectDamped Newton-Raphsonen
dc.subjectGlobal Newton-Raphsonen
dc.subjectInexact Newton-Raphsonen
dc.subjectLine-searchen
dc.subjectFEASTen
dc.subjectHigh performance computingen
dc.subjectHardware-orienteden
dc.subjectLarge-scaleen
dc.subjectParallel computingen
dc.subjectParallel efficiencyen
dc.subjectFinite element methoden
dc.subjectMixed formulationen
dc.subjectLBB stabilisationen
dc.subjectEqual-order finite elementsen
dc.subjectIrregular gridsen
dc.subjectSolid mechanicsen
dc.subjectStructural mechanicsen
dc.subjectElasticityen
dc.subjectElastostaticen
dc.subjectElastodynamicen
dc.subjectTransienten
dc.subjectIncompressible materialen
dc.subjectFinite deformationen
dc.subjectLarge deformationen
dc.subjectVolume lockingen
dc.subjectShear lockingen
dc.subject.ddc510-
dc.titleEfficient multilevel solvers and high performance computing techniques for the finite element simulation of large-scale elasticity problemsen
dc.typeText-
dc.contributor.refereeSuttmeier, F.-T.-
dc.date.accepted2010-03-11-
dc.type.publicationtypedoctoralThesis-
dc.identifier.urnurn:nbn:de:hbz:290-2003/26998-3-
dcterms.accessRightsopen access-
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



Items in Eldorado are protected by copyright, with all rights reserved, unless otherwise indicated.