Gradient recovery type a posteriori error estimate for finite element approximation on non-uniform meshes

In this paper, we derive gradient recovery type a posteriori error estimate for the finite element approximation of elliptic equations. We show that a posteriori error estimate provide both upper and lower bounds for the discretization error on the non-uniform meshes. Moreover, it is proved that a posteriori error estimate is also asymptotically exact on the uniform meshes if the solution is smooth enough. The numerical results demonstrating the theoretical results are also presented in this paper.     

Robust error estimates for the finite element approximation of elliptic optimal control problems

In this paper, we consider the finite element approximation of an elliptic optimal control problem. Based on an assumption on the adjoint state of the continuous problem with a small parameter, which represents a regularization of the bang-bang type control problem, we derive robust a priori error estimates for optimal control and state and a posteriori error estimate is also presented. Numerical experiments confirm our theoretical results.     

Analysis of recovery type a posteriori error estimators for mildly structured grids

Some recovery type error estimators for linear finite elements are analyzed under 0)$"> regular grids. Superconvergence of order is established for recovered gradients by three different methods. As a consequence, a posteriori error estimators based on those recovery methods are asymptotically exact.

     

A posteriori error control for finite element approximations of elliptic eigenvalue problems

We develop a new approach to a posteriori error estimation for Galerkin finite element approximations of symmetric and nonsymmetric elliptic eigenvalue problems. The idea is to embed the eigenvalue approximation into the general framework of Galerkin methods for nonlinear variational equations. In this context residual-based a posteriori error representations are available with explicitly given remainder terms. The careful evaluation of these error representations for the concrete situation of an eigenvalue problem results in a posteriori error estimates for the approximations of eigenvalues as well as eigenfunctions. These suggest local error indicators that are used in the mesh refinement process.     

A note on the residual type a posteriori error estimates for finite element eigenpairs of nonsymmetric elliptic eigenvalue problems

In this paper we study the residual type a posteriori error estimates for general elliptic (not necessarily symmetric) eigenvalue problems. We present estimates for approximations of semisimple eigenvalues and associated eigenvectors. In particular, we obtain the following new results: 1) An error representation formula which we use to reduce the analysis of the eigenvalue problem to the analysis of the associated source problem; 2) A local lower bound for the error of an approximate finite element eigenfunction in a neighborhood of a given mesh element T.     

A posteriori error analysis of the discontinuous finite element methods for first order hyperbolic problems

In this paper, we consider the a posteriori error analysis of discontinuous Galerkin finite element methods for the steady and nonsteady first order hyperbolic problems with inflow boundary conditions. We establish several residual-based a posteriori error estimators which provide global upper bounds and a local lower bound on the error. Further, for nonsteady problem, we construct a fully discrete discontinuous finite element scheme and derive the a posteriori error estimators which yield global upper bound on the error in time and space. Our a posteriori error analysis is based on the mesh-dependent a priori estimates for the first order hyperbolic problems. These a posteriori error analysis results can be applied to develop the adaptive discontinuous finite element methods.     

Residual-based a posteriori error estimates of nonconforming finite element method for elliptic problems with Dirac delta source terms

Two residual-based a posteriori error estimators of the nonconforming Crouzeix-Raviart element are derived for elliptic problems with Dirac delta source terms.One estimator is shown to be reliable and efficient,which yields global upper and lower bounds for the error in piecewise W1,p seminorm.The other one is proved to give a global upper bound of the error in Lp-norm.By taking the two estimators as refinement indicators,adaptive algorithms are suggested,which are experimentally shown to attain optimal convergence orders.     

Local a posteriori estimates for pointwise gradient errors in finite element methods for elliptic problems

