双曲型积分微分方程的非协调任意四边形H1Galerkin混合有限元方法.pdfVIP

第44卷第4期 郑州大学学报(理学版) V01．44No．4 2012年12月 Univ．(Nat．Sci．Ed．) DeC．2012 J．Zhengzhou DOI：10．3969／j．issn／1671-6841．2012．04．007 双曲型积分微分方程的非协调任意四边形 H1．Galerkin混合有限元方法 王海红1， 郭 城2 摘要：针对双曲型积分微分方程问题，研究了非协调任意四边形H1．Galerkjn混合有限元方法．在半离散格式下，利 用所选单元本身的特点，在不需要Ritz．Volterra投影的情况下得到了与传统协调混合有限元方法相同的误差估计． 关键词：H1．Galerkin混合有限元方法；非协调有限元；双曲型积分微分方程；误差估计 中图分类号：0242．21 文献标志码：A 文章编号：1671—6841(2012)04—0031—04 0 引言 双曲型积分微分方程在具有记忆性质材料的热传导、核反应动力学、生物学、力学以及松散介质中压力 等问题中有着广泛的应用．关于双曲型积分微分方程的有限元研究中，常用的方法有协调的Galerkin方 法¨≈J，传统的混合有限元方法【3j等．在常规的混合有限元方法研究中常涉及两个有限元逼近空间，且它们 需要满足所谓的Inf-Sup条件或LBB条件㈠』，这样使得混合元的构造或自由度的选取变得相当复杂．最近出 现了许多绕开条件的混合元方法，例如最小二乘法混合元方法∞“j，虽然该方法不要求位移函数和通量函 数的有限元满足LBB条件，然而位移函数的逼近却不能使用非协调有限元．为了克服上述缺点，[7]提出了 H1．Galerkin混合有限元方法，该方法显著的优点是两个有限元空间可以取成不同次数的多项式空间，而且 q 不必满足LBB条件，这使得该方法成为近年来受关注的研究课题之一【8J．然而遗憾的是这些研究只局限 于协调元的情况，而且一些投影是必不可少的分析工具．由于非协调元在降低计算量方面有显著的优势，最 考虑的是矩形元的情况，对于任意四边形的非协调混合有限元方法研究还比较少． 作者主要研究了非协调H1．Galerkin混合有限元方法在双曲型积分微分方程中的应用，其两个逼近空间 了与传统方法相同的误差估计． 有限元空间的选取 2，3，4，在K上构造有限元(K，P‘“，三‘n，江1，2)如下： 三1={秽1，勘2，秽3，秽4}， P1=span{札(戈，Y)}，秽。=v(A。)，i=1，2，3，4， 挚=∽’，∥，∥，∥，∥，∥，歹=l，2}， 收稿日期：2012—06—09 基金项目：国家自然科学基金资助项目，编61143002． 作者简介：王海红(1980一)，女，讲师，博士，主要从事有限元方法研究，E—mml：waih777@163．com 万方数据 32 郑州大学学报(理学版) 第44卷 P2=span{^l(戈，y)，i=1，2，3，4，妒(戈)，妒(，，)}Xspan{Ⅳi(茄，)，)，i=1，2，3，4，妒(x)，9(，，)}， Ⅳ1(菇，，，)=0．25(1一菇)(1～y)，Ⅳ2(戈，)，)=0．25(1+石)(1一)，)， ；』’)：(【￡华矗d；)／I文l，靠，：(【￡学矗面)／I玄I． 。K a酊 Jk d， 上述的插值j2；和P；是适定的¨2|． 设n是平面菇一Y上的一个有界凸多边形区域，厶是力的一个正则剖分且满足正则性假设¨4|，任意K∈ ^，设其四个顶点坐标为力=。，U，K满足正则性条件或拟一致假设条件正则剖分．对任意的单元K的中心点 2，3，4，而且存