Z—Z小片恢复技术超收敛性的研究.pdfVIP

第29卷第 1期 长春师范学院学报 (自然科学版) 2010年 2月 V01．29 No．1 JournalofChangehunNormalUniversity(Natural~ienee) Feb．2O1O Z—Z小片恢复技术超收敛性的研究 白春艳，裴晓雯，李明辉 (沈阳化工学院数理系，辽宁沈阳 110142) [摘 要]本文介绍了近年来发展起来的新的有限元的后处理技术一z—z小片恢复技术，并探讨了运 用z—z小片恢复技术对两点边值问题的有限元解的一阶导数进行后处理所产生的超收敛性和后验误 差估计．着重在高斯点利用最Jb--乘拟合技术对于k=1时的一次元和 k=2时的二次元的后处理公式 进行了推导，得出了后处理公式．并举例验证后处理公式对于两点边值问题一P(Ⅱ())+qu()=f ()，u(a)=u(b)=0，∈[a，b]的正确性。 [关键词】有限元；超收敛性；高斯点；后处理技术 [中图分类号]O241．2 [文献标识码]A [文章编号】1008—178X(2010)01—0007—03 有限元超收敛性的研究始于 l9世纪 7O年代 ，现有的研究工作基本遵循两种途径 ：一种是先找出有限元插 值逼近的超收敛点，再利用插值弱估计等手段导出有限元解本身所具有的超收敛性质；另一种是利用各种后处 理技术来导出经过后处理的有限元解的超收敛性．…Z—Z小片恢复技术是由Zienkiewice．o．c和 Zhu．J．z于 1992年发明的一种新的导数恢复技术，称为插值恢复的超收敛性．]由于这种技术的可行性及它在不规则网 格和普通几何体上的灵活性，它已被广泛应用到实际问题中．[3]但是，现有对有限元的研究都仅限于在理论上 证明对有限元解进行后处理确实能产生超收敛性或强超收敛性，而没有对后处理公式进行具体的推导和验证． 本文主要在高斯点利用最小二乘拟合技术[3]对后处理公式进行了推导，得出了后处理公式，并验证了对于一次 元，u^()经过后处理在节点兢处确实具有强超收敛性． 1 超收敛性的分析及后验误差估计 【】 对于两点边值问题 f一(a2() )一(al() )+ao()Ⅱ=f， (1) 1¨(o)：u(1)：0． (2) 其中，=(0，1)，假设 af(i=0，1，2)和厂都是严格光滑的，且V ∈，，a2()a0(a是一个常数)，(1)， (2)的等价方程就是找一u∈矾(，)，使得 (a2u，)+(alu，)+(aoH，)=(f，)，V ∈月6．()． (3) 设rh(0h寺)是 ，的一个子划分， ={}No，0=o…，v=1．记 = (一l，)，h= 一f—l，h=maxhf，S^={口∈H ()，f，∈P(，1)}，s2={秽∈础(，)，l，∈P(，1)}， s^和52都是在划分r^下，，上次数不超过r的连续分片多项式空间．设U^∈Is2是(3)的有限元解，即 (a2uh，口)+(alUh， )+(aoUh，)=(f，)，Vt7∈so． (4) 那么，从(3)中减去(4)得 (a2(u一u^)，)+(a1(u— ^)，)+(ao(u—u^)，)=0，VvE．s2． (5) 设U^，t，