Interpolation error-based a posteriori error estimation for two-point boundary value problems and parabolic equations in one space dimension

Summary. I derive a posteriori error estimates for two-point boundary value problems and parabolic equations in one dimension based on interpolation error estimates. The interpolation error estimates are obtained from an extension of the error formula for the Lagrange interpolating polynomial in the case of symmetrically-spaced interpolation points. From this formula pointwise and seminorm a priori estimates of the interpolation error are derived. The interpolant in conjunction with the a priori estimates is used to obtain asymptotically exact a posteriori error estimates of the interpolation error. These a posteriori error estimates are extended to linear two-point boundary problems and parabolic equations. Computational results demonstrate the convergence of a posteriori error estimates and their effectiveness when combined with an hp-adaptive code for solving parabolic systems. Received April 17, 2000 / Revised version received September 25, 2000 / Published online May 30, 2001     

Superconvergence analysis of a linearized three‐step backward differential formula finite element method for nonlinear Sobolev equation

A linearized three‐step backward differential formula (BDF) Galerkin finite element method (FEM) is developed for nonlinear Sobolev equation with bilinear element. Temporal error and spatial error are discussed through introducing a time‐discrete system. Solutions of the time‐discrete system are bounded in H²‐norm by the temporal error. Superconvergence results of order O(h² + τ³) in H¹‐norm for the original variable are deduced based on the spatial error. Some new tricks are utilized to get higher order of the temporal error and the spatial error. At last, two numerical examples are provided to support the theoretical analysis. Here, h is the subdivision parameter, and τ is the time step.     

On Mixed Error Estimates for Elliptic Obstacle Problems

We establish in this paper sharp error estimates of residual type for finite element approximation to elliptic obstacle problems. The estimates are of mixed nature, which are neither of a pure a priori form nor of a pure a posteriori form but instead they are combined by an a priori part and an a posteriori part. The key ingredient in our derivation for the mixed error estimates is the use of a new interpolator which enables us to eliminate inactive data from the error estimators. One application of our mixed error estimates is to construct a posteriori error indicators reliable and efficient up to higher order terms, and these indicators are useful in mesh-refinements and adaptive grid generations. In particular, by approximating the a priori part with some a posteriori quantities we can successfully track the free boundary for elliptic obstacle problems.     

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.     

Fast enclosure for the minimum norm least squares solution of the matrix equation AXB = C

Fast algorithms for enclosing the minimum norm least squares solution of the matrix equation AXB = C are proposed. To develop these algorithms, theories for obtaining error bounds of numerical solutions are established. The error bounds obtained by these algorithms are verified in the sense that all the possible rounding errors have been taken into account. Techniques for accelerating the enclosure and obtaining smaller error bounds are introduced. Numerical results show the properties of the proposed algorithms. Copyright © 2015 John Wiley & Sons, Ltd.     

A Posteriori Error Estimators and Adaptivity for Finite Element Approximation of the Non-Homogeneous Dirichlet Problem

Techniques are developed for a posteriori error analysis of the non-homogeneous Dirichlet problem for the Laplacian giving computable error bounds for the error measured in the energy norm. The techniques are based on the equilibrated residual method that has proved to be reliable and accurate for the treatment of problems with homogeneous Dirichlet data. It is shown how the equilibrated residual method must be modified to include the practically important case of non-homogeneous Dirichlet data. Explicit and implicit a posteriori error estimators are derived and shown to be efficient and reliable. Numerical examples are provided illustrating the theory.     

A posteriori error estimation for convection dominated problems on anisotropic meshes

A singularly perturbed convection–diffusion problem in two and three space dimensions is discretized using the streamline upwind Petrov Galerkin (SUPG) variant of the finite element method. The dominant convection frequently gives rise to solutions with layers; hence anisotropic finite elements can be applied advantageously. The main focus is on a posteriori energy norm error estimation that is robust in the perturbation parameter and with respect to the mesh anisotropy. A residual error estimator and a local problem error estimator are proposed and investigated. The analysis reveals that the upper error bound depends on the alignment of the anisotropies of the mesh and of the solution. Hence reliable error estimation is possible for suitable anisotropic meshes. The lower error bound depends on the problem data via a local mesh Peclet number. Thus efficient error estimation is achieved for small mesh Peclet numbers. Altogether, error estimation approaches for isotropic meshes are successfully extended to anisotropic elements. Several numerical experiments support the analysis. Copyright © 2003 John Wiley & Sons, Ltd.     

A posteriori error estimates for H 1-Galerkin mixed finite-element method for parabolic problems

The purpose of this article is to derive a posteriori error estimates for the H ¹-Galerkin mixed finite element method for parabolic problems. We study both semidiscrete and fully discrete a posteriori error analyses using standard energy argument. A fully discrete a posteriori error analysis based on the backward Euler method is analysed and upper bounds for the errors are derived. The estimators yield upper bounds for the errors which are global in space and time. Our analysis is based on residual approach and the estimators are free from edge residuals.     

Mixed FEM of higher-order for a frictional contact problem

This paper presents mixed finite element methods of higher-order for an idealized frictional contact problem in linear elasticity. The approach relies on a saddle point formulation where the frictional contact condition is captured by a Lagrange multiplier. The convergence of the mixed scheme is proven and some a priori estimates for the h- and p-method are derived. Furthermore, a posteriori error estimates are presented which rely on the estimation of the discretization error of an auxiliary problem and some further terms capturing the error in the friction and complementary conditions. Numerical results confirm the applicability of the a posteriori error estimates within h- and hp-adaptive schemes. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Global Error Bounds for <Emphasis Type="Italic">γ</Emphasis>-paraconvex Multifunctions

