半线性椭圆问题Petrov-Galerkin逼近及亏量迭代.pdfVIP

2014年 8月 计 算 数 学 第 36卷第 3期 Aug．，2014 MATHEMATICA NUMERICA SINICA Vo1．36，No．3 半线性椭圆问题 Petrov—Galerkin逼近 及亏量迭代术1) 司红颖 (商丘师范学院数学系，河南商丘476000) 陈绍春 (郑州大学数学系，郑州450052) 摘 要 本文考虑了二阶半线性椭圆问题的Petrov-Galerkin逼近格式，用双二次多项式空间作为形 函数空间，用双线性多项式空间作为试探函数空间，证明了此逼近格式与标准的二次有限元逼近 格式有同样的收敛阶．并且根据插值算子的逼近性质，进一步证明了半线性有限元解的亏量迭代 序列收敛到 Petrov—Galerkin解． 关键词：Petrov—Galerkin逼近；亏量迭代；插值算子 MR (2000)主题分类：65N30，65N15 1．引 言 有限元方法是求解各种微分方程的一种好的离散化方法．但有时用这种方法得到的有限 元解的精度并不高，特别是对一些奇异问题，因此有限元的高精度算法也变得越来越重要．目 前，有限元的高精度算法主要有超收敛，外推及校正等．但对有限元的校正至今研究的还比较 少，而有限元的校正法又显示了许多优点：l、它比差分法的校正简单，由于使用分部积分，降 低了应逼近的阶次，在实际计算时采用插值校正格式更是简单统一；2、校正法的证明虽然比 外推法复杂些，但用直接法求解时计算量却很少，而且数值例子表明它的精度更高，因此，有 限元的校正是很有前景的高精度算法．Frank等在 11【]中建立了椭圆边值问题的有限元亏量迭 代校正格式，对线性椭圆边值问题来说，在拟一致剖分条件下，文 (2-6]应用超收敛和渐进展开 式 已经讨论了这种格式的有效性．证明了给出的线性有限元解作为第一次迭代校正的原始逼 近解，能够获得与标准的二次有限元有相同的收敛阶．然而，如果剖分是分片拟一致，则在原 始光滑性条件下，标准二次有限元的逼近阶，比被首次迭代校正获得的逼近阶低．文 2【-3]中指 出经过多次迭代校正可以弥补损失的亏量．对两点边值问题以及二阶椭圆问题，文 [7-9]已经 证明了有限元解的亏量迭代校正收敛到 Petrov-Galerkin近似解．在本文中我们给出了半线性 椭圆边值问题的亏量迭代校正格式，并对文 [9]中插值算子的收敛性结论进行了改进，结论简 单且比 {9l中结果较优． 本文首先讨论了插值算子的收敛性质，其次考虑半线性椭圆边值问题，建立了Petrov— Galerkin逼近格式，最后证明了半线性有限元解的亏量迭代校正均匀的收敛到Petrov—Galerkin 解． 2013年 l1月 15日收到． )基金项 目：国家自然科学基金 (No；河南省 自然科学基金 (No．132300410272) 3期 司红颖 等：半线性椭圆问题 Petrov—Galerkin逼近及亏量迭代 317 2．插值算子的性质 令 T是矩形剖分，长为 2hl，宽为 2h2，对矩形 通过连接四边中点，得到四个子矩形 ， ， ， ， 如图： gt5(Xo，Yo) 乃 定义 = ＼／／(()。+( 令 hu是双二次拉格明 日捕值多项式，定义在 ’的节点 ai(i=1，．．．，9)上，Ilu是定义在四个 子矩形 ， ， ，T4上的分片双线性插值多项式．则有