Least-Squares FEM: Literature Review

During the last years the interest in least squares finite element methods (LSFEM) has grown continuously. Least squares finite element methods offer some advantages over the widely used Galerkin variational principle. One reason is the ability to cope with first order differential operators without special treatment as required by the Galerkin FEM. The other reason comes from the numerical point of view, where the LSFEM leads to symmetric positive definite matrices which can be solved very efficiently under some conditions. This report gives an overview about the recent literature which appeared in the field of least squares finite element methods and summarises the essential results and facts about the LSFEM.Während der letzten Jahre hat das Interesse an Least Squares Finite Element Methoden (LSFEM) stetig zugenommen. Least Squares Finite Element Methoden bieten einige Vorteile gegenüber dem etablierten Galerkin Variationsansatz. So können Differentialoperatoren erster Ordnung ohne besondere numerische Techniken, wie z.B. Stabilisierung, direkt behandelt werden. Ein anderer Grund für den Einsatz der LSFEM liegt in den entstehenden algebraischen Gleichungssystemen, die immer symmetrisch positiv definit sind und unter bestimmten Vorraussetzungen eine effiziente Lösung ermöglichen.Dieser Bericht gibt einen Überblick über die aktuelle Literatur zur LSFEM und faßt die entscheidenden Ergebnisse zusammen

Hp-spectral Methods for Structural Mechanics and Fluid Dynamics Problems

We consider the usage of higher order spectral element methods for the solution of problems in structures and fluid mechanics areas. In structures applications we study different beam theories, with mixed and displacement based formulations, consider the analysis of plates subject to external loadings, and large deformation analysis of beams with continuum based formulations. Higher order methods alleviate the problems of locking that have plagued finite element method applications to structures, and also provide for spectral accuracy of the solutions. For applications in computational fluid dynamics areas we consider the driven cavity problem with least squares based finite element methods. In the context of higher order methods, efficient techniques need to be devised for the solution of the resulting algebraic systems of equations and we explore the usage of element by element bi-orthogonal conjugate gradient solvers for solving problems effectively along with domain decomposition algorithms for fluid problems. In the context of least squares finite element methods we also explore the usage of Multigrid techniques to obtain faster convergence of the the solutions for the problems of interest. Applications of the traditional Lagrange based finite element methods with the Penalty finite element method are presented for modelling porous media flow problems. Finally, we explore applications to some CFD problems namely, the flow past a cylinder and forward facing step

Algebraic multigrid for stabilized finite element discretizations of the Navier Stokes equation

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Aeronautics and Astronautics, 2002.Includes bibliographical references (p. 141-152).A multilevel method for the solution of systems of equations generated by stabilized Finite Element discretizations of the Euler and Navier Stokes equations on generalized unstructured grids is described. The method is based on an elemental agglomeration multigrid which produces a hierarchical sequence of coarse subspaces. Linear combinations of the basis functions from a given space form the next subspace and the use of the Galerkin Coarse Grid Approximation (GCA) within an Algebraic Multigrid (AMG) context properly defines the hierarchical sequence. The multigrid coarse spaces constructed by the elemental agglomeration algorithm are based on a semi-coarsening scheme designed to reduce grid anisotropy. The multigrid transfer operators are induced by the graph of the coarse space mesh and proper consideration is given to the boundary conditions for an accurate representation of the coarse space operators. A generalized line implicit relaxation scheme is also described where the lines are constructed to follow the direction of strongest coupling. The solution algorithm is motivated by the decomposition of the system characteristics into acoustic and convective modes. Analysis of the application of elemental agglomeration AMG (AMGe) to stabilized numerical schemes shows that a characteristic length based rescaling of the numerical stabilization is necessary for a consistent multigrid representation.by Tolulope Olawale Okusanya.Ph.D

Simplex space-time meshes in thermally coupled two-phase flow simulations of mold filling

The quality of plastic parts produced through injection molding depends on many factors. Especially during the filling stage, defects such as weld lines, burrs, or insufficient filling can occur. Numerical methods need to be employed to improve product quality by means of predicting and simulating the injection molding process. In the current work, a highly viscous incompressible non-isothermal two-phase flow is simulated, which takes place during the cavity filling. The injected melt exhibits a shear-thinning behavior, which is described by the Carreau-WLF model. Besides that, a novel discretization method is used in the context of 4D simplex space-time grids [2]. This method allows for local temporal refinement in the vicinity of, e.g., the evolving front of the melt [10]. Utilizing such an adaptive refinement can lead to locally improved numerical accuracy while maintaining the highest possible computational efficiency in the remaining of the domain. For demonstration purposes, a set of 2D and 3D benchmark cases, that involve the filling of various cavities with a distributor, are presented.Comment: 14 pages, 11 Figures, 4 Table

Well-posedness and Robust Preconditioners for the Discretized Fluid-Structure Interaction Systems

In this paper we develop a family of preconditioners for the linear algebraic systems arising from the arbitrary Lagrangian-Eulerian discretization of some fluid-structure interaction models. After the time discretization, we formulate the fluid-structure interaction equations as saddle point problems and prove the uniform well-posedness. Then we discretize the space dimension by finite element methods and prove their uniform well-posedness by two different approaches under appropriate assumptions. The uniform well-posedness makes it possible to design robust preconditioners for the discretized fluid-structure interaction systems. Numerical examples are presented to show the robustness and efficiency of these preconditioners.Comment: 1. Added two preconditioners into the analysis and implementation 2. Rerun all the numerical tests 3. changed title, abstract and corrected lots of typos and inconsistencies 4. added reference

Simultaneous solution approaches for large optimization problems

AbstractIn this paper, efficient simultaneous strategies are presented for the optimization of practical problems involving PDE-models. In particular, reduced sequential quadratic programming methods for problems with only few influence variables and simultaneous quadratic programming iterations are discussed. As a result we obtain algorithms whose overall computational complexity is reduced considerably in comparison to a black-box approach

Potential-based Formulations of the Navier-Stokes Equations and their Application

Based on a Clebsch-like velocity representation and a combination of classical variational principles for the special cases of ideal and Stokes flow a novel discontinuous Lagrangian is constructed; it bypasses the known problems associated with non-physical solutions and recovers the classical Navier-Stokes equations together with the balance of inner energy in the limit when an emerging characteristic frequency parameter tends to infinity. Additionally, a generalized Clebsch transformation for viscous flow is established for the first time. Next, an exact first integral of the unsteady, three-dimensional, incompressible Navier-Stokes equations is derived; following which gauge freedoms are explored leading to favourable reductions in the complexity of the equation set and number of unknowns, enabling a self-adjoint variational principle for steady viscous flow to be constructed. Concurrently, appropriate commonly occurring physical and auxiliary boundary conditions are prescribed, including establishment of a first integral for the dynamic boundary condition at a free surface. Starting from this new formulation, three classical flow problems are considered, the results obtained being in total agreement with solutions in the open literature. A new least-squares finite element method based on the first integral of the steady two-dimensional, incompressible, Navier-Stokes equations is developed, with optimal convergence rates established theoretically. The method is analysed comprehensively, thoroughly validated and shown to be competitive when compared to a corresponding, standard, primitive-variable, finite element formulation. Implementation details are provided, and the well-known problem of mass conservation addressed and resolved via selective weighting. The attractive positive definiteness of the resulting linear systems enables employment of a customized scalable algebraic multigrid method for efficient error reduction. The solution of several engineering related problems from the fields of lubrication and film flow demonstrate the flexibility and efficiency of the proposed method, including the case of unsteady flow, while revealing new physical insights of interest in their own right