An accurate a posteriori error estimator for the Steklov eigenvalue problem and its applications

In this paper, a type of accurate a posteriori error estimator is proposed for the Steklov eigenvalue problem based on the complementary approach, which provides an asymptotic exact estimate for the approximate eigenpair. Besides, we design a type of cascadic adaptive finite element method for the Steklov eigenvalue problem based on the proposed a posteriori error estimator. In this new cascadic adaptive scheme, instead of solving the Steklov eigenvalue problem in each adaptive space directly, we only need to do some smoothing steps for linearized boundary value problems on a series of adaptive spaces and solve some Steklov eigenvalue problems on a low dimensional space. Furthermore, the proposed a posteriori error estimator provides the way to refine mesh and control the number of smoothing steps for the cascadic adaptive method. Some numerical examples are presented to validate the efficiency of the algorithm in this paper.

     

A POSTERIORI ERROR ESTIMATE OF OPTIMAL CONTROL PROBLEM OF PDE WITH INTEGRAL CONSTRAINT FOR STATE

In this paper, we study adaptive finite element discretization schemes for an optimal control problem governed by elliptic PDE with an integral constraint for the state. We derive the equivalent a posteriori error estimator for the finite element approximation, which particularly suits adaptive multi-meshes to capture different singularities of the control and the state. Numerical examples are presented to demonstrate the efficiency of a posteriori error estimator and to confirm the theoretical results.     

A posteriori error estimate for discontinuous Galerkin finite element method on polytopal mesh

In this work, we derive a posteriori error estimates for discontinuous Galerkin finite element method on polytopal mesh. We construct a reliable and efficient a posteriori error estimator on general polygonal or polyhedral meshes. An adaptive algorithm based on the error estimator and DG method is proposed to solve a variety of test problems. Numerical experiments are performed to illustrate the effectiveness of the algorithm.     

稳态Poisson-Nernst-Planck方程的残量型后验误差估计

针对稳态的Poisson-Nernst-Planck方程研究了一种残量型的后验误差估计子,对方程的两个解-浓度和电势,都分别给出了上界和下界估计.数值实验表明,基于这种后验误差估计子构造的自适应有限元算法对于稳态的Poisson-Nernst-Planck方程是有效的.     

Steklov eigenvalue problem and its applications An accurate a posteriori error estimator for the

In this paper, a type of accurate a posteriori error estimator is proposed for the Steklov eigenvalue problem based on the complementary approach, which provides an asymptotic exact estimate for the approximate eigenpair. Besides, we design a type of cascadic adaptive finite element method for the Steklov eigenvalue problem based on the proposed a posteriori error estimator. In this new cascadic adaptive scheme, instead of solving the Steklov eigenvalue problem in each adaptive space directly, we only need to do some smoothing steps for linearized boundary value problems on a series of adaptive spaces and solve some Steklov eigenvalue problems on a low dimensional space. Furthermore, the proposed a posteriori error estimator provides the way to refine meshes and control the number of smoothing steps for the cascadic adaptive method. Some numerical examples are presented to validate the efficiency of the algorithm in this paper.     

An adaptive finite element method for the sparse optimal control of fractional diffusion

We propose and analyze an a posteriori error estimator for a partial differential equation (PDE)-constrained optimization problem involving a nondifferentiable cost functional, fractional diffusion, and control-constraints. We realize fractional diffusion as the Dirichlet-to-Neumann map for a nonuniformly PDE and propose an equivalent optimal control problem with a local state equation. For such an equivalent problem, we design an a posteriori error estimator which can be defined as the sum of four contributions: two contributions related to the approximation of the state and adjoint equations and two contributions that account for the discretization of the control variable and its associated subgradient. The contributions related to the discretization of the state and adjoint equations rely on anisotropic error estimators in weighted Sobolev spaces. We prove that the proposed a posteriori error estimator is locally efficient and, under suitable assumptions, reliable. We design an adaptive scheme that yields, for the examples that we perform, optimal experimental rates of convergence.     

AN ADAPTIVE INVERSE ITERATION FEM FOR THE INHOMOGENEOUS DIELECTRIC WAVEGUIDES

We introduce an adaptive finite element method for computing electromagnetic guided waves in a closed, inhomogeneous, pillared three-dimensional waveguide at a given frequency based on the inverse iteration method. The problem is formulated as a generalized eigenvalue problems. By modifying the exact inverse iteration algorithm for the eigenvalue problem, we design a new adaptive inverse iteration finite element algorithm. Adaptive finite element methods based on a posteriori error estimate are known to be successful in resolving singularities of eigenfunctions which deteriorate the finite element convergence. We construct a posteriori error estimator for the electromagnetic guided waves problem. Numerical results are reported to illustrate the quasi-optimal performance of our adaptive inverse iteration finite element method.     

