Quasi-norm error bounds for the finite element approximation of a non-Newtonian flow

Summary. We consider the finite element approximation of a non-Newtonian flow, where the viscosity obeys a general law including the Carreau or power law. For sufficiently regular solutions we prove energy type error bounds for the velocity and pressure. These bounds improve on existing results in the literature. A key step in the analysis is to prove abstract error bounds initially in a quasi-norm, which naturally arises in degenerate problems of this type. Received May 25, 1993 / Revised version received January 11, 1994     

The artificial boundary conditions for incompressible materials on an unbounded domain

Summary. In this paper we consider the numerical simulations of the incompressible materials on an unbounded domain in . A series of artificial boundary conditions at a circular artificial boundary for solving incompressible materials on an unbounded domain is given. Then the original problem is reduced to a problem on a bounded domain, which be solved numerically by a mixed finite element method. The numerical example shows that our artificial boundary conditions are very effective. ReceivedJune 7, 1995 / Revised version received August 19, 1996     

On the quadratic finite element approximation to the obstacle problem

Summary. In this paper, we obtain the error bound for any , for the piecewise quadratic finite element approximation to the obstacle problem, without the hypothesis that the free boundary has finite length (see [3]). Received October 31, 2000 / Revised version received July 23, 2001 / Published online October 17, 2001 The project was supported by the National Natural Science Foundation of China     

Mixed finite element approximation of a degenerate elliptic problem

Summary. We present a mixed finite element approximation of an elliptic problem with degenerate coefficients, arising in the study of the electromagnetic field in a resonant structure with cylindrical symmetry. Optimal error bounds are derived. Received May 4, 1994 / Revised version received September 27, 1994     

A term by term stabilization algorithm for finite element solution of incompressible flow problems

This paper introduces a stabilization technique for Finite Element numerical solution of 2D and 3D incompressible flow problems. It may be applied to stabilize the discretization of the pressure gradient, and also of any individual operator term such as the convection, curl or divergence operators, with specific levels of numerical diffusion for each one of them. Its computational complexity is reduced with respect to usual (residual-based) stabilization techniques. We consider piecewise affine Finite Elements, for which we obtain optimal error bounds for steady Navier-Stokes and also for generalized Stokes equations (including convection). We include some numerical experiment in well known 2D test cases, that show its good performances. Received March 15, 1996 / Revised version received January 17, 1997     

On a decoupled algorithm for solving a finite element problem for the approximation of viscoelastic fluid flow

Summary. The aim of this work is to study a decoupled algorithm of a fixed point for solving a finite element (FE) problem for the approximation of viscoelastic fluid flow obeying an Oldroyd B differential model. The interest for this algorithm lies in its applications to numerical simulation and in the cost of computing. Furthermore it is easy to bring this algorithm into play. The unknowns are the viscoelastic part of the extra stress tensor, the velocity and the pressure. We suppose that the solution is sufficiently smooth and small. The approximation of stress, velocity and pressure are resp. discontinuous, continuous, continuous FE. Upwinding needed for convection of , is made by discontinuous FE. The method consists to solve alternatively a transport equation for the stress, and a Stokes like problem for velocity and pressure. Previously, results of existence of the solution for the approximate problem and error bounds have been obtained using fixed point techniques with coupled algorithm. In this paper we show that the mapping of the decoupled fixed point algorithm is locally (in a neighbourhood of ) contracting and we obtain existence, unicity (locally) of the solution of the approximate problem and error bounds. Received July 29, 1994 / Revised version received March 13, 1995     

Least-squares mixed finite element methods for non-selfadjoint elliptic problems: I. Error estimates

Summary. A least-squares mixed finite element method for general second-order non-selfadjoint elliptic problems in two- and three-dimensional domains is formulated and analyzed. The finite element spaces for the primary solution approximation and the flux approximation consist of piecewise polynomials of degree and respectively. The method is mildly nonconforming on the boundary. The cases and are studied. It is proved that the method is not subject to the LBB-condition. Optimal - and -error estimates are derived for regular finite element partitions. Numerical experiments, confirming the theoretical rates of convergence, are presented. Received October 15, 1993 / Revised version received August 2, 1994     

A mixed formulation for 3D magnetostatic problems: theoretical analysis and face-edge finite element approximation

Summary. A mixed field-based variational formulation for the solution of threedimensional magnetostatic problems is presented and analyzed. This method is based upon the minimization of a functional related to the error in the constitutive magnetic relationship, while constraints represented by Maxwell's equations are imposed by means of Lagrange multipliers. In this way, both the magnetic field and the magnetic induction field can be approximated by using the most appropriate family of vector finite elements, and boundary conditions can be imposed in a natural way. Moreover, this method is more suitable than classical approaches for the approximation of problems featuring strong discontinuities of the magnetic permeability, as is usually the case. A finite element discretization involving face and edge elements is also proposed, performing stability analysis and giving error estimates. Received January 23, 1998 / Revised version received July 23, 1998 / Published online September 24, 1999     

The approximation theory for the p-version finite element method and application to non-linear elliptic PDEs