We prove local a posteriori error estimates for pointwise gradient errors in finite element methods for a second-order linear elliptic model problem. First we split the local gradient error into a computable local residual term and a weaker global norm of the finite element error (the ``pollution term'). Using a mesh-dependent weight, the residual term is bounded in a sharply localized fashion. In specific situations the pollution term may also be bounded by computable residual estimators. On nonconvex polygonal and polyhedral domains in two and three space dimensions, we may choose estimators for the pollution term which do not employ specific knowledge of corner singularities and which are valid on domains with cracks. The finite element mesh is only required to be simplicial and shape-regular, so that highly graded and unstructured meshes are allowed.

     

Functional a posteriori error estimates for discontinuous Galerkin approximations of elliptic problems

In this article, we develop functional a posteriori error estimates for discontinuous Galerkin (DG) approximations of elliptic boundary‐value problems. These estimates are based on a certain projection of DG approximations to the respective energy space and functional a posteriori estimates for conforming approximations developed by S. Repin (see e.g., Math Comp 69 (2000) 481–500). On these grounds, we derive two‐sided guaranteed and computable bounds for the errors in “broken” energy norms. A series of numerical examples presented confirm the efficiency of the estimates. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Sharp upper global a posteriori error estimates for nonlinear elliptic variational problems

The paper is devoted to the problem of verification of accuracy of approximate solutions obtained in computer simulations. This problem is strongly related to a posteriori error estimates, giving computable bounds for computational errors and detecting zones in the solution domain where such errors are too large and certain mesh refinements should be performed. A mathematical model embracing nonlinear elliptic variational problems is considered in this work. Based on functional type estimates developed on an abstract level, we present a general technology for constructing computable sharp upper bounds for the global error for various particular classes of elliptic problems. Here the global error is understood as a suitable energy type difference between the true and computed solutions. The estimates obtained are completely independent of the numerical technique used to obtain approximate solutions, and are sharp in the sense that they can be, in principle, made as close to the true error as resources of the used computer allow. The latter can be achieved by suitably tuning the auxiliary parameter functions, involved in the proposed upper error bounds, in the course of the calculations.     

Energy norm a posteriori error estimates for mixed finite element methods

This paper deals with the a posteriori error analysis of mixed finite element methods for second order elliptic equations. It is shown that a reliable and efficient error estimator can be constructed using a postprocessed solution of the method. The analysis is performed in two different ways: under a saturation assumption and using a Helmholtz decomposition for vector fields.

     

Functional a posteriori error estimates for incremental models in elasto-plasticity

We consider incremental problem arising in elasto-plastic models with isotropic hardening. Our goal is to derive computable and guaranteed bounds of the difference between the exact solution and any function in the admissible (energy) class of the problem considered. Such estimates are obtained by an advanced version of the variational approach earlier used for linear boundary-value problems and nonlinear variational problems with convex functionals [24, 30]. They do no contain mesh-dependent constants and are valid for any conforming approximations regardless of the method used for their derivation. It is shown that the structure of error majorant reflects properties of the exact solution so that the majorant vanishes only if an approximate solution coincides with the exact one. Moreover, it possesses necessary continuity properties, so that any sequence of approximations converging to the exact solution in the energy space generates a sequence of positive numbers (explicitly computable by the majorant functional) that tends to zero.      

A posteriori error estimates with the finite element method of lines for a Sobolev equation

A posteriori error estimates for semidiscrete finite element methods for a nonlinear Sobolev equation are considered. The error estimates are obtained by solving local nonlinear or linear pseudo‐parabolic equations for corrections to the solution on each element. The ratios of these estimates and the true errors are proved to converge to 1, implying that the estimates can be used as indicators in adaptive schemes for the problem. Numerical results underline our theoretical results. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

A posteriori error estimates of edge residual type of finite element method for nonmonotone quasi‐linear elliptic problems

In this article, we study the edge residual‐based a posteriori error estimates of conforming linear finite element method for nonmonotone quasi‐linear elliptic problems. It is proven that edge residuals dominate a posteriori error estimates. Up to higher order perturbations, edge residuals can act as a posteriori error estimators. The global reliability and local efficiency bounds are established both in H ¹‐norm and L ²‐norm. Numerical experiments are provided to illustrate the performance of the proposed error estimators. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 813–837, 2014     

A posteriori error estimate for a modified weak Galerkin method solving elliptic problems

A residual‐type a posteriori error estimator is proposed and analyzed for a modified weak Galerkin finite element method solving second‐order elliptic problems. This estimator is proven to be both reliable and efficient because it provides computable upper and lower bounds on the actual error in a discrete H¹‐norm. Numerical experiments are given to illustrate the effectiveness of the this error estimator. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 381–398, 2017     

Local error estimates of the finite element method for an elliptic problem with a Dirac source term

The solutions of elliptic problems with a Dirac measure right‐hand side are not in dimension and therefore the convergence of the finite element solutions is suboptimal in the ‐norm. In this article, we address the numerical analysis of the finite element method for the Laplace equation with Dirac source term: we consider, in dimension 3, the Dirac measure along a curve and, in dimension 2, the punctual Dirac measure. The study of this problem is motivated by the use of the Dirac measure as a reduced model in physical problems, for which high accuracy of the finite element method at the singularity is not required. We show a quasioptimal convergence in the ‐norm, for on subdomains which exclude the singularity; in the particular case of Lagrange finite elements, an optimal convergence in ‐norm is shown on a family of quasiuniform meshes. Our results are obtained using local Nitsche and Schatz‐type error estimates, a weak version of Aubin‐Nitsche duality lemma and a discrete inf‐sup condition. These theoretical results are confirmed by numerical illustrations.     

A posteriori error estimate for finite volume element method of the parabolic equations

In this article, residual‐type a posteriori error estimates are studied for finite volume element (FVE) method of parabolic equations. Residual‐type a posteriori error estimator is constructed and the reliable and efficient bounds for the error estimator are established. Residual‐type a posteriori error estimator can be used to assess the accuracy of the FVE solutions in practical applications. Some numerical examples are provided to confirm the theoretical results. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 259–275, 2017     

Adaptive finite element algorithms for eigenvalue problems based on local averaging type a posteriori error estimates

The local averaging technique has become a popular tool in adaptive finite element methods for solving partial differential boundary value problems since it provides efficient a posteriori error estimates by a simple postprocessing. In this paper, the technique is introduced to solve a class of symmetric eigenvalue problems. Its efficiency and reliability are proved by both the theory and numerical experiments structured meshes as well as irregular meshes. Dedicated to Charles A. Micchelli on his 60th birthday Mathematics subject classifications (2000) 65N15, 65N25, 65N30, 65N50. Subsidized by the Special Funds for Major State Basic Research Projects, and also supported in part by the Chinese National Natural Science Foundation and the Knowledge Innovation Program of the Chinese Academy of Sciences.     

A posteriori error estimates for the Crank-Nicolson method for parabolic equations

We derive optimal order a posteriori error estimates for time discretizations by both the Crank-Nicolson and the Crank-Nicolson-Galerkin methods for linear and nonlinear parabolic equations. We examine both smooth and rough initial data. Our basic tool for deriving a posteriori estimates are second-order Crank-Nicolson reconstructions of the piecewise linear approximate solutions. These functions satisfy two fundamental properties: (i) they are explicitly computable and thus their difference to the numerical solution is controlled a posteriori, and (ii) they lead to optimal order residuals as well as to appropriate pointwise representations of the error equation of the same form as the underlying evolution equation. The resulting estimators are shown to be of optimal order by deriving upper and lower bounds for them depending only on the discretization parameters and the data of our problem. As a consequence we provide alternative proofs for known a priori rates of convergence for the Crank-Nicolson method.

     

A posteriori estimates for eigenvalue/vector approximations

We combine abstract eigenvalue/eigenvector estimates (from our earlier work) with a saturation assumption for finite element solution of associated stationary problem to obtain a posteriori estimates of the accuracy of finite element Rayleigh–Ritz approximations. Attention will be payed to the interplay between the accuracy estimate for the finite element method and a strategy for generating an adapted mesh. The obtained results use a preconditioned residuum of Neymeyr and extend his study of eigenvalue approximations with eigenvector estimates. We also prove that this eigenvalue estimator is equivalent to the global error. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)