一种治愈强激波数值不稳定性的混合方法-LSEC.PDF

　 ， ３６ ５ Ａｐｐｌｉｅｄ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｍｅｃｈａｎｉｃｓ　 网络出版时间：2015-07-03 10:04 　 ２０１５年５月１５ 日出版 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Ｖｏｌ．３６，Ｎｏ．５，Ｍａｙ．１５，２０１５ 网络出版地址：/kcms/detail/50.1060.O31004.018.html 文章编号：１０００⁃０８８７（２０１５）０５⁃０４８２⁃１２ ⓒ 应用数学和力学编委会，ＩＳＳＮ １０００⁃０８８７ 一种治愈强激波数值不稳定性的混合方法∗ 胡立军，　 袁　 礼 （中国科学院数学与系统科学研究院计算数学所； 科学与工程计算国家重点实验室（中国科学院），北京 １００１９０） 摘要：　 ＨＬＬＣ（Ｈａｒｔｅｎ⁃Ｌａｘ⁃Ｌｅｅｒ⁃ｃｏｎｔａｃｔ）格式是一种高分辨率格式，能够准确捕捉激波、接触间断和 稀疏波．但是使用ＨＬＬＣ格式计算多维问题时，在强激波附近会出现激波不稳定现象．ＦＯＲＣＥ（ｆｉｒｓｔ⁃ ｏｒｄｅｒ ｃｅｎｔｒｅｄ）格式在强激波附近表现出很好的稳定性，并且其数值耗散比ＨＬＬ（Ｈａｒｔｅｎ⁃Ｌａｘ⁃Ｌｅｅｒ） 格式小．分析了ＨＬＬＣ格式和ＦＯＲＣＥ格式在特定流动条件下的稳定性，构造了ＨＬＬＣ⁃ＦＯＲＣＥ混合 格式并且进一步结合开关函数来消除ＨＬＬＣ格式的激波不稳定现象．数值试验表明新构造的混合 格式不仅能够消除ＨＬＬＣ格式的激波不稳定现象，还最大程度地保留ＨＬＬＣ格式高分辨率的优点． 关　 键　 词：　 数值激波不稳定性；　 稳定性分析；　 ＨＬＬＣ⁃ＦＯＲＣＥ混合格式；　 ＭＵＳＣＬ重构 中图分类号：　 Ｏ３５４；Ｏ２４１．８２　 　 　 文献标志码：　 Ａ ｄｏｉ：１０．３８７９／ ｊ．ｉｓｓｎ．１０００⁃０８８７．２０１５．０５．００４ 引　 　 言 ＨＬＬＣ近似Ｒｉｅｍａｎｎ（黎曼）解法器是非常流行的激波捕捉方法，广泛应用于可压缩高速 流的计算．然而，使用ＨＬＬＣ格式计算多维问题时，不少研究人员都报道了它所存在的缺陷． ［１］ Ｐｅｅｒｙ等 发现，计算二维钝体绕流时，某些格式会产生“红斑”（ｃａｒｂｕｎｃｌｅ）现象，这其中就包 ［２］ 括了ＨＬＬＣ格式．Ｑｕｉｒｋ 指出：用ＨＬＬＣ格式计算某些多维流动算例时会导致非物理现象的产 生，其中包括膨胀激波、负内能、慢行激波、ｃａｒｂｕｎｃｌｅ现象、双Ｍａｃｈ（马赫）杆和奇偶失联现象． 这些现象统称为数值激波不稳定现象，探究这些现象产生的根源并找到相应的治愈方法，对于 流体力学数值计算研究具有深刻的意义． ［２］ 在探究根源方面，Ｑｕｉｒｋ 表明，如果压力的扰动影响到了密度，使用该格式计算会出现不 ［３］ 稳定现象．Ｌｉｏｕ 认为，数值格式的压力项反馈到了质量方程的数值通量，是造成格式不稳定 ［４］ ［３］ 性的根源．Ｘｕ 在Ｌｉｏｕ 的基础上更近了一步，他认为激波面切向耗散的消失也是影响格式 不稳定性的重要原因之一，并指出数值激波结构中的亚声速部分是不稳定现象产生的区域．此 ［５］ 外，Ｄｕｍｂｓｅｒ等 表明，数值激波的结构影响了格式的不稳定性，并且认为不稳定性源于超声 速区域． ∗ 收稿日期：　 ２０１４⁃０７⁃０４；修订日期：　 ２０１５⁃０１⁃１３ 基金项目：　 国家重点基础研究发展计划（９７３计划）（２０１０ＣＢ７３１５０５）；国家自然科学基金（面上项目） ；创新群体项目；国际（地区）合作与交流项目（１１２６１１６０４８６） 作者简介：　 胡立军（１９８５—），男，湖南人，博士生（Ｅ⁃ｍａｉｌ：ｈｕｌｉｊｕｎ＠ｌｓｅｃ．ｃｃ．ａｃ．ｃｎ）； 袁礼（１９６３—），男，重庆人，研究员，博士，博士生导师（通讯作者． Ｅ⁃ｍａｉｌ：ｌｙｕａｎ＠ｌｓｅｃ．ｃｃ．