Iterative schemes for mixed finite element methods with applications to elasticity and compressible flow problems

Summary Iterative schemes for mixed finite element methods are proposed and analyzed in two abstract formulations. The first one has applications to elliptic equations and incompressible fluid flow problems, while the second has applications to linear elasticity and compressible Stokes problems. These schemes are constructed through iteratively penalizing the mixed finite element scheme, of which iterated penalty method and augmented Lagrangian method are special cases. Convergence theorems are demonstrated in abstract formulations in Hilbert spaces, and applications to individual physical problems are considered as examples. Theoretical analysis and computational experiments both show that the proposed schemes have very fast convergence; a few iterations are normally enough to reduce the iterative error to a prescribed precision. Numerical examples with continuous and discontinuous coefficients are presented.     

An augmented mixed finite element method for 3D linear elasticity problems

In this paper we introduce and analyze a new augmented mixed finite element method for linear elasticity problems in 3D. Our approach is an extension of a technique developed recently for plane elasticity, which is based on the introduction of consistent terms of Galerkin least-squares type. We consider non-homogeneous and homogeneous Dirichlet boundary conditions and prove that the resulting augmented variational formulations lead to strongly coercive bilinear forms. In this way, the associated Galerkin schemes become well posed for arbitrary choices of the corresponding finite element subspaces. In particular, Raviart-Thomas spaces of order 0 for the stress tensor, continuous piecewise linear elements for the displacement, and piecewise constants for the rotation can be utilized. Moreover, we show that in this case the number of unknowns behaves approximately as 9.5 times the number of elements (tetrahedrons) of the triangulation, which is cheaper, by a factor of 3, than the classical PEERS in 3D. Several numerical results illustrating the good performance of the augmented schemes are provided.     

A posteriori error estimates for a finite element approximation of transmission problems with sign changing coefficients

We perform the a posteriori error analysis of residual type of transmission problem with sign changing coefficients. According to Bonnet-BenDhia et al. (2010) [9], if the contrast is large enough, the continuous problem can be transformed into a coercive one. We further show that a similar property holds for the discrete problem for any regular meshes, extending the framework from Bonnet-BenDhia et al. [9]. The reliability and efficiency of the proposed estimator are confirmed by some numerical tests.     

An approximation of incompressible miscible displacement in porous media by mixed finite element method and characteristics-mixed finite element method

An approximation scheme is defined for incompressible miscible displacement in porous media. This scheme is constructed by using two methods. Standard mixed finite element is used for the Darcy velocity equation. A characteristics-mixed finite element method is presented for the concentration equation. Characteristic approximation is applied to handle the convection part of the concentration equation, and a lowest-order mixed finite element spatial approximation is adopted to deal with the diffusion part. Thus, the scalar unknown concentration and the diffusive flux can be approximated simultaneously. In order to derive the optimal L²$L^2$-norm error estimates, a post-processing step is included in the approximation to the scalar unknown concentration. This scheme conserves mass globally; in fact, on the discrete level, fluid is transported along the approximate characteristics. Numerical experiments are presented finally to validate the theoretical analysis.     

Some equilibrium finite element methods for two-dimensional elasticity problems

Summary We consider some equilibrium finite element methods for two-dimensional elasticity problems. The stresses and the displacements are approximated by using piecewise linear functions. We establishL ₂-estimates of orderO(h ²) for both stresses and displacements.     

A posteriori error estimators for nonconforming finite element methods of the linear elasticity problem

In this work we derive and analyze a posteriori error estimators for low-order nonconforming finite element methods of the linear elasticity problem on both triangular and quadrilateral meshes, with hanging nodes allowed for local mesh refinement. First, it is shown that equilibrated Neumann data on interelement boundaries are simply given by the local weak residuals of the numerical solution. The first error estimator is then obtained by applying the equilibrated residual method with this set of Neumann data. From this implicit estimator we also derive two explicit error estimators, one of which is similar to the one proposed by Dörfler and Ainsworth (2005) [24] for the Stokes problem. It is established that all these error estimators are reliable and efficient in a robust way with respect to the Lamé constants. The main advantage of our error estimators is that they yield guaranteed, i.e., constant-free upper bounds for the energy-like error (up to higher order terms due to data oscillation) when a good estimate for the inf-sup constant is available, which is confirmed by some numerical results.     

