【摘 要】
:
多重网格法是求解偏微分方程大规模离散化方程的最为有效的方法,粗略地讲,它可分为几何多重网格法和代数多重网格(AMG)法。这里,我们将结合几何和代数两种途径来研究多重网格法,并称之为所谓的基于几何和分析的代数多重网格法,这是目前国际上代数多重网格法研究领域中新发展起来的方法。 本文分为两部分。第一部分,结合几何和代数多重网格法的特点,我们为两类典型的复杂有限元方程组,设计了具有很强的Robus
论文部分内容阅读
多重网格法是求解偏微分方程大规模离散化方程的最为有效的方法,粗略地讲,它可分为几何多重网格法和代数多重网格(AMG)法。这里,我们将结合几何和代数两种途径来研究多重网格法,并称之为所谓的基于几何和分析的代数多重网格法,这是目前国际上代数多重网格法研究领域中新发展起来的方法。 本文分为两部分。第一部分,结合几何和代数多重网格法的特点,我们为两类典型的复杂有限元方程组,设计了具有很强的Robust性和高效性的代数多重网格法。第一类是Rd,d=2,3中的高次Lagrange有限元方程组,其系数矩阵的具有较强稠密性:另一类是所谓的Criss-Cross网格下的线性有限元方程组的约化线性子系统,其相应的几何粗空间一般不具有嵌套性。通过对问题和有限元空间作深入、细致的分析,发现了许多重要的代数特征,形成了一些新的关于求解复杂有限元方程组的网格粗化技术和提升算子构造的代数方法,从本质上克服了粗网格层自由度难以控制等通常代数多重网格法的缺陷,该类AMG法还具有预处理(Setup)时间少、Robust性好和运算效率高等特性。进一步,通过引入新的证明方法,即利用所谓的Xu-Zikatanov恒等式等,我们从理论上严格证明了新算法的最优收敛(下降)率,数值试验验证了理论的正确性。另外,通过引入两套代数矩阵,对高次Lagrange有限元方程组,我们设计并分析了相应的基于AMG法的预条件共轭梯度法。这些主要的算法设计思想和理论分析方法,具有相当的普适性。 第二部分,我们针对两种应用问题,讨论和分析相应的代数多重网格法。第一种是晶格材料的离散模型。我们首先设计了一种基于AMG法的块预条件共轭梯度法,并就方形晶格模型,利用其近似连续模型,从理论上严格证明了其关于参数α是一致收敛性。接着又构造了对更广泛的晶格模型具有高效性和Robust性的AMG法和相应的APCG法,数值试验表明我们的算法对许多晶格模型,关于其规模和重要参数α是一致收敛的。第二种应用问题来源于辐射流体力学方程组,我们讨论其中的二维三温能量方程离散系统的代数多重网格法。我们针对二维三温能量方程的特殊性,建立了一种半粗化的代数多重网格法(SAMG)和以该SAMG为预条件子的Krylov子空间迭代法,并将其嵌入到能量方程与流体力学方程耦合后得到的应用程序中,通过与经典预条件GEMRES(m)和ORTHOMIN(m)迭代法作对比数值实验,表明我们的AMG方法具有高效性和很好的Robust性。
其他文献
对于非线性扩散占优Volterra型偏泛函微分方程的时间离散化,以往最常用的方法是Volterra泛函微分方程(VFDE)的隐式Euler方法及梯形方法,有时候也使用隐式中点法及Bellen和Zennaro所建议的二级Lobatto IIIC Runge-Kutta方法。为了避免出现Bellen和Zennaro所指出的“order failure”和“stability failure”现象,本文
分数微分方程(Fractional Differential Equations)在许多学科领域有广泛的应用,这些学科领域的许多数学模型都是用分数微分方程来描述的.近十年来,分数微分方程边值问题得到了迅速的发展.随着这一方向研究的进一步深入,研究内容和研究方法不断得到丰富,非线性分数微分方程边值问题获得了许多研究成果.本文主要研究非线性分数微分方程边值问题解的存在性.本文由四章组成,主要内容如下:
本文研究二维各向异性Delaunay三角化网格的生成与优化算法及其应用.当物体在不同方向上的变化速率不同时,例如树木向上生长与横向生长的速率是不同的,我们称其具有各向异性性质.各向异性问题在各向异性晶体材料、多孔介质、石油工程、粘性流模拟、图像处理等科学工程计算中有着广泛的应用.给出一个各向异性问题,进行有限元数值求解时,合适的网格剖分至关重要,选择符合各向异性特征的三角化网格是自然而然的想法.各
有限体积法由于能保持某些物理量(如质量、能量等)的局部守恒性,它已成为一种广泛应用于科学与工程计算领域的重要数值方法.有限体元(FVE)法是一种重要的有限体积法,虽然对它已有大量的研究工作,但仍存在许多需要研究的问题.ICF问题描述的是一种聚变等离子体流体力学问题,其中三温辐射热传导方程的数值求解在该问题的数值模拟中占有主要工作量,因此,对该类问题的有限体元格式及其快速算法进行研究,是一项具有重要
本文主要对Schro(?)dinger方程保持质量守恒的DDG(directdiscontinuousGalerkin)方法进行了研究.在第三章,对线性一维及二维Schro(?)dinger方程的守恒数值流进行了选取,对于其半离散格式和全离散格式进行了守恒性分析.其中,在空间离散上采用DDG方法,在时间离散上采用Crank-Nicolson格式.数值实验表明,对于k次元逼近能得到k+1阶精度,而且
代数多层网格(AMG)法和区域分解法(DDM)是国际上流行的两类求解大规模偏微分方程(PDEs)离散化系统的快速方法.目前,对于求解复杂PDEs离散化系统的AMG法和DDM,还有许多需要进一步研究的问题.本文针对两类具有广泛应用背景的PDEs离散化系统,研究其高效(并行)AMG法和DDM,主要工作如下:针对一种关于三温辐射扩散方程离散化系统的代数两层预条件子B01,测试并分析了相应的PGMRES法
玻色-爱因斯坦凝聚体(BEC)的实现,为研究非线性物质波的动力学行为提供了一个独一无二的平台。尤其是其中暗、亮和矢量孤子的成功观察使BEC中的非线性研究成为当前非线性学科和冷原子物理等交叉领域的研究热点之一。理论上,描述BEC物性是基于平均场近似下的Gross-Pitaevskii方程。该方程中有两个可宏观观测物理量:囚禁外部势阱和描述原子间相互作用的S-波散射长度。实验上,外部势阱可通过外加磁场
有限元渐近展式与外推的研究是有限元高精度理论的一个重要组成部分,是有限元超收敛理论研究的延续和深入,本文就几类基本的三角形网格,对线性元和二次元的渐近展式与外推作出了系统而深入的分析,发现了一些重要的具有一般性质的基本特征,在对被逼近函数光滑性要求比较高的条件下,获得了上述有限元的渐近展式与外推结果,所用方法主要是构造性的,包括有限元空间的正交分解方法,单元合并技术,能量嵌入方法,分离变量法以及F
泛函微分方程广泛出现于生物学、物理学、经济及社会学、控制论及工程技术等诸多领域。其算法理论的研究对推动这些科技领域的发展无疑非常重要。近年来,泛函微分方程,特别是其特例——延迟微分方程的算法理论的研究得到了很大发展,获得了丰硕成果,这些成果可参见Barwell,Bellen,Torelli,Zennaro,Spijker,Watanabe,Roth,in’t Hout,Baker,Paul,Kot
泛函微分方程(FDEs)在自动控制、生物学、医学、化学、人口学、经济学等众多领域有着广泛应用,其理论和算法研究具有无可置疑的重要性,近三十年来,Volterra泛函微分方程(VFDEs),特别是其重要子类——延迟微分方程(DDEs)的算法理论研究得到了众多学者的高度关注,取得了大量研究成果.例如在DDEs数值方法线性稳定性研究领域,Barwell、Watanabe、Zennaro、Spijker、