五阶FD-WENO格式与二阶Godunov格式MUSCL的数值测试与比较

来源 :湘潭大学 | 被引量 : 0次 | 上传用户:wudidewohaha
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
由于求解双曲守恒律组的高阶加权实质上无振荡有限差分格式(简记为FD-WENO)是最近十年才发展起来的,它与上世纪70-80年代发展起来的著名高阶Godunov格式比较,究竟优缺点各如何,哪种格式更好,这个问题目前在国内外文献中研究尚少,尚无定评。特别是对于辐射流体力学问题的数值计算,以上两种格式中究竟哪种格式的实际效果更好,这是惯性约束聚变(ICF)数值模拟研究中迫切需要回答的问题,为此,本文作者发展了用2阶Godunov格式MUSCL求解气体动力学Euler方程的应用软件,并通过同时用五阶FD-WENO格式(下文中简记为WENO5)和二阶MUSCL格式求解各种有代表性的流体力学问题,进行了大量的数值测试。例如,通过求解线性初边值问题、一维Burgers方程初边值问题、若干Riemann问题、一维激波相互碰撞问题及二维周期漩涡问题等对这两种格式进行了测试和定量比较,通过求解Richtmyer-Meshkov(RM)不稳定性问题和Rayleigh-Taylor(RT)不稳定性问题对这两种格式进行了测试和定性比较。我们发现对于简单的Sod Riemann问题或LaX Riemann问题,两种格式都易于算出令人满意的具有较高精度和较高分辨率的数值结果,相对来说,WEN05得到的数值结果粘性稍强,MUSCL得到的结果则在局部地方有微小的非物理振荡,但计算速度更快;对于PeakRiemann问题及、Woodward-Colella激波相互碰撞问题,WENO5所获数值结果在激波及稀疏波附近明显优于MUSCL,前者不仅分辨率较高,误差也较小;对于二维周期漩涡问题,采用81×81的一致网格用WEN05格式计算至时刻t=10的数值结果甚至比采用161×161的一致网格用MUSCL格式计算至时刻t=10的数值结果要好得多,但前者花费的计算时间更短;对于RM不稳定性问题和RT不稳定性问题,WENO5所得到的数值结果比MUSCL所得到的数值结果分辨率更高,图形较为符合实际情况。本文的研究结果表明,尽管对于简单的Riemann问题,二阶MUSCL具有一定优势,但对于既具有强激波又具有大面积复杂流动结构的问题,特别是对于ICF数值模拟问题,二阶格式一般难以提供具有可接受分辨率的数值结果,而高阶FD-WENO格式特别适合于求解这类问题。  
其他文献
德育是培养人品德的活动.要想将中等职业学校的德育落到实处,必须打造一支具有实效性的有战斗力的德育团队.德育教师队伍的建设要与整个教师队伍建设统一规划、部署,着重进行
K是一个有υ个点的有向完全图,G是一个不带孤立点的简单有向图.K的一个G-设计,记作(υ,k,1)-G-GD(其中k表示G的顶点个数,υ表示K的顶点数),是指一个二元组(x,B),其中X为K的点集,B为K的
内部控制是强化风险防范,规范管理行为,强调以预防为主的制度。目前事业单位的内部控制存在诸多弊端,内控体系不健全,本文对事业单位内部控制存在的问题进行简要分析,并就加强事业
期刊
When the fundamental frequency is shifting, it is hard for traditional repetitive controller to work at the resonant frequencies. In this paper, a novel adaptiv
本文研究利用外延基本元判定形式背景属性约简的方法,讨论基于外延和内涵基本元的上、下近似等问题。主要内容有:(1)讨论了形式背景的外延集和内涵集的基本构成单位,并由此提出
任何竞技体育运动项目,都有其自身独特的训练方式和特征.包括摔跤竞技运动不例外,其需要练习运用相应的基本功,而这里指的基本功即最基本的功底,与房屋地基的道理相同.若是摔
2000年左右,Lemaréchal,Oustry和Sagastizábal等人提出uv-分解理论,给出了研究非光滑凸函数的二阶性质的新方法.uv-分解理论的基本思想是将空间分解为两个正交的子空间的直和
期刊
Gorenstein投射模,内射模,平坦模已经被许多学者所研宄,近几十年来,Goren- stein同调代数成为了一个非常活跃的研宄领域,各种形式的Gorenstein同调模被引入和研宄,都是建立在特殊的
本文针对基于连分式的多元混合切触有理插值问题展开了深入的讨论,其内容主要包括基于连分式的一阶混合切触有理插值、基于连分式的高阶混合切触有理插值以及基于连分式的高阶