Hierarchical error estimates for the energy functional in obstacle problems

We present a hierarchical a posteriori error analysis for the minimum value of the energy functional in symmetric obstacle problems. The main result is that the error in the energy minimum is, up to oscillation terms, equivalent to an appropriate hierarchical estimator. The proof does not invoke any saturation assumption. We even show that small oscillation implies a related saturation assumption. In addition, we prove efficiency and reliability of an a posteriori estimate of the discretization error and thereby cast some light on the theoretical understanding of previous hierarchical estimators. Finally, we illustrate our theoretical results by numerical computations.     

一阶双曲问题间断有限元的后验误差分析

一阶双曲问题的有限元后验误差估计至今没有得到很好的解决.本文对d维区域上一阶双曲问题的k次间断有限元逼近提出了一种新的后验误差分析方法, 进而建立了间断有限元解在DG范数下(强于L₂范数)基于误差余量型的后验误差估计. 数值计算验证了本文理论分析的有效性. 本文方法也适用于其他变分问题有限元逼近的后验误差分析.     

A FE‐BE coupling for a fluid‐structure interaction problem: Hierarchical a posteriori error estimates

In this article, we establish a hierarchical a posteriori error estimate for a coupling of finite elements and boundary elements for a fluid‐structure interaction problem posed in two and three dimensions. These methods combine boundary elements for the exterior fluid and finite elements for the elastic structure. We consider two weak formulations, a nonsymmetric one and a symmetric one, which are both uniquely solvable. We present the reliability and efficiency of the error estimates. For the two dimensional case, we compute local error indicators which allow us to develop an adaptive mesh refinement strategy on triangles. For the three dimensional case, we use hexahedrons as elements. Numerical experiments underline our theoretical results. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2012     

High-order adaptive finite element-singly implicit Runge-Kutta methods for parabolic differential equations

We describe an adaptive mesh refinement finite element method-of-lines procedure for solving one-dimensional parabolic partial differential equations. Solutions are calculated using Galerkin's method with a piecewise hierarchical polynomial basis in space and singly implicit Runge-Kutta (SIRK) methods in time. A modified SIRK formulation eliminates a linear systems solution that is required by the traditional SIRK formulation and leads to a new reduced-order interpolation formula. Stability and temporal error estimation techniques allow acceptance of approximate solutions at intermediate stages, yielding increased efficiency when solving partial differential equations. A priori energy estimates of the local discretization error are obtained for a nonlinear scalar problem. A posteriori estimates of local spatial discretization errors, obtained by order variation, are used with the a priori error estimates to control the adaptive mesh refinement strategy. Computational results suggest convergence of the a posteriori error estimate to the exact discretization error and verify the utility of the adaptive technique.This research was partially supported by the U.S. Air Force Office of Scientific Research, Air Force Systems Command, USAF, under Grant Number AFOSR-90-0194; the U.S. Army Research Office under Contract Number DAAL 03-91-G-0215; by the National Science Foundation under Grant Number CDA-8805910; and by a grant from the Committee on Research, Tulane University.     

Estimation of nonbinary random response

The paper treats a new approach to reducing the dimension of factors which affect the non-binary response variable Y. This is relevant in analysis of a number of stochastic models, for instance, in biological and medical studies. The quality of Y estimation by means of a function in those factors is described by a specified error functional. It involves a penalty function to take into account the importance of the forecast for different response values. The joint distribution of factors and response variable is unknown. Thus it is quite natural to employ for statistical inference the estimates of the error functional constructed by prediction algorithm and cross-validation procedure. One of our main results provides the criterion of strong consistency of such estimates as the number of observations tends to infinity. Due to this result one can identify the significant factors. We introduce also the regularized versions of estimates and establish for them the central limit theorem (CLT). The statistical variant of our CLT permits to construct the approximate confidence intervals for unknown error functional.     

Analysis of the parareal approach based on discontinuous Galerkin method for time-dependent Stokes equations

