On the Convergence of a Multigrid Method for Linear Reaction-Diffusion Problems

In this note we consider discrete linear reaction-diffusion problems. For the discretization a standard conforming finite element method is used. For the approximate solution of the resulting discrete problem a multigrid method with a damped Jacobi or symmetric Gauss-Seidel smoother is applied. We analyze the convergence of the multigrid V- and W-cycle in the framework of the approximation- and smoothing property. The multigrid method is shown to be robust in the sense that the contraction number can be bounded by a constant smaller than one which does not depend on the mesh size or on the diffusion-reaction ratio. Received June 15, 2000     

An Algebraic Multigrid Method for Higher-order Finite Element Discretizations

In this paper, we will design and analyze a class of new algebraic multigrid methods for algebraic systems arising from the discretization of second order elliptic boundary value problems by high-order finite element methods. For a given sparse stiffness matrix from a quadratic or cubic Lagrangian finite element discretization, an algebraic approach is carefully designed to recover the stiffness matrix associated with the linear finite element disretization on the same underlying (but nevertheless unknown to the user) finite element grid. With any given classical algebraic multigrid solver for linear finite element stiffness matrix, a corresponding algebraic multigrid method can then be designed for the quadratic or higher order finite element stiffness matrix by combining with a standard smoother for the original system. This method is designed under the assumption that the sparse matrix to be solved is associated with a specific higher order, quadratic for example, finite element discretization on a finite element grid but the geometric data for the underlying grid is unknown. The resulting new algebraic multigrid method is shown, by numerical experiments, to be much more efficient than the classical algebraic multigrid method which is directly applied to the high-order finite element matrix. Some theoretical analysis is also provided for the convergence of the new method.     

A priori and a posteriori error analysis of an augmented mixed finite element method for incompressible fluid flows

In this paper we extend recent results on the a priori and a posteriori error analysis of an augmented mixed finite element method for the linear elasticity problem, to the case of incompressible fluid flows with symmetric stress tensor. Similarly as before, the present approach is based on the introduction of the Galerkin least-squares type terms arising from the constitutive and equilibrium equations, and from the relations defining the pressure in terms of the stress tensor and the rotation in terms of the displacement, all of them multiplied by stabilization parameters. We show that these parameters can be suitably chosen so that the resulting augmented variational formulation is defined by a strongly coercive bilinear form, whence the associated Galerkin scheme becomes well-posed for any choice of finite element subspaces. Next, we present a reliable and efficient residual-based a posteriori error estimator for the augmented mixed finite element scheme. Finally, several numerical results confirming the theoretical properties of this estimator, and illustrating the capability of the corresponding adaptive algorithm to localize the singularities and the large stress regions of the solution, are reported.     

Stable Evaluation of Jacobian Matrices on Highly Refined Finite Element Meshes

In the adaptive finite element method, the solution of a p.d.e. is approximated by finer and finer meshes, which are controlled from error estimators. So, starting from a given coarse mesh, some elements are subdivided a couple of times. We investigate the question of avoiding instabilities which limit this process from the fact that nodal coordinates of one element coincide in more and more leading digits. To overcome this problem we demonstrate a simple mechanism for red subdivision of triangles (and hanging nodes) and a more sophisticated technique for general quadrilaterals.     

Composite Finite Elements for Elliptic Boundary Value Problems with Discontinuous Coefficients

In this paper, we will introduce composite finite elements for solving elliptic boundary value problems with discontinuous coefficients. The focus is on problems where the geometry of the interfaces between the smooth regions of the coefficients is very complicated. On the other hand, efficient numerical methods such as, e.g., multigrid methods, wavelets, extrapolation, are based on a multi-scale discretization of the problem. In standard finite element methods, the grids have to resolve the structure of the discontinuous coefficients. Thus, straightforward coarse scale discretizations of problems with complicated coefficient jumps are not obvious. In this paper, we define composite finite elements for problems with discontinuous coefficients. These finite elements allow the coarsening of finite element spaces independently of the structure of the discontinuous coefficients. Thus, the multigrid method can be applied to solve the linear system on the fine scale. We focus on the construction of the composite finite elements and the efficient, hierarchical realization of the intergrid transfer operators. Finally, we present some numerical results for the multigrid method based on the composite finite elements (CFE–MG).     

Approximating Local Averages of Fluid Velocities: The Stokes Problem

As a first step to developing mathematical support for finite element approximation to the large eddies in fluid motion we consider herein the Stokes problem. We show that the local average of the usual approximate flow field u ^h over radius δ provides a very accurate approximation to the flow structures of O(δ) or greater. The extra accuracy appears for quadratic or higher velocity elements and degrades to the usual finite element accuracy as the averaging radius δ→h (the local meshwidth). We give both a priori and a posteriori error estimates incorporating this effect. Received December 3, 1999; revised October 16, 2000     

A Domain Decomposition Preconditioner for p-FEM Discretizations of Two-dimensional Elliptic Problems

