求解一般椭圆形方程不标记震荡项自适应方法.pdfVIP

2009年 l2月 高 等 学 校 计 算 数 学 学 报 第 31卷第 4期 求解一般椭圆形方程不标记 震荡项的自适应方法 邹青松 (中山大学数学与科学计算学院，广州 510275) ADAPTIVE M ETHoDSW ITHOUT M ARKING OSCILLATION FOR GENERAL ELLIPTICEQUATIONS ZouQingsong (SchoolofMathematicsandScientificComputing，ZhongshanUniversiY Guangzhou510275) Abstract In this note，we constructand analyze an adaptivefinite element method(AFEM)forsolvinggeneralellipticequations．Inthemarkingstepof ouradaptiveprocedure，weneednottomarkthegridstoberefinedaccordingto theoscillationterm．WeprovethatthisAFEM hasthesameconvergenceratesas thestandardAFEM proposedinl12，16，15I． Key words finiteelementmethods，ellipticequations，oscillation． AMS(2000)subjectclassifications 65N15，65N06 中图法分类号 024 士 l 置 国家 自然科学基金项 目资助，基金号 收稿 日期：2007—02—15． 320 · 邹青松：求解一般椭圆形方程不标记震荡项的自适应方法 第 4期 求解偏微分方程的自适应算法起始于 20世纪 70年代晚期，现在已经成为了科学和 工程中的通用工具．一般地，一个 自适应算法包括如下四个步骤： 方程求解— 误差估计—+单元标记— 网格加细 本文将不讨论第一步方程求解，因为迄今为止已有大量关于这方面的工作．第二步 的误差估计指的是后验误差估计，在这一步我们需要构造可计算的指示因子去指示精确 解和近似解的局部误差．在已过去的数十年中所构造的指示因子主要有三类： 基于残差的误差指示因子、基于复原技术的误差指示 因子 和 分层基误差指示因子． 本文所要讨论的是基于残差的误差指示因子．线性椭 圆型方程的残差指示因子的构造可 以追索到Babu~ka等人 1978年的工作 2【,3]．用适当局部测试函数去分析局部误差下界 的方法起始于 Verfurth[．关于基于残差的误差估计的综述有 Verfiirth [221、Ainsworth 和 Oden[．在第三步单元标记 中，我们通常必须使用一个所谓的标记策略 (或所谓的自 适应策略)，该策略对整个 自适应算法的收敛性和收敛阶非常重要．通常，标记策略都是 基于著名的误差等分布原理 6【，l1，141．这一原理的基本思想是，最恰当的网格应该是每一 单元上的局部误差相等 (或基本相等)．在第四步网格加细中，有很多的网格加细技术被 使用．这些技术包括最新节点对分技术 [20]，最长边对分技术 1【8】和红绿加细技术 [4】等． 自适应算法构造后，我们还需对它进行收敛性和收敛阶等进行分析．D6rfler分析了二 维Poisson方程的线性元的自适应算法的收敛性 (cf．，1【2】)．Nochetto及其合作者在一系列 文章中分析了线性椭圆型