An a Posteriori Error Estimator for Two-Body Contact Problems on Non-Matching Meshes

A posteriori error estimates for two-body contact problems are established. The discretization is based on mortar finite elements with dual Lagrange multipliers. To define locally the error estimator, Arnold–Winther elements for the stress and equilibrated fluxes for the surface traction are used. Using the Lagrange multiplier on the contact zone as Neumann boundary conditions, equilibrated fluxes can be locally computed. In terms of these fluxes, we define on each element a symmetric and globally H(div)-conforming approximation for the stress. Upper and lower bounds for the discretization error in the energy norm are provided. In contrast to many other approaches, the constant in the upper bound is, up to higher order terms, equal to one. Numerical examples illustrate the reliability and efficiency of the estimator. This work was supported in part by the Deutsche Forschungsgemeinschaft, SFB 404, B8.     

A Bi–Hyperbolic Finite Volume Method on Quadrilateral Meshes

A non-oscillatory, high resolution reconstruction method on quadrilateral meshes in two dimensions (2D) is presented. It is a two-dimensional extension of Marquina’s hyperbolic method. The generalization to quadrilateral meshes allows the method to simulate realistic flow problems in complex domains. An essential point in the construction of the method is a second order accurate approximation of gradients on an irregular, quadrilateral mesh. The resulting scheme is optimal in the sense that it is third order accurate and the reconstruction requires only nearest neighbour information. Numerical experiments are presented and the computational results are compared to experimental data.     

A Posteriori Error Estimates of Two-Grid Finite Element Methods for Nonlinear Elliptic Problems

In this article, we study the residual-based a posteriori error estimates of the two-grid finite element methods for the second order nonlinear elliptic boundary value problems. Computable upper and lower bounds on the error in the \(H^1\)-norm are established. Numerical experiments are also provided to illustrate the performance of the proposed estimators.     

Some a Posteriori Error Estimators for p-Laplacian Based on Residual Estimation or Gradient Recovery

In this paper, we first derive a posteriori error estimators of residual type for the finite element approximation of the p-Laplacian, and show that they are reliable, and efficient up to higher order terms. We then construct some a posteriori error estimators based on gradient recovery. We further compare the two types of a posteriori error estimators. It is found that there exist some relationships between the two types of estimators, which are similar to those held in the case of the Laplacian. It is shown that the a posteriori error estimators based on gradient recovery are equivalent to the discretization error in a quasi-norm provided the solution is sufficiently smooth and mesh is uniform. Under stronger conditions, superconvergnece properties have been established for the used gradient recovery operator, and then some of the gradient recovery based estimates are further shown to be asymptotically exact to the discretization error in a quasi-norm. Numerical results demonstrating these a posteriori estimates are also presented.     

Fully Computable a Posteriori Error Bounds for Hybridizable Discontinuous Galerkin Finite Element Approximations

We introduce an entirely new class of high-order methods for computational fluid dynamics based on the Gaussian process (GP) family of stochastic functions. Our approach is to use kernel-based GP prediction methods to interpolate/reconstruct high-order approximations for solving hyperbolic PDEs. We present a new high-order formulation to solve (magneto)hydrodynamic equations using the GP approach that furnishes an alternative to conventional polynomial-based approaches.     

基于Slicing的快速有限元网格体绘制算法

在Incremental Slicing算法的基础上,针对有限元网格的特点,提出了一种适合于多种有限元单元类型的快速体绘制方法,由于采用硬件加速的多边形ＲＧＢＡ填充方法混合显示数据场与几何体的外表面,优化了切面多边形的形成算法,并对于小单元进行了特殊处理,大大提高了有限元网格体绘制的速度与质量。     

A Posteriori Error Estimates for HDG Methods

We present the first a posteriori error analysis of the so-called hybridizable discontinuous Galerkin (HDG) methods for second-order elliptic problems. We show that the error in the flux can be controlled by only two terms. The first term captures the so-called data oscillation. The second solely depends on the difference between the trace of the scalar approximation and the corresponding numerical trace. Numerical experiments verifying the reliability and efficiency of the estimate in two-space dimensions are presented.     

A Posteriori Error Estimates for Weak Galerkin Finite Element Methods for Second Order Elliptic Problems

A residual type a posteriori error estimator is presented and analyzed for Weak Galerkin finite element methods for second order elliptic problems. The error estimator is proved to be efficient and reliable through two estimates, one from below and the other from above, in terms of an $H^1$ -equivalent norm for the exact error. Two numerical experiments are conducted to demonstrate the effectiveness of adaptive mesh refinement guided by this estimator.     

A Family of Finite Volume Schemes of Arbitrary Order on Rectangular Meshes

In this paper, we analyze vertex-centered finite volume method (FVM) of any order for elliptic equations on rectangular meshes. The novelty is a unified proof of the inf-sup condition, based on which, we show that the FVM approximation converges to the exact solution with the optimal rate in the energy norm. Furthermore, we discuss superconvergence property of the FVM solution. With the help of this superconvergence result, we find that the FVM solution also converges to the exact solution with the optimal rate in the $L^2$ -norm. Finally, we validate our theory with numerical experiments.     

A Posteriori Error Estimators for a Class of Variational Inequalities

In this paper, we present an a posteriori error analysis for the finite element approximation of a variational inequality. We derive a posteriori error estimators of residual type, which are shown to provide upper bounds on the discretization error for a class of variational inequalities provided the solutions are sufficiently regular. Furthermore we derive sharp a posteriori error estimators with both lower and upper error bounds for a subclass of the obstacle problem which are frequently met in many physical models. For sufficiently regular solutions, these estimates are shown to be equivalent to the discretization error in an energy type norm. Our numerical tests show that these sharp error estimators are both reliable and efficient in guiding mesh adaptivity for computing the free boundaries.     

