四阶抛物偏微分方程的H1Galerkin混合元方法及数值模拟.pdfVIP

2012年8月 计算数学 第34卷第3期 Aug．，2012 MATHEMATICANUMERICASINICA V01．34，No．3 方法及数值模拟术1) 刘 洋李宏 何斯日古楞 高 巍方志朝 (内蒙古大学数学科学学院，呼和浩特010021) 摘 要 到目前为止，Hl_Galerkin混合有限元方法研究的问题仅局限于二阶发展方程．然而对于高 阶发展方程，特别是重要的四阶发展方程问题的研究却没有出现．本文首次提出四阶发展方程的 日1一Calerkin混合有限元方法，为了给出理论分析的需要，我们考虑四阶抛物型发展方程．通过引 进三个适当的中间辅助变量，形成四个一阶方程组成的方程组系统，提出四阶抛物型方程的日1一 Galerkin混合有限元方法．得到了一维情形下的半离散和全离散格式的最优收敛阶误差估计和 多维情形的半离散格式误差估计，并采用迭代方法证明了全离散格式的稳定性．最后，通过数值例 子验证了提出算法的可行性．在一维情况下我们能够同时得到未知纯量函数、一阶导数、负二阶 导数和负三阶导数的最优逼近解，这一点是以往混合元方法所不能得到的． 关键词：四阶抛物偏微分方程；H1一Galerkin混合元方法；稳定性；最优阶误差估计 MR(2000)主题分类：65N30，65M60 1．引 言 考虑四阶抛物方程的初边值问题 U+ z，￡)，(X，t)∈Q×Z 让 ，o t)=0，(z，t)∈aQ×歹 (1．1) ，●●●J、●●●．．．L 让r@@．q 心归旷 一舭姒 ，∽功 z∈Q． 这里QC 函数，k为正常数． 四阶问题在很多方面得到了应用．由于四阶方程具有高阶空间导数的特点，数值求解有一 定困难．对于四阶问题的研究主要包括有限元方法、有限差分方法、有限体积元方法、混合 有限元方法[1-9】等．文献【2]提出了四阶椭圆问题的扩展混合有限元方法，同时可以得到三 个变量的误差估计结果．文[3】通过使用径向基函数的方法研究了四阶椭圆和抛物问题的混 合有限元方法．文【4】研究了四阶线性椭圆问题和四阶线性抛物问题的最优收敛阶误差估计． 文[5】利用时间间断而空间连续的时空有限元方法研究了四阶抛物型积分微分方程．文[7]研 的最优阶误差估计和导数逼近的超收敛结果．文【8]研究了四阶障碍问题的稳定混合元方法． 文【9】研究了在各向异性网络剖分下四阶特征值问题的双三次Hermite有限元逼近． +2011年8月1日收到． 万方数据 260 计算数学 混合有限元方法的研究备受国际计算数学界的关注，特别是比传统混合有限元更具有优 法．该方法克服了传统混合有限元必须满足LBB相容性条件的缺点，允许逼近有限元空间％ 和Wh具有不同次数的多项式，进而空间的选择比较自由．基于该方法的诸多优点，得到了足 够的重视，并将该方法相继应用到了抛物型和双曲型积分微分方程[11-13】、双曲和伪双曲型方 阻尼波动方程【20]等．但到目前为止，所有研究的问题都是具有二阶空间导数问题，而对于四 阶空间导数的问题研究还没有看到． 元格式．分析一维情况的半离散和全离散格式的最优收敛阶误差估计，采用迭代方法证明了 全离散格式的稳定性，并得到了多维情形的半离散格式误差估计的证叽最后通过数值实验的 结果验证了所提出算法是可行的． 为了研究问题的需要，首先给出一些记号．L2(Q)内积用(·，·)表示，令础={u∈H1v(o) 时，把w“，2记作日”，其范数简记为¨．‰．注意在本文中作估计时所有的C都是与空间网 格参数^和时间步长△t无关的正常数． 2．H1一Galerkin混合有限元格式 f札t+ku。。。。=，(z，t)，(z