A posteriori error estimators for optimal distributed control governed by the first‐order linear hyperbolic equation: DG method

In this article, We analyze the ‐version of the discontinuous Galerkin finite element method (DGFEM) for the distributed first‐order linear hyperbolic optimal control problems. We derive a posteriori error estimators on general finite element meshes which are sharp in the mesh‐width . These error estimators are shown to be useful in adaptive finite element approximation for the optimal control problems. For the DGFEM we admit very general irregular meshes. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

A posteriori error analysis of an augmented fully mixed formulation for the nonisothermal Oldroyd–Stokes problem

In this article, we consider an augmented fully mixed variational formulation that has been recently proposed for the nonisothermal Oldroyd–Stokes problem, and develop an a posteriori error analysis for the 2‐D and 3‐D versions of the associated mixed finite element scheme. More precisely, we derive two reliable and efficient residual‐based a posteriori error estimators for this problem on arbitrary (convex or nonconvex) polygonal and polyhedral regions. The reliability of the proposed estimators draws mainly upon the uniform ellipticity of the bilinear forms of the continuous formulation, suitable assumptions on the domain and the data, stable Helmholtz decompositions, and the local approximation properties of the Clément and Raviart–Thomas operators. On the other hand, inverse inequalities, the localization technique based on bubble functions, and known results from previous works are the main tools yielding the efficiency estimate. Finally, several numerical results confirming the properties of the a posteriori error estimators and illustrating the performance of the associated adaptive algorithms are reported.     

Adaptive finite element for semi-linear convection–diffusion problems

In this article a strategy of adaptive finite element for semi-linear problems, based on minimizing a residual-type estimator, is reported. We get an a posteriori error estimate which is asymptotically exact when the mesh size h tends to zero. By considering a model problem, the quality of this estimator is checked. It is numerically shown that without constraint on the mesh size h, the efficiency of the a posteriori error estimate can fail dramatically. This phenomenon is analysed and an algorithm which equidistributes the local estimators under the constraint h ⩽ h _max is proposed. This algorithm allows to improve the computed solution for semi-linear convection–diffusion problems, and can be used for stabilizing the Lagrange finite element method for linear convection–diffusion problems. This revised version was published online in June 2006 with corrections to the Cover Date.     

A posteriori finite element error estimators for indefinite elliptic boundary value problems ∗

A general construction technique is presented for a posteriori error estimators of finite element solutions of elliptic boundary value problems that satisfy a Gång inequality. The estimators are obtained by an element–by–element solution of ‘weak residual’ with or without considering element boundary residuals. There is no order restriction on the finite element spaces used for the approximate solution or the error estimation; that is, the design of the estimators is applicable in connection with either one of the h–p–, or hp– formulations of the finite element method. Under suitable assumptions it is shown that the estimators are bounded by constant multiples of the true error in a suitable norm. Some numerical results are given to demonstrate the effectiveness and efficiency of the approach.     

Mixed finite element methods in fluid structure systems

Summary. The eigenvalue problem describing the frequencies of a fluid vibrating in a rigid cavity or within moving boundaries is considered. Based on the method of Lagrange multipliers, a three field mixed formulation is introduced in order to avoid the spurious circulating modes. Stability and optimal error bounds are proved for two choices of finite element spaces. Received October 20, 1992 / Revised version received May 23, 1995     

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.

     

Residual Based a Posteriori Error Estimates for Convex Optimal Control Problems Governed by Stokes-Darcy Equations

In this paper, we derive a posteriori error estimates for finite element approximations of the optimal control problems governed by the Stokes-Darcy system. We obtain a posteriori error estimators for both the state and the control based on the residual of the finite element approximation. It is proved that the a posteriori error estimate provided in this paper is both reliable and efficient.     

Superconvergence of the - version of the finite element method in one dimension

In this paper, we investigate the superconvergence properties of the h-p version of the finite element method (FEM) for two-point boundary value problems. A postprocessing technique for the h-p finite element approximation is analyzed. The analysis shows that the postprocess improves the order of convergence. Furthermore, we obtain asymptotically exact a posteriori error estimators based on the postprocessing results. Numerical examples are included to illustrate the theoretical analysis.     

特征值问题协调/非协调有限元后验误差分析

