A posteriori error estimates for elliptic problems with Dirac delta source terms

The aim of this paper is to introduce residual type a posteriori error estimators for a Poisson problem with a Dirac delta source term, in L ^p norm and W^1,p seminorm. The estimators are proved to yield global upper and local lower bounds for the corresponding norms of the error. They are used to guide adaptive procedures, which are experimentally shown to lead to optimal orders of convergence.     

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.     

A posteriori error estimates for mixed finite element solutions of convex optimal control problems

In this paper, we present an a posteriori error analysis for mixed finite element approximation of convex optimal control problems. We derive a posteriori error estimates for the coupled state and control approximations under some assumptions which hold in many applications. Such estimates can be used to construct reliable adaptive mixed finite elements for the control problems.     

Adaptive multilevel finite element approximations of semilinear elliptic boundary value problems

Summary. In this paper we consider two aspects of the problem of designing efficient numerical methods for the approximation of semilinear boundary value problems. First we consider the use of two and multilevel algorithms for approximating the discrete solution. Secondly we consider adaptive mesh refinement based on feedback information from coarse level approximations. The algorithms are based on an a posteriori error estimate, where the error is estimated in terms of computable quantities only. The a posteriori error estimate is used for choosing appropriate spaces in the multilevel algorithms, mesh refinements, as a stopping criterion and finally it gives an estimate of the total error. Received April 8, 1997 / Revised version received July 27, 1998 / Published online September 24, 1999     

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.     

Efficient computable error bounds for discontinuous Galerkin approximations of elliptic problems

We present guaranteed and computable both sided error bounds for the discontinuous Galerkin (DG) approximations of elliptic problems. These estimates are derived in the full DG-norm on purely functional grounds by the analysis of the respective differential problem, and thus, are applicable to any qualified DG approximation. Based on the triangle inequality, the underlying approach has the following steps for a given DG approximation: (1) computing a conforming approximation in the energy space using the Oswald interpolation operator, and (2) application of the existing functional a posteriori error estimates to the conforming approximation. Various numerical examples with varying difficulty in computing the error bounds, from simple problems of polynomial-type analytic solution to problems with analytic solution having sharp peaks, or problems with jumps in the coefficients of the partial differential equation operator, are presented which confirm the efficiency and the robustness of the estimates.     

A comparison of a posteriori error estimates for biharmonic problems solved by the FEM

The classical a posteriori error estimates are mostly oriented to the use in the finite element h$h$-methods while the contemporary higher-order hp$hp$-methods usually require new approaches in a posteriori error estimation. These methods hold a very important position among adaptive numerical procedures for solving ordinary as well as partial differential equations arising from various technical applications.     

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 posteriori error estimators for the first-order least-squares finite element method

In this paper, we propose a posteriori error estimators for certain quantities of interest for a first-order least-squares finite element method. In particular, we propose an a posteriori error estimator for when one is interested in where . Our a posteriori error estimators are obtained by assigning proper weight (in terms of local mesh size h_T) to the terms of the least-squares functional. An a posteriori error analysis yields reliable and efficient estimates based on residuals. Numerical examples are presented to show the effectivity of our error estimators.     

A unifying theory of a posteriori finite element error control

Summary Residual-based a posteriori error estimates are derived within a unified setting for lowest-order conforming, nonconforming, and mixed finite element schemes. The various residuals are identified for all techniques and problems as the operator norm |||| of a linear functional of the formin the variable of a Sobolev space V. The main assumption is that the first-order finite element space is included in the kernel Ker of . As a consequence, any residual estimator that is a computable bound of |||| can be used within the proposed frame without further analysis for nonconforming or mixed FE schemes. Applications are given for the Laplace, Stokes, and Navier-Lamè equations.Supported by the DFG Research Center Matheon Mathematics for key technologies in Berlin.     

A posteriori error estimates for mixed FEM in elasticity

A residue based reliable and efficient error estimator is established for finite element solutions of mixed boundary value problems in linear, planar elasticity. The proof of the reliability of the estimator is based on Helmholtz type decompositions of the error in the stress variable and a duality argument for the error in the displacements. The efficiency follows from inverse estimates. The constants in both estimates are independent of the Lamé constant , and so locking phenomena for are properly indicated. The analysis justifies a new adaptive algorithm for automatic mesh–refinement. Received July 17, 1997     

