非线性Sobolev方程的非协调H~1-Galerkin混合有限元方法

利用带约束的非协调旋转Q_1元和零阶R-T元对一类非线性Sobolev方程构造了总体自由度较小的非协调H~1-Galerkin混合元格式,基于单元插值算子的特殊性质,利用积分恒等式和平均值技巧,分别在半离散和全离散格式下得到了与以往文献中使用协调元方法完全相同的超逼近和超收敛结果.     

一个改进的Reissner-Mindlin矩形元

本文提出了一个改进的Reissner-Mindlin矩形非协调元方法：旋度用连续双线性元逼近，横向位移用旋转矩形非协调元逼近，而作为中间变量的剪切力用增广的分片常数元逼近，我们证明：该方法具有关于板厚一致稳定性和一致最优收敛性。     

变分不等式的不完全双二次非协调板元逼近

本文证明了不完全双二次非协调板元离散的 Poincaré型不等式成立,利用它建立了类似于文的引理,从而由此讨论了边界固支薄板弯曲障碍问题和带平均曲率型约束问题的不完全双二次非协调板元逼近,得到了离散解对真解的收敛性并讨论了离散问题的迭代解。     

关于不完全双二次非协调板元的收敛性

多年来,工程界普遍认为Irons的分片检验准则是检验非协调元收敛性的一个充要条件。作者在[3,4]中曾对三类四边形无证明了非协调元可以不通过分片检验仍然收敛,可见分片检验并非必要。最近,吴茂庆在[5]中给出了一个八个自由度的不完全双二次矩形板元,其形状函数由矩形四个角点上的函数值与四边中点上的法向导数值确定.这是一个非协调元,形状函数及其一阶偏导数在相邻单元的共同边界上不连续,有点象Morley元.[5]称此非协调元不通过分片检验,但却收敛,并给出收敛速度的一个估计:     

广义神经传播方程的非协调有限元分析

讨论了带约束的旋转Q_1元对广义神经传播方程的应用.利用Bramble-Hilbert引理及插值技巧,在不需要传统的Ritz投影的和任何修正格式情况下导出了相应的最优误差估计和超逼近结果.     

板的低阶间断与连续有限元法统一分析

基于Reissner-Mindlin板问题的间断Galerkin有限元逼近, 建立了一个对挠度空间和角位移空间取连续或间断元都适用的低阶有限元离散格式. 取剪切力空间为分片常数元, 挠度空间和角位移空间无论取间断元还是连续元, 格式都是一致稳定的, 并给出了H¹范数估计及L²范数估计. 作为应用,对几类低阶有限元空间讨论. 结果表明, 该格式对常见的低阶有限元空间都适用, 并且若至少有一个元连续时, 该格式需要的空间比[1,2]中的都要简单.       

电报方程的H~1-Galerkin非协调混合有限元逼近

在几乎均匀矩形剖分下取双线性Q_(11)元和类Wilson元为逼近空间,研究了一类电报方程的H~1-Galerkin非协调混合有限元方法.利用单元的特殊性质,积分恒等式和平均值技巧,在不需要验证LBB相容性条件及抛弃传统的Ritz投影的情形下,得到了半离散和全离散格式下原始变量及流量分别在H~1模和H(div,Ω)模意义下的超逼近性质.进一步地,借助插值后处理技术,导出了相应的整体超收敛结果.     

任意维空间Stokes算子特征值的确切下界

本文研究Stokes算子特征值的确切下界.主要思想有二:其一是在散度算子的核空间中,该算子是正定的;其二是扩充非协调旋转Q_1元标准的插值算子有可交换和质量守恒性质,因此可以证明插值误差常数及相应的后处理常数可以显式表示出来且与空间维数无关.于是,可以利用扩充非协调旋转Q_1元产生的真实特征值的渐近下界,通过一个简单的后处理,在任意正则张量积网格上得到真实特征值的确切下界.     

Stokes问题的变网格非协调有限元法

众所周知,由于LBB条件的限制,用非协调元格式求解速度—压力型的Stokes问题具有构造简单,计算经济和误差阶匹配等优点而在实际计算中经常被采用。用非协调格式处理Stokes问题首先是由Crouzeix-Raviart提出来的,他们采用分片线性三中点三角元这一非协调元作为速度逼近空间,用分片常数有限元空间作为压力逼近空间(即C—R格式),得     

非定常线性化Navier-Stokes方程的非协调流线扩散有限元法分析

对非定常线性化Navier-Stokes方程提出了非协调流线扩散有限元方法．用向后Euler格式离散时间,用流线扩散法处理扩散项带来的非稳定性．速度采用不连续的分片线性逼近,压力采用分片常数逼近．得到了离散解的存在唯一性以及在一定范数意义下离散解的稳定性和误差估计．     

Two lower order nonconforming quadrilateral elements for the Reissner-Mindlin plate

