自适应间断有限元方法求解三维可压缩欧拉方程

来源 :中国科学技术大学 | 被引量 : 3次 | 上传用户:juhaixu
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
偏微分方程(组)广泛应用于大气物理、天体物理、燃烧与爆炸理论、航空与航天、惯性约束聚变、石油勘探等众多领域.本文研究的欧拉方程组和辐射扩散方程组具有很强的物理意义,是辐射流体力学研究的主要内容.如何高效快捷的求解该类非线性问题成为广大科研工作者探索的一个重要课题.   自适应间断有限元方法和预处理的Jacobian-free Newton-Krylov(JFNK)方法具有计算量小,易于处理大规模、多尺度的高度非线性问题的特点.因此,这两类方法被广泛应用于流体力学和传热学领域.本文通过设计合理的算法,保证自适应间断有限元方法和预处理的JFNK方法具有高效求解该类非线性问题的特点,利用数值算例验证了我们方法的有效性.   第一章首先回顾了各种间断有限方法,其中主要包括:Runge-Kutta间断有限元方法,Lax-Wendroff时间离散的间断有限元方法,局部间断有限元方法.其次,介绍了三类自适应方法的发展和基本思想.最后,介绍了非精确Newton方法和Krylov子空间方法,并且详细阐述了GMRES方法、BiCGSTAB方法和TFQMR方法,我们分析了这三类方法之间的联系,并比较了其优缺点.   第二章应用自适应LWDG方法求解三维双曲守恒律方程组,文中给出一种自适应策略,其中均衡折中策略适用于非相容四面体网格.将二维情形下的后验误差指示子推广到三维双曲守恒律方程组中,与传统的二阶RKDG方法相比,该方法具有计算量小和精度高的特点.数值实验证明了方法的有效性.   第三章利用自适应间断有限元方法在非相容四而体网格上求解三维可压缩欧拉方程.通过选取不同物理量的梯度作为网格加密或粗化的条件,比较了在三维非相容网格情形下用密度,能量,压力和焓作为误差指示子的效果.数值实验表明选用不同物理量作误差指示子的效率.   第四章研究了相容四面体网格下重构的自适应间断有限元方法.通过对DG解在单元中心处进行泰勒展开,并且采用间断有限元方法计算本单元以及相邻单元的平均值和一阶导数值,有效地减少了求解三维可压缩Euler方程的计算量.与二阶的自适应RKDG方法相比,重构的自适应间断有限元方法具有计算量小的特点.   第五章采用局部时间步长的自适应间断有限元方法求解三维可压缩欧拉方程.我们将改进的Osher NC流通量应用到间断有限元方法数值流量的计算中,与统一步长的自适应间断有限元方法相比,该方法具有计算量小的特点.   最后一章针对大规模非线性标量方程组的离散求解问题,研究JFNK方法.该方法最重要的优点在于不需要形成和存储非线性方程的Jacobi矩阵.为了应用预处理技术求解一类非平衡辐射扩散与物质温度热传导耦合方程组,我们设计了两类用于构造预处理的线性化方法.数值结果表明这两种预处理方法能够改进JFNK方法的收敛性.
其他文献
随着科技的不断发展,人脸识别已经成为一种重要的生物识别手段。在人脸识别领域,由于人脸图像的维数相当高,直接在原图上进行处理,将加大算法的复杂度,并且对计算机的硬件性能也是
随着科学技术与经济的迅速腾飞,我国化工行业也有着良好的发展势头。而在化工生产的产业链中化工分析是关键的组成部分,一旦缺失化工分析的过程,必然会使化工生产过程的复杂性大
本文通过对荣华二采区10
学位
(e)p正则化问题在变量选择、信号处理、压缩传感、数据挖掘、金融最优化等许多领域有广泛的应用背景.对该问题的理论与算法的研究是目前国际优化界关注的一个热点.本文研究求
轨道方法对表示论中许多困难(遗留)问题的突破有重要作用。在余伴随轨道研究中,有一个基本的结果:对一般线性李代数gl(G),它的子代数的对偶空间可以实现为gl(G)的某个子空间。本文
生活,就是在不断地继续过去,创造未来,我们继承着传统的精髓,却也要发扬个性的独特,于是我们听到这样一句话“Do It Yourself”,是时候—— DIY,英文Do It Yourself的缩写,自
请下载后查看,本文暂不支持在线获取查看简介。 Please download to view, this article does not support online access to view profile.
期刊
日本京都先进通信研究实验室 (TAR)的研究人员提出一种新的“混沌反射镜”,可能会形成新的自由空间光通信方法的基础。这种 7面形的腔式反射镜正由该室的自适应通信研究实验
A series of unconfined compression tests (UCTs) were conducted to investigate the effects of content of reactive magnesia (MgO) and carbonation time on the engi