A Zienkiewicz-type finite element applied to fourth-order problems

This paper deals with convergence analysis and applications of a Zienkiewicz-type (Z-type) triangular element, applied to fourth-order partial differential equations. For the biharmonic problem we prove the order of convergence by comparison to a suitable modified Hermite triangular finite element. This method is more natural and it could be applied to the corresponding fourth-order eigenvalue problem. We also propose a simple postprocessing method which improves the order of convergence of finite element eigenpairs. Thus, an a posteriori analysis is presented by means of different triangular elements. Some computational aspects are discussed and numerical examples are given.     

Mixed hp finite element methods for problems in elasticity and Stokes flow

Summary. We consider the mixed formulation for the elasticity problem and the limiting Stokes problem in , . We derive a set of sufficient conditions under which families of mixed finite element spaces are simultaneously stable with respect to the mesh size and, subject to a maximum loss of , with respect to the polynomial degree . We obtain asymptotic rates of convergence that are optimal up to in the displacement/velocity and up to in the "pressure", with arbitrary (both rates being optimal with respect to ). Several choices of elements are discussed with reference to properties desirable in the context of the -version. Received March 4, 1994 / Revised version received February 12, 1995     

Streamline diffusion finite element method for quasilinear elliptic problems

Summary We extend the analysis of the streamline diffusion finite element method to quasilinear elliptic problems of second order. An existence theorem and error estimates are given in the case of branches of nonsingular solutions following a recent abstract approach in [12, 13, 26].     

Finite element approximation of multi-scale elliptic problems using patches of elements

In this paper we present a method for the numerical solution of elliptic problems with multi-scale data using multiple levels of not necessarily nested grids. The method consists in calculating successive corrections to the solution in patches whose discretizations are not necessarily conforming. This paper provides proofs of the results published earlier (see C. R. Acad. Sci. Paris, Ser. I 337 (2003) 679–684), gives a generalization of the latter to more than two domains and contains extensive numerical illustrations. New results including the spectral analysis of the iteration operator and a numerical method to evaluate the constant of the strengthened Cauchy-Buniakowski-Schwarz inequality are presented. Supported by CTI Project 6437.1 IWS-IW.     

Estimation of the effect of numerical integration in finite element eigenvalue approximation

Summary Finite element approximations of the eigenpairs of differential operators are computed as eigenpairs of matrices whose elements involve integrals which must be evaluated by numerical integration. The effect of this numerical integration on the eigenvalue and eigenfunction error is estimated. Specifically, for 2nd order selfadjoint eigenvalue problems we show that finite element approximations with quadrature satisfy the well-known estimates for approximations without quadrature, provided the quadrature rules have appropriate degrees of precision.The work of this author was partially supported by the National Science Foundation under Grant DMS-84-10324     

The combined effect of numerical integration and approximation of the boundary in the finite element method for eigenvalue problems

