A posteriori error estimates of discontinuous Galerkin methods for non-standard Volterra integro-differential equations

^** Email: jingtang{at}lsec.cc.ac.cn^*** Email: hermann{at}math.mun.ca In this paper we establish a posteriori error estimates forthe discontinuous Galerkin (DG) method applied to linear, semilinearand non-standard (non-linear) Volterra integro-differentialequations. We also present an analysis of the DG method withquadrature for the memory term. Numerical experiments basedon three integro-differential equations are used to illustratevarious aspects of the error analysis.     

Guaranteed error bounds for finite element approximations of noncoercive elliptic problems and their applications

In this paper, we discuss with guaranteed a priori and a posteriori error estimates of finite element approximations for not necessarily coercive linear second order Dirichlet problems. Here, ‘guaranteed’ means we can get the error bounds in which all constants included are explicitly given or represented as a numerically computable form. Using the invertibility condition of concerning elliptic operator, guaranteed a priori and a posteriori error estimates are formulated. This kind of estimates plays essential and important roles in the numerical verification of solutions for nonlinear elliptic problems. Several numerical examples that confirm the actual effectiveness of the method are presented.     

A stabilized finite-element method for the Stokes problem including element and edge residuals

^** Email: raraya{at}ing-mat.udec.cl^*** Corresponding author. Email: gbarrene{at}ing-mat.udec.cl^**** Email: valentin{at}lncc.br A new stabilized finite-element method is presented for theStokes problem. The method is of a Douglas–Wang type,and includes a positive jump term controlling the residual ofthe Cauchy stress tensor on the internal edges of the triangulation.A priori error estimates are obtained in the natural norms ofthe unknowns and an a posteriori error estimator is proposed,analysed and tested through numerical experiments.     

Linear estimates of accuracy for approximate solutions of inverse problems

In this paper, we explore the question of which non-linear inverse problems, which are solved by a selected regularization method, may have so-called linear a priori accuracy estimates – that is, the accuracy of corresponding approximate solutions linearly depends on the error level of the data. In particular, we prove that if such a linear estimate exists, then the inverse problem under consideration is well posed, according to Tikhonov. For linear inverse problems, we find that the existence of linear estimates lead to, under some assumptions, the well-posedness (according to Tikhonov) on the whole space of solutions. Moreover, we consider a method for solving inverse problems with guaranteed linear estimates, called the residual method on the correctness set (RMCS). The linear a priori estimates of absolute and relative accuracy for the RMCS are presented, as well as analogous a posteriori estimates. A numerical illustration of obtaining linear a priori estimates for appropriate parametric sets of solutions using RMCS is given in comparison with Tikhonov regularization. The a posteriori estimates are calculated on these parametric sets as well.     

A posteriori error estimates for finite element approximations to the wave equation with discontinuous coefficients

We derive residual‐based a posteriori error estimates of finite element method for linear wave equation with discontinuous coefficients in a two‐dimensional convex polygonal domain. A posteriori error estimates for both the space‐discrete case and for implicit fully discrete scheme are discussed in L^∞(L²) norm. The main ingredients used in deriving a posteriori estimates are new Clément type interpolation estimates in conjunction with appropriate adaption of the elliptic reconstruction technique of continuous and discrete solutions. We use only an energy argument to establish a posteriori error estimates with optimal order convergence in the L^∞(L²) norm.     

A posteriori error estimates of edge residual type of finite element method for nonmonotone quasi‐linear elliptic problems

In this article, we study the edge residual‐based a posteriori error estimates of conforming linear finite element method for nonmonotone quasi‐linear elliptic problems. It is proven that edge residuals dominate a posteriori error estimates. Up to higher order perturbations, edge residuals can act as a posteriori error estimators. The global reliability and local efficiency bounds are established both in H ¹‐norm and L ²‐norm. Numerical experiments are provided to illustrate the performance of the proposed error estimators. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 813–837, 2014     

Functional a posteriori error estimates for the reaction-convection-diffusion problem

In this paper, a general form of functional type a posteriori error estimates for linear reaction-convection-diffusion problems is presented. It is derived by purely functional arguments without attracting specific properties of the approximation method. The estimate provides a guaranteed upper bound of the difference between the exact solution and any conforming approximation from the energy functional class. It is also proved that the derived error majorants give computable quantities, which are equivalent to the error evaluated in the energy and combined primal-dual norms. Bibliography: 14 titles. Published in Zapiski Nauchnykh Seminarov POMI, Vol. 348, 2007, pp. 127–146.     

A convergent adaptive finite element method for an optimal design problem

The optimal design problem for maximal torsion stiffness of an infinite bar of given geometry and unknown distribution of two materials of prescribed amounts is one model example in topology optimisation. It eventually leads to a degenerate convex minimisation problem. The numerical analysis is therefore delicate for possibly multiple primal variables u but unique derivatives σ : = DW(D u). Even fine a posteriori error estimates still suffer from the reliability-efficiency gap. However, it motivates a simple edge-based adaptive mesh-refining algorithm (AFEM) that is not a priori guaranteed to refine everywhere. Its convergence proof is therefore based on energy estimates and some refined convexity control. Numerical experiments illustrate even nearly optimal convergence rates of the proposed AFEM. Supported by the DFG Research Center MATHEON “Mathematics for key technologies” in Berlin.     