This paper analyzes a parareal approach based on discontinuous Galerkin (DG) method for the time-dependent Stokes equations. A class of primal discontinuous Galerkin methods, namely variations of interior penalty methods, are adopted for the spatial discretization in the parareal algorithm (we call it parareal DG algorithm). We study three discontinuous Galerkin methods for the time-dependent Stokes equations, and the optimal continuous in time error estimates for the velocities and pressure are derived. Based on these error estimates, the proposed parareal DG algorithm is proved to be unconditionally stable and bounded by the error of discontinuous Galerkin discretization after a finite number of iterations. Finally, some numerical experiments are conducted which confirm our theoretical results, meanwhile, the efficiency of the parareal DG algorithm can be seen through a parallel experiment.     

Efficiency of a posteriori BEM--error estimates for first-kind integral equations on quasi--uniform meshes

In the numerical treatment of integral equations of the first kind using boundary element methods (BEM), the author and E. P. Stephan have derived a posteriori error estimates as tools for both reliable computation and self-adaptive mesh refinement. So far, efficiency of those a posteriori error estimates has been indicated by numerical examples in model situations only. This work affirms efficiency by proving the reverse inequality. Based on best approximation, on inverse inequalities and on stability of the discretization, and complementary to our previous work, an abstract approach yields a converse estimate. This estimate proves efficiency of an a posteriori error estimate in the BEM on quasi--uniform meshes for Symm's integral equation, for a hypersingular equation, and for a transmission problem.

     

On Mixed Error Estimates for Elliptic Obstacle Problems

We establish in this paper sharp error estimates of residual type for finite element approximation to elliptic obstacle problems. The estimates are of mixed nature, which are neither of a pure a priori form nor of a pure a posteriori form but instead they are combined by an a priori part and an a posteriori part. The key ingredient in our derivation for the mixed error estimates is the use of a new interpolator which enables us to eliminate inactive data from the error estimators. One application of our mixed error estimates is to construct a posteriori error indicators reliable and efficient up to higher order terms, and these indicators are useful in mesh-refinements and adaptive grid generations. In particular, by approximating the a priori part with some a posteriori quantities we can successfully track the free boundary for elliptic obstacle problems.     

Euler-type schemes for weakly coupled forward-backward stochastic differential equations and optimal convergence analysis

We introduce a new Euler-type scheme and its iterative algorithm for solving weakly coupled forward-backward stochastic differential equations (FBSDEs). Although the schemes share some common features with the ones proposed by C. Bender and J. Zhang [Ann. Appl. Probab., 2008, 18: 143–177], less computational work is needed for our method. For both our schemes and the ones proposed by Bender and Zhang, we rigorously obtain first-order error estimates, which improve the half-order error estimates of Bender and Zhang. Moreover, numerical tests are given to demonstrate the first-order accuracy of the schemes.     

AN ANISOTROPIC NONCONFORMING FINITE ELEMENT WITH SOME SUPERCONVERGENCE RESULTS

The main aim of this paper is to study the error estimates of a nonconforming finite element with some superconvergence results under anisotropic meshes. The anisotropic interpolation error and consistency error estimates are obtained by using some novel approaches and techniques, respectively. Furthermore, the superclose and a superconvergence estimate on the central points of elements are also obtained without the regularity assumption and quasi-uniform assumption requirement on the meshes. Finally, a numerical test is carried out, which coincides with our theoretical analysis.     

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     

Composite hierachical linear quantile regression

Multilevel （hierarchical） modeling is a generalization of linear and generalized linear modeling in which regression coefficients are modeled through a model, whose parameters are also estimated from data. Multilevel model fails to fit well typically by the use of the EM algorithm once one of level error variance （like Cauchy distribution） tends to infinity. This paper proposes a composite multilevel to combine the nested structure of multilevel data and the robustness of the composite quantile regression, which greatly improves the efficiency and precision of the estimation. The new approach, which is based on the Gauss-Seidel iteration and takes a full advantage of the composite quantile regression and multilevel models, still works well when the error variance tends to infinity, We show that even the error distribution is normal, the MSE of the estimation of composite multilevel quantile regression models nearly equals to mean regression. When the error distribution is not normal, our method still enjoys great advantages in terms of estimation efficiency.     

Some numerical quadrature schemes of a non-conforming quadrilateral finite element

