非匹配网格上Stokes-Darcy 问题非协调元离散的两水平和多水平方法

本文讨论了非匹配网格上Stokes-Darcy 问题的两种低阶非协调元方法, 给出了误差估计, 对耦合的非协调元离散问题, 通过粗网格求得的界面条件, 我们提出了一个解耦的两水平算法. 并且我们将两水平方法推广到多水平情形, 其只需在一个很粗的网格上解一耦合问题, 然后在逐步加细的网格上求解解耦的问题, 理论分析和数值试验都说明方法的高效性.     

带间断系数椭圆问题的P1非协调四边形元的加性 Schwarz 方法

本文讨论了带间断系数的二阶椭网问题的P1非协调四边形元的加性Schwarz方法.通过分析加性Schwaxz预处理后系统的特征值分布,我们证明了除少数小特征值外,其余所有特征值都有正的关于间断系数和网格尺寸拟一致的上下界.数值试验验证了我们的结论.     

非匹配网格上Stokes-Darcy模型的非协调元方法及其预条件技术

本文讨论了非匹配网格上Stokes-Darcy模型的两种低阶非协调元方法,证明了离散问题的适定性并得到了最优的误差估计.对离散出来的非对称不定线性方程组,我们提出了几种有效的预条件子,证明了预条件子的最优性.最后,数值试验验证了我们的理论结果.     

混合有限元非嵌套网格型预处理法

1 引言 众所周知,二阶椭圆型问题混合有限元离散以后的矩阵是不定的,所以对混合法很难形成一种有效的区域分解法,在文[9]、[10]、[11]中提出了一些混合有限元方法的区域分解法,但在实际计算中有很多局限性。最近Chen对混合有限元法提出一种全新的解释并把它应用到多重网格法中,他的基本思想是混合有限元离散的代数系统实际上等价于某个非协调有限元离散的代数系统,这样可把一个不定问题转化为一个正定问题,本文将基于这种思想考虑混合有限元的区域分解法。 若按传统的Dryia-widlund两水平加性Schwarz方法,要求两层网格间具有嵌套关系,这样在应用中将带来很大的不便。本文将不要求粗网格嵌入细网格中,减少两层网格间的     

非嵌套有限元的W循环多重网格方法

本文讨论了下述情形：１非嵌套网格；２曲边有限元；３非协调元；４拟协调元；５有限元的型函数有特殊性质，都能导致非嵌套的有限元空间．对一个包括上述情形的问题给出了非嵌套有限元的Ｗ循环多重网格方法，并证明了它的收敛性。     

粘弹性方程的非协调变网格有限元方法

讨论了粘弹性方程的Crouzeix-Raviart型非协调变网格有限元方法,在不需要引入传统分析中Riesz投影的情况下得到了最优误差估计.     

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

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

基于非协调元的加法型Schwarz交替法：——弱重迭情形

基于非协调元的加法型Ｓｃｈｗａｒｚ交替法──弱重迭情形黄建国（上海交通大学应用数学系）ＡＤＤＩＴＩＶＥＳＣＨＷＡＲＺＡＬＴＥＲＮＡＴＩＮＧＭＥＴＨＯＤＦＯＲＮＯＮＣＯＮＦＯＲＭＩＮＧＦＩＮＩＴＥＥＬＥＭＥＮＴ──ＣＡＳＥＯＦＷＥＡＫＯＶＥＲＬＡＰ￥Ｈ...     

各向异性网格下抛物方程一个新的非协调混合元收敛性分析

本文将 Crouzeix-Raviart 型非协调线性三角形元应用到抛物方程,建立了一个新的混合元格式.在抛弃传统有限元分析的必要工具 Ritz 投影算子的前提下,直接利用单元的插值性质和导数转移技巧, 分别得到了各向异性剖分下关于原始变量u 的H^-1-模和积分意义下L²-模以及通量p=-▽u 在L²-模下的最优阶误差估计.数值结果与我们的理论分析是相吻合的.     

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

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

MULTIGRID METHODS FOR MORLEY ELEMENT ON NONNESTED MESHES

1.IntroductionWeconsidersomemultigridalgorithmsforthebiharmonicequationdiscretizedbyMoneyelementonnonnestedmeshes.TOdefineamultigridalgorithm,certainintergridtransferoperatorhastobeconstructed.Throughtakingtheaveragesofthenodalvariables,weconstructanintergridtransferoperatorforMoneyelementonnonnestedmeshesthatsatisfiesacertainstableapproximationpropertywhichplaysakeyroleinmultigridmethodsfornonconformingplateelementsonnonnestedmeshes.Theso--calledregularity-approximaticnassurnptionisestablis…     

CONVERGENCE ANALYSIS OF MORLEY ELEMENT ON ANISOTROPIC MESHES

The main aim of this paper is to study tile convergence of a nonconforming triangular plate element-Morley element under anisotropic meshes. By a novel approach, an explicit bound for the interpolation error is derived for arbitrary triangular meshes （which even need not satisfy the maximal angle condition and the coordinate system condition ）, the optimal consistency error is obtained for a family of anisotropically graded finite element meshes.     