Duality Based A Posteriori Error Estimator for the Dirichlet Problem

The present work is devoted to the a posteriori error estimation for 2nd order elliptic problems with Dirichlet boundary conditions. Using the duality technique we derive a reliable and efficient a posteriori error estimator that measures the error in the energy norm. All the derivations are done on continuous level, and the estimator can be used in assessing the error of any approximate solution which belongs to the Sobolev space H¹, independently of the discretization method chosen. In particular, we make no use of the Galerkin orthogonality, which enables us to implement the estimator for measuring the error of the fictitious domain/penalty finite element method. The estimator is easily computable, and the only constant present in the estimator is the one from Friedrichs' inequality; the constant depends solely on the domain geometry, and the estimator is quite non‐sensitive to the error in the constant evaluation. Finally, we show how accurately the estimator captures the local error distribution, thus, creating a base for a justified adaptivity of an approximation.     

A posterior error analysis for the nonconforming discretization of Stokes eigenvalue problem

In this paper, we present a posteriori error estimator for the nonconforming finite element approximation, including using Crouzeix–Raviart element and extended Crouzeix–Raviart element, of the Stokes eigenvalue problem. With the technique of Helmholtz decomposition, we first give out a posteriori error estimator and prove it as the global upper bound and local lower bound of the approximation error. Then, by deleting a jump term in the indicator, another simpler but equivalent indicator is obtained. Some numerical experiments are provided to verify our analysis.     

Residual error estimate for BEM-based FEM on polygonal meshes

A conforming finite element method on polygonal meshes is reviewed which handles hanging nodes naturally. Trial functions are defined to fulfil the homogeneous PDE locally and they are treated by means of local boundary integral equations. Using a quasi-interpolation operator of Clément type a residual-based error estimate is obtained. This a posteriori estimator can be used to rate the accuracy of the approximation over polygonal elements or it can be applied to an adaptive BEM-based FEM. The numerical experiments confirm our results and show optimal convergence for the adaptive strategy on general meshes.     

Adaptive Hybridized Interior Penalty Discontinuous Galerkin Methods for H(curl)-Elliptic Problems

We develop and analyze an adaptive hybridized Interior Penalty Discontinuous Galerkin (IPDG-H) method for H(curl)-elliptic boundary value problems in 2D or 3D arising from a semi-discretization of the eddy currents equations. The method can be derived from a mixed formulation of the given boundary value problem and involves a Lagrange multiplier that is an approximation of the tangential traces of the primal variable on the interfaces of the underlying triangulation of the computational domain. It is shown that the IPDG-H technique can be equivalently formulated and thus implemented as a mortar method. The mesh adaptation is based on a residual-type a posteriori error estimator consisting of element and face residuals. Within a unified framework for adaptive finite element methods, we prove the reliability of the estimator up to a consistency error. The performance of the adaptive symmetric IPDG-H method is documented by numerical results for representative test examples in 2D.     

Superconvergence analysis and a posteriori error estimation of a finite element method for an optimal control problem governed by integral equations

In this paper, we discuss the numerical simulation for a class of constrained optimal control problems governed by integral equations. The Galerkin method is used for the approximation of the problem. A priori error estimates and a superconvergence analysis for the approximation scheme are presented. Based on the results of the superconvergence analysis, a recovery type a posteriori error estimator is provided, which can be used for adaptive mesh refinement. The research project is supported by the National Basic Research Program under the Grant 2005CB321701 and the National Natural Science Foundation of China under the Grant 10771211.     

Low cost a posteriori error estimators for an augmented mixed FEM in linear elasticity

We consider an augmented mixed finite element method applied to the linear elasticity problem and derive a posteriori error estimators that are simpler and easier to implement than the ones available in the literature. In the case of homogeneous Dirichlet boundary conditions, the new a posteriori error estimator is reliable and locally efficient, whereas for non-homogeneous Dirichlet boundary conditions, we derive an a posteriori error estimator that is reliable and satisfies a quasi-efficiency bound. Numerical experiments illustrate the performance of the corresponding adaptive algorithms and support the theoretical results.     

An adaptive least‐squares mixed finite element method for the Signorini problem

