A POSTERIORI ERROR ESTIMATE FOR BOUNDARY CONTROL PROBLEMS GOVERNED BY THE PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS

In this paper, we discuss the a posteriori error estimate of the finite element approximation for the boundary control problems governed by the parabolic partial differential equations. Three different a posteriori error estimators are provided for the parabolic boundary control problems with the observations of the distributed state, the boundary state and the final state. It is proven that these estimators are reliable bounds of the finite element approximation errors, which can be used as the indicators of the mesh refinement in adaptive finite element methods.     

SUPERCONVERGENCE OF LEAST-SQUARES MIXED FINITE ELEMENT FOR SECOND-ORDER ELLIPTIC PROBLEMS

In this paper the least-squares mixed finite element is considered for solving secondorder elliptic problems in two dimensional domains. The primary solution u and the flux er are approximated using finite element spaces consisting of piecewise polynomials of degree k and r respectively. Based on interpolation operators and an auxiliary projection,superconvergent H^1-error estimates of both the primary solution approximation uh and the flux approximation σh are obtained under the standard quasi-uniform assumption on finite element partition. The superconvergence indicates an accuracy of O(h^r 2) for the least-squares mixed finite element approximation if Raviart-Thomas or Brezzi-DouglasFortin-Marini elements of order r are employed with optimal error estimate of O(h^r l).     

Strong uniqueness of best approximations in spaces of bounded linear operators

The present paper is concerned with problems of the strong uniqueness of thebest approximation and the characterization of a uniqueness element in operator spaces.Some results on the strong uniqueness of the best approximation operator from RS-setsare proved and the uniqueness element of a sun in the compact operator space from c_0to c_0 is characterized by the strict Kolmogorov's condition. Some recent results due toLewicki and others are extended and improved.     

美式债券期权定价问题的有限元方法

The aim of this paper is to investigate the finite element methods for pricing the American put option on bonds. Based on a new variational inequality equation for the option pricing problems, both semidiscrete and fully discretized finite element approximation schemes are established. It is proved that the finite element methods are stable and convergent under L2 and H^1 norms.     

Numerical Simulation for Convection-Diffusion Problem with Periodic Micro-Structure

In this article, the convection dominated convection-diffusion problems with the periodic micro-structure are discussed. A two-scale finite element scheme based on the homogenization technique for this kind of problems is provided. The error estimates between the exact solution and the approximation solution, of the homogenized equation or the two-scale finite element scheme are analyzed. It is shown that the scheme provided in this article is convergent for any fixed diffusion coefficient 5, and it may be convergent independent of δ under some conditions. The numerical results demonstrating the theoretical results are presented in this article.     

SOME ARGUMENTS FOR RECOVERING THE FINITE ELEMENT~*

The dual argument is well known for recoving the optimal L_2-error of the finite element method in elliptic context. This argument, however, will lose efficacy in hyperbolic case. An expansion argument and an approximation argument are presented in this paper to recover the optimal L_2-error of finite element methods for hyperbolic problems. In particular, a second order error estimate in L_2-norm for the standard linear finite element method of hyperbolic problems is obtained if the exact solution is smooth and the finite element mesh is almost uniform, and some superconvergence estimates are also established for leas smooth solution.     

HOW TO PROVE THE DISCRETE RELIABILITY FOR NONCONFORMING FINITE ELEMENT METHODS

Optimal convergence rates of adaptive finite element methods are well understood in terms of the axioms of adaptivity.One key ingredient is the discrete reliability of a residualbased a posteriori error estimator,which controls the error of two discrete finite element solutions based on two nested triangulations.In the error analysis of nonconforming finite element methods,like the Crouzeix-Raviart or Morley finite element schemes,the difference of the piecewise derivatives of discontinuous approximations to the distributional gradients of global Sobolev functions plays a dominant role and is the object of this paper.The nonconforming interpolation operator,which comes natural with the definition of the aforementioned nonconforming finite element in the sense of Ciarlet,allows for stability and approximation properties that enable direct proofs of the reliability for the residual that monitors the equilibrium condition.The novel approach of this paper is the suggestion of a right-inverse of this interpolation operator in conforming piecewise polynomials to design a nonconforming approximation of a given coarse-grid approximation on a refined triangulation.The results of this paper allow for simple proofs of the discrete reliability in any space dimension and multiply connected domains on general shape-regular triangulations beyond newest-vertex bisection of simplices.Particular attention is on optimal constants in some standard discrete estimates listed in the appendices.     

