On the Singular Neumann Problem in Linear Elasticity

The Neumann problem of linear elasticity is singular with a kernel formed by the rigid motions of the body. There are several tricks that are commonly used to obtain a non-singular linear system. However, they often cause reduced accuracy or lead to poor convergence of the iterative solvers. In this paper, different well-posed formulations of the problem are studied through discretization by the finite element method, and preconditioning strategies based on operator preconditioning are discussed. For each formulation we derive preconditioners that are independent of the discretization parameter. Preconditioners that are robust with respect to the first Lam\'e constant are constructed for the pure displacement formulations, while a preconditioner that is robust in both Lam\'e constants is constructed for the mixed formulation. It is shown that, for convergence in the first Sobolev norm, it is crucial to respect the orthogonality constraint derived from the continuous problem. Based on this observation a modification to the conjugate gradient method is proposed that achieves optimal error convergence of the computed solution

An algebraic multigrid method for Q2−Q1Q_2-Q_1 mixed discretizations of the Navier-Stokes equations

Algebraic multigrid (AMG) preconditioners are considered for discretized systems of partial differential equations (PDEs) where unknowns associated with different physical quantities are not necessarily co-located at mesh points. Specifically, we investigate a $Q_2-Q_1$ mixed finite element discretization of the incompressible Navier-Stokes equations where the number of velocity nodes is much greater than the number of pressure nodes. Consequently, some velocity degrees-of-freedom (dofs) are defined at spatial locations where there are no corresponding pressure dofs. Thus, AMG approaches leveraging this co-located structure are not applicable. This paper instead proposes an automatic AMG coarsening that mimics certain pressure/velocity dof relationships of the $Q_2-Q_1$ discretization. The main idea is to first automatically define coarse pressures in a somewhat standard AMG fashion and then to carefully (but automatically) choose coarse velocity unknowns so that the spatial location relationship between pressure and velocity dofs resembles that on the finest grid. To define coefficients within the inter-grid transfers, an energy minimization AMG (EMIN-AMG) is utilized. EMIN-AMG is not tied to specific coarsening schemes and grid transfer sparsity patterns, and so it is applicable to the proposed coarsening. Numerical results highlighting solver performance are given on Stokes and incompressible Navier-Stokes problems.Comment: Submitted to a journa

Grid generation for the solution of partial differential equations

A general survey of grid generators is presented with a concern for understanding why grids are necessary, how they are applied, and how they are generated. After an examination of the need for meshes, the overall applications setting is established with a categorization of the various connectivity patterns. This is split between structured grids and unstructured meshes. Altogether, the categorization establishes the foundation upon which grid generation techniques are developed. The two primary categories are algebraic techniques and partial differential equation techniques. These are each split into basic parts, and accordingly are individually examined in some detail. In the process, the interrelations between the various parts are accented. From the established background in the primary techniques, consideration is shifted to the topic of interactive grid generation and then to adaptive meshes. The setting for adaptivity is established with a suitable means to monitor severe solution behavior. Adaptive grids are considered first and are followed by adaptive triangular meshes. Then the consideration shifts to the temporal coupling between grid generators and PDE-solvers. To conclude, a reflection upon the discussion, herein, is given

Geometry–aware finite element framework for multi–physics simulations: an algorithmic and software-centric perspective

In finite element simulations, the handling of geometrical objects and their discrete representation is a critical aspect in both serial and parallel scientific software environments. The development of codes targeting such envinronments is subject to great development effort and man-hours invested. In this thesis we approach these issues from three fronts. First, stable and efficient techniques for the transfer of discrete fields between non matching volume or surface meshes are an essential ingredient for the discretization and numerical solution of coupled multi-physics and multi-scale problems. In particular L2-projections allows for the transfer of discrete fields between unstructured meshes, both in the volume and on the surface. We present an algorithm for parallelizing the assembly of the L2-transfer operator for unstructured meshes which are arbitrarily distributed among different processes. The algorithm requires no a priori information on the geometrical relationship between the different meshes. Second, the geometric representation is often a limiting factor which imposes a trade-off between how accurately the shape is described, and what methods can be employed for solving a system of differential equations. Parametric finite-elements and bijective mappings between polygons or polyhedra allow us to flexibly construct finite element discretizations with arbitrary resolutions without sacrificing the accuracy of the shape description. Such flexibility allows employing state-of-the-art techniques, such as geometric multigrid methods, on meshes with almost any shape.t, the way numerical techniques are represented in software libraries and approached from a development perspective, affect both usability and maintainability of such libraries. Completely separating the intent of high-level routines from the actual implementation and technologies allows for portable and maintainable performance. We provide an overview on current trends in the development of scientific software and showcase our open-source library utopia

