可压缩湍流建模 — 故障排除指南

可压缩湍流建模 — 故障排除指南

故障排除

1. 混合层增长率过大

原因：可压缩性补正无效。标准k-epsilon或k-omega是用不可压缩流校准的，因此无法预测可压缩性效应（增长率抑制）。

对策：

• 启用Sarkar或Zeman可压缩性补正
• 在Fluent中：Viscous Model → Options → Compressibility Effects 勾选
• 确认对流马赫数 $M_c = (U_1 - U_2)/(a_1 + a_2)$。若 $M_c > 0.3$ 则必须进行补正

2. 衝撃波附近湍流量的非物理增加

原因：衝撃波产生的压力梯度被湍流生成项 $P_k$ 过度评估。生产限制器未设置。

对策：

• 启用SST k-omega的生产限制器：$P_k = \min(P_k, C_{lim} \cdot \bar{\rho} \beta^* k \omega)$，$C_{lim} = 10$
• 在Fluent中使用"Production Limiter"或"Kato-Launder modification"
• 在衝撃波附近细化网格，减少衝撃波的数值厚度

3. 基于密度的求解器发散

原因和对策：

• 初始场不适当：整个区域初始化为M=0的静止场会导致与超音速入口的不匹配，衝撃波在初期形成引起发散。→ 用等熵关系设置适当的初始场
• CFL数过大：初期从 CFL = 0.5～1.0 开始，逐步提高
• 边界条件不匹配：在超音速出口使用压力出口可能导致背压与流场矛盾。→ 在超音速出口使用外推边界条件

4. DES计算中出现GID（网格诱导分离）

原因：DES的RANS→LES切换在边界层内发生，应该用RANS求解的区域因LES网格分辨率不足而被求解。

对策：

• 使用DDES（延迟DES）。用 $f_d$ 屏障函数保护边界层为RANS区域
• 切换至IDDES（改进型DDES）。包含WMLES功能，更稳定
• 重新检查网格设计：检查壁面平行方向的网格宽高比是否过大