A parallel solver for adaptive finite element discretizations

A parallel solver for the adaptive finite element analysis is presented. The primary aim of this work has been to establish an efficient parallel computational procedure which requires only local computations to update the solution of the system of equations arising from the finite element discretization after a local mesh-adaptation step. For this reason a set of algorithms has been developed (two-level domain decomposition, recursive hierarchical mesh-refinement, selective solution-update of linear systems of equations) which operate upon general and easily available partitioning, meshing and linear systems solving algorithms. AMS subject classification 15A23, 65N50, 65N60     

Parallel finite element splitting-up method for parabolic problems

An efficient method for solving parabolic systems is presented. The proposed method is based on the splitting-up principle in which the problem is reduced to a series of independent 1D problems. This enables it to be used with parallel processors. We can solve multidimensional problems by applying only the 1D method and consequently avoid the difficulties in constructing a finite element space for multidimensional problems. The method is suitable for general domains as well as rectangular domains. Every 1D subproblem is solved by applying cubic B-splines. Several numerical examples are presented.     

Parallel triangularization of substructured finite element problems

Much of the computational effort of the finite element process involves the solution of a system of linear equations. The coefficient matrix of this system, known as the global stiffness matrix, is symmetric, positive definite, and generally sparse. An important technique for reducing the time required to solve this system is substructuring or matrix partitioning. Substructuring is based on the idea of dividing a structure into pieces, each of which can then be analyzed relatively independently. As a result of this division, each point in the finite element discretization is either interior to a substructure or on a boundary between substructures. Contributions to the global stiffness matrix from connections between boundary points form the K_bb matrix. This paper focuses on the triangularization of a general K_bb matrix on a parallel machine.     

An iterative adaptive finite element method for elliptic eigenvalue problems

We consider the task of resolving accurately the n$n$th eigenpair of a generalized eigenproblem rooted in some elliptic partial differential equation (PDE), using an adaptive finite element method (FEM). Conventional adaptive FEM algorithms call a generalized eigensolver after each mesh refinement step. This is not practical in our situation since the generalized eigensolver needs to calculate n$n$ eigenpairs after each mesh refinement step, it can switch the order of eigenpairs, and for repeated eigenvalues it can return an arbitrary linear combination of eigenfunctions from the corresponding eigenspace. In order to circumvent these problems, we propose a novel adaptive algorithm that only calls a generalized eigensolver once at the beginning of the computation, and then employs an iterative method to pursue a selected eigenvalue–eigenfunction pair on a sequence of locally refined meshes. Both Picard’s and Newton’s variants of the iterative method are presented. The underlying partial differential equation (PDE) is discretized with higher-order finite elements (hp$hp$-FEM) but the algorithm also works for standard low-order FEM. The method is described and accompanied with theoretical analysis and numerical examples. Instructions on how to reproduce the results are provided.     

Mesh and solver co‐adaptation in finite element methods for anisotropic problems

Mesh generation and algebraic solver are two important aspects of the finite element methodology. In this article, we are concerned with the joint adaptation of the anisotropic triangular mesh and the iterative algebraic solver. Using generic numerical examples pertaining to the accurate and efficient finite element solution of some anisotropic problems, we hereby demonstrate that the processes of geometric mesh adaptation and the algebraic solver construction should be adapted simultaneously. We also propose some techniques applicable to the co‐adaptation of both anisotropic meshes and linear solvers. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

An adaptive least-squares mixed finite element method for nonlinear parabolic problems

A least-squares mixed finite element method for nonlinear parabolic problems is investigated in terms of computational efficiency. An a posteriori error estimator, which is needed in an adaptive refinement algorithm, was composed with the least-squares functional, and a posteriori errors were effectively estimated.     

An adaptive finite element method with asymptotic saturation for eigenvalue problems

This paper discusses adaptive finite element methods for the solution of elliptic eigenvalue problems associated with partial differential operators. An adaptive method based on nodal-patch refinement leads to an asymptotic error reduction property for the computed sequence of simple eigenvalues and eigenfunctions. This justifies the use of the proven saturation property for a class of reliable and efficient hierarchical a posteriori error estimators. Numerical experiments confirm that the saturation property is present even for very coarse meshes for many examples; in other cases the smallness assumption on the initial mesh may be severe.     

Convergence of goal‐oriented adaptive finite element methods for nonsymmetric problems