Schnelle Löser für partielle Differentialgleichungen

The workshop Schnelle Löser für partielle Differentialgleichungen, organised by Randolph E. Bank (La Jolla), Wolfgang Hackbusch(Leipzig), Gabriel Wittum (Heidelberg) was held May 22nd - May 28th, 2005. This meeting was well attended by 47 participants with broad geographic representation from 9 countries and 3 continents. This workshop was a nice blend of researchers with various backgrounds

Multilevel Artificial Neural Network Training for Spatially Correlated Learning

Multigrid modeling algorithms are a technique used to accelerate relaxation models running on a hierarchy of similar graphlike structures. We introduce and demonstrate a new method for training neural networks which uses multilevel methods. Using an objective function derived from a graph-distance metric, we perform orthogonally-constrained optimization to find optimal prolongation and restriction maps between graphs. We compare and contrast several methods for performing this numerical optimization, and additionally present some new theoretical results on upper bounds of this type of objective function. Once calculated, these optimal maps between graphs form the core of Multiscale Artificial Neural Network (MsANN) training, a new procedure we present which simultaneously trains a hierarchy of neural network models of varying spatial resolution. Parameter information is passed between members of this hierarchy according to standard coarsening and refinement schedules from the multiscale modelling literature. In our machine learning experiments, these models are able to learn faster than default training, achieving a comparable level of error in an order of magnitude fewer training examples.Comment: Manuscript (24 pages) and Supplementary Material (4 pages). Updated January 2019 to reflect new formulation of MsANN structure and new training procedur

Towards a robust Terra

In this work mantle convection simulation with Terra is investigated from a numerical point of view, theoretical analysis as well as practical tests are performed. The stability criteria for the numerical formulation of the physical model will be made clear. For the incompressible case and the Terra specific treatment of the anelastic approximation, two inf-sup stable grid modifications are presented, which are both compatible with hanging nodes. For the Q1hQ12h element pair a simple numeric test is introduced to prove the stability for any given grid. For the Q1h Pdisc 12h element pair and 1-regular refinements with hangig nodes an existing general proof can be adopted. The influence of the slip boundary condition is found to be destabilizing. For the incompressible case a cure can be adopted from the literature. The necessary conditions for the expansion of the stability results to the anelastic approximation will be pointed out. A numerical framework is developed in order to measure the effect of different numerical approaches to improve the handling of strongly varying viscosity. The framework is applied to investigate how block smoothers with different block sizes, combination of different block smoothers, different prolongation schemes and semi coarsening influence the multigrid performance. A regression-test framework for Terra will be briefly introduced

On Multilevel Methods Based on Non-Nested Meshes

