CORC  > 力学研究所  > 中国科学院力学研究所
题名欧拉方程的高效高精度时间及空间算法研究
作者唐荣
答辩日期2019-05-26
文献子类硕士
授予单位中国科学院大学
授予地点北京
导师崔凯
关键词欧拉方程 共享权重 Weno 格式 多步加权 Crweno 格式 最后一级隐式 Runge-kutta 方法
其他题名The Study of High Efficiency and High Accuracy Numerical Methods for Euler Equations
学位专业流体力学
英文摘要

复杂流动的大规模、长时间数值模拟,除了要求一定的精度和分辨率之外,还要求
其计算过程具有较强的稳定性及较高的计算效率,即需要发展高效高精度的时间及空间
算法。加权基本无振荡(Weighted essentially non-oscillatory,简称 WENO)格式因其具
有优秀的高精度激波捕捉能力而成为当下应用和发展最为广泛的流动空间离散格式;而
多级 Runge-Kutta 方法具有精度高、稳定性好、存储量低等特点,作为经典的常微分方
程数值解法,常被用于计算流体力学领域的时间方向离散推进计算。本文以离散 WENO
空间格式以及 Runge-Kutta 时间推进方法为基础,发展了计算效率更高,而且数值耗散
误差小、稳定性强的相应方法。主要工作包括:
1)发展了一种针对求解流动 Euler 方程的共享权重 WENO 方法,该方法采用
Steger-Warming 分裂后的能量通量来计算 WENO 格式中的子模板权值,所有通量分量
皆利用该权值进行加权组合,因而极大地提高了 Euler 方程的计算效率。对于共享权重
WENO 方法,加权变量的选取起着十分重要的作用,不当的变量选取会导致格式的稳
定性不足或分辨率降低。数值实验结果表明本文所发展的格式具有较好的鲁棒性以及较
低的耗散误差。
2)进一步发展了多步加权 CRWENO 格式。多步 CRWENO 格式在光滑区域具有紧
致格式的谱分辨性质,同时在间断点附近具有较高的精度和分辨率,但与多步 WENO
格式类似,最初的多步 CRWENO 格式同样存在计算量大的问题。本文将最近曾方军等
人提出的关于改进多步 WENO 格式的思想和方法推广到多步 CRWENO 格式,在保证
精度不降低的情况下极大地提高了多步 CRWENO 的计算效率,并通过数值算例对格式
性能进行了比较和验证。
3)提出和发展了最后一级为隐式的高精度 Runge-Kutta 格式。多级 Runge-Kutta 方
法具有精度高、稳定性好、存储量低等特点,但由于稳定性的要求,显式 Runge-Kutta
方法的时间步长受到一定的限制,为了增大时间步长,提高方法的稳定性,我们提出发
展最后一级为隐式的 Runge-Kutta 方法。与显式方法相比,最后一级为隐式的
Runge-Kutta 方法能够在更大的时间步长下推进并保持计算稳定和较低的耗散误差,因
此在对刚性问题或大规模长时间数值模拟时有助于提高计算效率,且由于只有一级需要
进行隐式迭代求解,因此比对角隐式方法具有更快的计算速度。结合稳定性及耗散色散
误差分析确定格式系数,本文给出了最后一级为隐式的三级三阶和四级四阶
Runge-Kutta 方法,并与空间激波捕捉格式相结合,应用于流体力学的数值模拟研究,
验证了方法的有效性。

语种中文
内容类型学位论文
源URL[http://dspace.imech.ac.cn/handle/311007/79085]  
专题中国科学院力学研究所
力学研究所_高温气体动力学国家重点实验室
推荐引用方式
GB/T 7714
唐荣. 欧拉方程的高效高精度时间及空间算法研究[D]. 北京. 中国科学院大学. 2019.
个性服务
查看访问统计
相关权益政策
暂无数据
收藏/分享
所有评论 (0)
暂无评论
 

除非特别说明,本系统中所有内容都受版权保护,并保留所有权利。


©版权所有 ©2017 CSpace - Powered by CSpace