Stokes方程的一个新的非协调四边形单元格式

对于Stokes方程给出了一个新的非协调四边形单元格式.新单元具有构造简单,自由度较少等优势.特别指出的是,该单元在矩形网格下,还是一个Locking-free元,可用于平面弹性问题.尽管该单元不含协调部分,其相容误差估计较困难,通过采用新的技巧和方法得到了最优误差估计.     

A new a priori error estimate of nonconforming finite element methods

This paper is devoted to a new error analysis of nonconforming finite element methods.Compared with the classic error analysis in literature,only weak continuity,the F-E-M-Test for nonconforming finite element spaces,and basic Hm regularity for exact solutions of 2m-th order elliptic problems under consideration are assumed.The analysis is motivated by ideas from a posteriori error estimates and projection average operators.One main ingredient is a novel decomposition for some key average terms on(n.1)-dimensional faces by introducing a piecewise constant projection,which defines the generalization to more general nonconforming finite elements of the results in literature.The analysis and results herein are conjectured to apply for all nonconforming finite elements in literature.     

A quadrilateral nonconforming finite element for linear elasticity problem

In this paper, a four-parameter quadrilateral nonconforming finite element with DSP (double set parameters) is presented. Then we discuss the quadrilateral nonconforming finite element approximation to the linear elastic equations with pure displacement boundary. The optimal convergence rate of the method is established in the broken energy and -norms, and in particular, the convergence is uniform with respect to the Lamé parameter . Also the performance of the scheme does not deteriorate as the material becomes nearly incompressible. Lastly, a numerical test is carried out, which coincides with our theoretical analysis. The research is supported by NSF of China (No. 10471133) and the project of Creative Engineering of Henan Province of China.     

A mixed finite element for the stokes problem using quadrilateral elements

The object of this paper is to complete the results obtained in [3] by showing that the new mixed finite element that we have constructed in [3] also works for quadrilateral elements and to compare this method with the standard finite volume method. Estimates of optimal order are derived for both the new mixed finite element and an associated finite volume method.     

DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD FOR A FORWARD-BACKWARD HEAT EQUATION

A space-time finite element method,discontinuous in time but continuous in space, is studied to solve the nonlinear forward-backward heat equation. A linearized technique is introduced in order to obtain the error estimates of the approximate solutions. And the numerical simulations are given.     

Mixed finite element methods for general quadrilateral grids

We study a new mixed finite element of lowest order for general quadrilateral grids which gives optimal order error in the H(div)-norm. This new element is designed so that the H(div)-projection Π_h satisfies ∇ · Π_h = P_hdiv. A rigorous optimal order error estimate is carried out by proving a modified version of the Bramble-Hilbert lemma for vector variables. We show that a local H(div)-projection reproducing certain polynomials suffices to yield an optimal L²-error estimate for the velocity and hence our approach also provides an improved error estimate for original Raviart-Thomas element of lowest order. Numerical experiments are presented to verify our theory.     

Numerical integration over polygons using an eight-node quadrilateral spline finite element

In this paper, a cubature formula over polygons is proposed and analysed. It is based on an eight-node quadrilateral spline finite element [C.-J. Li, R.-H. Wang, A new 8-node quadrilateral spline finite element, J. Comp. Appl. Math. 195 (2006) 54–65] and is exact for quadratic polynomials on arbitrary convex quadrangulations and for cubic polynomials on rectangular partitions. The convergence of sequences of the above cubatures is proved for continuous integrand functions and error bounds are derived. Some numerical examples are given, by comparisons with other known cubatures.     

Technical Note: Reduction of rational integrals related to linear and convex quadrilateral finite elements

We reduce a rational function of bivariate nth degree polynomial numerator with a linear denominator to a simple bivariate polynomial of degree (n ? 1) and a rational function of a single variate nth degree polynomial numerator with the same bivariate linear denominator. This has very greatly contributed to the evaluation of (n + 1)(n + 2)/2 rational integrals in bivariates to mere (n + 1) rational integral of a single variate and an integration of simple polynomial in bivariates. Thus the effort of integration is reduced several times and leads to simple analytical expressions in terms of the nodal coordinates. In order to illustrate the numerical process two examples are considered. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 759–770, 2002; Published online in Wiley InterScience (www.interscience.wiley.com); DOI 10.1002/num.10026.     

A priori and a posteriori analysis of non-conforming finite elements with face penalty for advection-diffusion equations

^** Corresponding author. Email: l.elalaoui{at}imperial.ac.uk^*** Email: ern{at}cermics.enpc.fr^**** Email: erik.burman{at}epfl.ch We analyse a non-conforming finite-element method to approximateadvection–diffusion–reaction equations. The methodis stabilized by penalizing the jumps of the solution and thoseof its advective derivative across mesh interfaces. The a priorierror analysis leads to (quasi-)optimal estimates in the meshsize (sub-optimal by order in the L²-norm and optimal in thebroken graph norm for quasi-uniform meshes) keeping the Pécletnumber fixed. Then, we investigate a residual a posteriori errorestimator for the method. The estimator is semi-robust in thesense that it yields lower and upper bounds of the error whichdiffer by a factor equal at most to the square root of the Pécletnumber. Finally, to illustrate the theory we present numericalresults including adaptively generated meshes.     