In this article, we develop convergence theory for a class of goal‐oriented adaptive finite element algorithms for second‐order nonsymmetric linear elliptic equations. In particular, we establish contraction results for a method of this type for Dirichlet problems involving the elliptic operator with A Lipschitz, symmetric positive definite, with b divergence‐free, and with . We first describe the problem class and review some standard facts concerning conforming finite element discretization and error‐estimate‐driven adaptive finite element methods (AFEM). We then describe a goal‐oriented variation of standard AFEM. Following the recent work of Mommer and Stevenson for symmetric problems, we establish contraction and convergence of the goal‐oriented method in the sense of the goal function. Our analysis approach is signficantly different from that of Mommer and Stevenson, combining the recent contraction frameworks developed by Cascon, Kreuzer, Nochetto, and Siebert; by Nochetto, Siebert, and Veeser; and by Holst, Tsogtgerel, and Zhu. We include numerical results, demonstrating performance of our method with standard goal‐oriented strategies on a convection problem. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 479–509, 2016     

Feedback and adaptive finite element solution of one-dimensional boundary value problems

Summary This paper examines the concepts of feedback and adaptivity for the Finite Element Method. The model problem concernsC ⁰ elements of arbitrary, fixed degree for a one-dimensional two-point boundary value problem. Three different feedback methods are introduced and a detailed analysis of their adaptivity is given.Dedicated to F.L. Bauer on the occasion of his 60th birthdayThis research was partially supported by the Office of Naval Research under grant number N00014-77-C-0623     

Phase field modeling of complex crack patterns in multi-physics problems

Recently developed continuum phase field models for brittle fracture show excellent modeling capability in situations with complex crack topologies including branching in the small and large strain applications. This work presents a generalization towards fully coupled multi-physics problems at large strains. A modular concept is outlined for the linking of the diffusive crack modeling with complex multi field material response, where the focus is put on the model problem of finite thermo-elasticity. This concerns a generalization of crack driving forces from the energetic definitions towards stress-based criteria, the constitutive modeling of degradation of non-mechanical fluxes on generated crack faces. Particular assumptions are made on the generation of convective heat exchanges approximating surface load integrals of the sharp crack approach by distinct volume integrals. The coupling effect is also shown in generation of cracks due to thermally induced stress states. We finally demonstrate the performance of the phase field formulation of fracture at large strains by means of representative numerical examples. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Parallel multigrid method for finite element simulations of complex flow problems on locally refined meshes

We present a parallel multigrid solver on locally refined meshes for solving very complex three‐dimensional flow problems. Besides describing the parallel implementation in detail, we prove the smoothing property of the suggested iteration for a simple model problem. For demonstration of the efficiency and feasibility of the solver, we show a chemically reactive flow simulation for a Methane burner using detailed chemical reaction modeling. Further, we give the results of an ocean flow simulation. All described methods are implemented in the finite element toolbox Gascoigne. Copyright © 2010 John Wiley & Sons, Ltd.     

Convergence and quasi-optimality of adaptive nonconforming finite element methods for some nonsymmetric and indefinite problems

Recently an adaptive nonconforming finite element method (ANFEM) has been developed by Carstensen and Hoppe (in Numer Math 103:251–266, 2006). In this paper, we extend the result to some nonsymmetric and indefinite problems. The main tools in our analysis are a posteriori error estimators and a quasi-orthogonality property. In this case, we need to overcome two main difficulties: one stems from the nonconformity of the finite element space, the other is how to handle the effect of a nonsymmetric and indefinite bilinear form. An appropriate adaptive nonconforming finite element method featuring a marking strategy based on the comparison of the a posteriori error estimator and a volume term is proposed for the lowest order Crouzeix–Raviart element. It is shown that the ANFEM is a contraction for the sum of the energy error and a scaled volume term between two consecutive adaptive loops. Moreover, quasi-optimality in the sense of quasi-optimal algorithmic complexity can be shown for the ANFEM. The results of numerical experiments confirm the theoretical findings.     

Multi-mesh adaptive finite element algorithms for constrained optimal control problems governed by semi-linear parabolic equations

In this paper, we derive a posteriori error estimators for the constrained optimal control problems governed by semi-linear parabolic equations under some assumptions. Then we use them to construct reliable and efficient multi-mesh adaptive finite element algorithms for the optimal control problems. Some numerical experiments are presented to illustrate the theoretical results.     

