激波捕捉方法的数值稳定性研究

激波捕捉方法的数值稳定性研究

论文摘要

尽管CFD方法和计算技术都取得了长足的进展,高超声速流动的数值模拟仍然存在巨大的挑战。目前应用于高超声速流动模拟的激波捕捉方法通常会遭遇激波不稳定问题,如著名的“痈”(carbuncle)等激波异常现象。对于简单的流动问题,准确的流场结构一般是可以预知的,激波异常流场的识别是相对容易的。然而,实际的数值计算往往涉及复杂的几何外形和复杂的流动现象,如果产生激波不稳定等异常现象,几乎无法对准确流场与异常流场做区分。这一困难大大降低了现有CFD方法的可靠性,影响了它们在高超声速流动中的应用。因此,研究激波捕捉方法的数值稳定性对于提升现有CFD方法的可靠性,实现高超声速流动的准确预测具有十分重要的意义。综合应用线性扰动分析与数值试验相结合的研究方法,针对典型的Godunov型激波捕捉格式进行了数值耗散性分析,阐明了不同类型激波捕捉方法对线性退化波的数值耗散特性。研究发现,熵波与剪切波的数值黏性对激波稳定性具有不同的影响作用,与熵波相关的数值耗散并不能够有效地抑制数值激波不稳定,而与剪切波相关的数值耗散则能够有效地稳定激波。采用数值试验的研究方法,针对大量的激波捕捉格式进行了激波稳定性的试验研究,研究发现导致激波不稳定的扰动来源于数值激波结构内部,如果激波前后质量流量能够保持一致,那么激波捕捉格式则是稳定的。与线性扰动分析相结合,本文进一步阐明了数值激波不稳定的发生机理。研究发现导致激波不稳定的扰动产生于激波结构内部并随着熵波与声波向激波下游传播,继而导致激波后出现非物理的扰动误差,引发激波前后质量流量的不一致,最终导致激波不稳定的发生。从热力学第二定律出发,综合运用熵产生分析与数值试验相结合的研究方法,对激波捕捉方法的数值熵产生与激波稳定性之间的关系进行系统性的研究,研究揭示了数值激波不稳定的内在机理。研究发现数值激波不稳定问题是由于激波结构内部不适量的熵产生所引起的。如果激波捕捉方法能够保证激波结构内部产生足够的熵,那么数值格式捕捉到的激波将是稳定的。基于研究结论,本文提出了治愈低耗散激波捕捉方法数值激波不稳定问题的熵控制法,该方法可以应用于多种激波捕捉格式,具备一般性。不同于已有的激波不稳定修正方法,熵控制法并不依赖于额外的数值黏性抑制激波不稳定,从而不影响原有格式的求解精度,非常适合于高超声速复杂流动问题的计算。为了满足全速域流动对数值方法的要求,本文针对激波捕捉方法在低马赫不可压缩极限下的数值特性展开了理论分析与数值试验研究。采用全马赫修正法将适用于可压缩流动的激波捕捉方法拓展到低马赫不可压缩流域。结合抑制数值激波不稳定的熵控制法,本文提出了一种构建全速域格式的简单框架,该全速域修正框架具有一般性,适用于多种激波捕捉格式。采用一系列多种复杂的验证算例验证了典型全速域格式的求解精度与鲁棒性,试验结果表明全速域激波捕捉方法能够准确稳定地模拟全速域范围内多种复杂的流动问题。