This paper generalizes two nonconforming rectangular elements of the Reissner-Mindlin plate to the quadrilateral mesh. The first quadrilateral element uses the usual conforming bilinear element to approximate both components of the rotation, and the modified nonconforming rotated Q ₁ element enriched with the intersected term on each element to approximate the displacement, whereas the second one uses the enriched modified nonconforming rotated Q ₁ element to approximate both the rotation and the displacement. Both elements employ a more complicated shear force space to overcome the shear force locking, which will be described in detail in the introduction. We prove that both methods converge at optimal rates uniformly in the plate thickness t and the mesh distortion parameter in both the H ¹-and the L ²-norms, and consequently they are locking free. This work was supported by the National Natural Science Foundation of China (Grant No. 10601003) and National Excellent Doctoral Dissertation of China (Grant No. 200718)     

Two lower order nonconforming rectangular elements for the Reissner-Mindlin plate

In this paper, we propose two lower order nonconforming rectangular elements for the Reissner-Mindlin plate. The first one uses the conforming bilinear element to approximate both components of the rotation, and the modified nonconforming rotated element to approximate the displacement, whereas the second one uses the modified nonconforming rotated element to approximate both the rotation and the displacement. Both elements employ a projection operator to overcome the shear force locking. We prove that both methods converge at optimal rates uniformly in the plate thickness in both the - and -norms, and consequently they are locking free.

     

New mixed finite elements for plane elasticity and Stokes equations

We consider mixed finite elements for the plane elasticity system and the Stokes equation. For the unmodified Hellinger-Reissner formulation of elasticity in which the stress and displacement fields are the primary unknowns, we derive two new nonconforming mixed finite elements of triangle type. Both elements use piecewise rigid motions to approximate the displacement and piecewise polynomial functions to approximate the stress, where no vertex degrees of freedom are involved. The two stress finite element ...     

Quadratic nonconforming finite element method for 3D Stokes equations on cuboid meshes

In this paper, the quadratic nonconforming brick element(MSLK element) introduced in [10] is used for the 3D Stokes equations. The instability for the mixed element pair MSLK-P_1 is analyzed, where the vector-valued MSLK element approximates the velocity and the piecewise P_1 element approximates the pressure. As a cure, we adopt the piecewise P_1 macroelement to discretize the pressure instead of the standard piecewise P_1 element on cuboid meshes. This new pair is stable and the optimal error estimate is achieved. Numerical examples verify our theoretical analysis.     

A superconvergent nonconforming mixed finite element method for the Navier–Stokes equations

The superconvergence for a nonconforming mixed finite element approximation of the Navier–Stokes equations is analyzed in this article. The velocity field is approximated by the constrained nonconforming rotated Q₁ (CNRQ₁) element, and the pressure is approximated by the piecewise constant functions. Under some regularity assumptions, the superconvergence estimates for both the velocity in broken H¹‐norm and the pressure in L²‐norm are obtained. Some numerical examples are presented to demonstrate our theoretical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 646–660, 2016     

非线性抛物方程混合有限元方法的高精度分析

采用双线性元及零阶Raviart-Thomas元（Q₁₁+Q₁₀×Q₀₁）对非线性抛物方程讨论了一种H¹-Galerkin混合有限元方法.提出一个线性化的二阶格式,利用数学归纳法有技巧的导出了原始变量u在H¹（Ω）模意义下及流量p=▽u在L²（Ω）模意义下的O（h²+τ²）阶超逼近性质.引入一个有关初始点的时间离散方程,并利用其得到了▽ ·在L²（Ω）模意义下的O（h²+τ²）阶的超逼近结果.同时利用插值后处理技巧得到整体超收敛.最后,数值算例结果验证了理论分析（其中,h是剖分参数,τ是时间步长）.     

平面线弹性纯位移边值问题低阶非协调矩形元的Robust多重网格方法

1引 言 对于各向同性,均匀介质的平面线弹性问题,当Lamé常数λ→∞(泊松率v→0.5)时,即对于几乎不可压介质,通常的协调有限元格式的解往往不再收敛到原问题的解,或者达不到最优收敛阶,这就是所谓的闭锁现象(见[3],[7],[8]及[10]).究其原因,在通常的有限元分析中,其误差估计的系数与λ有关,当λ→∞时,该系数将趋于无穷大.因此为克服闭锁现象就需要构造特殊的有限元格式,使得当λ→∞时,有限元逼近解仍然收敛到原问题的解.     

Nonconforming elements in least-squares mixed finite element methods

In this paper we analyze the finite element discretization for the first-order system least squares mixed model for the second-order elliptic problem by means of using nonconforming and conforming elements to approximate displacement and stress, respectively. Moreover, on arbitrary regular quadrilaterals, we propose new variants of both the rotated nonconforming element and the lowest-order Raviart-Thomas element.

     

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.