四阶强阻尼波动方程新混合元模式的高精度分析.pdfVIP

应 用 数 学 MATHEMATICA APPLICATA 2015，28(2)：368—377 四阶强阻尼波动方程 新混合元模式的高精度分析 王萍莉 ，石东伟 ，王芬玲 (1．许 昌学院数学与统计学院，河南 许 昌461000；2．河南科技学院数学系，河南 新乡453000) 摘要：本文针对 四阶强阻尼波动方程研究一种新混合元逼近格式．基 于双线性元 Q 。 及其梯度空间Q0×Q。。的高精度分析 ，并借助于插值后处理技术，在半离散和全离 散格式下，分别导 出原始变量 M在 H 模和 中间变量 在L。模意义下相应的超逼近 性质及超收敛结果． 关键词 ：四阶强阻尼波动方程 ；新混合元格式；双线性元；超逼近及超收敛 中图分类号 ：0212．6 AMS(2000)主题分类 ：65N30；65N15 文献标识码 ：A 文章编号 ：1001—9847(2015)02—0368—10 1．引言 本文考虑如下半线性强阻尼波动方程 的初值问题 r ～ ·( “+V 十 )一 ，()， (z，Y，)∈力×(0，T3， u(x，Y，￡)一0， (z，，￡)∈an×(0，T2， (1．1) Iu(x，Y，O)一甜o(z，)， (z，Y，0)===1(，)， (z，)∈n， 其中 c R 是z—Y平面上具有 Lipschitz连续边界a 的有界凸区域 ，0T∞， 与 · 分别表示梯度与散度算子，U。(，)，乱(z，)为已知函数，厂(“)关于 满足Lipschitz连续，即 存在常数 L 0，使得 lf(u1)一f(u2)I≤ LIl一 “2I，V l，M2∈R． 问题 (1．1)是一类含有时间和空间混合偏导数的双曲型方程，并有一些数值方法 的研究， 例如 ，文r1]利用 Ritz—Volterra投影的技巧 ，考虑了一维和二维情况下问题 (1．1)的H Galer— kin混合元方法的半离散和全离散格式 ，给出了两个变量在 L 模和 H 模意义下相应的最优 阶误差估计．文[2]研究了该方程的最低阶的Raviart—Thomas三角形混合元控制体积法，得到 了两个变量在L 模和 H(div；n)模意义下的最优阶误差估计．文[3]将文[13的工作进一步深 化 ，研究了问题 (1．1)的 H Galerkin混合元方法的超收敛性 ，直接采用插值算子，在半离散 和全离散格式下均得到了比文[1]高一阶的超逼近结果．文E4]研究了方程(1．1)的非协调元方 法，并在半离散和全离散格式下给出了 “在H 模意义下的超逼近性质． * 收稿 日期 ：2014—05—28 基金项 目：国家 自然科学基金 (1l101381 作者简介：王萍莉，女，汉族，河南人，讲师，研究方向：有限元方法及应用 第 2期 王萍莉等 ：四阶强阻尼波动方程新混合元模式的高精度分析 369 众所周知，传统混合有限元方法所涉及的两个有限元逼近空间通常要求满足所谓的B-B 条件．实际上，要构造稳定的混合元格式并不是一件容易的事，因此，出现了一些不需要这一限 制条件的有限元格式E]．最近 ，文EloJ针对二阶椭圆问题首先给出了一种新的混合元格式，该 格式只需要试探空间对满足一个简单的包含关系，则 B-B条件 自然成立．基于这一 明显优势， 随后，该方法被广泛地用于解决其它问题[]1-14]，并得到了满意的收敛结果．但到 目前为止，该方 法对方程 (1．1)的应用还未见报道． 本文的主要 目的是基于文ElO3的思想利用双线性元 Q 及其梯度空间 Q0×Q】o对问题 (1．1)建立一种新的混合元格式．借助于高精度分析和插值后处理技术，在半离散 ／全离散格式 下，导 出原始变量 在H 和 中间变量 在L。模意义下具有二阶精度 O(h。)／O(h+rz)的超 逼近结果和整体超收敛结果(h是关于 剖分参数，r是时间步长)． 2．单元构造及混合元格式 设 是 R 中的一个矩形区域 ，其边 a力平行坐标轴， 是n满足正则假设 的矩形单元剖 分．任意K ∈n ，