Maximum-norm error analysis of a non-monotone FEM for a singularly perturbed reaction-diffusion problem

A non-monotone FEM discretization of a singularly perturbed one-dimensional reaction-diffusion problem whose solution exhibits strong layers is considered. The method is shown to be maximum-norm stable although it is not inverse monotone. Both a priori and a posteriori error bounds in the maximum norm are derived. The a priori result can be used to deduce uniform convergence of various layer-adapted meshes proposed in the literature. Numerical experiments complement the theoretical results. AMS subject classification (2000)  65L10, 65L50, 65L60     

Maximum norm a posteriori error estimates for a 1D singularly perturbed semilinear reaction-diffusion problem

^{A singularly perturbed semilinear two-point boundary-value problemis discretized on arbitrary non-uniform meshes. We present second-ordermaximum norm a posteriori error estimates that hold true uniformlyin the small parameter. Their application to monitor-functionequidistribution and a posteriori mesh refinement are discussed.Numerical results are presented that support our theoreticalestimates.}     

Robust a posteriori error estimators for a singularly perturbed reaction-diffusion equation

We derive robust a posteriori error estimators for a singularly perturbed reaction-diffusion equation. Here, robust means that the estimators yield global upper and local lower bounds on the error measured in the energy norm such that the ratio of the upper and lower bounds is bounded from below and from above by constants which do neither depend on any meshsize nor on the perturbation parameter. The estimators are based either on the evaluation of local residuals or on the solution of discrete local Dirichlet or Neumann problems. Received June 5, 1996     

Robust a posteriori error estimates for conforming discretizations of a singularly perturbed reaction-diffusion problem on anisotropic meshes

In this paper, we derive robust a posteriori error estimates for conforming approximations to a singularly perturbed reaction-diffusion problem on anisotropic meshes, since the solution in general exhibits anisotropic features, e.g., strong boundary or interior layers. Based on the anisotropy of the mesh elements, we improve the a posteriori error estimates developed by Cheddadi et al., which are reliable and efficient on isotropic meshes but fail on anisotropic ones. Without the assumption that the mesh is shape-regular, the resulting mesh-dependent error estimator is shown to be reliable, efficient and robust with respect to the reaction coefficient, as long as the anisotropic mesh sufficiently reflects the anisotropy of the solution. We present our results in the framework of the vertex-centered finite volume method but their nature is general for any conforming one, like the piecewise linear finite element one. Our estimates are based on the usual H(div)-conforming, locally conservative flux reconstruction in the lowest-order Raviart-Thomas space on a dual mesh associated with the original anisotropic simplex one. Numerical experiments in 2D confirm that our estimates are reliable, efficient and robust on anisotropic meshes.     

A posteriori error estimation for arbitrary order FEM applied to singularly perturbed one-dimensional reaction-diffusion problems

FEM discretizations of arbitrary order r are considered for a singularly perturbed one-dimensional reaction-diffusion problem whose solution exhibits strong layers. A posteriori error bounds of interpolation type are derived in the maximum norm. An adaptive algorithm is devised to resolve the boundary layers. Numerical experiments complement our theoretical results.     

Maximum norm a posteriori error estimate for a 3d singularly perturbed semilinear reaction-diffusion problem

A singularly perturbed semilinear reaction-diffusion problem in the unit cube, is discretized on arbitrary nonuniform tensor-product meshes. We establish a second-order maximum norm a posteriori error estimate that holds true uniformly in the small diffusion parameter. No mesh aspect ratio condition is imposed. This result is obtained by combining (i) sharp bounds on the Green’s function of the continuous differential operator in the Sobolev W ^1,1 and W ^2,1 norms and (ii) a special representation of the residual in terms of an arbitrary current mesh and the current computed solution. Numerical results on a priori chosen meshes are presented that support our theoretical estimate.     

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     

Uniform error estimates in the finite element method for a singularly perturbed reaction-diffusion problem

Consider the problem with homogeneous Neumann boundary condition in a bounded smooth domain in . The whole range is treated. The Galerkin finite element method is used on a globally quasi-uniform mesh of size ; the mesh is fixed and independent of .

A precise analysis of how the error at each point depends on and is presented. As an application, first order error estimates in , which are uniform with respect to , are given.

     

Functional a posteriori error estimates for elliptic problems in exterior domains

