标准k-ω模型(Wilcox) — 故障排除指南
更丰富的内容可在 k-omega-standard.html 中查看。
故障排除
标准k-ω有特有的问题吗?
与SST有许多共同之处,但标准k-ω有其特有的问题。
1. 自由流敏感性(1988版)
改变入口的ω值后结果发生了很大变化。
原因:1988版的k-ω对自由流的 $\omega$ 边界值敏感。$\omega_{inlet}$ 改变10倍时,壁面Cf可能改变20%以上。
解决方案:
- 切换到2006版(通过交叉扩散项改进)
- 切换到SST(根本解决)
- 如果必须使用1988版,进行 $\omega$ 敏感性分析并评估结果的不确定性
2. 无转换的全湍流假设
在低雷诺数流中,层流区域没有被计算。
原因:k-ω假设全湍流。即使使用Low-Re修正,自然转换的预测也不充分。
解决方案:
- 使用γ-Reθ转换模型(Transition SST)
- 手动指定层流区域(将湍流粘性固定为零)
- 层流-湍流区域分离网格创建
3. 尾流区域的过度扩散
圆柱体尾流中的卡门涡快速消散。
原因:定常RANS的局限性。涡粘性在尾流区域过大,导致涡扩散。k-ω和SST都有同样的问题。
解决方案:
4. ω的初始值设置
我不知道如何设置ω的初始值。
按以下步骤估计。
1. 假设湍流强度 $I$ 和长度尺度 $l_t$
2. $k = \frac{3}{2}(U I)^2$
3. $\omega = \frac{k^{0.5}}{C_\mu^{0.25} l_t}$ 或 $\omega = \frac{\rho k}{\mu \cdot (\mu_t/\mu)}$
粘性比 $\mu_t/\mu$ 在外部流中为1-10,在内部流中为10-100。ω过小($\mu_t$ 过大)会导致初始发散,过大($\mu_t$ 过小)会导致湍流效应较弱。
标准k-ω的问题主要源于自由流敏感性。
正是如此。SST被广泛使用的最大原因是消除了这种自由流敏感性。
壁面上的ω→∞问题——初始值设置的烦恼
在k-ω模型实现中,初学者常遇到"理论上壁面ω为无穷大"的问题。由于数值计算无法处理无穷大,因此需要在壁面第一个单元处根据解析解设置ω的近似值。这个值过大会导致数值不稳定,过小则边界层求解不正确。OpenFOAM的omegaWallFunction实际上是一个非常敏感的实现,不同版本的行为不同。"为壁面ω初始化而苦恼一整天"是k-ω用户之间常听到的故事。
相关主题
更多细节
错误