Error analysis of discontinuous Galerkin finite element method for optimal control problem governed by the transport equation

This article discusses a priori and a posteriori error estimates of discontinuous Galerkin finite element method for optimal control problem governed by the transport equation. We use variational discretization concept to discretize the control variable and discontinuous piecewise linear finite elements to approximate the state and costate variable. Based on the error estimates of discontinuous Galerkin finite element method for the transport equation, we get a priori and a posteriori error estimates for the transport equation optimal control problem. Finally, two numerical experiments are carried out to confirm the theoretical analysis.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1493–1512, 2017     

A posteriori error estimates for the Crank-Nicolson method for parabolic equations

We derive optimal order a posteriori error estimates for time discretizations by both the Crank-Nicolson and the Crank-Nicolson-Galerkin methods for linear and nonlinear parabolic equations. We examine both smooth and rough initial data. Our basic tool for deriving a posteriori estimates are second-order Crank-Nicolson reconstructions of the piecewise linear approximate solutions. These functions satisfy two fundamental properties: (i) they are explicitly computable and thus their difference to the numerical solution is controlled a posteriori, and (ii) they lead to optimal order residuals as well as to appropriate pointwise representations of the error equation of the same form as the underlying evolution equation. The resulting estimators are shown to be of optimal order by deriving upper and lower bounds for them depending only on the discretization parameters and the data of our problem. As a consequence we provide alternative proofs for known a priori rates of convergence for the Crank-Nicolson method.

     

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     

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)     

A discontinuous Galerkin method for the Rosenau equation

A priori error estimates for the Rosenau equation, which is a K-dV like Rosenau equation modelled to describe the dynamics of dense discrete systems, have been studied by one of the authors. But since a priori error bounds contain the unknown solution and its derivatives, it is not effective to control error bounds with only a given step size. Thus we need to estimate a posteriori errors in order to control accuracy of approximate solutions using variable step sizes. A posteriori error estimates of the Rosenau equation are obtained by a discontinuous Galerkin method and the stability analysis is discussed for the dual problem. Numerical results on a posteriori error and wave propagation are given, which are obtained by using various spatial and temporal meshes controlled automatically by a posteriori error.     

A reduced local C0 discontinuous galerkin method for Kirchhoff plates

We propose and analyze a reduced local C⁰ discontinuous Galerkin (reduced LCDG) method with minimal penalization for Kirchhoff plate bending problems. The resulting linear system of the method can be solved efficiently by Gaussian elimination. Based on the key observation that the reduced LCDG method can be viewed as the localization of Hellan–Herrmann–Johnson method, we are inspired to use the techniques for dealing with the latter method to derive the well‐posedness and a priori error estimates of the reduced LCDG method. With the help of Zienkiewicz–Guzmán–Neilan element space, the a posteriori error analysis is also developed. Some numerical results are provided to demonstrate the theoretical estimates.© 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1902–1930, 2014     

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     

A posteriori error estimation of the p‐Laplace problem

Two a posteriori error estimates are discussed for the p‐Laplace problem. Up to errors in their numerical computation, they provide a guaranteed upper bound for the W^1,p‐seminorm and a weighted W^1,2‐seminorm of u – u_h. The first, sharper a posteriori estimate is based on the numerical solution of local interface problems. The standard residual‐based error estimate is addressed with emphasis on involved constants, determined as local eigenvalues. Numerical examples that illustrate the performance of these estimators can be found in [3].     

A posteriori error estimation for nonlinear variational problems

We introduce a new modus operandi for a posteriori error estimation for nonlinear (and linear) variational problems based on the duality theory of the calculus of variations. We derive what we call duality error estimates and show that they yield computable a posteriori error estimates without directly solving the dual problem.     

A mortar element method for elliptic problems with discontinuous coefficients

This paper proposes a mortar finite element method for solvingthe two-dimensional second-order elliptic problem with jumpsin coefficients across the interface between two subregions.Non-matching finite element grids are allowed on the interface,so independent triangulations can be used in different subregions.Explicitly realizable mortar conditions are introduced to couplethe individual discretizations. The same optimal L²-norm andenergy-norm error estimates as for regular problems are achievedwhen the interface is of arbitrary shape but smooth, thoughthe regularity of the true solution is low in the whole physicaldomain.     

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

定常Navier-Stokes方程流函数形式两重网格算法的残量型后验误差估计

运用七种两重网格协调元方法得出了不可压Navier-Stokes方程流函数形式的残量型后验误差估计．对比标准有限元方法的后验误差估计,两重网格算法的后验误差估计多了一些额外项(三线性项)．说明了这些额外项在误差估计中对研究离散解渐近性的重要性,推出了对于最优网格尺寸,这些额外项的收敛阶不高于标准离散解的收敛阶．