几类典型PDEs高次有限元方程的快速算法研究

来源 :湘潭大学 | 被引量 : 0次 | 上传用户:woyaodefeixiang
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
H(D)(D=grad,curl,div)型椭圆偏微分方程和Maxwell鞍点问题是几类典型的微分方程组(PDEs).高次有限元方法是求解这几类偏微分方程的重要离散化方法,由于这些离散系统系数矩阵的条件数较强地依赖于网格规模、跳系数分布及有限元的次数(如HBk分层基下的高次有限元)等,因此研究其快速求解算法非常必要.本文利用基于辅助空间的预条件子方法、代数多层网格(AMG)法和非重叠区域分解法(DDM),比较系统地研究了上述几类典型PDEs的高次有限元离散系统的高效预条件子及相应的快速求解算法,获得以下主要结果.针对含跳系数的H(grad)型椭圆问题在HBk分层基下的高次有限元方程,给出了一种基于块磨光方法的并行两水平预条件子算法(TLB-p),它本质性地将高次元预条件子构造问题转化为相应的线性元预条件子构造问题.通过将线性元预条件子取为经典并行AMG预条件子,获得了求解高次元方程的第一种预条件子(TLB-AMG-p),数值实验结果表明,相应的并行PCG法的迭代次数基本不依赖于网格规模,弱依赖于高次元的次数以及系数的跳幅.接着,设计了一种基于非重叠DDM的线性元预条件子,与已有的非重叠DDM预条件子相比,它具有粗空间简单和计算复杂度低等优点,理论和数值实验结果表明该预条件子对应的预条件系统的有效条件数是渐近最优的.基于该线性元预条件子,我们得到了另一种求解高次元方程的预条件子(TLB-DDM-p),数值实验结果表明,相应的PCG法也是高效的.针对H(curl)型椭圆问题的两类高次棱有限元离散系统,利用高次元空间的稳定性分解理论和上述构造的H(grad)型高次元预条件子,本质性地将H(curl)型高次元预条件子构造问题转化为相应的线性元预条件子构造问题.进一步,通过将线性元预条件子取为一种基于辅助空间的预条件子,得到了求解H(curl)型高次棱元方程的第一种预条件子,数值实验结果表明,不论是对光滑系数还是对有无浮动子区域的跳系数情形,相应的并行PCG法的迭代次数都基本不依赖于网格规模,弱依赖于系数跳幅,且具有很好的算法可扩展性.接着,设计了一种基于非重叠DDM的线性元预条件子,它具有与上述H(grad)型线性元DDM预条件子同样的优点,特别对常系数的情形,证明了该预条件系统的条件数是渐近最优的.数值实验验证了理论的正确性,同时表明对大跳系数的情形,该预条件子也是高效的.针对Maxwell鞍点问题零阶项系数γ=0的情形,利用正则化思想和辅助空间预条件子方法,为高次棱元鞍点系统设计了一种新的Uzawa算法.特别针对光滑系数下的线性棱元鞍点系统,证明了该Uzawa法的收敛率与网格规模无关.数值实验结果表明,不论是对于有无浮动子区域及有无内交叉点的跳系数情形,该 Uzawa 法的迭代次数基本不依赖于系数跳幅及网格规模,且比常用的 Uzawa 法具有更强的健壮性和更高的运算效率.针对H(div)型椭圆问题的两类高次有限元离散系统,建立了一般高次元空间的稳定性分解理论,它本质性地将H(div)型高次元预条件子构造问题转化为相应的线性元预条件子和H(curl)型高次元预条件子的构造问题.进一步,利用上述构造的H(curl)型高次元预条件子和基于辅助空间的线性元预条件子构造方法,获得了一种H(div)型高次元预条件子,理论上证明了相应的预条件系统的条件数不依赖于网格规模.数值实验验证了理论的正确性,同时也表明了算法的高效性和健壮性.
其他文献
目的:检测NUSAP1在膀胱尿路上皮癌中的表达及临床意义。方法:采用免疫组化方法检测NUSAP1在84例膀胱尿路上皮癌组织及10例正常膀胱上皮组织中的表达情况,并分析其表达与临床
继2017年首次被写入政府工作报告后,人工智能今年又出现在李克强总理所作的政府工作报告当中,在人工智能技术当中,语音识别是一项重要的应用。两会期间,全国人大代表、科大讯飞股
报纸
当下传统的信息检索不能充分理解问题的语义,给出的查询结果往往形似而神不似,且是杂乱的知识需要人为筛选。尤其是进行医疗咨询,医学知识的复杂性和严谨性使得获取该领域相
本文的工作分两部分。第一部分主要研究了一类基于局部多项式逼近空间的单位分解法的最优误差估计。近年来,无网格方法被大量应用到科学与工程计算中。与经典有限元方法相比
本文试图对欧盟《单一数字市场版权指令》第17条关于网络内容分享服务商的平台内容使用新规定进行分析,从法经济学和法社会学的角度论证:该条款是各方利益针对当前网络版权产
在3月18日的俄罗斯总统大选中,俄罗斯现任总统普京以76.67%的得票率强势赢得连任,将于5月7日举行就职仪式,开启第四个总统任期。普京当天在竞选总部对媒体表示,中国是俄罗斯的战略
二氧化钛(TiO2)作为n型宽带隙无机半导体材料,具有无毒性、耐光化学腐蚀、无二次污染且优异的光催化活性,在光催化、气体传感器、空气净化、自清洁材料及太阳能电池等领域得到
<正> 目前,在井矿盐区以硫酸钠型盐岩卤为原料的制盐厂(矿)的盐产量约占井矿盐区盐产量的一半。这些企业在转换经营机制,调整产品结构、以盐为主,多种经营中首先瞄准的目标是
目的:(1)对健康志愿者足踝部进行CT扫描,通过三维重建及有限元前处理,建立正常成年人单足静止站立状态下的有限元模型,验证模型的有效性及真实性,为踝关节损伤的生物力学研究提供可靠模型。(2)对距骨颈骨折后不同内固定方法进行三维有限元分析评价,探究不同内固定方法的力学稳定性及骨折断面最佳的生物力学环境,为临床内固定物的选用提供可靠的力学论据。(3)距骨骨折是足部的严重损伤,而距骨颈骨折又是距骨骨折中
目前物流信息化建设已经成为企业热门话题之一,众多制造型企业开始投资建设属于企业自己的物流管理信息化系统,但是只有少数企业实施物流信息化取得了成功,究其根本原因在于