A posteriori error estimator for eigenvalue problems by mixed finite element method

In this paper,a residual type of a posteriori error estimator for the general second order elliptic eigenpair approximation by the mixed finite element method is derived and analyzed,based on a type of superconvergence result of the eigenfunction approximation.Its efficiency and reliability are proved by both theoretical analysis and numerical experiments.     

SUPERCONVERGENCE ANALYSIS OF FINITE ELEMENT METHODS FOR OPTIMAL CONTROL PROBLEMS OF THE STATIONARY BENARD TYPE

In this paper, we consider the finite element approximation of the distributed optimal control problems of the stationary Benard type under the pointwise control constraint. The states and the co-states are approximated by polynomial functions of lowest-order mixed finite element space or piecewise linear functions and the control is approximated by piecewise constant functions. We give the superconvergence analysis for the control; it is proved that the approximation has a second-order rate of convergence. We further give the superconvergence analysis for the states and the co-states. Then we derive error estimates in L^∞-norm and optimal error estimates in L^2-norm.     

A posterior error analysis for the nonconforming discretization of Stokes eigenvalue problem

In this paper, we present a posteriori error estimator for the nonconforming finite element approximation, including using Crouzeix–Raviart element and extended Crouzeix–Raviart element, of the Stokes eigenvalue problem. With the technique of Helmholtz decomposition, we first give out a posteriori error estimator and prove it as the global upper bound and local lower bound of the approximation error. Then, by deleting a jump term in the indicator, another simpler but equivalent indicator is obtained. Some numerical experiments are provided to verify our analysis.     

Stability analysis of frame structures by bubble-function method

In the stability analysis of frame structures, the results by conventional finite element method (FEM) in which one member is taken as one element are sometimes unavailable. This paper took a new basic function system with bubble functions as the shape function of a bar element to develop a bubble function finite element method (BFEM), in which the bending and the geometric stiffness matrices were derived from the principle of virtual work. Bubble functions are finite element modes that are located entirely within a single element and are zero on boundaries of the element, but are nonzero at the other points. BFEM is as concise as conventional bar FEM but has better accuracy, and is adaptable to the buckling analysis of all kinds of frame structures. The use of bubble functions significantly improves the convergence of finite element analysis, and efficiently reduces the computation cost for the buckling analysis of frame structures. Numerical results show that using bubble functions in finite element for the stability analysis of structures is very efficient, especially for high-rise and large-scale frame structures.     

关于MZ1元,MZ2元和MB1元的误差估计

本文给出MZ1元、MZ2元和MB1元的收敛阶,并且证明MZ1元与MB1元等价。     

Rail-bridge coupling element of unequal lengths for analysing train-track-bridge interaction systems

This paper presents a rail-bridge coupling element of unequal lengths, in which the length of a bridge element is longer than that of a rail element, to investigate the dynamic problem of train-track-bridge interaction systems. The equation of motion in matrix form is given for a train-track-bridge interaction system with the proposed element. The first two numerical examples with two types of bridge models are chosen to illustrate the application of the proposed element. The results show that, for the same length of rail element, (1) the dynamic responses of train, track and bridge obtained by the proposed element are almost identical to those obtained by the rail-bridge coupling element of equal length, and (2) compared with the rail-bridge coupling element of equal length, the proposed element can help to save computer time. Furthermore, the influence of the length of rail element on the dynamic responses of rail is significant. However, the influence of the length of rail element on the dynamic responses of bridge is insignificant. Therefore, the proposed element with a shorter rail element and a longer bridge element may be adopted to study the dynamic responses of a train-track-bridge interaction system. The last numerical example is to investigate the effects of two types of track models on the dynamic responses of vehicle, rail and bridge. The results show that: (1) there are differences of the dynamic responses of vehicle, rail and bridge based on the single-layer and double-layer track models, (2) the maximum differences increase with the increase of the mass of sleeper, (3) the double-layer track model is more accurate.     

