两层网格算法和多重网格算法及其应用

来源 :浙江大学 | 被引量 : 0次 | 上传用户:gmgan
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
在现代科学计算中,寻找和设计简单高效的数值算法变得非常重要。在求解椭圆边值问题的方法中,由Xu[112]提出的两层网格方法是一种非常实用的有限元离散方法。借助于两层网格子空间,对于复杂的问题(例如非线性或非对称不定问题)只需花很小的代价在粗网格空间求解原问题,而在细网格空间上求解相对简单的问题。因此两层网格离散方法可以减少很多计算量。   对于求解由有限元方法离散椭圆问题得到的线性方程,多重网格方法是一种最优的数值方法。其收敛速度与网格的步长无关且总体工作量为O(N),其中N表示问题的未知量个数。多重网格方法已成功地应用于流体计算、结构力学计算和高度非线性的半导体器件模拟计算等领域。   本文在前人工作的基础上继续讨论了两层网格离散方法和多重网格方法。主要内容共分为四章。在第一章中,首先介绍了有限元离散方法,并通过两个简单问题描述了两层网格离散方法。然后简单介绍了多重网格的发展历史和现状,并列举了常用的光滑子。   在第二章中,针对非线性椭圆问题,基于修正牛顿迭代法给出了新的两层网格算法。在细网格空间上求解非线性问题简化成在粗网格空间求解两个问题(一个线性和一个非线性)和在细网格空间上求解两个具有相同系数矩阵的线性问题。理论分析表明粗网格步长H和细网格步长h在H1(Ω)范数下可选取为H=O(h(λ/7)),在L2(Ω)范数下可选取为H=O(h(1/4))。因此在粗网格空间非常小的情况下,仍能得到问题的最优解。基于Bai等[9]提出的HSS迭代法,针对非对称和非线性椭圆方程分别给出了新的两层网格算法。在细网格空间上求解原问题简化成在粗网格空间上求解原问题和在细网格空间上求解对称正定问题和带位移的反对称问题。收敛性分析表明,在一定的条件下,粗网格空间相对很小的情况下仍能得到最优精确解。数值例子表明新的算法具有很好的数值效果。   第三章主要讨论了HSS迭代法做光滑子时多重网格算法求解非对称线性椭圆问题。基于HSS迭代方法,给出了修改的加性和乘性光滑子。利用扰动分析的方法,在粗网格步长充分小的条件下,证明了基于这些新的光滑子的多重网格算法的收敛性。数值实验表明在这些新的光滑子作用下,多重网格算法具有最优复杂性,与阻尼Jacobi光滑子和Gauss-Seidel光滑子相比,HSS光滑子对粗网格步长的依赖性更小。   第四章主要讨论了自适应多重网格方法求解复系数二阶椭圆问题。近年来,随着电磁场问题研究的不断深入,产生了一些重要的复系数问题,例如时谐散射和辐射问题。以及由完美匹配层(PML)方法应用到散射问题和共振问题所产生的复系数问题。虽然这些实际问题比椭圆问题复杂,但是研究复系数椭圆问题是研究这些问题的基础[60]。在粗网格步长充分小的条件下,给出了基于Gauss-Seidel光滑子的自适应多重网格算法的收敛性证明。数值实验表明该算法对于求解复系数椭圆问题的有效性。
其他文献
期刊
本文主要研究了两类特征值优化问题:一类是关于区域密度分布的约束优化问题,一类是关于边界控制的约束优化问题。针对这两类不同的问题,本论文提出了贪婪算法、单调算法以及边界
乔多(Ambrogio ou Angiolotto di Bondone Giotto 1266—1330)可说是基督教者圣法朗梭阿大西士(Saint Fran□ois d’Assise)的历史画家。他一生重要的壁画分布在三所教堂中,
在这篇论文里,我们研究了几类平面含参数的分段光滑系统的滑动分支现象。通过运用分段光滑系统的Filippov理论和微分方程定性理论,我们系统地研究了两类分段光滑线性二次系统
分支控制是指通过设计一个控制器去改变非线性系统的特性来实现一些理想的动力学行为。传统的控制方法包括推延分支的产生和在合适的位置引入新的分支等。此外,分支控制的具体
最优化理论是应用数学的一个重要组成部分,它在实际生活中有着广泛的应用。研究最优化的方法很多,17世纪牛顿和莱布尼茨在他们所创建的微积分中,形成变分法。1947年Dantizig
城市用电量预测在城市规划建设和城市供电系统调度中具有重要的作用。准确、合理的预测城市用电量不仅能直接指导发电厂的生产,更能为电厂间的优化调度提供可靠的技术支持。本
近些年,非线性弹性理论和拟共形映射的发展促使微分形式椭圆方程的研究取得了极大的进展,已经从最初的Laplace方程扩展到了A-调和方程。Hodge-Dirac算子的发展来源于理论物理学
曲线曲面是计算机辅助几何设计(CAGD)中的主要研究对象,CAGD中的大多数操作都是以曲线曲面为基础的.而不论是参数曲线或者数据点的逼近问题,还是形状可控的曲线曲面的设计问题,它
现代医用CT希望减少X射线的使用量,从而减少对病人的伤害,但是同时又希望得到质量较高的图像,而这就需要好的算法.当前针对这个问题的研究主要是从预处理与后处理两个方面来