解双曲型守恒律的修正加权群速度控制格式.pdfVIP

LHD 2011 年度夏季学术研讨会 133 解双曲型守恒律的修正加权群速度控制格式 1 1 2 2 何志伟 ，李新亮 ，傅德薰 ，马延文 （1 中国科学院力学研究所高温气体动力学国家重点实验室 （筹），北京海淀区 100190） (2 中国科学院力学研究所非线性力学重点实验室，北京海淀区 100190) 摘要 基于群速度控制理论，本文提出了一类新的格式——修正加权群速度控制格式。通过常见的 测试案例——包括：激波/密度波干扰问题，双马赫反射问题，瑞利-泰勒不稳定问题，并和经典的 WENO 格式的数值结果相比，新格式的数值结果显示了对于小尺度结构很好的分辨率和良好的激波捕 捉能力。因此，新格式适用于复杂流动的数值模拟。 关键词 高精度；双曲型守恒律；群速度控制；WENO 1. 引言 高精度有限差分格式在复杂流动的数值模 关函数来实现判定。另一类方法是修改 WENO 拟方面，如湍流的 DNS 和 LES，起着很重要的 格式的加权方式从而实现在光滑处具有较小的 角色。对于可压缩湍流，当马赫数足够大时， 数值粘性[6-8,10,11]。 流场中就会出现激波（激波簇），这些激波 为了解释和克服在数值结果中间断处产生 （激波簇）给湍流的DNS 和 LES 带来了很大的 数值震荡的现象，文献[9]提出了群速度控制 挑战。为了克服激波附近的数值震荡，数值格 (GVC)理论。该理论认为间断处产生数值震荡 式要提供足够的数值粘性，但是过大的数值粘 的原因是由于波包的群速度不一致造成的。本 性则会耗散掉小尺度结构，从而破坏数值模拟 文根据群速度控制理论，系统的提出了如何定 的结果。所以，设计具有高分辨率和低耗散的 量地实现慢格式和混合格式，如何通过加权实 高精度激波捕捉格式对于湍流DNS 和 LES 是非 现群速度控制(GVC)。为了克服 GVC 理论的缺 常重要的。 点，借助于 WENO 加权理论，本文提出了一类 WENO 格式[1,2]是目前最为成功的数值格 名为修正加权群速度控制格式(WGVC-M)。大量 式。在复杂流动的数值模拟中有着广泛的运 的数值案例的测试表明，和同阶的经典的 WENO 用。WENO 格式已被相当多的数值结果证明具有 格式[1,2]相比，WGVC-M 显示了非常好的对于 非常好的激波捕捉能力。但是大量的数值结果 小尺度结果的分辨率和良好的激波捕捉能力。 也表明经典的 WENO 格式除了计算量很大的缺 本文结构如下：第 2 部分将系统的解释群 点外，对于小尺度结构的耗散也过大，因而不 速度控制理论。第 3 部分将系统的解释如何构 是可压缩湍流数值模拟最为理想的数值格式。 造修正加权群速度控制格式，包括：如何定量 最近，为了克服 WENO 格式耗散过大的缺点， 地实现慢格式和混和格式；如何通过加权实现 很多研究者做了各种改进。一类方法是混合格 群速度控制，如何构造修正加权群速度控制格 式。即间断处使用 WENO 格式，而流场光滑处 式(WGVC-M)。第 4 部分将给出新的格式对于一 使用一类低耗散格式[3,4]或者是精度最优的 些经典的数值格式测试案例的数值结果。结论 格式[5]。而间断处和光滑处则是通过设计开 放在第5 部分。 2. 群速度控制(GVC)理论 134 LHD 2011 年度夏季学术研讨会 文献[9]提出了一种群速度控制(GVC)理 du (t )