Robust hierarchical a posteriori error estimators for stabilized convection–diffusion problems

We construct a hierarchical a posteriori error estimator for a stabilized finite element discretization of convection‐diffusion equations with height Péclet number. The error estimator is derived without the saturation assumption and without any comparison with the classical residual estimator. Besides, it is robust, such that the equivalence between the norm of the exact error and the error estimator is independent of the meshsize or the diffusivity parameter. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2012     

Robust a Posteriori Error Estimation for a Singularly Perturbed Reaction–Diffusion Equation on Anisotropic Tetrahedral Meshes

We consider a singularly perturbed reaction–diffusion problem and derive and rigorously analyse an a posteriori residual error estimator that can be applied to anisotropic finite element meshes. The quotient of the upper and lower error bounds is the so-called matching function which depends on the anisotropy (of the mesh and the solution) but not on the small perturbation parameter. This matching function measures how well the anisotropic finite element mesh corresponds to the anisotropic problem. Provided this correspondence is sufficiently good, the matching function is O(1). Hence one obtains tight error bounds, i.e. the error estimator is reliable and efficient as well as robust with respect to the small perturbation parameter. A numerical example supports the anisotropic error analysis.     

A posteriori error analysis for a boundary element method with nonconforming domain decomposition

We present and analyze an a posteriori error estimator based on mesh refinement for the solution of the hypersingular boundary integral equation governing the Laplacian in three dimensions. The discretization under consideration is a nonconforming domain decomposition method based on the Nitsche technique. Assuming a saturation property, we establish quasireliability and efficiency of the error estimator in comparison with the error in a natural (nonconforming) norm. Numerical experiments with uniform and adaptively refined meshes confirm our theoretical results. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 947–963, 2014     

A POSTERIORI ERROR ESTIMATES FOR A MODIFIED WEAK GALERKIN FINITE ELEMENT APPROXIMATION OF SECOND ORDER ELLIPTIC PROBLEMS WITH DG NORM

In this paper,we derive a residual based a posteriori error estimator for a modified weak Galerkin formulation of second order elliptic problems.We prove that the error estimator used for interior penalty discontinuous Galerkin methods still gives both upper and lower bounds for the modified weak Galerkin method,though they have essentially different bilinear forms.More precisely,we prove its reliability and efficiency for the actual error measured in the standard DG norm.We further provide an improved a priori error estimate under minimal regularity assumptions on the exact solution.Numerical results are presented to verify the theoretical analysis.     

Convergence of adaptive 3D BEM for weakly singular integral equations based on isotropic mesh‐refinement

We consider the adaptive lowest‐order boundary element method based on isotropic mesh refinement for the weakly‐singular integral equation for the three‐dimensional Laplacian. The proposed scheme resolves both, possible singularities of the solution as well as of the given data. The implementation thus only deals with discrete integral operators, that is, matrices. We prove that the usual adaptive mesh‐refining algorithm drives the corresponding error estimator to zero. Under an appropriate saturation assumption which is observed empirically, the sequence of discrete solutions thus tends to the exact solution within the energy norm. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

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 estimator for a mixed finite element method for Reissner-Mindlin plate

We present an a posteriori error estimator for a mixed finite element method for the Reissner-Mindlin plate model. The finite element method we deal with, was analyzed by Durán and Liberman in 1992 and can also be seen as a particular example of the general family analyzed by Brezzi, Fortin, and Stenberg in 1991. The estimator is based on the evaluation of the residual of the finite element solution. We show that the estimator yields local lower and global upper bounds of the error in the numerical solution in a natural norm for the problem, which includes the norms of the terms corresponding to the deflection and the rotation and a dual norm for the shearing force. The estimates are valid uniformly with respect to the plate thickness.

     

A Posteriori Error Estimators of Gradient Recovery Type for Fem of a Model Optimal Control Problem

In this paper, we derive an a posteriori error estimator of gradient recovery type for a model optimal control problem. We show that the a posteriori error estimator is equivalent to the discretization error in a norm of energy type on general meshes. Furthermore, when the solution of the control problem is smooth and the meshes are uniform, it is shown to be asymptotically exact.     

A Posteriori Error Analysis for Mixed Finite Element Solution of the Two-Dimensional Stationary Stokes Problem

In this paper we present a posteriori error estimator in a suitable norm of mixed finite element solution for two-dimensional stationary Stokes problem. The estimator is optimal in the sense that, up to multiplicative constant, the upper and lower bounds of the error are the same. The constants are independent of the mesh and the true solution of the problem.     

A posteriori estimators for nonlinear elliptic partial differential equations

Many works have reported results concerning the mathematical analysis of the performance of a posteriori error estimators for the approximation error of finite element discrete solutions to linear elliptic partial differential equations. For each estimator there is a set of restrictions defined in such a way that the analysis of its performance is made possible. Usually, the available estimators may be classified into two types, i.e., the implicit estimators (based on the solution of a local problem) and the explicit estimators (based on some suitable norm of the residual in a dual space). Regarding the performance, an estimator is called asymptotically exact if it is a higher-order perturbation of a norm of the exact error. Nowadays, one may say that there is a larger understanding about the behavior of estimators for linear problems than for nonlinear problems. The situation is even worse when the nonlinearities involve the highest derivatives occurring in the PDE being considered (strongly nonlinear PDEs). In this work we establish conditions under which those estimators, originally developed for linear problems, may be used for strongly nonlinear problems, and how that could be done. We also show that, under some suitable hypothesis, the estimators will be asymptotically exact, whenever they are asymptotically exact for linear problems. Those results allow anyone to use the knowledge about estimators developed for linear problems in order to build new reliable and robust estimators for nonlinear problems.     

