求解Helmholtz方程的组合紧致差分格式及其快速预处理算子

来源 :华南师范大学 | 被引量 : 0次 | 上传用户:ustczl
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文第一次成功地将组合紧致差分(CCD)格式应用于离散一维和二维Helmholtz方程,并显式给出了对应的用矩阵表示的非对称CCD线性系统和借助合适的边界条件实现了理论上的高阶收敛精度.特别地,对于带有Neumann边界条件的Helmholtz方程,还提出了一组新的高阶CCD边界格式用于提高边界的逼近精度.总体而言,本文所给出的CCD格式对于一般边界情况至少达到了六阶精度.   对于一维Helmholtz方程,先从理论上讨论了CCD离散线性系统满足非奇异性的充分必要条件,然后提出并分析了一个新的块下三角预处理算子用于迭代求解所考虑的CCD线性系统.同时,注意到CCD线性系统所带有的鞍点结构,本文还考虑了一类基于Schilders’分解的约束预处理算子作为比较.数值结果表明本文提出的块下三角预处理算子要明显快于现有的约束预处理算子.更加重要的是,基于块下三角预处理算子的GMRES迭代算法在数值上实现了与网格规模无关的收敛速度,对块下三角预处理算子的特征值分析也部分地证实了这一数值现象.对于规模为n的CCD线性系统,块下三角预处理算子的预处理方程求解仅需()(n)的计算量.   对于二维Helmholtz方程,借助Kronecker积和一维CCD线性系统简洁地给出了二维CCD线性系统的矩阵表示,同时,为了保证离散系统的非奇异性,四个CCD边界公式用于替换四个角点的原边界方程.为了构造针对二维CCD线性系统的块下三角类型预处理算子,一种创新的矩阵拉伸技巧用于重组原系数矩阵.同样,本文还构造了基于Schilders’分解的约束预处理算子用于比较.数值结果显示本文提出的块下三角类型预处理算子比现有的约束预处理算子要更加有效,这与理论上的特征值分析是吻合的.数值上也观察到了与网格规模无关的收敛速度,这一特点对于求解超大规模问题特别重要.   最后,本文得到了一些结论并给出了部分后续研究的建议.
其他文献
本博士后报告从数学角度研究了气体动力学中几种含有亚音速流及跨音速激波的特殊流动模式的唯一性。这些流动模式包括:   ●三维无限长扩张管道内定常亚音速可压缩位势流;
本文主要分两部分。第一部分系统研究了一阶线性ODE系统的单值同构方法,研究了Painlevé方程的对应的一阶线性ODE系统的单值同构形变,导出Lax对,阐明了Painlevé方程与单值同构
本文研究4维球面S4到CPn的常曲率弱Lagrangian极小浸入(ρ):S4→CPn,其诱导度量ds2具有常曲率c.文中证明了存在一个整数s≥1使得c=4/[S(S+3)].浸入(ρ)被两个四元齐次多项式f
数学模型对描述种群增长起着重要的作用,近年来,在对描述生物种群增长的模型的研究也已经取得了大量的结果。这就意味着由随机微分方程描述的生物模型对现实世界中应用的重要性
无网格Galerkin方法(EFG)是近年来迅速兴起的一种数值方法。该方法构造形函数时只需要具体的节点信息,而不要网格,因此可显著减少因网格畸变带来的困难。在涉及大变形、自由
经典的Besov空间和Triebel-Lizorkin空间在偏微分方程的研究中起了非常重要的作用。J.Bourgain,T. Tao,C.E.Kenig,T.Kato等人将它们运用到非线性发展方程的研究中,获得了令人瞩
函数空间上的算子理论已成为人们研究的热点,由于研究的载体是函数空间,所以这些常见的算子必是由某些函数导出的,从而我们需要探讨这些算子的性质与它们的诱导函数有怎样的
图G的一个正常k-全染色是指一个映射φ:V(G)∪E(G)→{1,2,…,k},使得V(G)∪E(G)中任意两个相邻或关联的元素染不同的颜色.图G的全色数x"(G)是指图G有一个k-全染色的最小k值.令Cφ
本文主要研究了算子的q-数值域,这是对算子的数值域的结论的推广。我们首先给出了算子的q-数值域及q-数值域半径的定义,介绍了算子的q-数值域的基本性质,并给出算子的q-数值域的
我国是世界上种植转基因农产品规模比较大的国家,但主要是在棉花而不在食品上。我国进口的转基因农产品主要集中在大豆上,因为我们自身大豆的产量是1300万~1400万t,但是我们需