Summary. The paper deals with the finite element analysis of second order elliptic eigenvalue problems when the approximate domains are not subdomains of the original domain and when at the same time numerical integration is used for computing the involved bilinear forms. The considerations are restricted to piecewise linear approximations. The optimum rate of convergence for approximate eigenvalues is obtained provided that a quadrature formula of first degree of precision is used. In the case of a simple exact eigenvalue the optimum rate of convergence for approximate eigenfunctions in the -norm is proved while in the -norm an almost optimum rate of convergence (i.e. near to is achieved. In both cases a quadrature formula of first degree of precision is used. Quadrature formulas with degree of precision equal to zero are also analyzed and in the case when the exact eigenfunctions belong only to the convergence without the rate of convergence is proved. In the case of a multiple exact eigenvalue the approximate eigenfunctions are compard (in contrast to standard considerations) with linear combinations of exact eigenfunctions with coefficients not depending on the mesh parameter . Received September 18, 1993 / Revised version received September 26, 1994     

Continuous-discontinuous finite element method for convection-diffusion problems with characteristic layers

We study convergence properties of a numerical method for convection-diffusion problems with characteristic layers on a layer-adapted mesh. The method couples standard Galerkin with an h-version of the nonsymmetric discontinuous Galerkin finite element method with bilinear elements. In an associated norm, we derive the error estimate as well as the supercloseness result that are uniform in the perturbation parameter. Applying a post-processing operator for the discontinuous Galerkin method, we construct a new numerical solution with enhanced convergence properties.     

Two-grid finite volume element method for linear and nonlinear elliptic problems

Two-grid finite volume element discretization techniques, based on two linear conforming finite element spaces on one coarse and one fine grid, are presented for the two-dimensional second-order non-selfadjoint and indefinite linear elliptic problems and the two-dimensional second-order nonlinear elliptic problems. With the proposed techniques, solving the non-selfadjoint and indefinite elliptic problem on the fine space is reduced into solving a symmetric and positive definite elliptic problem on the fine space and solving the non-selfadjoint and indefinite elliptic problem on a much smaller space; solving a nonlinear elliptic problem on the fine space is reduced into solving a linear problem on the fine space and solving the nonlinear elliptic problem on a much smaller space. Convergence estimates are derived to justify the efficiency of the proposed two-grid algorithms. A set of numerical examples are presented to confirm the estimates. The work is supported by the National Natural Science Foundation of China (Grant No: 10601045).     

Error estimates of the lumped mass finite element method for semilinear elliptic problems

We are concerned with the semilinear elliptic problems. We first investigate the L²-error estimate for the lumped mass finite element method. We then use the cascadic multigrid method to solve the corresponding discrete problem. On the basis of the finite element error estimates, we prove the optimality of the proposed multigrid method. We also report some numerical results to support the theory.     

A posteriori error estimates for mixed finite element approximations of elliptic problems

We derive residual based a posteriori error estimates of the flux in L ²-norm for a general class of mixed methods for elliptic problems. The estimate is applicable to standard mixed methods such as the Raviart–Thomas–Nedelec and Brezzi–Douglas–Marini elements, as well as stabilized methods such as the Galerkin-Least squares method. The element residual in the estimate employs an elementwise computable postprocessed approximation of the displacement which gives optimal order.     

Existence and approximation results for nonlinear mixed problems: Application to incompressible finite elasticity

Summary This paper considers the problems of minimizing Gateaux-differentiable functionals over subsets of real Banach spaces defined by a non-linear equality constraint. The existence of a Lagrange multiplier is proved, together with approximation results on the constrained subset, provided a nonlinear compatibility condition, generalizing the classical inf-sup condition, is satisfied. These ideas are applied to equilibrium problems in incompressible finite elasticity and lead to convergence results for these problems.     

Nonconforming mixed finite element approximation to the stationary Navier-Stokes equations on anisotropic meshes

The main aim of this paper is to study the error estimates of a rectangular nonconforming finite element for the stationary Navier-Stokes equations under anisotropic meshes. That is, the nonconforming rectangular element is taken as approximation space for the velocity and the piecewise constant element for the pressure. The convergence analysis is presented and the optimal error estimates both in a broken H¹-norm for the velocity and in an L²-norm for the pressure are derived on anisotropic meshes.     

Dual-mixed finite element approximation of Stokes and nonlinear Stokes problems using trace-free velocity gradients

In this work a finite element method for a dual-mixed approximation of Stokes and nonlinear Stokes problems is studied. The dual-mixed structure, which yields a twofold saddle point problem, arises in a formulation of this problem through the introduction of unknown variables with relevant physical meaning. The method approximates the velocity, its gradient, and the total stress tensor, but avoids the explicit computation of the pressure, which can be recovered through a simple postprocessing technique. This method improves an existing approach for these problems and uses Raviart-Thomas elements and discontinuous piecewise polynomials for approximating the unknowns. Existence, uniqueness, and error results for the method are given, and numerical experiments that exhibit the reduced computational cost of this approach are presented.     

Biquadratic finite volume element methods based on optimal stress points for parabolic problems

In this paper, the semi-discrete and full discrete biquadratic finite volume element schemes based on optimal stress points for a class of parabolic problems are presented. Optimal order error estimates in H¹ and L² norms are derived. In addition, the superconvergences of numerical gradients at optimal stress points are also discussed. A numerical experiment confirms some results of theoretical analysis.