Global a posteriori error estimates for convection–reaction–diffusion problems

In this note we propose a nonstandard technique for constructing global a posteriori error estimates for the stationary convection–reaction–diffusion equation. In order to estimate the approximation error in appropriate weighted energy norms, which measures the overall quality of the approximations, the underlying bilinear form is decomposed into several terms which can be directly computed or easily estimated from above using elementary tools of functional analysis. Several auxiliary parameters are introduced to construct such a splitting and tune the resulting upper error bound. It is demonstrated how these parameters can be chosen in some natural and convenient way for computations so that the weighted energy norm of the error is almost recovered, which shows that the estimates proposed are, in fact, quasi-sharp. The presented methodology is completely independent of numerical techniques used to compute approximate solutions. In particular, it is applicable to approximations which fail to satisfy the Galerkin orthogonality, e.g., due to an inconsistent stabilization, flux limiting, low-order quadrature rules, round-off and iteration errors etc. Moreover, the only constant that appears in the proposed error estimates is of global nature and comes from the Friedrichs–Poincaré inequality.     

Estimateur hiérarchique robuste pour un problème de perturbation singulière

In this Note, we show that a modified and simplified version of the estimator of Bank–Weiser can be used to define a robust a posteriori error estimator for singularly perturbed problem. We prove without comparison with a residual estimator or saturation assumption, the equivalence of the estimator with the error in the energy norm and the robusteness with respect to the diffusion coefficient. To cite this article: B. Achchab et al., C. R. Acad. Sci. Paris, Ser. I 336 (2003).     

A solution to anisotropic suboptimal filtering problem by convex optimization

This paper considers a robust filtering problem for a linear discrete time invariant system with measured and estimated outputs. The system is exposed to random disturbances with imprecisely known distributions generated by an unknown stable shaping filter from the Gaussian white noise. The stochastic uncertainty of the input disturbance is measured by the mean anisotropy functional. The estimation error is quantified by the anisotropic norm which is a stochastic analogue of the H _∞ norm. A sufficient condition for an estimator to exist and ensure that the error is less than a given threshold value is derived in form of a convex inequality on the determinant of a positive definite matrix and two linear matrix inequalities. The suboptimal problem setting results to a set of the estimators ensuring the anisotropic norm of the error to be strictly bounded thereby providing some additional degree of freedom to impose some additional constraints on the estimator performance specification.     

A Posteriori Error Estimates and a Local Refinement Strategy for a Finite Element Method to Solve Structural-Acoustic Vibration Problems

This paper deals with an adaptive technique to compute structural-acoustic vibration modes. It is based on an a posteriori error estimator for a finite element method free of spurious or circulation nonzero-frequency modes. The estimator is shown to be equivalent, up to higher order terms, to the approximate eigenfunction error, measured in a useful norm; moreover, the equivalence constants are independent of the corresponding eigenvalue, the physical parameters, and the mesh size. This a posteriori error estimator yields global upper and local lower bounds for the error and, thus, it may be used to design adaptive algorithms. We propose a local refinement strategy based on this estimator and present a numerical test to assess the efficiency of this technique.     

A posteriori error estimators for mixed finite element methods in linear elasticity

Summary. Three a posteriori error estimators for PEERS and BDMS elements in linear elasticity are presented: one residual error estimator and two estimators based on the solution of auxiliary local problems with different boundary conditions. All of them are reliable and efficient with respect to the standard norm and furthermore robust for nearly incompressible materials.Correspondence to: R. Verfürth     

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.     

Minimum Hellinger distance estimation for supercritical Galton–Watson processes

This paper studies the asymptotic behavior of the minimum Hellinger distance estimator of the underlying parameter in a supercritical branching process whose offspring distribution is known to belong to a parametric family. This estimator is shown to be asymptotically normal, efficient at the true model and robust against gross errors. These extend the results of Beran (Ann. Statist. 5, 445–463 (1977)) from an i.i.d., continuous setup to a dependent, discrete setup.     

Analysis of a subdomain‐based error estimator for finite element approximations of elliptic problems

In this article we analyze a subdomain residual error estimator for finite element approximations of elliptic problems. It is obtained by solving local problems on patches of elements in weighted spaces and provides an upper bound on the energy norm of the error when the local problems are solved in sufficiently enriched discrete spaces. A guaranteed lower bound on the error is also derived by a simple postprocess of the solutions to the local problems. Numerical tests show very good effectivity indices for both the upper and lower bounds and a strong reliability of this estimator even for coarse meshes. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 165–192, 2004     

Robust a posteriori error estimation for the nonconforming Fortin-Soulie finite element approximation

We obtain a computable a posteriori error bound on the broken energy norm of the error in the Fortin-Soulie finite element approximation of a linear second order elliptic problem with variable permeability. This bound is shown to be efficient in the sense that it also provides a lower bound for the broken energy norm of the error up to a constant and higher order data oscillation terms. The estimator is completely free of unknown constants and provides a guaranteed numerical bound on the error.

     

不可压Navier-Stokes方程迎风格式及后验误差估计

1引言不可压Navier-Stokes方程作为流体力学的基本方程,其数值计算一直是科学与工程计算关心的问题．本文考虑定常问题: -ε△u (u·▽)u ▽p = f x∈Ω,▽·u=0 x∈Ω, (1) u =0 x∈(?)Ω．这里ε=1／Re是Reynolds数的倒数,u=(u1,u2,…,ud)为待求流速场,p是待求压力场,f=(f1,f2,…,fd)是给定的体力．Ωv(?) Rd(d=2,3)是有界区域,且具有分片Lipschitz连续边界(?)Ω．