In this paper, a uniformly elliptic second order boundary value problem in 2-D discretized by the p-version of the finite element method is considered. An inexact Dirichlet-Dirichlet domain decomposition pre-conditioner for the system of linear algebraic equations is investigated. Two solvers for the problem in the sub-domains, a pre-conditioner for the Schur-complement and an extension operator operating from the edges of the elements into the interior are proposed as ingredients for the inexact DD-pre-conditioner. In the main part of the paper, several numerical experiments on a parallel computer are given.     

Arbitrary-level hanging nodes and automatic adaptivity in the hp-FEM

In this paper we present a new automatic adaptivity algorithm for the hp-FEM which is based on arbitrary-level hanging nodes and local element projections. The method is very simple to implement compared to other existing hp-adaptive strategies, while its performance is comparable or superior. This is demonstrated on several numerical examples which include the L-shape domain problem, a problem with internal layer, and the Girkmann problem of linear elasticity. With appropriate simplifications, the proposed technique can be applied to standard lower-order and spectral finite element methods.     

The Streamline–Diffusion Method for Conforming and Nonconforming Finite Elements of Lowest Order Applied to Convection–Diffusion Problems

We consider the streamline-diffusion finite element method with finite elements of lowest order for solving convection-diffusion problems. Our investigations cover both conforming and nonconforming finite element approximations on triangular and quadrilateral meshes. Although the considered finite elements are of the same interpolation order their stability and approximation properties are quite different. We give a detailed overview on the stability and the convergence properties in the L ²- and in the streamline–diffusion norm. Numerical experiments show that often the theoretical predictions on the convergence properties are sharp. Received December 7, 1999; revised October 5, 2000     

Algebraic Multigrid Based on Computational Molecules, 1: Scalar Elliptic Problems

We consider the problem of splitting a symmetric positive definite (SPD) stiffness matrix A arising from finite element discretization into a sum of edge matrices thereby assuming that A is given as the sum of symmetric positive semidefinite (SPSD) element matrices. We give necessary and sufficient conditions for the existence of an exact splitting into SPSD edge matrices and address the problem of best positive (nonnegative) approximation. Based on this disassembling process we present a new concept of ``strong' and ``weak' connections (edges), which provides a basis for selecting the coarse-grid nodes in algebraic multigrid methods. Furthermore, we examine the utilization of computational molecules (small collections of edge matrices) for deriving interpolation rules. The reproduction of edge matrices on coarse levels offers the opportunity to combine classical coarsening algorithms with effective (energy minimizing) interpolation principles yielding a flexible and robust new variant of AMG.     

Semi-Unstructured Grids

We present a novel automatic grid generator for the finite element discretization of partial differential equations in 3D. The grids constructed by this grid generator are composed of a pure tensor product grid in the interior of the domain and an unstructured grid which is only contained in boundary cells. The unstructured component consists of tetrahedra, each of which satisfies a maximal interior angle condition. By suitable constructing the boundary cells, the number of types of boundary subcells is reduced to 12 types. Since this grid generator constructs large structured grids in the interior and small unstructured grids near the boundary, the resulting semi-unstructured grids have similar properties as structured tensor product grids. Some appealing properties of this method are computational efficiency and natural construction of coarse grids for multilevel algorithms. Numerical results and an analysis of the discretization error are presented. Received July 17, 2000; revised October 27, 2000     

Weighted Error Estimates for Finite Element Solutions of Variational Inequalities

In this note the studies begun in Blum and Suttmeier (1999) on adaptive finite element discretisations for nonlinear problems described by variational inequalities are continued. Similar to the concept proposed, e.g., in Becker and Rannacher (1996) for variational equalities, weighted a posteriori estimates for controlling arbitrary functionals of the discretisation error are constructed by using a duality argument. Numerical results for the obstacle problem demonstrate the derived error bounds to be reliable and, used for an adaptive grid refinement strategy, to produce economical meshes. Received September 6, 1999; revised February 8, 2000     

An Iterative Substructuring Method for div-stable Finite Element Approximations of the Oseen Problem

We apply an iterative substructuring algorithm with transmission conditions of Robin–Robin type to the discretized Oseen problem appearing as a linearized variant of the incompressible Navier–Stokes equations. Here we consider finite element approximations using velocity/pressure pairs which satisfy the Babuška–Brezzi stability condition. After proving well-posedness and strong convergence of the method, we derive an a-posteriori error estimate which controls convergence of the discrete subdomain solutions to the global discrete solution by measuring the jumps of the velocities at the interface. Additionally we obtain information how to design a parameter of the Robin interface condition which essentially influences the convergence speed. Numerical experiments confirm the theoretical results and the applicability of the method. Received February 18, 2000; revised February 21, 2001     

Anisotropic interpolation error estimates for isoparametric quadrilateral finite elements

Anisotropic local interpolation error estimates are derived for quadrilateral and hexahedral Lagrangian finite elements with straight edges. These elements are allowed to have diameters with different asymptotic behaviour in different space directions. The case of affine elements (parallel-epipeds) with arbitrarily high degree of the shape functions is considered first. Then, a careful examination of the multi-linear map leads to estimates for certain classes of more general, isoparametric elements. As an application, the Galerkin finite element method for a reaction diffusion problem in a polygonal domain is considered. The boundary layers are resolved using anisotropic trapezoidal elements.     

