一种健壮的低耗散通量分裂格式-应用数学和力学.PDFVIP

　 ， ４０ ２ Ａｐｐｌｉｅｄ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｍｅｃｈａｎｉｃｓ　 　 ２０１９年２月１ 日出版 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Ｖｏｌ．４０，Ｎｏ．２，Ｆｅｂ．１，２０１９ 文章编号：１０００⁃０８８７（２０１９）０２⁃０１５０⁃１７ ⓒ 应用数学和力学编委会，ＩＳＳＮ １０００⁃０８８７ 一种健壮的低耗散通量分裂格式∗ １ ２ ３ 胡立军 ，　 袁　 礼 ，　 翟　 健 （１． 衡阳师范学院数学与统计学院，湖南 衡阳４２１００２； ２． 中国科学院数学与系统科学研究院，北京 １００１９０； ３． 曙光信息产业（北京）有限公司，北京 １００１９３） 摘要：　 随着计算流体力学的快速发展，设计精确、高效并且健壮的数值格式变得尤为重要．通过对 ３种流行的通量分裂方法（ＡＵＳＭ、Ｚｈａ⁃Ｂｉｌｇｅｎ和Ｔｏｒｏ⁃Ｖáｚｑｕｅｚ）的对流通量和压力通量进行特征分 析，构造了一种简单、低耗散并且健壮的通量分裂格式（命名为Ｒ⁃ＺＢ格式）．采用Ｚｈａ⁃Ｂｉｌｇｅｎ分裂方 法将Ｅｕｌｅｒ方程的通量分裂成对流通量和压力通量，其中对流通量采用迎风方法来计算，压力通量 采用低耗散的ＨＬＬ格式来计算，从而克服了原始的ＨＬＬ格式不能精确分辨接触间断的缺点．数值 实验表明，该文给出的Ｒ⁃ＺＢ格式不仅保留了原始Ｚｈａ⁃Ｂｉｌｇｅｎ格式简单高效、能够精确分辨接触间 断等优点，而且具有更好的健壮性，在计算二维问题时不会出现数值激波不稳定现象． 关　 键　 词：　 Ｅｕｌｅｒ方程；　 Ｚｈａ⁃Ｂｉｌｇｅｎ分裂；　 ＨＬＬ；　 Ｒ⁃ＺＢ；　 低耗散；　 数值激波不稳定性 中图分类号：　 Ｏ３５４；Ｏ２４１．８２　 　 　 文献标志码：　 Ａ ＤＯＩ：１０．２１６５６／ １０００⁃０８８７．３９０１３２ 引　 　 言 近几十年来，尽管流体力学数值方法的研究取得了巨大的进步，设计出了许多高精度、高 ［１］ 效率的数值格式，但是在流体力学的数值模拟中仍然还有许多问题需要解决 ．中心格式、通 量差分裂（ＦＤＳ）格式、通量向量分裂（ＦＶＳ）格式和ＡＵＳＭ（ａｄｖｅｃｔｉｏｎ ｕｐｓｔｒｅａｍ ｓｐｌｉｔｔｉｎｇ ｍｅｔｈｏｄ） 类格式广泛地应用于计算流体力学中．其中，中心格式不能精确地分辨中间特征场并且强烈地 ［２］ 依赖与问题相关的参数 ，这极大地限制了它的应用．基于Ｇｏｄｕｎｏｖ方法的通量差分裂格式需 要求解两个相邻状态之间的局部Ｒｉｅｍａｎｎ（黎曼）问题．尽管Ｇｏｄｕｎｏｖ 和Ｖａｎ Ｌｅｅｒ 系统地阐述 了Ｅｕｌｅｒ方程精确Ｒｉｅｍａｎｎ解法器的构造过程，但是大多数情况下计算精确的Ｒｉｅｍａｎｎ解会很 困难，并且精确Ｒｉｅｍａｎｎ解的数值黏性过小，在计算某些算例时，会产生数值震荡．因此，设计 精确、健壮的近似Ｒｉｅｍａｎｎ解法器成为了流体力学数值方法研究的热点问题．著名的近似Ｒｉｅ⁃ ［３］ ［４］ ［５］ ｍａｎｎ解法器包括Ｒｏｅ 、ＨＬＬ （Ｈａｒｔｅｎ，Ｌａｘ ａｎｄ Ｖａｎ Ｌｅｅｒ）、ＨＬＬＣ （Ｈａｒｔｅｎ，Ｌａｘ ａｎｄ Ｖａｎ Ｌｅｅｒ ［６］ ｗｉｔｈ ｃｏｎｔａｃｔ）和Ｏｓｈｅｒ 等．其中，两波近似的ＨＬＬ格式非常高效和健壮，但是不能捕捉接触间 断极大地限制了它的应用．Ｒｏｅ、ＨＬＬＣ和Ｏｓｈｅｒ格式能够精确地分辨激波和接触间断，从而广 泛应用于流体力学的数值模拟中．但使用它们计算某些多维问题时，会导致非物理现象的出