superconvergence.docVIP

三角形有限元的超收敛性*  陈传淼单位：陈传淼(湖南师范大学计算研究所，长沙 410081 ??????? 1998-11-02收稿*数学天元基金重点资助项目(批准号 摘要　　基于三角形上的两类正交展开，对二阶椭圆问题研究了任意m次三角形有限元解(对偶数m)及平均梯度(对奇数m)在对称点上的超收敛性. 除此之外，再没有其他与方程系数无关的超收敛点.关键词　　超收敛　三角形元　任意次 　　设在平面域Ω上三角形剖分是拟一致的，记为m次三角形有限元子空间.　二阶椭圆方程的解u及在Ω的边界Γ上有u=v=0.　三角形有限元(Ritz投影)满足 (1)这里假定双线性型 用Wk,p(Ω)记Sobolev空间，范数为‖u‖k，p，Ω.　　　单元正交分析(EOA)的思想［1，2］是构造u的一个较好的展开uI∈，使得uI超接近于uh，由(1)式，这等价于要求R=u-uI近似正交于: A(uh-uI,v)=A(u-uI,v)=O(hm+2)‖v‖2,1,Ω，　m≥2，(2)这里网格范数‖v‖2,1,Ω=∑e‖v‖2,1,e.　若(2)式成立并取v=gh∈(离散Green函数，见下面引理1)，它导致在Ω中一致的超收敛估计 (uh-uI)(z)=O(hm+2lnh), D(uh-uI)(z)=O(hm+1lnh), m≥2.　(3)于是，R及其梯度DR的根正好分别是uh与Duh的超收敛点.　故单元正交分析由两部分组成，即用单元上的某种正交展开来构造所需的逼近uI,并用单元之间的消除技术证明(2)式.　　　最基本的工具是在区间e=(-h,h)上的L-型与M-型正交展开.　仍用u(t)记u(x)=u(ht),t∈(-1,1)，用D记对x求导，用记对t求导.　显然，iu=hiDiu=O(hi).　引进Legendre正交多项式lj(t)=i(t2-1)j/(2j)!!(注意lj(±1)=(±1)j).　对t积分后，得M-型多项式［3］：M0=1，M1=t,…,Mj+1=j-1(t2-1)j/(2j)!!.　若i-j=0或±2, 则有(Mi,Mj)≠0，否则(Mi,Mj)=0.　显然对n≥2有Mn(±1)=0.　　　为构造前述所需的展开，在过去20年中人们一直基于这样一种性质，即有限元解在每单元的顶点上有好的精度，于是常常使用M-型展开.　为了处理一般方程(组)，提出了单元消除技术，并研究了多种单元的整体超收敛性，如任意次一维元与矩形元，三角形线元与二次元［4～6］等.　由于奇次三角形元在每个单元的顶点没有超收敛, 奇次一维元与矩形元L2投影也是如此［7］.　为此，必须引进某些新思想和L-型展开.　　　文献［8～10］对高次三角形元得到了局部对称点上的内部超收敛性.　本文将用单元正交分析导出某些深刻结果.　　　定理　设u∈Wm+1,2(Ω)∩Wm+2，∞(Ωλ)，Ωλ=｛x∈Ω,dist(x,Ω)＞λ＞0｝，且uh∈是u的m次三角形有限元投影.　设对某个l＞0，下列负范数估计成立：对偶数m≥2有‖u-uh‖-l,Ω=O(hm+2)，或对奇数m≥1有O(hm+1).　设剖分在Ωλ上是均匀的，用Th记所有单元顶点与边中点的集合，则在z∈Th∩Ω2λ上， (u-uh)(z)=O(hm+2lnh)　(偶数m≥2),(4)或者平均梯度 (u-uh)(z)=O(hm+1lnh)　 (奇数m≥1).(5)除此之外，uh及Duh在所有单元中再没有其他与A的系数aij无关的超收敛点.　　　本文结果能推广到二阶椭圆组.　为了研究内部超收敛性，利用适当光滑的截断函数ψ∈在Ωλ中ψ≡1，可构造及其局部的Ritz投影对某个l＞0, 有［11］ ‖uh-uh‖1,∞，Ω2λ≤C(h2m‖u‖m+1,Ωλ+‖u-uh‖-l,Ω), 由此，uh的内部超收敛性能从局部投影uh的超收敛导出.　不失一般性，下面只要讨论u∈Wm+2,∞0(Ω′)就够了，这里Ω′=Ωλ，且Γ′=Ω′.　 1　在三角形上的M-型分解 　　记参考三角形E=｛(s,t)｜-1＜s,t＜1,s+t＜0｝，3个顶点是z1=(-1，-1)，z2=(1,-1)，z3=(-1,1)，且3边是E=S0+S1+S2，这里S0=｛s=-1,-1≤t≤1｝，S1=｛t=-1,-1≤s≤1｝，S2=｛t=-s,-1≤s≤1｝.　任意n次多项式u∈Pn包含Nn=(n+1)(n+2)/2项，其中3n项由u在E上的值确定，而其他Nn-2=Nn-3n项由u在E内的值确定.　对n≥3次齐次多项式｛sn,sn-1t,…,stn-1,tn｝，显然sn与tn由u在E上的值确定.　因为线性组合sn-1t-(-1)ns