有限元法中大单元的构造

在通常的有限元法中,单元内的插值多项式的阶数固定不变,通过加密剖分网格来提高精度.大单元法则剖分的网格固定不变而通过增加单元内逼近级数的项数来提高精度. 本文提出采用两套变量的办法来构造大单元,即单元内采用一套变量,单元的边界上采用另一套变量,然后用杂交罚函数法把两者联系起来.这种方法能适用于任何椭圆型方程,任意几何形状区域以及任何复杂的边界条件.本文用严密的数学方法证明了:在一般情况下,这种方法的精度比通常的有限元法和文[7]的大单元法高得多.即在达到相同的精度时,本文方法所需要的自由度(即未知数数目)比上述两种方法少得多.     

A quadrilateral plate bending element based on deformation modes

In this paper, a quadrilateral element is proposed for the analysis of thin plate bending. This element is non-conforming and consists of four-nodes and twelve degrees of freedom. A third-order field for the element displacement is written in terms of the deformation modes. Moreover, the rotational fields are obtained by utilizing the first-order Jacobean matrix. All interpolation functions are explicitly found by the presented formulation. The stiffness matrix of the element is then computed by using these functions. Finally, the accuracy of the suggested element is evaluated by solving some thin plate bending structures. Numerical findings reveal the new quadrilateral element MKQ12 is robust and accurate for analysis of thin plates.     

半导体器件瞬态模拟的对称正定混合元方法

提出具有对称正定特性的混合元格式求解非稳态半导体器件瞬态模拟问题。提出一个最小二乘混合元方法、一个新的具有分裂和对称正定性质的混合元格式和一个解经典混合元方程的对称正定失窃工格式求解电场位势和电场强度方程;提出一个最小二乘混合元格式求解关于电子与空穴浓度的非稳态对流扩散方程,浓度函数和流函数被同时求解;采用标准的有限元方法求解热传导方程。建立了误差分析理论。     

基于二重网格的定常Navier-Stokes方程的局部和并行有限元算法

对二维定常的不可压缩的Navier-Stokes方程的局部和并行算法进行了研究．给出的算法是多重网格和区域分解相结合的算法,它是基于两个有限元空间:粗网格上的函数空间和子区域的细网格上的函数空间．局部算法是在粗网格上求一个非线性问题,然后在细网格上求一个线性问题,并舍掉内部边界附近的误差相对较大的解．最后,基于局部算法,通过有重叠的区域分解而构造了并行算法,并且做了算法的误差分析,得到了比标准有限元方法更好的误差估计,也对算法做了数值试验,数值结果通过比较验证了本算法的高效性和合理性．     

九参数广义协调元的收敛性

众知周所,解板弯曲问题的Zienkiewicz不协调三次元只对特殊的单元剖分才收敛.但由于这种元采用单元顶点的函数值及二个一阶导数值作为节点参数,计算简单,总体自由度少,所以相继出现一些对Zienkiewicz元的改进形式,使之对任意剖分均收敛,如拟协调元,TRUNC元,Specht元.对这些元的分析见[7—9].最近龙     

三维热传导型半导体问题的交替方向特征有限元方法

提出交替方向特征有限元方法,对电场位势方程采用混合元格式,对电子,空穴浓度方程采用交替方向特征有限元格式,对温度方程提出交替方向格式.应用向量积计算及先验估计理论和技巧,得到最佳的L2误差估计.     

A posteriori finite element error estimators for indefinite elliptic boundary value problems ∗

A general construction technique is presented for a posteriori error estimators of finite element solutions of elliptic boundary value problems that satisfy a Gång inequality. The estimators are obtained by an element–by–element solution of ‘weak residual’ with or without considering element boundary residuals. There is no order restriction on the finite element spaces used for the approximate solution or the error estimation; that is, the design of the estimators is applicable in connection with either one of the h–p–, or hp– formulations of the finite element method. Under suitable assumptions it is shown that the estimators are bounded by constant multiples of the true error in a suitable norm. Some numerical results are given to demonstrate the effectiveness and efficiency of the approach.