Superconvergence and extrapolation of non-conforming low order finite elements applied to the Poisson equation

It is well-known that on uniform meshes the piecewise linearconforming finite element solution of the Poisson equation approximatesthe interpolant to a higher order than the solution itself.In this paper, this type of superclose property is studied forthe canonical interpolant defined by the nodal functionals ofseveral non-conforming finite elements of lowest order. By givingexplicit examples we show that some non-conforming finite elementsdo not admit the superclose property. In particular, we discusstwo non-conforming finite elements which satisfy the supercloseproperty. Moreover, applying a postprocessing technique, wecan also state a superconvergence property for the discretizationerror of the postprocessed discrete solution to the solutionitself. Finally, we show that an extrapolation technique leadsto a further improvement of the accuracy of the finite elementsolution.     

四边形单元h－收敛误差估计的探讨*

本文证明了四边形单元的h-收敛性,给出了相应的引理和定理,讨论了误差估计问题,为四边形单元h-收敛自适应有限元分析准备了基础。     

二阶常微初值问题的时间间断最小二乘有限元法

采用时间间断最小二乘线性有限元方法求解二阶常微分方程初值问题.利用回收技巧及离散Gronwall引理证明了方法的稳定性.通过引入有限元空间上的范数,给出了方法在该范数意义下丰满的误差估计.数值实验验证了理论分析结果.     

Gradient recovery type a posteriori error estimate for finite element approximation on non-uniform meshes

In this paper, we derive gradient recovery type a posteriori error estimate for the finite element approximation of elliptic equations. We show that a posteriori error estimate provide both upper and lower bounds for the discretization error on the non-uniform meshes. Moreover, it is proved that a posteriori error estimate is also asymptotically exact on the uniform meshes if the solution is smooth enough. The numerical results demonstrating the theoretical results are also presented in this paper.     

Some details of the Galerkin finite element method

The incorporation of the Galerkin technique in the finite element method has removed the constraint of finding a variational formulation for many problems of mathematical physics. The method has been successfully applied to many areas and has received wide acceptance. However, in the process of transplanting the concept from the Galerkin method for the entire domain to the Galerkin finite element method, some formal details have been overlooked or glossed over in the literature. This paper considers some of these details, including a possible reason for integration by parts and the contribution of interelement discontinuity terms.     

A posteriori error control for finite element approximations of elliptic eigenvalue problems

We develop a new approach to a posteriori error estimation for Galerkin finite element approximations of symmetric and nonsymmetric elliptic eigenvalue problems. The idea is to embed the eigenvalue approximation into the general framework of Galerkin methods for nonlinear variational equations. In this context residual-based a posteriori error representations are available with explicitly given remainder terms. The careful evaluation of these error representations for the concrete situation of an eigenvalue problem results in a posteriori error estimates for the approximations of eigenvalues as well as eigenfunctions. These suggest local error indicators that are used in the mesh refinement process.     

半线性抛物方程的时空有限元方法

讨论半线性抛物方程的连续Galerkin时空有限元方法,利用有限元方法和有限差分方法相结合的技巧,证明了弱解的存在唯一性,给出了时间最大模,空间L~2模,即L~∞(L~2)模误差估计.并给出数值算例证明了连续时空有限元方法对于半线性抛物方程的有效性.     

热传导型半导体瞬态问题特征变网格有限元法及其分析

热传导型半导体器件的瞬时状态由四个方程的非线性偏微分方程组的初边值问题所决定,其中电子位势方程是椭圆型的,电子和空穴浓度方程是对流扩散型的,温度方程为热传导型的。本文提出解这类问题的特征变网格有限元法,并进行了理论分析,在一定条件下,得到了某种意义下的最佳L^2误差估计结果。     

A non-conforming finite element method with anisotropic mesh grading for the Stokes problem in domains with edges

The solution of the Stokes problem in three-dimensional domainswith edges has anisotropic singular behaviour which is treatednumerically by using anisotropic finite element meshes. Thevelocity is approximated by Crouzeix–Raviart (nonconformingP₁ ) elements and the pressure by piecewise constants. Thismethod is stable for general meshes (without minimal or maximalangle condition). Denoting by N_e the number of elements in themesh, the interpolation and consistency errors are of the optimalorder h N_e^–1/3 which is proved for tensor product meshes.As a by-product, we analyse also nonconforming prismatic elementswith P₁ [oplus ] span {x₃²} as the local space for the velocitywhere x₃ is the direction of the edge.     

Adaptive finite element algorithms for eigenvalue problems based on local averaging type a posteriori error estimates

The local averaging technique has become a popular tool in adaptive finite element methods for solving partial differential boundary value problems since it provides efficient a posteriori error estimates by a simple postprocessing. In this paper, the technique is introduced to solve a class of symmetric eigenvalue problems. Its efficiency and reliability are proved by both the theory and numerical experiments structured meshes as well as irregular meshes. Dedicated to Charles A. Micchelli on his 60th birthday Mathematics subject classifications (2000) 65N15, 65N25, 65N30, 65N50. Subsidized by the Special Funds for Major State Basic Research Projects, and also supported in part by the Chinese National Natural Science Foundation and the Knowledge Innovation Program of the Chinese Academy of Sciences.