A MIXED VIRTUAL ELEMENT METHOD FOR THE BOUSSINESQ PROBLEM ON POLYGONAL MESHES

In this work we introduce and analyze a mixed virtual element method(mixed-VEM)for the two-dimensional stationary Boussinesq problem.The continuous formulation is based on the introduction of a pseudostress tensor depending nonlinearly on the velocity,which allows to obtain an equivalent model in which the main unknowns are given by the aforementioned pseudostress tensor,the velocity and the temperature,whereas the pressure is computed via a postprocessing formula.In addition,an augmented approach together with a fixed point strategy is used to analyze the well-posedness of the resulting continuous formulation.Regarding the discrete problem,we follow the approach employed in a previous work dealing with the Navier-Stokes equations,and couple it with a VEM for the convection-diffusion equation modelling the temperature.More precisely,we use a mixed-VEM for the scheme associated with the fluid equations in such a way that the pseudostress and the velocity are approximated on virtual element subspaces of H(div)and H¹,respectively,whereas a VEM is proposed to approximate the temperature on a virtual element subspace of H¹.In this way,we make use of the L²-orthogonal projectors onto suitable polynomial spaces,which allows the explicit integration of the terms that appear in the bilinear and trilinear forms involved in the scheme for the fluid equations.On the other hand,in order to manipulate the bilinear form associated to the heat equations,we define a suitable projector onto a space of polynomials to deal with the fact that the diffusion tensor,which represents the thermal conductivity,is variable.Next,the corresponding solvability analysis is performed using again appropriate fixed-point arguments.Further,Strang-type estimates are applied to derive the a priori error estimates for the components of the virtual element solution as well as for the fully computable projections of them and the postprocessed pressure.The corresponding rates of convergence are also established.Finally,several numerical examples illustrating the performance of the mixed-VEM scheme and confirming these theoretical rates are presented.     

MODIFIED MORLEY ELEMENT METHOD FOR A FOURTH ORDER ELLIPTIC SINGULAR PERTURBATION PROBLEM

This paper proposes a modified Morley element method for a fourth order ellipticsingular perturbation problem. The method also uses Morley element or rectangle Morleyelement, but linear or bilinear approximation of finite element functions is used in the lowerpart of the bilinear form. It is shown that the modified method converges uniformly in theperturbation parameter.     

CONVERGENCE ANALYSIS FOR A NONCONFORMING MEMBRANE ELEMENT ON ANISOTROPIC MESHES

Regular assumption of finite element meshes is a basic condition of most analysis offinite element approximations both for conventional conforming elements and nonconform-ing elements.The aim of this paper is to present a novel approach of dealing with theapproximation of a four-degree nonconforming finite element for the second order ellipticproblems on the anisotropic meshes.The optimal error estimates of energy norm and L~2-norm without the regular assumption or quasi-uniform assumption are obtained based onsome new special features of this element discovered herein.Numerical results are givento demonstrate validity of our theoretical analysis.     

各向异性网格下Bogner-Fox-Schmit元的超收敛

1引言 有限元超收敛的研究已有三十多年的历史,至今为止已取得了丰富的成果,可见[3],[18],[10],[6],[5]以及[17].1981年,陈传淼(见[2]345-372页)考虑了四阶板问题有限元解的超收敛性,得到了高一阶的超收敛结果.1995年,林群和罗平[8]用积分恒等式技巧再次研究这个问题,在均匀矩形网格的条件下,得到了更好的结论,有限元解与有限元插值函数之间的误差在H2范数下,有高二阶的超收敛.     

Hamilton-Jacobi方程的单调有限元格式

A monotone finite element scheme is obtained by applying the finite element method to the viscosity equation of the Hamilton-Jacobi equation on unstructured meshes. Under some constraints, we show that this scheme is monotone and its numerical solution converges to the viscosity solution of the Hamilton-Jacobi equa-tion. Numerical examples test the stability and the convergence of this scheme.     

双参数Morley有限元多重网格法

1、引言 多重网格方法是求解偏微分方程的高效快速算法,在实际中得到广泛应用．[2][6]中考察了Morley元的多重网格方法,并用于双调和方程问题。     

THE TWO-LEVEL STABILIZED FINITE ELEMENT METHOD BASED ON MULTISCALE ENRICHMENT FOR THE STOKES EIGENVALUE PROBLEM

In this paper, we first propose a new stabilized finite element method for the Stokes eigenvalue problem. This new method is based on multiscale enrichment, and is derived from the Stokes eigenvalue problem itself. The convergence of this new stabilized method is proved and the optimal priori error estimates for the eigenfunctions and eigenvalues are also obtained. Moreover, we combine this new stabilized finite element method with the two-level method to give a new two-level stabilized finite e...     

关于Morley元的误差估计

§1.引言 解薄板弯曲问题的三角形Morley元是六十年代出现的一种非协调元,它的形函数是完整的二次多项式,节点参数是单元顶点上的三个函数值及三边中点上的法向导数值.由于板弯曲问题的常应变是二次多项式,所以这是一个参数最少的非协调板元.由