Convergence on layer-adapted meshes and anisotropic interpolation error estimates of non-standard higher order finite elements

For a general class of finite element spaces based on local polynomial spaces E with Pp⊂E⊂Qp we construct a vertex-edge-cell and point-value oriented interpolation operators that fulfil anisotropic interpolation error estimates.Using these estimates we prove ε-uniform convergence of order p for the Galerkin FEM and the LPSFEM for a singularly perturbed convection-diffusion problem with characteristic boundary layers.     

An aspect ratio agglomeration multigrid for unstructured grids

A robust aspect ratio‐based agglomeration algorithm to generate high quality of coarse grids for unstructured and hybrid grids is proposed in this paper. The algorithm focuses on multigrid techniques for the numerical solution of Euler and Navier–Stokes equations, which conform to cell‐centered finite volume special discretization scheme, combines vertex‐based isotropic agglomeration and cell‐based directional agglomeration to yield large increases in convergence rates. Aspect ratio is used as fusing weight to capture the degree of cell convexity and give an indication of cell stretching. Agglomeration front queue is established to propagate inward from the boundaries, which stores isotropic vertex and also high‐stretched cell marked with different flag according to aspect ratio. We conduct the present method to solve Euler and Navier–Stokes equations on unstructured and hybrid grids and compare the results with single grid as well as MGridGen, which shows that the present method is efficient in reducing computational time for large‐scale system equations. Copyright © 2013 John Wiley & Sons, Ltd.     

An accurate gradient and Hessian reconstruction method for cell‐centered finite volume discretizations on general unstructured grids

In this paper, a novel reconstruction of the gradient and Hessian tensors on an arbitrary unstructured grid, developed for implementation in a cell‐centered finite volume framework, is presented. The reconstruction, based on the application of Gauss' theorem, provides a fully second‐order accurate estimate of the gradient, along with a first‐order estimate of the Hessian tensor. The reconstruction is implemented through the construction of coefficient matrices for the gradient components and independent components of the Hessian tensor, resulting in a linear system for the gradient and Hessian fields, which may be solved to an arbitrary precision by employing one of the many methods available for the efficient inversion of large sparse matrices. Numerical experiments are conducted to demonstrate the accuracy, robustness, and computational efficiency of the reconstruction by comparison with other common methods. Copyright © 2009 John Wiley & Sons, Ltd.     

Hybrid mesh adaptation applied to industrial numerical combustion

This paper presents an anisotropic mesh adaptation method applied to industrial combustion problems. The method is based on a measure of the distance between two Riemannian metrics called metric non‐conformity. This measure, which can be used to build a cost function to adapt meshes comprising several types of mesh elements, provides the basis for a generic mesh adaptation approach applicable to various types of physical problems governed by partial differential equations. The approach is shown to be applicable to industrial combustion problems, through the specification of a target metric computed as the intersection of several Hessian matrices reconstructed from the main variables of the governing equations. Numerical results show that the approach is cost effective in that it can drastically improve the prediction of temperature and species distributions in the flame region of a combustor while reducing computational cost. The results can be used as a basis for pollutant prediction models. Copyright © 2011 John Wiley & Sons, Ltd.     

Robust stabilization of uncertain linear dynamical systems with time-varying delay

In this paper, the problem of the robust stabilization for a class of uncertain linear dynamical systems with time-varying delay is considered. By making use of an algebraic Riccati equation, we derive some sufficient conditions for robust stability of time-varying delay dynamical systems with unstructured or structured uncertainties. In our approach, the only restriction on the delay functionh(t) is the knowledge of its upper boundh ^–. Some analytical methods are employed to investigate these stability conditions. Since these conditions are independent of the delay, our results are also applicable to systems with perturbed time delay. Finally, a numerical example is given to illustrate the use of the sufficient conditions developed in this paper.     

On the adjacencies of triangular meshes based on skeleton-regular partitions

For any 2D triangulation τ, the 1-skeleton mesh of τ is the wireframe mesh defined by the edges of τ, while that for any 3D triangulation τ, the 1-skeleton and the 2-skeleton meshes, respectively, correspond to the wireframe mesh formed by the edges of τ and the “surface” mesh defined by the triangular faces of τ. A skeleton-regular partition of a triangle or a tetrahedra, is a partition that globally applied over each element of a conforming mesh (where the intersection of adjacent elements is a vertex or a common face, or a common edge) produce both a refined conforming mesh and refined and conforming skeleton meshes. Such a partition divides all the edges (and all the faces) of an individual element in the same number of edges (faces). We prove that sequences of meshes constructed by applying a skeleton-regular partition over each element of the preceding mesh have an associated set of difference equations which relate the number of elements, faces, edges and vertices of the nth and (n−1)th meshes. By using these constitutive difference equations we prove that asymptotically the average number of adjacencies over these meshes (number of triangles by node and number of tetrahedra by vertex) is constant when n goes to infinity. We relate these results with the non-degeneracy properties of longest-edge based partitions in 2D and include empirical results which support the conjecture that analogous results hold in 3D.     

Domain decomposition preconditioners of Neumann-Neumann type for hp-approximations on boundary layer meshes in three dimensions

We develop and analyse Neumann–Neumann methods for hpfinite-element approximations of scalar elliptic problems ongeometrically refined boundary layer meshes in three dimensions.These are meshes that are highly anisotropic where the aspectratio typically grows exponentially with the polynomial degree.The condition number of our preconditioners is shown to be independentof the aspect ratio of the mesh and of potentially large jumpsof the coefficients. In addition, it only grows polylogarithmicallywith the polynomial degree, as in the case of p approximationson shape-regular meshes. This work generalizes our previousone on two-dimensional problems in Toselli & Vasseur (2003a,submitted to Numerische Mathematik, 2003c to appear in Comput.Methods Appl. Mech. Engng.) and the estimates derived here canbe employed to prove condition number bounds for certain typesof FETI methods.     

Locking-free DGFEM for elasticity problems in polygons

The h-version of the discontinuous Galerkin finite element method(h-DGFEM) for nearly incompressible linear elasticity problemsin polygons is analysed. It is proved that the scheme is robust(locking-free) with respect to volume locking, even in the absenceof H²-regularity of the solution. Furthermore, it is shown thatan appropriate choice of the finite element meshes leads torobust and optimal algebraic convergence rates of the DGFEMeven if the exact solutions do not belong to H².     

Superconvergence of a Nonconforming Finite Element Approximation to Viscoelasticity Type Equations on Anisotropic Meshes

The main aim of this paper is to study the approximation to viscoelasticity type equations with a Crouzeix-Raviart type nonconforming finite element on the anisotropic meshes. The superclose property of the exact solution and the optimal error estimate of its derivative with respect to time are derived by using some novel techniques. Moreover, employing a postprocessing technique, the global superconvergence property for the discretization error of the postprocessed discrete solution to the solution itself is studied.