抛物型方程一个新的非协调混合元超收敛性分析和外推.pdf

2013年 11月 计 算 数 学 第 35卷第 4期 NOV．，2013 MATHEMATICA NUMERICA SINICA V_01．35．No．4 抛物型方程一个新的非协调混合元 超收敛性分析及外推冰1) 石东洋 一 (郑州大学数学与统计学院，郑州450001) 张亚东 ， =，～一 (许昌学院数学与统计学院，河南许昌461000) ， 、 ， 摘 要 本文研究了抛物型方程在新混合元格式下的非协调混合有限元方法．在抛弃传统有限元分析 的必要工具 一Ritz投影算子的前提下，直接利用单元的插值性质，运用高精度分析和对时间t的 导数转移技巧，借助于插值后处理技术，分别导出了关于原始变量 乱的 日一模和通量 =一Vu 在 一模下的O(h)阶超逼近性质和整体超收敛．进一步，通过构造合适的辅助问题，运用 Richardson外推格式，得到了具有更高精度 O(h)阶的外推结果．最后，给出了一些数值结果验 证了理论分析的正确性． 关键词：抛物型方程；非协调元；新混合元格式；超收敛；外推 MR (2000)主题分类：65N15，65N30 1-引 言 考虑如下抛物型方程的初边值问题 在QX(0，T】中， 在 Q×(0，T]上， (1．1) 在 Q中， 其中QCR 是具有 Lipschitz连续边界的有界凸多边形区域，aQ为 2【的边界，T∈(0，O0)为 一 定值，X=(，Y)，f(X，t)及 uo(X)为已知的光滑函数． 标准的有限元空间对逼近解的光滑性要求很高，这给实际应用带来了很多困难，于是产生 了混合有限元方法．利用混合有限元方法有很多优点，例如在计算多孔介质流时，通常需要计 算速度，如果用通常的有限元法，只能先求出压力，然后求导得到速度，这样做精度将降低，而 利用混合有限元法求解，可同时求出压力和速度，提高了离散解的精度．此外，许多问题本身 自然的Galerkin逼近就只能采用混合有限元方法，例如 Stokes问题的有限元逼近．其所涉及 的两个逼近空间通常需要满足所谓的inf-sup条件或LBB条件，空间选择受限制．对抛物型方 程 (1．1)的混合有限元方法有许多研究 [1--8J，其中，5【】利用Raviart—Thomas元讨论了L2-模 2012年 10月 7日收到． )基金项 目： 国家 自然科学基金 11271340；高等学校博士学科点专项基金 (2009410111006)；