自适应间断Galerkin 有限元的多水平方法.pdfVIP

中国科学: 数学 2012 年 第42 卷 第5 期: 409 428 自适应间断 有限元的多水平方法 卢培培 石钟慈 许学军∗ 中国科学院数学与系统科学研究院, 北京 100190 E-mail: lupeipei@lsec.cc.ac.cn, shi@lsec.cc.ac.cn, xxj@lsec.cc.ac.cn 收稿日期: 2011-12-19; 接受日期: 2012-03-13; * 通信作者 国家自然科学基金(批准号: 和国家重点基础研究专项经费 (批准号: 2011CB309701) 资助项目 摘要 本文讨论在自适应网格上间断 Galerkin 有限元离散系统的局部多水平算法. 对于光滑系数和 间断系数情形, 利用 Schwarz 理论分析了算法的收敛性. 理论和数值试验均说明算法的收敛率与网格 层数以及网格尺寸无关. 对强间断系数情形算法是拟最优的, 即收敛率仅与网格层数有关. 关键词 间断有限元 局部多水平算法 间断系数 悬点 主题分类 65F10, 65N30 引言 间断Galerkin 有限元方法可以灵活地求解很多实际问题, 也能处理有悬点的网格. 近年来, 此方 法已受到越来越多的数值分析家的重视. 对二阶椭圆问题, Arnold 等在文献 [1] 中提出了一个一般性 的框架, 并在此框架下分析了几类间断Galerkin 有限元方法的稳定性和收敛性. 在 [2] 中, Hoppe 证明 了自适应内部惩罚间断Galerkin (IPDG) 有限元方法的收敛性, 其中自适应网格是根据后验误差估计 子 (详见 [3]) 利用最新顶点法加密而得. 然而在离散代数方程求解方面却没有太多的研究. 在 [4] 中, 作者在网格一致加密的情况下对扩散问题提出了一个一致有界的多重网格算法. 文献 [5] 针对间断系 数问题提出了一个多水平预条件子. 文献 [6] 讨论了二阶椭圆问题内部惩罚间断 Galerkin (IPDG) 有 限元方法的多水平预条件子, 证明了当网格正则并且每边至多有一个悬点时条件数一致有界. 在文献 [7–9] 的框架下, 重叠型区域分解和多重网格均可看成一类子空间校正方法. 本文在此框 架下从空间多水平分裂的角度出发, 通过构造一个间断 Galerkin 有限元空间的能量稳定的空间分解, 设计了相应的局部乘性和加性算法, 并利用Schwarz 理论对算法的收敛性进行了分析. 本文第2 节简要介绍了模型问题. 第3 节给出了抽象的局部多水平算法及其理论. 第4 节在光滑 系数情形下, 分别对两种网格剖分, 逐一验证第3 节中的假设, 证明了算法的一致收敛性. 在第5 节考 虑间断系数情形, 证明了算法的拟一致收敛性. 最后一节给出了几个数值算例验证理论分析的结果. 模型问题 本文采用标准的Lebesgue 和 Sobolev 空间记号 (详见 [10]). 记 (·, ·) 为L2 (Ω) 上的内积, ∥ · ∥1Ω , | · |1Ω 分别是Sobolev 空间H 1 (Ω) 上的全范和半范. 用A B 表示A CB , 其中C 是一个只和网格 英文引用格式: Lu P P, Shi Z C, Xu X J. Local multilevel methods for adaptive discontinuous Galerkin finite element methods (in Chinese). Sci Sin Math, 2012, 42(5): 409–428, doi: 10.1360/012011-1004 卢培培等: 自适应间断Galerkin 有限元的多水平方法 的正则性有关的正常数. A ≈ B 即A B A. 为方便起见, 下文中的C 表示仅和网格的正则性有关 的正常数, 它们的值并不一定相同. 设求解区域Ω ⊂ (d = 2, 3) 是一个有界的Lipschitz 多边形(多面体) 区域, 考虑下述二阶椭圆 边值问题   −div(ρ()∇u) = f, 在Ω 内,