抛物问题各向异性有限元的超收敛分析.pdfVIP

维普资讯 应 用 数 学 MATHEMATICA APPLICATA 2007，20(4)：659～665 抛物问题各 向异性有限元的超收敛分析 石东洋，高新慧 (郑州大学数学系，河南 郑州 450052) 摘要：本文研究具有各向异性特征的双二次元对二阶抛物方程的逼近．通过积分恒等 式和插值后处理技术，在各向异性网格下得到 了相应的超逼近和超收敛结果． 关键词：二阶抛物方程；各向异性元；超逼近及超收敛 ． 中图分类号：0242．21 AMS(2000)主题分类：65N15；65N30 文献标识码：A 文章编号：1001—9847(2007)04—0659—07 1．引言 众所周知 ，关于抛物方程解的有限元逼近及高精度分析已有一些研究[1]．但这些传统的 有限元方法不仅要求网格剖分满足正则性条件[6]，即^ ／ ≤ c，VK ∈ ，这里 为区域n 上的剖分，h 和p 分别表示单元K 的直径和内切于K 的最大 圆的直径 ，c为一个与h无关的 正常数，而且把 Ritz投影与修正格式作为必备的工具．这在很大程度上限制了有限元方法的 应用范围．近年来出现了有关各 向异性单元的研究E7-1o]，并着重分析了二阶椭 圆问题的误差 估计．然而关于抛物问题在各向异性剖分下高精度分析的研究到 目前为止还不多见．本文讨论 具有各 向异性特征 的双二次元对二阶抛物方程 的逼近 ，利用积分恒等式技巧，在不需要 Ritz 投影或任何修正格式情况下得到了各向异性网格下的超逼近结果．在此基础上，利用插值后处 理技术，给出了整体超收敛估计．由于本文的显著特色是摆脱了以往文献对网格剖分正则性的 要求以及把 Ritz投影或修正格式作为必备工具的限制 ，从而扩大了其应用范围． 2．抛物问题及单元构造 考虑模型问题[1]： r“，一△“一_厂， 在nx (0，T]内， “一0， 在anX (0，T]上， (1) lu(X，0)一‰(x)， 在 n内， 其中n(=二R。为有界凸区域 ，_厂∈L。(n)，an为n 的边界．与(1)等价 的变分问题为 ：求 “∈ H (n)，使 V ∈H3(n)， f( ，)+ ( “， )一 (_厂，)， 1u(X，0)一 (x)． * 收稿 日期 ：2006—10—20 基金项 目：国家 自然科学基金资助项 目10671184) 维普资讯 660 应 用 数 学 2007 为方便起见 ，设 为X—Y平面上边平行于坐标轴的区域 ， 为0上的一族矩形剖分 VK ∈丁^，记K 的中心点为 (xK，YK)，平行于z，Y轴的边长分别是 2h，2h 且不妨设h 》h，h — m axdiam(K)，其四顶点分别为z (zK—hx,K一^)，Z2(zK+hz，K一^)，Z3(zK+^z， K∈ YK+hy)及 Z (xK—h，YK+hy)，四边为 z一—ZiZ—~l(i一 1，2，3，4，Z 一Z)．记 R一 [一1， 1]×[一1，1]为 平面上的参考元，四顶点坐标依次为三(一1，一1)，三。(1，一1)，三。(1，1)及 三(一1，1)，其四边为2一 (一1,2，3，4三 一三)．那么存在R到K的仿射变换 — K： {x一=h^~++xK．’ 记Qz(R)一span{1，，，， ，矿， ，。，矿}．那么，存在唯一插值算子疗：Hz(R)一 Q