Constructing Order Two Superconvergent WG Finite Elements on Rectangular Meshes

In this paper, we introduce a stabilizer free weak Galerkin (SFWG) finite element method for second order elliptic problems on rectangular meshes. With a special weak Gradient space, an order two superconvergence for the SFWG finite element solution is obtained, in both $L^2$ and $H^1$ norms. A local post-process lifts such a $P_k$ weak Galerkin solution to an optimal order $P_{k+2}$ solution. The numerical results confirm the theory.     

A Stabilized Low Order Finite Element Method for Three Dimensional Elasticity Problems

We introduce a low order finite element method for three dimensional elasticity problems. We extend Kouhia-Stenberg element [12] by using two nonconforming components and one conforming component, adding stabilizing terms to the associated bilinear form to ensure the discrete Korn's inequality. Using the second Strang's lemma, we show that our scheme has optimal convergence rates in $L^2$ and piecewise $H^1$-norms even when Poisson ratio $\nu$ approaches $1/2$. Even though some efforts have been made to design a low order method for three dimensional problems in [11,16], their method uses some higher degree basis functions. Our scheme is the first true low order method. We provide three numerical examples which support our analysis. We compute two examples having analytic solutions. We observe the optimal $L^2$ and $H^1$ errors for many different choices of Poisson ratios including the nearly incompressible cases. In the last example, we simulate the driven cavity problem. Our scheme shows non-locking phenomena for the driven cavity problems also.     

二阶椭圆问题基于泡函数的简化的稳定化二阶混合有限元格式

本文研究二阶椭圆方程基于泡函数的稳定化的二阶混合有限元格式,通过消去泡函数导出一种自由度很少的简化的稳定化的二阶混合有限元格式, 误差分析表明消去泡函数的简化格式与带有泡函数的格式具有相同的精度而可以节省6N_p个自由度(其中N_p三角形剖分中的顶点数目).       

Mortar Finite Volume Method with Adini Element for Biharmonic Problem

In this paper, we construct and analyse a mortar finite volume method for the discretization for the biharmonic problem in R2. This method is based on the mortar-type Adini nonconforming finite element spaces. The optimal order H2-seminorm error estimate between the exact solution and the mortar Adini finite volume solution of the biharmonic equation is established.     

求解Brinkman方程的修正弱有限元方法

本文针对Brinkman方程引入了一种修正弱Galerkin(MWG)有限元方法.我们通过具有两个离散弱梯度算子的变分形式来逼近模型. 在MWG方法中, 分别用次数为$k$和$k-1$的不连续分段多项式来近似速度函数$u$和压力函数$p$. MWG方法的主要思想是用内部函数的平均值代替边界函数. 因此, 与WG方法相比, MWG方法在不降低准确性的同时, 具有更少的自由度, 对于任意次数不超过$k-1$ 的多项式,MWG方法均可以满足稳定性条件. MWG 方法具有高度的灵活性, 它允许在具有一定形状正则性的任意多边形或多面体上使用不连续函数. 针对$H^1$和$L^22$范数下的速度和压力近似解, 建立了最优阶误差估计. 数值算例表明了该方法的准确性, 收敛性和稳定性.     

A robust nonconforming -element

Finite element methods for some elliptic fourth order singular perturbation problems are discussed. We show that if such problems are discretized by the nonconforming Morley method, in a regime close to second order elliptic equations, then the error deteriorates. In fact, a counterexample is given to show that the Morley method diverges for the reduced second order equation. As an alternative to the Morley element we propose to use a nonconforming -element which is -conforming. We show that the new finite element method converges in the energy norm uniformly in the perturbation parameter.

     

A Modified Weak Galerkin Finite Element Method for Sobolev Equation

For Sobolev equation, we present a new numerical scheme based on a modified weak Galerkin finite element method, in which differential operators are approximated by weak forms through the usual integration by parts. In particular, the numerical method allows the use of discontinuous finite element functions and arbitrary shape of element. Optimal order error estimates in discrete $H^1$ and $L^2$ norms are established for the corresponding modified weak Galerkin finite element solutions. Finally, some numerical results are given to verify theoretical results.     

二个非协调膜元的收敛问题

为了用有限元方法求(0-2)的近似解,定义有限元空间.记P(K),?K∈?_h,为单元K上多项式组成的有限维空间,v_h∈P(K)可由K上的节点参数或其它类型的参数(例如函数的导数或函数本身在单元K上积分的数值)唯一决定.     

A new quadratic nonconforming finite element on rectangles

A new quadratic nonconforming finite element on rectangles (or parallelograms) is introduced. The nonconforming element consists of P₂ ⊕ Span{x²y,xy²} on a rectangle and eight degrees of freedom. Our element is essentially of seven degrees of freedom since the degree of freedom associated with the integration on rectangle is essentially of bubble‐function nature. Global basis functions are constructed for both Dirichlet and Neumann type of problems; accordingly the corresponding dimensions are counted. The local and global interpolation operators are defined. Error estimates of optimal order are derived in both broken energy and L²(Ω) norms for second‐order of elliptic problems. Brief numerical results are also shown to confirm the optimality of the presented quadratic nonconforming element. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

A Conforming Quadratic Polygonal Element and Its Application to Stokes Equations

