四阶抛物方程H1-Galerkin混合有限元方法的超逼近和最优误差估计.pdfVIP

2014年 11月 计 算 数 学 第 36卷第 4期 NOV．，2014 MATHEMATICA NUMERICA SINICA Vo1．36，No．4 四阶抛物方程 H1．Galerkin混合有限元方 法的超逼近及最优误差估计术1) 石东洋 (郑州大学数学与统计学院，郑州450001) 史艳华。) 王芬玲 (许昌学院数学与统计学院，河南许昌461000) 摘 要 本文基于双线性元及零阶 Raviart—Thomas元 ( )对四阶抛物方程建立了半离散和向后 欧拉全离散 日一Galerkin混合有限元格式．利用积分恒等式技巧和单元的特殊构造，证明了关于 上述两元的两个新的重要性质．进而导出了这两种格式下相关变量的最优误差估计和超逼近性质． 关键词：四阶抛物方程；H1Galerkin混合有限元方法；半离散和全离散；误差估计及超逼近 MR (2000)主题分类：65N15，65N30 1．引 言 X u=f( ， ， 其中QCR 是一个有界凸多角形区域，具有 Lipschitz边界 aQ，X=(X，)，f(X，t)和uo(X) 是已知的光滑函数，k为正常数． 四阶问题有深刻的物理背景，它可以用于描述梁的静态运动或刚体运动 _lJ，一维弹塑性 杆的纵振动问题和二维反平面剪切 问题 2【]及人 口问题中的弥散和增长 3[J等．4【,5]讨论了一 类四阶抛物方程解的存在唯一性．[6】研究了四阶椭圆和抛物问题的混合协调有限元方法，导 出了最优收敛阶．[7]对四阶椭圆问题构造了三个 非协调元，证明了其分别具有一阶和二阶 收敛性．[8]对于位移障碍下四阶变分不等式问题给出了非协调有限元方法的一般估计式．[9] 建立 了双调和方程混合元的一种新格式，通过对流函数和涡函数选取合适的逼近多项式空间 得到了与Ciarlet—Raviart格式相同的收敛阶．[10】讨论了问题 (1．1)的H -Galerkin混合元方 法，利用 Riesz投影导出了一维情形半离散格式和向后欧拉全离散格式下的 。和 日 模最优 收敛阶估计及二维情形半离散格式下的误差估计式，但没有达到最优，更未涉及超逼近研究． 2013年9月8日收到． )基金项 目：国家 自然科学基金 11271340；河南省教育厅资助基金 (14Al10009)；许 昌学院青年骨干教师项 目． 。)通讯作者：史艳华，Email：syhsdq@163．com． 364 计 算 数 学 H 一Galerkin混合元方法最早由Pani于 1998年提出并应用于抛物方程 _llJ．与传统的混 合元方法相比，它具有不需要满足 LBB相容性条件，两个逼近有限元空间可以自由选择等优 点．随后，Pani又将其应用于抛物积分微分方程 [12]和双曲方程 1【3]．[14—16]又分别将该方法 用于正则长波方程，双曲积分微分方程和伪双曲方程．[17]和 [18】采用与 1【1—16]不同的技巧， 即用插值代替传统的Riesz投影或Ritz—Volterra投影，分别研究了Sobolev方程和四阶强阻 尼波动方程，前者得到了与传统有限元分析相 同的误差估计，后者得到比传统结果高一阶的超 逼近和超收敛结果． 本文主要借助双线性元和零阶R-T元对问题 (1．1)进行了半离散和向后欧拉全离散 日一 Galerkin混合元方法逼近．首先利用积分恒等式技巧和单元的插值性质证明了两个新的重要 引理 (