本文研究对称椭圆特征值问题的有限元后验误差估计,包括协调元和非协调元,具有下列特色:(1)对协调/非协调元建立了有限元特征函数uh的误差与相应的边值问题有限元解的误差在局部能量模意义下的恒等关系式,该边值问题的右端为有限元特征值λh与uh的乘积,有限元解恰好为uh.从而边值问题有限元解在能量模意义下的局部后验误差指示子,包括残差型和重构型后验误差指示子,成为有限元特征函数在能量模意义下的局部后验误差指示子.(2)讨论了协调有限元特征函数的基于插值后处理的梯度重构型后验误差估计,对有限元特征函数的导数得到了最大模意义下的渐近准确局部后验误差指示子.     

On the Convergence of Adaptive Finite Element Methods

State of the art simulations in computational mechanics aim reliability and efficiency via adaptive finite element methods (AFEMs) with a posteriori error control. The a priori convergence of finite element methods is justified by the density property of the sequence of finite element spaces which essentially assumes a quasi‐uniform mesh‐refining. The advantage is guaranteed convergence for a large class of data and solutions; the disadvantage is a global mesh refinement everywhere accompanied by large computational costs. AFEMs automatically refine exclusively wherever the refinement indication suggests to do so and so violate the density property on purpose. Then, the a priori convergence of AFEMs is not guaranteed automatically and, in fact, crucially depends on algorithmic details. The advantage of AFEMs is a more effective mesh accompanied by smaller computational costs in many practical examples; the disadvantage is that the desirable error reduction property is not always guaranteed a priori. Efficient error estimators can justify a numerical approximation a posteriori and so achieve reliability. But it is not clear from the start that the adaptive mesh‐refinement will generate an accurate solution at all. This paper discusses particular versions of an AFEMs and their analyses for error reduction, energy reduction, and convergence results for linear and nonlinear problems. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Recovery type a posteriori estimates and superconvergence for nonconforming FEM of eigenvalue problems

The main goal of this paper is to present recovery type a posteriori error estimators and superconvergence for the nonconforming finite element eigenvalue approximation of self-adjoint elliptic equations by projection methods. Based on the superconvergence results of nonconforming finite element for the eigenfunction we derive superconvergence and recovery type a posteriori error estimates of the eigenvalue. The results are based on some regularity assumption for the elliptic problem and are applicable to the lowest order nonconforming finite element approximations of self-adjoint elliptic eigenvalue problems with quasi-regular partitions. Therefore, the results of this paper can be employed to provide useful a posteriori error estimators in practical computing under unstructured meshes.     

hp–Adaptive composite discontinuous Galerkin methods for elliptic problems on complicated domains

In this article, we develop the a posteriori error estimation of hp–version discontinuous Galerkin composite finite element methods for the discretization of second‐order elliptic partial differential equations. This class of methods allows for the approximation of problems posed on computational domains which may contain a huge number of local geometrical features, or microstructures. Although standard numerical methods can be devised for such problems, the computational effort may be extremely high, as the minimal number of elements needed to represent the underlying domain can be very large. In contrast, the minimal dimension of the underlying composite finite element space is independent of the number of geometric features. Computable bounds on the error measured in terms of a natural (mesh‐dependent) energy norm are derived. Numerical experiments highlighting the practical application of the proposed estimators within an automatic hp–adaptive refinement procedure will be presented. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1342–1367, 2014     

Stabilized finite element method for Navier-Stokes equations with physical boundary conditions

This paper deals with the numerical approximation of the 2D and 3D Navier-Stokes equations, satisfying nonstandard boundary conditions. This lays on the finite element discretisation of the corresponding Stokes problem, which is achieved through a three-fields stabilized mixed formulation. A priori and a posteriori error bounds are established for the nonlinear problem, ascertaining the convergence of the method. Finally, numerical tests are presented, including mesh refinement via error indicators.

     

A type of new posteriori error estimators for stokes problems

1. IntroductionIn the numerical approximation of PDE, it is often very importals to detect regionswhere the accuracy of the numerical solution is degraded by local singularities of the solutionof the continuous problem such as the singularity near the re-entrant corller. An obviousremedy is to refine the discretization in the critical regions, i.e., to place more gridpointswhere the solution is less regular. The question is how to identify these regions automdticallyand how to determine a goo…     

A posteriori error analysis for the normal compliance problem