In this paper, we construct an $H^1$-conforming quadratic finite element on convex polygonal meshes using the generalized barycentric coordinates. The element has optimal approximation rates. Using this quadratic element, two stable discretizations for the Stokes equations are developed, which can be viewed as the extensions of the $P_2$-$P_0$ and the $Q_2$-(discontinuous)$P_1$ elements, respectively, to polygonal meshes. Numerical results are presented, which support our theoretical claims.     

定常Navier-Stokes方程的一个二阶非协调三角型混合元格式（英文）

In this paper, a nonconforming triangular mixed finite element scheme with second order convergence behavior is proposed for the stationary Navier-Stokes equations.The new nonconforming triangular element is taken as approximation space for the velocity and the linear element for the pressure. The convergence analysis is presented and optimal error estimates of both broken H~1-norm and L~2-norm for velocity as well as the L~2-norm for the pressure are derived.     

A Second Order Nonconforming Triangular Mixed Finite Element Scheme for the Stationary Navier-Stokes Equations

In this paper,a nonconforming triangular mixed finite element scheme with second order convergence behavior is proposed for the stationary Navier-Stokes equations.The new nonconforming triangular element is taken as approximation space for the velocity and the linear element for the pressure.The convergence analysis is presented and optimal error estimates of both broken H1-norm and L2-norm for velocity as well as the L2-norm for the pressure are derived.     

位移障碍下四阶变分不等式问题的非协调有限元一般误差估计式

1.引 言 关于二阶变分不等式问题的非协调有限元逼近已有大量研究[1-5].但是,对于四阶变分不等式的研究相对而言较少[6-7].[8,9,10]给出了位移障碍问题的非协调有限元,包括C0元(如Zienkiewicz元及Adini元)和非C0元(如Morley元及De Veubeke元)逼近的理论分析及最优误差估计.经过仔细分析发现,其成功的关键技巧是充分利用上述单元的一个     

A quadrilateral Morley element for biharmonic equations

In this paper, we propose a Morley-type finite element for quadrilateral meshes to solve biharmonic problems. For each quadrilateral $Q$ , the finite element space is defined by the span of $P_2(Q)$ plus two functions in $P_3(Q)$ . Each of the cubic polynomials is the product of a pair of equations of opposite edges and the equation of the bimedian between them. The degrees of freedom consist of the values at vertices and integrals of normal derivatives over edges. Optimal orders of convergence are proved both in discrete $H^2$ and $H^1$ seminorms. Several numerical tests confirm the convergence analysis.     

PRECONDITIONING MIXED ELEMENT METHOD FOR NONSELFADJOINT AND INDEFINITE SECOND ORDER ELLIPTIC PROBLEMS*

In this paper, an equivalence between mixed element method and nonconforming element method for nonselfad joint and indefinite second order elliptic problems is established without using any bubble functions. It is proved that the H~1-condition number of preconditioned operator B_h~(-1)A_h is uniformly bounded and its B_h-singular values cluster in a positive finite interval, where A_h is the equivalent nonconforming element discretization of nonselfad joint and indefinite second order elliptic operator A, B_h is usual noncon forming element discretization of selfadjoint and positive definite second order elliptic operator B. Finally a simple V-cycle multigrid implementation of B_h~(-1) is given.     

二次有限元解的渐近展开与外推

在Poisson方程的求解域Ω存在一致的三角剖分,并且相邻两初始单元构成平行四边形的假设下,证明了若Poisson方程的解u属于H6(Ω),那么二次有限元的误差有h4的渐近展开.基于误差的渐近展开,可以利用h4-Richardson外推进一步提高数值解的精度阶,并且能够得到一个后验误差估计.最后,一个数值算例验证了理论分析.     

A new rotated nonconforming pyramid element

In this paper, a new nonparametric nonconforming pyramid finite element is introduced. This element takes the five face mean values as the degrees of the freedom and the finite element space is a subspace of P₂. Different from the other nonparametric elements, the basis functions of this new element can be expressed explicitly without solving linear systems locally, which can be achieved by introducing a new reference pyramid. To evaluate the integration, a class of new quadrature formulae with only two/three equally weighted points on pyramid are constructed. We present the error estimation in the presence of quadrature formulae. Numerical results are shown to confirm the optimality of the convergence order for the second order elliptic problems.     

各向异性网格下特征值问题的有限元分析

主要目的是在各向异性网格下研究二阶椭圆特征值问题的两类非协调有限元—类Wilson矩形元和Carey三角形元—的收敛性分析.通过新的技巧和方法,得到了与传统有限元网格剖分下相同的特征对的最优误差估计.推广了已有的结果.     

二阶椭圆问题的各向异性非协调Crouzeix-Raviart型有限元方法

对满足最大角条件和坐标系条件的二维区域中的各向异性一般三角形网格,研究了二阶椭圆问题的非协调Crouzeix-Raviart型线性三角形有限元逼近,得到了最优的能量模和L2-模误差估计结果．     

A simple nonconforming quadrilateral finite element

We introduce and analyze a simple nonconforming quadrilateral finite element and then we derive optimal a priori error estimates for arbitrary regular quadrilaterals. The idea of extension to some non-conforming elements for three-dimensional hexahedrons is also presented.