Numerical quadrature schemes of a non-conforming finite element method for general second order elliptic problems in two dimensional (2-D) and three dimensional (3-D) space are discussed in this paper. We present and analyze some optimal numerical quadrature schemes. One of the schemes contains only three sampling points, which greatly improves the efficiency of numerical computations. The optimal error estimates are derived by using some traditional approaches and techniques. Lastly, some numerical results are provided to verify our theoretical analysis.     

Convergence and superconvergence analysis of a new quadratic Hermite-type triangular element on anisotropic meshes

A new quadratic Hermite-type triangular finite element is conceived to solve a class of two-dimensional second-order elliptic boundary value problems. Its error estimates on anisotropic meshes are developed. Furthermore, we verify that some conditions set to the meshes contribute to the proof of its superconvergence properties, which can improve the approximation results. Numerical examples are given to confirm our theoretical analysis.     

FINITE ELEMENT METHODS FOR SOBOLEV EQUATIONS

A new high-order time-stepping finite element method based upon the high-order numerical integration formula is formulated for Sobolev equations, whose computations consist of an iteration procedure coupled with a system of two elliptic equations. The optimal and superconvergence error estimates for this new method are derived both in space and in time. Also, a class of new error estimates of convergence and superconvergence for the time-continuous finite element method is demonstrated in which there are no time derivatives of the exact solution involved, such that these estimates can be bounded by the norms of the known data. Moreover, some useful a-posteriori error estimators are given on the basis of the superconvergence estimates.     

Finite element method for a nonsmooth elliptic equation

In this paper, we study the finite element method for a non-smooth elliptic equation. Error analysis is presented, including a priori and a posteriori error estimates as well as superconvergence analysis. We also propose two algorithms for solving the underlying equation. Numerical experiments are employed to confirm our error estimations and the efficiency of our algorithms.     

导弹最大射程Bayes估计的分析和改进

本文用验前数据的质量因子及估计的相对均方误差分析了导弹最大射程的一类Ｂａｙｅｓ估计的性能，对不同的质量因子，给出了最佳验前数据量的一种近似公式。针对这类Ｂａｙｅｓ估计的冒进问题，本文对它们进行了改进并得到了一类新的估计。最后，通过ＭｏｎｔｅＣａｒｌｏ法比较了这些估计的相对均方误差，验证了新估计的优良性。     

Preconditioning and a posteriori error estimates using h‐ and p‐hierarchical finite elements with rectangular supports

We show some of the properties of the algebraic multilevel iterative methods when the hierarchical bases of finite elements (FEs) with rectangular supports are used for solving the elliptic boundary value problems. In particular, we study two types of hierarchies; the so‐called h‐ and p‐hierarchical FE spaces on a two‐dimensional domain. We compute uniform estimates of the strengthened Cauchy–Bunyakowski–Schwarz inequality constants for these spaces. Moreover, for diagonal blocks of the stiffness matrices corresponding to the fine spaces, the optimal preconditioning matrices can be found, which have tri‐ or five‐diagonal forms for h‐ or p‐refinements, respectively, after a certain reordering of the elements. As another use of the hierarchical bases, the a posteriori error estimates can be computed. We compare them in test examples for h‐ and p‐hierarchical FEs with rectangular supports. Copyright © 2008 John Wiley & Sons, Ltd.     

On Goal Oriented Time Adaptivity using Local Error Estimates

We consider adaptive time discretization methods for ordinary differential equations where one aims to control the error in a quantity of interest of the form , with j : ℝ^d → ℝ. In this setting we propose a new timestep controller based on local error estimates of the quantity of interest. The new method converges when the tolerance goes to zero. We experimentally compare the new scheme with the classic norm-based time-adaptivity based on local error estimates as well as the dual-weighted residual (DWR) method. The results show significantly lower efficiency for the DWR method. The local error based schemes are similarly efficient, with the new scheme showing significant improvement in some cases. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A posteriori error estimates for mixed FEM in elasticity

A residue based reliable and efficient error estimator is established for finite element solutions of mixed boundary value problems in linear, planar elasticity. The proof of the reliability of the estimator is based on Helmholtz type decompositions of the error in the stress variable and a duality argument for the error in the displacements. The efficiency follows from inverse estimates. The constants in both estimates are independent of the Lamé constant , and so locking phenomena for are properly indicated. The analysis justifies a new adaptive algorithm for automatic mesh–refinement. Received July 17, 1997