This thesis is concerned with multilevel methods for the efficient solution of partial differential equations in the field of scientific computing. Further, emphasis is put on an extensive study of the information transfer between finite element spaces associated with non-nested meshes. For the discretization of complicated geometries with a finite element method, unstructured meshes are often beneficial as they can easily be adjusted to the shape of the computational domain. Such meshes, and thus the corresponding discrete function spaces, do not allow for straightforward multilevel hierarchies that could be exploited to construct fast solvers. In the present thesis, we present a class of "semi-geometric" multilevel iterations, which are based on hierarchies of independent, non-nested meshes. This is realized by a variational approach such that the images of suitable prolongation operators in the next (finer) space recursively determine the coarse level spaces. The semi-geometric concept is of very general nature compared with other methods relying on geometric considerations. This is reflected in the relatively loose relations of the employed meshes to each other. The specific benefit of the approach based on non-nested meshes is the flexibility in the choice of the coarse meshes, which can, for instance, be generated independently by standard methods. The resolution of the boundaries of the actual computational domain in the constructed coarse level spaces is a characteristic feature of the devised class of methods. The flexible applicability and the efficiency of the presented solution methods is demonstrated in a series of numerical experiments. We also explain the practical implementation of the semi-geometric ideas and concrete transfer concepts between non-nested meshes. Moreover, an extension to a semi-geometric monotone multigrid method for the solution of variational inequalities is discussed. We carry out the analysis of the convergence and preconditioning properties, respectively, in the framework of the theory of subspace correction methods. Our technical considerations yield a quasi-optimal result, which we prove for general, shape regular meshes by local arguments. The relevant properties of the operators for the prolongation between non-nested finite element spaces are the H1-stability and an L2-approximation property as well as the locality of the transfer. This thesis is a contribution to the development of fast solvers for equations on complicated geometries with focus on geometric techniques (as opposed to algebraic ones). Connections to other approaches are carefully elaborated. In addition, we examine the actual information transfer between non-nested finite element spaces. In a novel study, we combine theoretical, practical and experimental considerations. A thourough investigation of the qualitative properties and a quantitative analysis of the differences of individual transfer concepts to each other lead to new results on the information transfer as such. Finally, by the introduction of a generalized projection operator, the pseudo-L2-projection, we obtain a significantly better approximation of the actual L2-orthogonal projection than other approaches from the literature.Nicht-geschachtelte Gitter in Multilevel-Verfahren Diese Arbeit beschäftigt sich mit Multilevel-Verfahren zur effizienten Lösung von Partiellen Differentialgleichungen im Bereich des Wissenschaftlichen Rechnens. Dabei liegt ein weiterer Schwerpunkt auf der eingehenden Untersuchung des Informationsaustauschs zwischen Finite-Elemente-Räumen zu nicht-geschachtelten Gittern. Zur Diskretisierung von komplizierten Geometrien mit einer Finite-Elemente-Methode sind unstrukturierte Gitter oft von Vorteil, weil sie der Form des Rechengebiets einfacher angepasst werden können. Solche Gitter, und somit die zugehörigen diskreten Funktionenräume, besitzen im Allgemeinen keine leicht zugängliche Multilevel-Struktur, die sich zur Konstruktion schneller Löser ausnutzen ließe. In der vorliegenden Arbeit stellen wir eine Klasse "semi-geometrischer" Multilevel-Iterationen vor, die auf Hierarchien voneinander unabhängiger, nicht-geschachtelter Gitter beruhen. Dabei bestimmen in einem variationellen Ansatz rekursiv die Bilder geeigneter Prolongationsoperatoren im jeweils folgenden (feineren) Raum die Grobgitterräume. Das semi-geometrische Konzept ist sehr allgemeiner Natur verglichen mit anderen Verfahren, die auf geometrischen Überlegungen beruhen. Dies zeigt sich in der verhältnismäßig losen Beziehung der verwendeten Gitter zueinander. Der konkrete Nutzen des Ansatzes mit nicht-geschachtelten Gittern ist die Flexibilität der Wahl der Grobgitter. Diese können beispielsweise unabhängig mit Standardverfahren generiert werden. Die Auflösung des Randes des tatsächlichen Rechengebiets in den konstruierten Grobgitterräumen ist eine Eigenschaft der entwickelten Verfahrensklasse. Die flexible Einsetzbarkeit und die Effizienz der vorgestellten Lösungsverfahren zeigt sich in einer Reihe von numerischen Experimenten. Dazu geben wir Hinweise zur praktischen Umsetzung der semi-geometrischen Ideen und konkreter Transfer-Konzepte zwischen nicht-geschachtelten Gittern. Darüber hinaus wird eine Erweiterung zu einem semi-geometrischen monotonen Mehrgitterverfahren zur Lösung von Variationsungleichungen untersucht. Wir führen die Analysis der Konvergenz- bzw. Vorkonditionierungseigenschaften im Rahmen der Theorie der Teilraumkorrekturmethoden durch. Unsere technische Ausarbeitung liefert ein quasi-optimales Resultat, das wir mithilfe lokaler Argumente für allgemeine, shape-reguläre Gitterfamilien beweisen. Als relevante Eigenschaften der Operatoren zur Prolongation zwischen nicht-geschachtelten Finite-Elemente-Räumen erweisen sich die H1-Stabilität und eine L2-Approximationseigenschaft sowie die Lokalität des Transfers. Diese Arbeit ist ein Beitrag zur Entwicklung schneller Löser für Gleichungen auf komplizierten Gebieten mit Schwerpunkt auf geometrischen Techniken (im Unterschied zu algebraischen). Verbindungen zu anderen Ansätzen werden sorgfältig aufgezeigt. Daneben untersuchen wir den Informationsaustausch zwischen nicht-geschachtelten Finite-Elemente-Räumen als solchen. In einer neuartigen Studie verbinden wir theoretische, praktische und experimentelle Überlegungen. Eine sorgfältige Prüfung der qualitativen Eigenschaften sowie eine quantitative Analyse der Unterschiede verschiedener Transfer-Konzepte zueinander führen zu neuen Ergebnissen bezüglich des Informationsaustauschs selbst. Schließlich erreichen wir durch die Einführung eines verallgemeinerten Projektionsoperators, der Pseudo-L2-Projektion, eine deutlich bessere Approximation der eigentlichen L2-orthogonalen Projektion als andere Ansätze aus der Literatur