Residual type a posteriori error estimates for elliptic obstacle problems

Summary. A posteriori error estimators of residual type are derived for piecewise linear finite element approximations to elliptic obstacle problems. An instrumental ingredient is a new interpolation operator which requires minimal regularity, exhibits optimal approximation properties and preserves positivity. Both upper and lower bounds are proved and their optimality is explored with several examples. Sharp a priori bounds for the a posteriori estimators are given, and extensions of the results to double obstacle problems are briefly discussed. Received June 19, 1998 / Published online December 6, 1999     

A posteriori error estimates for hp finite element solutions of convex optimal control problems

In this paper, we present a posteriori error analysis for hp finite element approximation of convex optimal control problems. We derive a new quasi-interpolation operator of Clément type and a new quasi-interpolation operator of Scott-Zhang type that preserves homogeneous boundary condition. The Scott-Zhang type quasi-interpolation is suitable for an application in bounding the errors in L²-norm. Then hp a posteriori error estimators are obtained for the coupled state and control approximations. Such estimators can be used to construct reliable adaptive finite elements for the control problems.     

A posteriori error estimates for nonconforming finite element methods

Summary. Computable a posteriori error bounds for a large class of nonconforming finite element methods are provided for a model Poisson-problem in two and three space dimensions. Besides a refined residual-based a posteriori error estimate, an averaging estimator is established and an -estimate is included. The a posteriori error estimates are reliable and efficient; the proof of reliability relies on a Helmholtz decomposition. Received March 4, 1997 / Revised version received September 4, 2001 / Published online December 18, 2001     

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).     

Guaranteed and robust discontinuous Galerkin a posteriori error estimates for convection-diffusion-reaction problems

We propose and study a posteriori error estimates for convection-diffusion-reaction problems with inhomogeneous and anisotropic diffusion approximated by weighted interior-penalty discontinuous Galerkin methods. Our twofold objective is to derive estimates without undetermined constants and to analyze carefully the robustness of the estimates in singularly perturbed regimes due to dominant convection or reaction. We first derive locally computable estimates for the error measured in the energy (semi)norm. These estimates are evaluated using -conforming diffusive and convective flux reconstructions, thereby extending the previous work on pure diffusion problems. The resulting estimates are semi-robust in the sense that local lower error bounds can be derived using suitable cutoff functions of the local Péclet and Damköhler numbers. Fully robust estimates are obtained for the error measured in an augmented norm consisting of the energy (semi)norm, a dual norm of the skew-symmetric part of the differential operator, and a suitable contribution of the interelement jumps of the discrete solution. Numerical experiments are presented to illustrate the theoretical results.     

A posteriori error estimation of residual type for anisotropic diffusion-convection-reaction problems

This paper presents a robust a posteriori residual error estimator for diffusion-convection-reaction problems with anisotropic diffusion, approximated by a SUPG finite element method on isotropic or anisotropic meshes in R^d, d=2 or 3. The equivalence between the energy norm of the error and the residual error estimator is proved. Numerical tests confirm the theoretical results.     

A posteriori error estimators in the finite element method

Summary In the paper we develop a structured approach to thea posteriori estimation of the error in the approximation obtained via the finite element method. This aids the classification of existing estimators as well as allowing new estimators to be proposed for new situations. A class of abstract estimators for finite elements of orderp>1 in ,n=2, 3 based on exploiting the superconvergence phenomenon are analyzed.     

A posteriori dual-mixed adaptive finite element error control for Lamé and Stokes equations

A unified and robust mathematical model for compressible and incompressible linear elasticity can be obtained by rephrasing the Herrmann formulation within the Hellinger-Reissner principle. This quasi-optimally converging extension of PEERS (Plane Elasticity Element with Reduced Symmetry) is called Dual-Mixed Hybrid formulation (DMH). Explicit residual-based a posteriori error estimates for DMH are introduced and are mathematically shown to be locking-free, reliable, and efficient. The estimator serves as a refinement indicator in an adaptive algorithm for effective automatic mesh generation. Numerical evidence supports that the adaptive scheme leads to optimal convergence for Lamé and Stokes benchmark problems with singularities.     

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.