A Posteriori Analysis of a Positive Streamwise Invariant Discretization of a Convection-Diffusion Equation

We consider the finite element discretization of a convection-diffusion equation, where the convection term is handled via a fluctuation splitting algorithm. We prove a posteriori error estimates which allow us to perform mesh adaptivity in order to optimize the discretization of these equations. Numerical results confirm the interest of such an approach.     

A Posteriori Error Estimates for Parabolic Variational Inequalities

We study a posteriori error estimates in the energy norm for some parabolic obstacle problems discretized with a Euler implicit time scheme combined with a finite element spatial approximation. We discuss the reliability and efficiency of the error indicators, as well as their localization properties. Apart from the obstacle resolution, the error indicators vanish in the so-called full contact set. The case when the obstacle is piecewise affine is studied before the general case. Numerical examples are given.     

A Posteriori Error Estimation for the Dirichlet Problem with Account of the Error in the Approximation of Boundary Conditions

H ¹, independently of the discretization method chosen. In particular, our error estimator can be applied also to problems and discretizations where the Galerkin orthogonality is not available. We will present different strategies for the evaluation of the error estimator. Only one constant appears in its definition which is the one from Friedrichs' inequality; that 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. Received April 15, 2002; revised March 10, 2003 Published online: June 23, 2003     

Star-Based a Posteriori Error Estimates for Elliptic Problems

We give an a posteriori error estimator for low order nonconforming finite element approximations of diffusion-reaction and Stokes problems, which relies on the solution of local problems on stars. It is proved to be equivalent to the energy error up to a data oscillation, without requiring Helmholtz decomposition of the error nor saturation assumption. Numerical experiments illustrate the good behavior and efficiency of this estimator for generic elliptic problems.     

Monotone Finite Volume Schemes for Diffusion Equation with Imperfect Interface on Distorted Meshes

In this paper we consider the numerical solution of stiff problems in which the eigenvalues are separated into two clusters, one containing the “stiff”, or fast, components and one containing the slow components, that is, there is a gap in their eigenvalue spectrum. By using exponential fitting techniques we develop a class of explicit Runge–Kutta methods, that we call stability fitted methods, for which the stability domain has two regions, one close to the origin and the other one fitting the large eigenvalues. We obtain the size of their stability regions as a function of the order and the fitting conditions. We also obtain conditions that the coefficients of these methods must satisfy to have a given stiff order for the Prothero–Robinson test equation. Finally, we construct an embedded pair of stability fitted methods of orders 2 and 1 and show its performance by means of several numerical experiments.     

Energy Norm shape A Posteriori Error Estimation for Mixed Discontinuous Galerkin Approximations of the Stokes Problem

In this paper, we develop the a posteriori error estimation of mixed discontinuous Galerkin finite element approximations of the Stokes problem. In particular, we derive computable upper bounds on the error, measured in terms of a natural (mesh-dependent) energy norm. This is done by rewriting the underlying method in a non-consistent form using appropriate lifting operators, and by employing a decomposition result for the discontinuous spaces. A series of numerical experiments highlighting the performance of the proposed a posteriori error estimator on adaptively refined meshes are presented.Paul Houston - Funded by the EPSRC (Grant GR/R76615). Thomas P. Wihler - Funded by the Swiss National Science Foundation (Grant PBEZ2-102321). This revised version was published online in July 2005 with corrected volume and issue numbers.     

Finite Volume Discretization and Multilevel Methods in Flow Problems

This article is intended as a preliminary report on the implementation of a finite volume multilevel scheme for the discretization of the incompressible Navier–Stokes equations. As is well known the use of staggered grids (e.g. MAC grids, Perić et al. Comput. Fluids, 16(4), 389–403, (1988)) is a serious impediment for the implementation of multilevel schemes in the context of finite differences. This difficulty is circumvented here by the use of a colocated finite volume discretization (Faure et al. (2004a) Submitted, Perić et al. Comput. Fluids, 16(4), 389–403, (1988)), for which the algebra of multilevel methods is much simpler than in the context of MAC type finite differences. The general ideas and the numerical simulations are presented in this article in the simplified context of a two-dimensional Burgers equations; the two-, and three-dimensional Navier–Stokes equations introducing new difficulties related to the incompressibility condition and the time discretization, will be considered elsewhere (see Faure et al. (2004a) Submitted and Faure et al. (2004b), in preparation).     

Reduced Basis Approximation and a Posteriori Error Estimation for Affinely Parametrized Elliptic Coercive Partial Differential Equations

In this paper we consider (hierarchical, La-grange)reduced basis approximation anda posteriori error estimation for linear functional outputs of affinely parametrized elliptic coercive partial differential equa-tions. The essential ingredients are (primal-dual)Galer-kin projection onto a low-dimensional space associated with a smooth “parametric manifold” - dimension re-duction; efficient and effective greedy sampling meth-ods for identification of optimal and numerically stable approximations - rapid convergence;a posteriori er-ror estimation procedures - rigorous and sharp bounds for the linear-functional outputs of interest; and Offine-Online computational decomposition strategies - min-imummarginal cost for high performance in the real-time/embedded (e.g., parameter-estimation, control)and many-query (e.g., design optimization, multi-model/ scale)contexts. We present illustrative results for heat conduction and convection-diffusion,inviscid flow, and linear elasticity; outputs include transport rates, added mass,and stress intensity factors. This work was supported by DARPA/AFOSR Grants FA9550-05-1-0114 and FA-9550-07-1-0425,the Singapore-MIT Alliance,the Pappalardo MIT Mechanical Engineering Graduate Monograph Fund,and the Progetto Roberto Rocca Politecnico di Milano-MIT.We acknowledge many helpful discussions with Professor Yvon Maday of University Paris6.