多尺度问题的守恒型区域分解法

来源 :中国科学院大学 | 被引量 : 0次 | 上传用户:lowsong1
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
在现代科学与工程计算中,区域分解算法已成为求解偏微分方程的重要工具.本文主要工作是设计了多尺度问题的有限体积法和混合元方法对应的离散代数系统的区域分解算法,包括算法的设计,理论分析以及数值模拟.具体可分为以下三个部分:  第一部分首先给出了Sobolev空间及向量有限元函数等一些预备知识,并以多尺度问题为例,简单介绍了有限体积法和混合元方法的离散格式和误差估计理论.其次介绍了两层重叠型区域分解算法的基本框架以及P-GMRES和P-MINRES这两种线性方程组迭代解法.  第二部分研究了多尺度问题有限体积法的两层区域分解算法.首先基于延拓算子,限制算子和校正算子构造了新的粗校正矩阵,并利用该矩阵以及一层的重叠型Schwarz算法构造了细尺度有限体代数系统的预条件子.然后在对偶粗网格上利用一些已有的多尺度方法具体构造了不同种的延拓,限制和校正算子.数值实验表明这些新的粗校正矩阵比文献中已有的算法更稳定,更适合处理复杂形状的间断系数.算法中的后处理迭代技巧可使得迭代解在每个粗网格单元上质量守恒,方便了速度场的局部重构.数值实验验证了几类粗校正矩阵的有效性.  第三部分设计了多尺度问题混合元方法的两层区域分解算法.原始鞍点系统经过等价变换后可得到一个新的鞍点系统,该系统可以用一个2×2的块对角矩阵的逆作为它的最优预条件子.该对角矩阵的第一块求逆对应于一个多尺度H(div)系统的求解,大规模计算下直接法求解是不现实的.针对这个多尺度H(div)系统设计了新的粗空间,证明了粗分解的能量稳定性.我们推导了加权的Helmholtz分解并借此证明了局部分解的稳定性,最后给出了预优H(div)系统的条件数估计,进而得到了预优鞍点系统的正负特征值分布区间的估计.数值实验验证了理论的正确性和有效性.
其他文献
组合测试是一种重要的软件测试方法。它以覆盖表作为测试用例集,能够以少量的测试用例有效地检测由软件各参数之间的交互作用所触发的失效。  如何得到尽可能小规模的覆盖表
学位
《语文课程标准》明确指出,写作教学应引导学生关注现实,热爱生活,表达真情实感.我们应该把作文看成是一种自然的生命状态,一种朴实的生活行为.著名教育家陶行知先生也说:“
期刊
国共合作时期,中国共产党先后派周恩来、恽代英、肖楚女、熊雄、聂荣臻等同志到黄埔军校有效地开展党的思想政治工作。为我党培养了大批优秀的革命人才。 建立党代表和政治
学位
学位
本文共分为五个章节,将致力于给出含小参数的强不定问题解的存在与集中现象的分析刻画。首先,我们将介绍含小参数的强不定问题的背景来源及本文主要结果。在第二章中,我们将给出