An adaptive finite element algorithm with reliable and efficient error control for linear parabolic problems

An efficient and reliable a posteriori error estimate is derived for linear parabolic equations which does not depend on any regularity assumption on the underlying elliptic operator. An adaptive algorithm with variable time-step sizes and space meshes is proposed and studied which, at each time step, delays the mesh coarsening until the final iteration of the adaptive procedure, allowing only mesh and time-step size refinements before. It is proved that at each time step the adaptive algorithm is able to reduce the error indicators (and thus the error) below any given tolerance within a finite number of iteration steps. The key ingredient in the analysis is a new coarsening strategy. Numerical results are presented to show the competitive behavior of the proposed adaptive algorithm.

     

On finite element methods for plasticity problems

Summary We prove an error estimate for an incremental finite element method for obtaining approximations to the stresses in an elastic-perfectly plastic body. We also comment on the limit load problem.     

Error estimators and the adaptive finite element method on large strain deformation problems

We generalize the well‐known residual‐based error estimator in linear elasticity to the case of non‐linear deformation problems based on large strain and demonstrate its use in adaptive mesh control. Copyright © 2009 John Wiley & Sons, Ltd.     

Adaptive finite element procedures for elastoplastic problems at finite strains

A major difficulty in the context of adaptive analysis of geometrically nonlinear problems is to provide a robust remeshing procedure that accounts both for the error caused by the spatial discretization and for the error due to the time discretization. For stability problems, such as strain localization and necking, it is essential to provide a step–size control in order to get a robust algorithm for the solution of the boundary value problem. For this purpose we developed an easy to implement step–size control algorithm. In addition we will consider possible a posteriori error indicators for the spatial error distribution of elastoplastic problems at finite strains. This indicator is adopted for a density–function–based adaptive remeshing procedure. Both error indicators are combined for the adaptive analysis in time and space. The performance of the proposed method is documented by means of representative numerical examples.     

Parallel interior-point solver for structured quadratic programs: Application to financial planning problems

Many practical large-scale optimization problems are not only sparse, but also display some form of block-structure such as primal or dual block angular structure. Often these structures are nested: each block of the coarse top level structure is block-structured itself. Problems with these characteristics appear frequently in stochastic programming but also in other areas such as telecommunication network modelling. We present a linear algebra library tailored for problems with such structure that is used inside an interior point solver for convex quadratic programming problems. Due to its object-oriented design it can be used to exploit virtually any nested block structure arising in practical problems, eliminating the need for highly specialised linear algebra modules needing to be written for every type of problem separately. Through a careful implementation we achieve almost automatic parallelisation of the linear algebra. The efficiency of the approach is illustrated on several problems arising in the financial planning, namely in the asset and liability management. The problems are modelled as multistage decision processes and by nature lead to nested block-structured problems. By taking the variance of the random variables into account the problems become non-separable quadratic programs. A reformulation of the problem is proposed which reduces density of matrices involved and by these means significantly simplifies its solution by an interior point method. The object-oriented parallel solver achieves high efficiency by careful exploitation of the block sparsity of these problems. As a result a problem with over 50 million decision variables is solved in just over 2 hours on a parallel computer with 16 processors. The approach is by nature scalable and the parallel implementation achieves nearly perfect speed-ups on a range of problems. Supported by the Engineering and Physical Sciences Research Council of UK, EPSRC grant GR/R99683/01     

An adaptive finite element method for singular parabolic equations

Summary. In this paper we consider additive Schwarz-type iteration methods for saddle point problems as smoothers in a multigrid method. Each iteration step of the additive Schwarz method requires the solutions of several small local saddle point problems. This method can be viewed as an additive version of a (multiplicative) Vanka-type iteration, well-known as a smoother for multigrid methods in computational fluid dynamics. It is shown that, under suitable conditions, the iteration can be interpreted as a symmetric inexact Uzawa method. In the case of symmetric saddle point problems the smoothing property, an important part in a multigrid convergence proof, is analyzed for symmetric inexact Uzawa methods including the special case of the additive Schwarz-type iterations. As an example the theory is applied to the Crouzeix-Raviart mixed finite element for the Stokes equations and some numerical experiments are presented. Mathematics Subject Classification (1991):65N22, 65F10, 65N30Supported by the Austrian Science Foundation (FWF) under the grant SFB F013}\and Walter Zulehner