A Cascadic Multigrid Algorithm for Semilinear Indefinite Elliptic Problems

We propose a cascadic multigrid algorithm for a semilinear indefinite elliptic problem. We use a standard finite element discretization with piecewise linear finite elements. The arising nonlinear equations are solved by a cascadic organization of Newton's method with frozen derivative on a sequence of nested grids. This gives a simple version of a multigrid method without projections on coarser grids. The cascadic multigrid algorithm starts on a comparatively coarse grid where the number of unknowns is small enough to obtain an approximate solution within sufficiently high precision without substantial computational effort. On each finer grid we perform exactly one Newton step taking the approximate solution from the coarsest grid as initial guess. The linear Newton systems are solved iteratively by a Jacobi-type iteration with special parameters using the approximate solution from the previous grid as initial guess. We prove that for a sufficiently fine initial grid and for a sufficiently good start approximation the algorithm yields an approximate solution within the discretization error on the finest grid and that the method has multigrid complexity with logarithmic multiplier. Received February 1999, revised July 13, 1999     

BDDC and FETI-DP preconditioners for spectral element discretizations

Two of the most recent and important nonoverlapping domain decomposition methods, the BDDC method (Balancing Domain Decomposition by Constraints) and the FETI-DP method (Dual-Primal Finite Element Tearing and Interconnecting) are here extended to spectral element discretizations of second-order elliptic problems. In spite of the more severe ill-conditioning of the spectral element discrete systems, compared with low-order finite elements and finite differences, these methods retain their good properties of scalability, quasi-optimality and independence on the discontinuities of the elliptic operator coefficients across subdomain interfaces.     

A Refined Finite Element Convergence Theory for Highly Indefinite Helmholtz Problems

It is well known that standard h-version finite element discretisations using lowest order elements for Helmholtz' equation suffer from the following stability condition: ``The mesh width h of the finite element mesh has to satisfy k <sup>2</sup> h≲1', where k denotes the wave number. This condition rules out the reliable numerical solution of Helmholtz equation in three dimensions for large wave numbers k≳50. In our paper, we will present a refined finite element theory for highly indefinite Helmholtz problems where the stability of the discretisation can be checked through an ``almost invariance' condition. As an application, we will consider a one-dimensional finite element space for the Helmholtz equation and apply our theory to prove stability under the weakened condition hk≲1 and optimal convergence estimates. Dedicated to Prof. Dr. Ivo Babuška on the occasion of his 80th birthday.     

Nitsche Type Mortaring for Singularly Perturbed Reaction-diffusion Problems

The paper is concerned with the Nitsche mortaring in the framework of domain decomposition where non-matching meshes and weak continuity of the finite element approximation at the interface are admitted. The approach is applied to singularly perturbed reaction-diffusion problems in 2D. Non-matching meshes of triangles being anisotropic in the boundary layers are applied. Some properties as well as error estimates of the Nitsche mortar finite element schemes are proved. In particular, using a suitable degree of anisotropy of triangles in the boundary layers of a rectangle, convergence rates as known for the conforming finite element method are derived. Numerical examples illustrate the approach and the results.     

Domain decomposition methods for the neutron diffusion problem

The neutronic simulation of a nuclear reactor core is performed using the neutron transport equation, and leads to an eigenvalue problem in the steady-state case. Among the deterministic resolution methods, simplified transport (S_PN$S{P}_N$) or diffusion approximations are often used. The MINOS solver developed at CEA Saclay uses a mixed dual finite element method for the resolution of these problems, and has shown his efficiency. In order to take into account the heterogeneities of the geometry, a very fine mesh is generally required, and leads to expensive calculations for industrial applications. In order to take advantage of parallel computers, and to reduce the computing time and the local memory requirement, we propose here two domain decomposition methods based on the MINOS solver. The first approach is a component mode synthesis method on overlapping subdomains: several eigenmodes solutions of a local problem on each subdomain are taken as basis functions used for the resolution of the global problem on the whole domain. The second approach is an iterative method based on a non-overlapping domain decomposition with Robin interface conditions. At each iteration, we solve the problem on each subdomain with the interface conditions given by the solutions on the adjacent subdomains estimated at the previous iteration. Numerical results on parallel computers are presented for the diffusion model on realistic 2D and 3D cores.     

Superconvergence of a Mixed Covolume Method for Elliptic Problems

We consider a mixed covolume method for a system of first order partial differential equations resulting from the mixed formulation of a general self-adjoint elliptic problem with a variable full diffusion tensor. The system can be used to model the transport of a contaminant carried by a flow. We use the lowest order Raviart-Thomas mixed finite element space. We show the first order convergence in L ² norm and the superconvergence in certain discrete norms both for the pressure and velocity. Finally some numerical examples illustrating the error behavior of the scheme are provided. Supported by the National Natural Science Foundation of China under grant No. 10071044 and the Research Fund of Doctoral Program of High Education by State Education Ministry of China.