论文目录

  • 摘要
  • Abstract
  • 第一章 绪论
  •   1.1 研究背景与意义
  •   1.2 国内外研究现状及发展动态分析
  •     1.2.1 激波捕捉方法的研究进展
  •     1.2.2 数值激波稳定性问题的研究进展
  •   1.3 本文工作
  • 第二章 控制方程与数值离散方法
  •   2.1 控制方程与有限体积法
  •   2.2 Godunov方法与迎风差分
  •   2.3 近似黎曼解算器
  •     2.3.1 Roe型格式
  •     2.3.2 HLL型格式
  •     2.3.3 矢通量分裂格式
  •   2.4 本章小结
  • 第三章 数值激波不稳定的机理研究——数值黏性
  •   3.1 Godunov型格式的数值耗散性研究
  •     3.1.1 两种耗散可控的HLL型格式
  •     3.1.2 Godunov型格式的数值耗散性分析
  •   3.2 Godunov型格式的数值激波不稳定特性
  •     3.2.1 数值试验设置
  •     3.2.2 一维数值激波不稳定特性
  •     3.2.3 多维数值激波不稳定特性
  •   3.3 Godunov型格式的数值激波稳定性分析
  •     3.3.1 数值激波不稳定机理的试验研究
  •     3.3.2 数值黏性与激波不稳定的关系
  •     3.3.3 “痈”现象的启发式解释
  •   3.4 抑制数值激波不稳定的剪切波黏性法
  •     3.4.1 一种激波稳定的Godunov型格式
  •   3.5 数值试验结果
  •     3.5.1 Quirk的奇偶网格扰动问题
  •     3.5.2 激波衍射问题
  •     3.5.3 双马赫反射问题
  •     3.5.4 高超声速钝头体绕流问题
  •     3.5.5 层流平板边界层问题
  •   3.6 本章小结
  • 第四章 数值激波不稳定的机理研究——熵产生
  •   4.1 熵产生与激波稳定性的关系
  •   4.2 Godunov型格式的熵产生分析
  •     4.2.1 HLL型格式的熵产生分析
  •     4.2.2 两种熵可控的Godunov型格式
  •   4.3 熵产生与激波稳定性关系的定量分析
  •     4.3.1 一维数值激波不稳定
  •     4.3.2 多维数值激波不稳定
  •   4.4 熵控制法的线性扰动分析
  •     4.4.1 一维数值激波不稳定
  •     4.4.2 多维数值激波不稳定
  •   4.5 抑制数值激波不稳定的熵控制法
  •     4.5.1 熵控制Godunov型格式
  •     4.5.2 矩阵稳定性分析
  •   4.6 数值试验结果
  •     4.6.1 Quirk的奇偶网格扰动问题
  •     4.6.2 激波衍射问题
  •     4.6.3 双马赫反射问题
  •     4.6.4 高超声速钝头体绕流问题
  •     4.6.5 层流平板边界层问题
  •     4.6.6 Elling的物理“痈”问题
  •   4.7 本章小结
  • 第五章 全速域激波捕捉方法研究
  •   5.1 迎风格式的低马赫数精度问题
  •     5.1.1 迎风格式低马赫奇异现象
  •     5.1.2 Euler方程的低马赫不可压缩极限
  •     5.1.3 迎风格式的不可压缩极限
  •   5.2 低马赫精度问题的修正
  •     5.2.1 通用的修正方法
  •     5.2.2 HLL型格式的低马赫修正
  •     5.2.3 渐近分析
  •   5.3 全速域激波捕捉格式构建的一般性框架
  •     5.3.1 全速域Godunov型格式
  •     5.3.2 数值通量的性质
  •   5.4 数值试验结果
  •     5.4.1 改进的Sod激波管问题
  •     5.4.2 双稀疏波问题
  •     5.4.3 Gresho涡
  •     5.4.4 带凸体的通道流动
  •     5.4.5 无黏低马赫NACA0012翼型绕流问题
  •     5.4.6 低马赫NACA0012翼型湍流绕流问题
  •     5.4.7 RAE2822跨声速翼型
  •     5.4.8 双马赫反射问题
  •     5.4.9 高超声速钝头体无粘绕流问题
  •     5.4.10 高超声速钝头体黏性绕流问题
  •   5.5 本章小结
  • 第六章 结论与展望
  •   6.1 论文的主要成果与结论
  •   6.2 论文的主要创新之处
  •   6.3 未来研究工作展望
  • 致谢
  • 参考文献
  • 作者在学期间取得的学术成果
  • 附录A 两种熵修正的Roe格式
  • 附录B 激波不稳定的线性扰动分析
  •   B.1 激波纵向质量流量扰动误差的线性分析
  •   B.2 激波横向质量流量扰动误差的线性分析
  • 附录C HLLEM-ρ格式的熵产生分析
  • 附录D 熵控制法的线性扰动分析
  • 附录E AM-HLLC格式的渐进分析
  •   E.1 AM-HLLC格式的渐进方程
  •   E.2 无量纲分析
  • 文章来源

    类型: 博士论文

    作者: 谢文佳

    导师: 李桦

    关键词: 有限体积法,激波捕捉格式,激波不稳定,痈现象,高超声速,全速域格式

    来源: 国防科技大学

    年度: 2019

    分类: 基础科学

    专业: 力学

    单位: 国防科技大学

    分类号: O354.5

    DOI: 10.27052/d.cnki.gzjgu.2019.000051

    总页数: 178

    文件大小: 12277k

    下载量: 3

    相关论文文献

    • [1].基于纹影法的高压柴油喷雾反射激波特性[J]. 内燃机学报 2020(02)
    • [2].高压柴油射流激波特性影响研究[J]. 现代车用动力 2020(01)
    • [3].激波汇聚效应对球形气泡演化影响的数值研究[J]. 物理学报 2017(06)
    • [4].圆弧激波缓聚点火及后续燃烧传播[J]. 实验流体力学 2019(05)
    • [5].激波的传播与干扰[J]. 力学进展 2016(00)
    • [6].气体体积粘性对二维环形激波聚焦的影响[J]. 计算物理 2017(04)
    • [7].垂直无碰撞激波的离子加速机制[J]. 地球物理学报 2008(04)
    • [8].真空管道超级列车激波簇结构研究[J]. 机械工程学报 2020(02)
    • [9].受限通道中三维激波与侧壁面作用的数值模拟研究[J]. 热科学与技术 2018(04)
    • [10].柱状汇聚激波冲击重气柱不稳定性的数值研究[J]. 中国科学技术大学学报 2014(06)
    • [11].进口斜激波、膨胀波干扰下等直隔离段内的激波串特性[J]. 航空学报 2010(09)
    • [12].基于特征线理论的辨识法在激波装配中的应用[J]. 气体物理 2018(01)
    • [13].柱形汇聚激波冲击球形重气体界面的实验研究[J]. 力学学报 2012(03)
    • [14].管内激波串振荡和壁面脉动压力特性[J]. 实验流体力学 2010(05)
    • [15].暂冲式激波聚焦及起爆爆震的实验与数值模拟研究[J]. 实验流体力学 2016(01)
    • [16].基于气动/弹道耦合的激波针外形优化研究[J]. 弹道学报 2020(02)
    • [17].强激波冲击多晶碳/氦界面诱导的微观Richtmyer-Meshkov不稳定性研究[J]. 中国科学:物理学 力学 天文学 2020(10)
    • [18].具有分流叶片跨声速离心叶轮激波及能量损失分析[J]. 航空动力学报 2016(12)
    • [19].导流块深度对楔形空腔内激波聚焦过程的影响[J]. 航空动力学报 2017(10)
    • [20].激波辐射对红外成像探测影响的抑制方法[J]. 红外与激光工程 2013(02)
    • [21].密切弯曲激波乘波体技术的应用及有效性分析[J]. 推进技术 2018(02)
    • [22].高压燃油喷雾诱导激波产生机理及频率特性[J]. 内燃机学报 2018(05)
    • [23].行星际慢激波的观测与研究[J]. 云南师范大学学报(自然科学版) 2016(03)
    • [24].激波聚焦及起爆爆震波的研究进展[J]. 推进技术 2015(10)
    • [25].同频干涉激波清灰技术在余热锅炉上的应用[J]. 设备管理与维修 2014(05)
    • [26].绕射激波和反射激波冲击下的Richtmyer-Meshkov不稳定性[J]. 中国科学:物理学 力学 天文学 2014(11)
    • [27].地球弓激波的旋转非对称性[J]. 地球物理学报 2010(04)
    • [28].一类非线性奇摄动问题的激波性态[J]. 工程数学学报 2008(02)
    • [29].两类涡轮激波结构及演化规律研究[J]. 北京理工大学学报 2015(06)
    • [30].激波与可运动颗粒群相互作用反射与透射机理实验研究[J]. 应用力学学报 2010(02)

    标签:;  ;  ;  ;  ;  ;  

    激波捕捉方法的数值稳定性研究
    下载Doc文档

    猜你喜欢