In this work, a contact problem between a linear elastic material and a deformable obstacle is numerically analyzed. The contact is modeled using the well-known normal compliance contact condition. The weak formulation leads to a nonlinear variational equation which is approximated by using the finite element method. A priori error estimates are recalled. Then, we define an a posteriori error estimator of residual type to evaluate the accuracy of the finite element approximation of the problem. Upper and lower bounds of the discretization error are proved for this estimator.     

A posteriori error estimation for the Stokes–Darcy coupled problem on anisotropic discretization

This paper presents an a posteriori error analysis for the stationary Stokes–Darcy coupled problem approximated by finite element methods on anisotropic meshes in or 3. Korn's inequality for piecewise linear vector fields on anisotropic meshes is established and is applied to non‐conforming finite element method. Then the existence and uniqueness of the approximation solution are deduced for non‐conforming case. With the obtained finite element solutions, the error estimators are constructed and based on the residual of model equations plus the stabilization terms. The lower error bound is proved 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 useful tool. With its help, the upper error bound is shown by means of the corresponding anisotropic interpolation estimates and a special Helmholtz decomposition in both media. Copyright © 2016 John Wiley & Sons, Ltd.     

Estimate of the Constant in Two Strengthened C.B.S. Inequalities for F.E.M. Systems of 2D Elasticity: Application to Multilevel Methods and a Posteriori Error Estimators

The constant γ in the strengthened Cauchy-Buniakowski-Schwarz (C.B.S.) inequality plays a crucial role in the convergence rate of multilevel iterative methods as well as in the efficiency of a posteriori error estimators, that is in the framework of finite element approximations of SPD problems. We consider the approximation of the 2D elasticity problem by the Courant element. Concerning multilevel convergence rate, that is the γ corresponding to nested general triangular meshes of size h and 2h, we have proved that γ²≤ 3/4$ uniformly on the mesh and the Poisson ratio. Concerning error estimator, that is the γ corresponding to quadratic and linear approximations on the same mesh, numerical computations have shown that the exact γ for a reference element deteriorates that is goes to one, when the Poisson ratio tends to 1/2     

A posteriori FE error control for p-Laplacian by gradient recovery in quasi-norm

A posteriori error estimators based on quasi-norm gradient recovery are established for the finite element approximation of the p-Laplacian on unstructured meshes. The new a posteriori error estimators provide both upper and lower bounds in the quasi-norm for the discretization error. The main tools for the proofs of reliability are approximation error estimates for a local approximation operator in the quasi-norm.

     

Numerical analysis of the primal problem of elastoplasticity with hardening

Summary. The finite element method is a reasonable and frequently utilised tool for the spatial discretization within one time-step in an elastoplastic evolution problem. In this paper, we analyse the finite element discretization and prove a priori and a posteriori error estimates for variational inequalities corresponding to the primal formulation of (Hencky) plasticity. The finite element method of lowest order consists in minimising a convex function on a subspace of continuous piecewise linear resp. piecewise constant trial functions. An a priori error estimate is established for the fully-discrete method which shows linear convergence as the mesh-size tends to zero, provided the exact displacement field u is smooth. Near the boundary of the plastic domain, which is unknown a priori, it is most likely that u is non-smooth. In this situation, automatic mesh-refinement strategies are believed to improve the quality of the finite element approximation. We suggest such an adaptive algorithm on the basis of a computable a posteriori error estimate. This estimate is reliable and efficient in the sense that the quotient of the error by the estimate and its inverse are bounded from above. The constants depend on the hardening involved and become larger for decreasing hardening. Received May 7, 1997 / Revised version received August 31, 1998     

A posteriori error estimates in quantities of interest for the finite element heterogeneous multiscale method

We present an “a posteriori” error analysis in quantities of interest for elliptic homogenization problems discretized by the finite element heterogeneous multiscale method. The multiscale method is based on a macro‐to‐micro formulation, where the macroscopic physical problem is discretized in a macroscopic finite element space, and the missing macroscopic data are recovered on‐the‐fly using the solutions of corresponding microscopic problems. We propose a new framework that allows to follow the concept of the (single‐scale) dual‐weighted residual method at the macroscopic level in order to derive a posteriori error estimates in quantities of interests for multiscale problems. Local error indicators, derived in the macroscopic domain, can be used for adaptive goal‐oriented mesh refinement. These error indicators rely only on available macroscopic and microscopic solutions. We further provide a detailed analysis of the data approximation error, including the quadrature errors. Numerical experiments confirm the efficiency of the adaptive method and the effectivity of our error estimates in the quantities of interest. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013