分辨率优化的混合WENO格式.PDF

第３５卷 第４期 计　 　 算　 　 物　 　 理 Ｖｏｌ．３５，Ｎｏ．４ ２０１８年７月 　 ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＣＯＭＰＵＴＡＴＩＯＮＡＬ ＰＨＹＳＩＣＳ　 Ｊｕｌ． ，２０１８ 文章编号：１００１⁃２４６Ｘ（２０１８）０４⁃０３９７⁃０８ 分辨率优化的混合ＷＥＮＯ格式 １ １ １，∗ ２，３ 郭　 元 ，　 田　 奇 ，　 梁　 贤 ，　 李新亮 （１ 北方民族大学数学与信息科学学院，宁夏 银川　 ７５００２１；　 ２ 中国科学院大学工程科学学院，北京　 １０００４９） ３．中国科学院力学研究所高温气体动力学国家重点实验室，北京　 １００１９０） 摘　 要：为提高有限差分格式的分辨率，利用傅里叶分析对ＷＥＮＯ格式进行色散及耗散优化，并给出优化的线性 权重．用优化后的ＷＥＮＯ格式与保单调格式（ＭＰ）进行加权混合，得到新的加权混合ＷＥＮＯ格式（Ｈ⁃ＷＥＮＯ）．通过 一维激波管问题、Ｓｈｕ⁃Ｏｓｈｅｒ 问题及二维双Ｍａｃｈ反射问题及Ｒ⁃Ｔ不稳定性问题对格式进行数值测试．结果显示，新 格式具有强健的激波捕捉能力和对小尺度波结构的高分辨率，与原ＷＥＮＯ格式相比改进明显． 关键词：混合格式；数值耗散；色散；保单调；优化 中图分类号：Ｏ３５４５ 文献标识码：Ａ ＤＯＩ：１０．１９５９６／ ｊ．ｃｎｋｉ．１００１⁃２４６ｘ．７６６４ ０　 引言 ［１］ 随着直接数值模拟（ＤＮＳ）及大涡模拟（ＬＥＳ）等复杂流场精细模拟 的广泛开展，高精度、高分辨率数值 ［２］ 方法越来越受到重视 ．对于可压缩湍流，流场中既包含了多尺度脉动又包含间断． 由于可压缩湍流精细模 拟的需求驱动，高分辨率、低耗散的高精度激波捕捉格式研究成为了计算流体力学研究的热点方向之一． 近年来以加权基本无振荡格式（ＷｅｉｇｈｔｅｄＥｓｓｅｎｔｉａｌＮｏｎ⁃ＯｓｃｉｌｌａｔｉｏｎＳｃｈｅｍｅ，ＷＥＮＯ）为代表的高精度激波 ［３］ 捕捉格式得到了广泛地应用 ．该格式采用子模板加权思想，在光滑区逼近高精度高分辨率的线性格式，而 在间断区利用较低精度的子模板格式尽量抑制数值振荡．Ｊｉａｎｇ 和Ｓｈｕ 分析并改进了Ｌｉｕ等人提出的ＷＥＮＯ ［４］ ［５］ 格式 ，并使用了一种新的光滑度量因子计算方法，提升了计算效率及分辨率 ．该格式具有非常好的鲁棒 性，在可压缩湍流的精细模拟中发挥了重要作用．ＷＥＮＯ⁃ＪＳ格式的不足是数值耗散仍有些偏大，近年来学者 对其进行了各种改进，Ｈｅｎｒｉｃｋ等人指出ＷＥＮＯ⁃ＪＳ格式的光滑度量因子未能在极值点处提高ＷＥＮＯ格式的 ［６］ ［７］ 精度 ，提出一映射函数对格式进行优化，提升了格式极值点处的精度． Ｂｏｒｇｅｓ 等人 设计出一种新的 ＷＥＮＯ光滑度量因子，通过修改 ＷＥＮＯ 格式的权重，构造了ＷＥＮＯ⁃Ｚ 格式． 与经典的ＷＥＮＯ 格式相比， ＷＥＮＯ⁃Ｚ格式能够更快地得到ＷＥＮＯ格式的理想权重．通过对ＷＥＮＯ⁃ＪＳ增加一组额外的迎风基架点并且 ［８］ 选择模板的最优权重，Ｍａｒｔｉｎ等人发展了一种ＷＥＮＯ⁃ＳＹＭＢＯ格式 ，使得ＷＥＮＯ格式的分辨率和宽带耗散 性质得到了改善．除了优化之外，混合格式（Ｈｙｂｒｉｄ Ｓｃｈｅｍｅ）也是减少格式耗散的一种方法 ［９－１０］ ，该方法基 于激波识别器，利用切换或者加权方法，将两种格式