Asymptotically exact functional error estimators based on superconvergent gradient recovery

The use of dual/adjoint problems for approximating functionals of solutions of PDEs with great accuracy or to merely drive a goal-oriented adaptive refinement scheme has become well-accepted, and it continues to be an active area of research. The traditional approach involves dual residual weighting (DRW). In this work we present two new functional error estimators and give conditions under which we can expect them to be asymptotically exact. The first is of DRW type and is derived for meshes in which most triangles satisfy an -approximate parallelogram property. The second functional estimator involves dual error estimate weighting (DEW) using any superconvergent gradient recovery technique for the primal and dual solutions. Several experiments are done which demonstrate the asymptotic exactness of a DEW estimator which uses a gradient recovery scheme proposed by Bank and Xu, and the effectiveness of refinement done with respect to the corresponding local error indicators. Resubmitted to Numerische Mathematik, June 30, 2005, with changes suggested by referees.     

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

A posteriori error estimations of residual type for Signorini's problem

This paper presents an a posteriori error analysis for the linear finite element approximation of the Signorini problem in two space dimensions. A posteriori estimations of residual type are defined and upper and lower bounds of the discretization error are obtained. We perform several numerical experiments in order to compare the convergence of the terms in the error estimator with the discretization error.     

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 mixed finite element method for solving the nonstationary stokes equation

Summary This paper deals with a mixed finite element method for approximating a fourth order initial value problem arising from the nonstationary Stokes problem. For piecewise linear shape functions error estimates are given with convergence rates similar to the elliptic case. Some numerical computations will illustrate the theoretical results.     

A discontinuous Galerkin method on refined meshes for the two-dimensional time-harmonic Maxwell equations in composite materials

In this paper, a discontinuous Galerkin method for the two-dimensional time-harmonic Maxwell equations in composite materials is presented. The divergence constraint is taken into account by a regularized variational formulation and the tangential and normal jumps of the discrete solution at the element interfaces are penalized. Due to an appropriate mesh refinement near exterior and interior corners, the singular behaviour of the electromagnetic field is taken into account. Optimal error estimates in a discrete energy norm and in the L²$L^2$-norm are proved in the case where the exact solution is singular.     

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.     

Numerical verifications of solutions for elliptic equations in nonconvex polygonal domains

Summary In this paper, methods for numerical verifications of solutions for elliptic equations in nonconvex polygonal domains are studied. In order to verify solutions using computer, it is necessary to determine some constants which appear in a priori error estimations. We propose some methods for determination of these constants. In numerical examples, calculating these constants for anL-shaped domain, we verify the solution of a nonlinear elliptic equation.     

Interpolation error estimates in W 1,p for degenerate Q 1 isoparametric elements

Optimal order error estimates in H ¹, for the Q ₁ isoparametric interpolation were obtained in Acosta and Durán (SIAM J Numer Anal37, 18–36, 1999) for a very general class of degenerate convex quadrilateral elements. In this work we show that the same conlusions are valid in W ^1,p for 1≤ p < 3 and we give a counterexample for the case p ≥ 3, showing that the result cannot be generalized for more regular functions. Despite this fact, we show that optimal order error estimates are valid for any p ≥ 1, keeping the interior angles of the element bounded away from 0 and π, independently of the aspect ratio. We also show that the restriction on the maximum angle is sharp for p ≥ 3.     

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.