In this paper, error bounds for γ-paraconvex multifunctions are considered. A Robinson-Ursescu type Theorem is given in normed spaces. Some results on the existence of global error bounds are presented. Perturbation error bounds are also studied.     

An a posteriori residual error estimator for the finite element method on anisotropic tetrahedral meshes

Summary. A new a posteriori residual error estimator is defined and rigorously analysed for anisotropic tetrahedral finite element meshes. All considerations carry over to anisotropic triangular meshes with minor changes only. The lower error bound is obtained by means of bubble functions and the corresponding anisotropic inverse inequalities. In order to prove the upper error bound, it is vital that an anisotropic mesh corresponds to the anisotropic function under consideration. To measure this correspondence, a so-called matching function is defined, and its discussion shows it to be a useful tool. With its help anisotropic interpolation estimates and subsequently the upper error bound are proven. Additionally it is pointed out how to treat Robin boundary conditions in a posteriori error analysis on isotropic and anisotropic meshes. A numerical example supports the anisotropic error analysis. Received April 6, 1999 / Revised version received July 2, 1999 / Published online June 8, 2000     

Reliable,goal‐oriented postprocessing for FE‐discretizations

In this article we present strategies to improve the quality of adaptive FE‐approximations measured in terms of linear functionals. The ideas are based on the so called dual‐weighted‐residual (DWR) approach to a posteriori error control for FE‐schemes. In more details, we exploit those parts of an underlying error representation, which are completely computable, to improve the FE‐solution. Furthermore, the remaining parts of the error identity can be estimated by well‐established a posteriori energy estimates yielding reliable error bounds for the postprocessed values. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

A Posteriori Error Estimators of Gradient Recovery Type for Elliptic Obstacle Problems

In this paper, we present a posteriori error estimates of gradient recovery type for elliptic obstacle problems. The a posteriori error estimates provide both lower and upper error bounds. It is shown to be equivalent to the discretization error in an energy type norm for general meshes. Furthermore, when the solution is smooth and the mesh is uniform, it is shown to be asymptotically exact. Some numerical results which demonstrate the theoretical results are also reported in this paper.     

On the p version of the finite element method in the presence of numerical integration

We analyze the error in the p version of the finite element method when the effect of the quadrature error is taken into account. We extend some results by Banerjee and Suri [Math. Comput. 59 , 1–20 (1992)] on the H¹-norm error to the case of the error in the L₂ norm. We investigate three sources of quadrature error that can occur: the error due to the numerical integration of the right-hand side, that due to nonconstant coefficients, and that due to the presence of mapped elements. Presented are various theoretical and computational examples regarding the sharpness of our results. In addition, we make a note on the use of numerical quadrature in conjunction with p-adaptive procedures and on the necessity of overintegration in the h version with linear elements, when the L₂ norm is of interest. © 1993 John Wiley & Sons, Inc.     

A posteriori error estimates of constrained optimal control problem governed by convection diffusion equations

In this paper, we study a posteriori error estimates of the edge stabilization Galerkin method for the constrained optimal control problem governed by convection-dominated diffusion equations. The residual-type a posteriori error estimators yield both upper and lower bounds for control u measured in L ²-norm and for state y and costate p measured in energy norm. Two numerical examples are presented to illustrate the effectiveness of the error estimators provided in this paper.      

hp-Adaptive Finite Element Methods for Variational Inequalities

In this work, we combine an hp–adaptive strategy with a posteriori error estimates for variational inequalities, which are given by contact problems. The a posteriori error estimates are obtained using a general approach based on the saddle point formulation of contact problems and making use of a yposteriori error estimates for variational equations. Error estimates are presented for obstacle problems and Signorini problems with friction. Numerical experiments confirm the reliability of the error estimates for finite elements of higher order. The use of the hp–adaptive strategy leads to meshes with the same characteristics as geometric meshes and to exponential convergence. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stable modification of explicitLU factors for simplex updates

This paper presents a method for modifying anLU factorization of a matrixA when a Simplex update is performed onA. At each stage of the algorithm, two options are allowed, and in one of these a row permutation ofL is made so that theLU structure is always maintained. The expanding file of operators used in other methods is not required. The result, of an error analysis is presented in which the error is bounded by the growth in the partial sums ofLU. Various pivotal strategies which minimize a growth bound are discussed and an example is given in which the worst case for growth occurs repeatedly. Numerical results from experiments on ill-conditioned problems show no tendency for error growth.     

Minimizing Maximum Weighted Error for Imprecise Computation Tasks

We consider the problem of preemptively scheduling a set of imprecise computation tasks on m ≥ 1 identical processors, with each task T_i having two weights, w_i and w′_i. Two performance metrics are considered: (1) the maximum w′-weighted error; (2) the total w-weighted error subject to the constraint that the maximum w′-weighted error is minimized. For the problem of minimizing the maximum w′-weighted error, we give an algorithm which runs in O(n³ log²n) time for multiprocessors and O(n²) time for a single processor. For the problem of minimizing the total w-weighted error subject to the constraint that the maximum w′-weighted error is minimized, we give an algorithm which also has the same time complexity.     

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     

Quadrature of smooth stochastic processes

Summary The integral of a stochastic process is estimated by means of classical quadrature formulae. In contrast to certain conventional methods, knowledge of covariances is not required, and no regularity conditions are assumed. Explicit error representations and error bounds with respect to theL ^p -norm are established.