We present and analyze a least squares formulation for contact problems in linear elasticity which employs both, displacements and stresses, as independent variables. As a consequence, we obtain stability and high accuracy of our discretization also in the incompressible limit. Moreover, our formulation gives rise to a reliable and efficient a posteriori error estimator. To incorporate the contact constraints, the first‐order system least squares functional is augmented by a contact boundary functional which implements the associated complementarity condition. The bilinear form related to the augmented functional is shown to be coercive and therefore constitutes an upper bound, up to a constant, for the error in displacements and stresses in . This implies the reliability of the functional to be used as an a posteriori error estimator in an adaptive framework. The efficiency of the use of the functional as an a posteriori error estimator is monitored by the local proportion of the boundary functional term with respect to the overall functional. Computational results using standard conforming linear finite elements for the displacement approximation combined with lowest‐order Raviart‐Thomas elements for the stress tensor show the effectiveness of our approach in an adaptive framework for two‐dimensional and three‐dimensional Hertzian contact problems. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 276–289, 2017     

An adaptive strategy for elliptic problems including a posteriori controlled boundary approximation

We derive a posteriori error estimates for the approximation of linear elliptic problems on domains with piecewise smooth boundary. The numerical solution is assumed to be defined on a Finite Element mesh, whose boundary vertices are located on the boundary of the continuous problem. No assumption is made on a geometrically fitting shape.

A posteriori error estimates are given in the energy norm and the -norm, and efficiency of the adaptive algorithm is proved in the case of a saturated boundary approximation. Furthermore, a strategy is presented to compute the effect of the non-discretized part of the domain on the error starting from a coarse mesh. This especially implies that parts of the domain, where the measured error is small, stay non-discretized. The presented algorithm includes a stable path following to supply a sufficient polygonal approximation of the boundary, the reliable computation of the a posteriori estimates and a mesh adaptation based on Delaunay techniques. Numerical examples illustrate that errors outside the initial discretization will be detected.

     

An adaptive weak Galerkin finite element method with hierarchical bases for the elliptic problem

Based on a posteriori error estimator with hierarchical bases, an adaptive weak Galerkin finite element method (WGFEM) is proposed for the elliptic problem with mixed boundary conditions. For the posteriori error estimator, we are only required to solve a linear algebraic system with diagonal entries corresponding to the degree of freedoms, which significantly reduces the computational cost. The upper and lower bounds of the error estimator are shown to addresses the reliability and efficiency of the adaptive approach. Numerical simulations are provided to demonstrate the effectiveness and robustness of the proposed method.     

An a posteriori error estimator for anisotropic refinement

Summary. Besides an algorithm for local refinement, an a posteriori error estimator is the basic tool of every adaptive finite element method. Using information generated by such an error estimator the refinement of the grid is controlled. For 2nd order elliptic problems we present an error estimator for anisotropically refined grids, like -d cuboidal and 3-d prismatic grids, that gives correct information about the size of the error; additionally it generates information about the direction into which some element has to be refined to reduce the error in a proper way. Numerical examples are presented for 2-d rectangular and 3-d prismatic grids. Received March 15, 1994 / Revised version received June 3, 1994     

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.     

The backward euler anisotropic a posteriori error analysis for parabolic integro‐differential equations

We derive two optimal a posteriori error estimators for an implicit fully discrete approximation to the solutions of linear integro‐differential equations of the parabolic type. A continuous, piecewise linear finite element space is used for the space discretization and the time discretization is based on an implicit backward Euler method. The a posteriori error indicator corresponding to space discretization is derived using the anisotropic interpolation estimates in conjunction with a Zienkiewicz‐Zhu error estimator to approach the error gradient. The error due to time discretization is derived using continuous, piecewise linear polynomial in time. We use the linear approximation of the Volterra integral term to estimate the quadrature error in the second estimator. Numerical experiments are performed on the isotropic mesh to validate the derived results.© 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1309–1330, 2016     

Residual-based a posteriori error estimate for a mixed Reißner-Mindlin plate finite element method

Reliable and efficient residual-based a posteriori error estimates are established for the stabilised locking-free finite element methods for the Reissner-Mindlin plate model. The error is estimated by a computable error estimator from above and below up to multiplicative constants that do neither depend on the mesh-size nor on the plate's thickness and are uniform for a wide range of stabilisation parameter. The error is controlled in norms that are known to converge to zero in a quasi-optimal way. An adaptive algorithm is suggested and run for improving the convergence rates in three numerical examples for thicknesses 0.1, .001 and .001.