[论文解读] Dispersion-dissipation condition for finite difference schemes

Xiao Hu, V. K. Tritschler+2 more|arXiv (Cornell University)|Apr 23, 2012

Computational Fluid Dynamics and Aerodynamics参考文献 10被引用 26

一句话总结

本文提出了一种色散-耗散条件，用于确定有限差分格式中抑制虚假高频波所需的最小数值耗散。通过采用修正波数分析将群速度与耗散率关联，该条件优化了线性和非线性格式中的色散与耗散特性——在提升WENO-CU6-M2格式以实现对湍流流动的高精度解耦模拟方面得到验证，包括自相似能量衰减和湍流转捩过程。

ABSTRACT

A general dispersion-dissipation condition for finite difference schemes is derived by analyzing the numerical dispersion and dissipation of explicit finite-difference schemes. The proper dissipation required to damp spurious high-wavenumber waves in the solution is determined from a physically motivated relation between group velocity and dissipation rate. The application to a previously developed low-dissipation weighted essentially non-oscillatory scheme (WENO-CU6-M2) [X. Y. Hu and N. A. Adams; Scale separation for implicit large eddy simulation, J. Comput. Phys. 230 (2011) 7240-7249] demonstrates that this condition can serve as general guideline for optimizing the dispersion and dissipation of linear and non-linear finite-difference schemes. Moreover, the improved WENO-CU6-M2 which satisfies the dispersion-dissipation condition can be used for under-resolved simulations. We demonstrate this capability by considering transition to turbulence and self-similar energy decay of the three-dimensional Taylor-Green vortex. Simulations of the inviscid and the viscous Taylor-Green vortex at Reynolds numbers ranging from $Re=400$ to $Re=3000$ show a significant improvement over the classical dynamic Smagorinsky model and demonstrate competitiveness with state-of-the-art implicit LES models, while preserving shock-capturing properties.

研究动机与目标

• 建立有限差分格式中色散与耗散之间的定量关联，以控制虚假高频波的传播。
• 解决长期存在的难题：确定能够有效抑制非物理波但不损害已解析尺度精度的最优耗散水平。
• 为线性和非线性格式中的色散与耗散优化提供通用指导，尤其适用于解耦模拟。
• 通过满足色散-耗散条件改进WENO-CU6-M2格式，实现对湍流流动的高精度大涡模拟。

提出的方法

• 基于线性有限差分格式的修正波数，从线性对流方程推导出色散-耗散条件。
• 通过相位误差（|dξ̃R/dξ − 1|）与耗散率（−ξ̃I）的比值定义该条件，确保对临界分辨波的有效衰减。
• 通过修改其重构过程，将该条件应用于WENO-CU6-M2格式，以调整耗散而不改变其高分辨率特性。
• 采用谱分析方法评估不同波数下的色散与耗散特性，尤其关注奈奎斯特截止频率附近区域。
• 在基准问题上验证改进后的格式：Sod激波管、双激波相互作用以及Taylor-Green涡流。
• 将结果与DNS数据及成熟的LES模型（动态Smagorinsky模型、ALDM）进行对比，评估其在能量衰减与耗散率预测方面的精度。

实验结果

研究问题

• RQ1 有限差分格式中，为有效抑制虚假高频波，所需的最小数值耗散量是多少？
• RQ2 如何在保持已解析尺度精度的同时，平衡色散与耗散以抑制非物理振荡？
• RQ3 该色散-耗散条件能否应用于非线性、激波捕捉型格式（如WENO-CU6-M2）？
• RQ4 改进后的WENO-CU6-M2格式在解耦湍流流动模拟中的表现，相较于经典LES模型如何？
• RQ5 满足色散-耗散条件在多大程度上提升了Taylor-Green涡流中自相似能量衰减的预测精度？

主要发现

• 满足色散-耗散条件的改进WENO-CU6-M2格式，能够准确捕捉三维无粘Taylor-Green涡流的自相似能量衰减，其衰减速率与理论预测的t−1.2一致。
• 在雷诺数Re=400至Re=3000的黏性Taylor-Green涡流模拟中，改进格式在耗散率预测方面显著优于经典动态Smagorinsky模型。
• 在Re=400和Re=800条件下，该格式的耗散率与最先进的隐式LES模型——自适应局部重构模型（ALDM）相当。
• Kolmogorov惯性区的k−5/3谱特性在截止波数之前得到良好保持，表明已解析尺度中的数值耗散极小。
• 该格式在保持精确激波捕捉能力的同时，有效减少了虚假波的传播，Sod激波管与双激波相互作用问题的验证结果证实了这一点。
• 色散-耗散条件为线性和非线性格式中的数值耗散提供了系统性优化方法，为高保真湍流流动模拟提供了通用框架。