Approximation theoretic results are obtained for approximation using continuous piecewise polynomials of degree p on meshes of triangular and quadrilateral elements. Estimates for the rate of convergence in Sobolev spaces , are given. The results are applied to estimate the rate of convergence when the p-version finite element method is used to approximate the -Laplacian. It is shown that the rate of convergence of the p-version is always at least that of the h-version (measured in terms of number of degrees of freedom used). If the solution is very smooth then the p-version attains an exponential rate of convergence. If the solution has certain types of singularity, the rate of convergence of the p-version is twice that of the h-version. The analysis generalises the work of Babuska and others to the case . In addition, the approximation theoretic results find immediate application for some types of spectral and spectral element methods. Received August 2, 1995 / Revised version received January 26, 1998     

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     

The finite element approximation of a coupled reaction-diffusion problem with non-Lipschitz nonlinearities

Summary. A coupled semilinear elliptic problem modelling an irreversible, isothermal chemical reaction is introduced, and discretised using the usual piecewise linear Galerkin finite element approximation. An interesting feature of the problem is that a reaction order of less than one gives rise to a "dead core" region. Initially, one reactant is assumed to be acting as a catalyst and is kept constant. It is shown that error bounds previously obtained for a scheme involving numerical integration can be improved upon by considering a quadratic regularisation of the nonlinear term. This technique is then applied to the full coupled problem, and optimal and error bounds are proved in the absence of quadrature. For a scheme involving numerical integration, bounds similar to those obtained for the catalyst problem are shown to hold. Received May 25, 1993 / Revised version received July 5, 1994     

An asymptotically optimal finite element scheme for the arch problem

Received November 22, 1993     

Superconvergence analysis and error expansion for the Wilson nonconforming finite element

Summary. In this paper the Wilson nonconforming finite element is considered for solving a class of two-dimensional second-order elliptic boundary value problems. Superconvergence estimates and error expansions are obtained for both uniform and non-uniform rectangular meshes. A new lower bound of the error shows that the usual error estimates are optimal. Finally a discussion on the error behaviour in negative norms shows that there is generally no improvement in the order by going to weaker norms. Received July 5, 1993     

Error bounds for Romberg quadrature

Denote by the error of a Romberg quadrature rule applied to the function f. We determine approximately the constants in the bounds of the types and for all classical Romberg rules. By a comparison with the corresponding constants of the Gaussian rule we give the statement “The Gaussian quadrature rule is better than the Romberg method” a precise meaning. Received September 10, 1997 / Revised version received February 16, 1998     

Multiscale finite element for problems with highly oscillatory coefficients

Summary. In this paper, we study a multiscale finite element method for solving a class of elliptic problems with finite number of well separated scales. The method is designed to efficiently capture the large scale behavior of the solution without resolving all small scale features. This is accomplished by constructing the multiscale finite element base functions that are adaptive to the local property of the differential operator. The construction of the base functions is fully decoupled from element to element; thus the method is perfectly parallel and is naturally adapted to massively parallel computers. We present the convergence analysis of the method along with the results of our numerical experiments. Some generalizations of the multiscale finite element method are also discussed. Received April 17, 1998 / Revised version received March 25, 2000 / Published online June 7, 2001     

A minimal stabilisation procedure for mixed finite element methods

Summary. Stabilisation methods are often used to circumvent the difficulties associated with the stability of mixed finite element methods. Stabilisation however also means an excessive amount of dissipation or the loss of nice conservation properties. It would thus be desirable to reduce these disadvantages to a minimum. We present a general framework, not restricted to mixed methods, that permits to introduce a minimal stabilising term and hence a minimal perturbation with respect to the original problem. To do so, we rely on the fact that some part of the problem is stable and should not be modified. Sections 2 and 3 present the method in an abstract framework. Section 4 and 5 present two classes of stabilisations for the inf-sup condition in mixed problems. We present many examples, most arising from the discretisation of flow problems. Section 6 presents examples in which the stabilising terms is introduced to cure coercivity problems. Received August 9, 1999 / Revised version received May 19, 2000 / Published online March 20, 2001     

Mixed finite element for the linear plate problem: the Hermann-Miyoshi model revisited

Summary. In this paper we study the relationship between the Hermann-Miyoshi and the Ciarlet-Raviart formulations of the first biharmonic problem. This study will be based on a decomposition principle which will leads us to a new convergence analysis explaining some discrepancies between numerical results obtained with the first formulation on certain meshes and some theoretical convergence results. Received May 24, 1994 / Revised version received August 11, 1995     

Adaptive finite element methods for elliptic equations with non-smooth coefficients

Summary. We consider a second-order elliptic equation with discontinuous or anisotropic coefficients in a bounded two- or three dimensional domain, and its finite-element discretization. The aim of this paper is to prove some a priori and a posteriori error estimates in an appropriate norm, which are independent of the variation of the coefficients. Received February 5, 1999 / Published online March 16, 2000     

Convergence of a finite element method for non-parametric mean curvature flow

Summary. Convergence for the spatial discretization by linear finite elements of the non-parametric mean curvature flow is proved under natural regularity assumptions on the continuous solution. Asymptotic convergence is also obtained for the time derivative which is proportional to mean curvature. An existence result for the continuous problem in adequate spaces is included. Received September 30, 1993