两层网格算法和多重网格算法的应用.pdf

摘 要 在现代科学计算中，寻找和设计简单高效的数值算法变得非常重要。在求 解椭圆边值问题的方法中，由Xu【112】提出的两层网格方法是一种非常实用的 有限元离散方法。借助于两层网格子空间，对于复杂的问题(例如非线性或非 对称不定问题1只需花很小的代价在粗网格空间求解原问题，而在细网格空间 上求解相对简单的问题。因此两层网格离散方法可以减少很多计算量。 对于求解由有限元方法离散椭圆问题得到的线性方程，多重网格方法是一 种最优的数值方法。其收敛速度与网格的步长无关且总体工作量为0(Ⅳ)，其 中Ⅳ表示问题的未知量个数。多重网格方法已成功地应用于流体计算、结构力 学计算和高度非线性的半导体器件模拟计算等领域。 本文在前人工作的基础上继续讨论了两层网格离散方法和多重网格方法。 主要内容共分为四章。在第一章中，首先介绍了有限元离散方法，并通过两个 简单问题描述了两层网格离散方法。然后简单介绍了多重网格的发展历史和现 状，并列举了常用的光滑子。 在第二章中，针对非线性椭圆问题，基于修正牛顿迭代法给出了新的两层 网格算法。在细网格空间上求解非线性问题简化成在粗网格空间求解两个问 题(一个线性和一个非线性)和在细网格空间上求解两个具有相同系数矩阵的线 性问题。理论分析表明粗网格步长日和细网格步长h在日1(Q)范数下可选取 小的情况下，仍能得到问题的最优解。基于Bai等|9】提出的HSS迭代法，针对 非对称和非线性椭圆方程分别给出了新的两层网格算法。在细网格空间上求解 原问题简化成在粗网格空间上求解原问题和在细网格空间上求解对称正定问题 和带位移的反对称问题。收敛性分析表明，在一定的条件下，粗网格空间相对 很小的情况下仍能得到最优精确解。数值例子表明新的算法具有很好的数值效 果。 第三章主要讨论了HSS迭代法做光滑子时多重网格算法求解非对称线性椭 圆问题。基于HSS迭代方法，给出了修改的加性和乘性光滑子。利用扰动分析 的方法，在粗网格步长充分小的条件下，证明了基于这些新的光滑子的多重网 格算法的收敛性。数值实验表明在这些新的光滑子作用下，多重网格算法具有 ii 两层网格算法和多重网格算法及其应用 粗网格步长的依赖性更小。 第四章主要讨论了自适应多重网格方法求解复系数二阶椭圆问题。近年 来，随着电磁场问题研究的不断深入，产生了一些重要的复系数问题，例如时 谐散射和辐射问题。以及由完美匹配层(PML)方法应用到散射问题和共振问题 所产生的复系数问题。虽然这些实际问题比椭圆问题复杂，但是研究复系数椭 圆问题是研究这些问题的基础f60】。在粗网格步长充分小的条件下，给出了基 于Gauss．Seidel光滑子的自适应多重网格算法的收敛性证明。数值实验表明该 算法对于求解复系数椭圆问题的有效性。 关键词：两层网格方法，多重网格方法，HSS迭代方法，光滑子，自适应有 限元。 Abstract Inmodernscientific is to and computing，itveryimportantinvestigatesimple efficientnumerical levelmethod akindof proposedby algorithm．Two Xu【112】is methodfor value efficientfiniteelementdiscretizationsolvingellipticboundary ontwo asnonlinear subspaces，forcomplicatedproblem(such problems．Based and andindefinite