This paper is concerned with the derivation of computable and guaranteed upper bounds of the difference between the exact and approximate solutions of an exterior domain boundary value problem for a linear elliptic equation. Our analysis is based upon purely functional argumentation and does not attract specific properties of an approximation method. Therefore, the estimates derived in the paper at hand are applicable to any approximate solution that belongs to the corresponding energy space. Such estimates (also called error majorants of functional type) were derived earlier for problems in bounded domains of RN. Bibliography: 4 titles. Illustrations: 1 figure.     

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.     

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     

A posteriori error estimation for a singularly perturbed Volterra integro-differential equation

Numerical Algorithms - A singularly perturbed Volterra integro-differential equation with an integrable singularity in the integral term is considered. The upwind difference method is used to...     

Two-sided a posteriori error estimates for linear elliptic problems with mixed boundary conditions

The paper is devoted to 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 consisting of a linear elliptic (reaction-diffusion) equation with a mixed Dirichlet/Neumann/Robin boundary condition is considered in this work. On the base of this model, we present simple technologies for straightforward constructing computable upper and lower bounds for the error, which is understood as the difference between the exact solution of the model and its approximation measured in the corresponding energy norm. The estimates obtained are completely independent of the numerical technique used to obtain approximate solutions and are “flexible” in the sense that they can be, in principle, made as close to the true error as the resources of the used computer allow. This work was supported by the Academy Research Fellowship No. 208628 from the Academy of Finland.     

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 mixed finite element approximations of elliptic problems

We derive residual based a posteriori error estimates of the flux in L ²-norm for a general class of mixed methods for elliptic problems. The estimate is applicable to standard mixed methods such as the Raviart–Thomas–Nedelec and Brezzi–Douglas–Marini elements, as well as stabilized methods such as the Galerkin-Least squares method. The element residual in the estimate employs an elementwise computable postprocessed approximation of the displacement which gives optimal order.     

Superconvergence error estimates of discontinuous Galerkin time stepping for singularly perturbed parabolic problems

Numerical Algorithms - A parabolic convection-diffusion-reaction problem is discretized by the non-symmetric interior penalty Galerkin (NIPG) method in space and discontinuous Galerkin (DG) method...     

Functional a posteriori estimates for elliptic optimal control problems

In this note, new results on functional type a posteriori estimates for elliptic optimal control problems with control constraints are presented. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A posteriori error estimation for a fully discrete discontinuous Galerkin approximation to a kind of singularly perturbed problems

The time variable in the semi-discrete problem is still continuous. In order to obtain an expected numerical solution, discretization the time variable from the semi-discrete form (full discretization) is needed. For a kind of non-stationarily singular perturbation problem in 1D, a fully discrete discontinuous Galerkin (DG) method is considered. That is to say, space variable is discretized with a primal DG method with penalty, and time variable is done using the backward Euler method. By virtue of duality arguments, inverse estimation of finite element method and interpolation theory, we present a residual-type a posteriori error indicator, which is usually used for adaptivity.     

Pointwise a posteriori error control for elliptic obstacle problems

We consider a finite element method for the elliptic obstacle problem over polyhedral domains in ^d, which enforces the unilateral constraint solely at the nodes. We derive novel optimal upper and lower a posteriori error bounds in the maximum norm irrespective of mesh fineness and the regularity of the obstacle, which is just assumed to be Hölder continuous. They exhibit optimal order and localization to the non-contact set. We illustrate these results with simulations in 2d and 3d showing the impact of localization in mesh grading within the contact set along with quasi-optimal meshes. Partially supported by NSF Grant DMS-9971450 and NSF/DAAD Grant INT-9910086.Partially suported by DAAD/NSF grant ``Projektbezogene Förderung des Wissenschaftleraustauschs in den Natur-, Ingenieur- und den Sozialwissenschaften mit der NSF'.Partially supported by DAAD/NSF grant ``Projektbezogene Förderung des Wissenschaftleraustauschs in den Natur-, Ingenieur- und den Sozialwissenschaften mit der NSF', and by the TMR network ``Viscosity solutions and their Applications', Italian M.I.U.R. projects ``Scientific Computing: Innovative Models and Numerical Methods' and ``Symmetries, Geometric Structures, Evolution and Memory in PDEs'.Mathematics Subject Classification (1991):65N15, 65N30, 35J85