标准k-ω模型（Wilcox）
标准k-ω模型（Wilcox）的理论基础
概要
老师！标准k-ω模型与SST k-ω有什么区别？
标准k-ω模型由Wilcox (1988, 2006)开发，是SST的基础模型。使用乱流比消散率（Specific Dissipation Rate）$\omega = \varepsilon/(C_\mu k)$ 作为输运变量，壁面附近的处理比k-ε模型更自然。
也就是说对壁面处理更强？
是的。ω的壁面边界条件可以解析定义（Dirichlet条件），不像ε在壁面处有奇异性。可以不用壁函数直接解析粘性底层。但是对自由流边界条件的敏感性很高，这是个重大弱点。
支配方程
请告诉我方程。
Wilcox (2006)改进版k-ω方程：
k方程：
ω方程：
涡粘性：$\mu_t = \rho k / \omega$
模型常数：$\alpha = 13/25$、$\beta = \beta_0 f_\beta$、$\beta^* = 9/100$、$\sigma = 1/2$、$\sigma^* = 3/5$、$\sigma_d = 1/8$（仅当$\nabla k \cdot \nabla \omega > 0$）。
$\sigma_d$的交叉扩散项是有条件的，这很有趣。
这是Wilcox 2006版的重要改进。这大大改善了自由流敏感性问题。1988年版没有这一项，自由流的$\omega$值会使结果产生巨大变化。
自由流敏感性问题
自由流敏感性具体是什么问题？
自由流（远离壁面的区域）中$\omega$的边界值改变，会导致壁面附近的解也改变。这在物理上是不合理的。Menter开发k-ε混合的SST的主要动机就是为了解决这个问题。
2006年版的交叉扩散项大大改善了这个问题，但没有完全解决。因此在工业CFD中，SST k-ω远比标准k-ω广泛使用。
Wilcox花50年持续更新模型
标准k-ω模型自1988年David C. Wilcox提出以来，经历了1998年版、2006年版、2008年版的持续改进。特别是自由流敏感性（入口边界ω值过于敏感的问题）长期以来是众所周知的弱点，Wilcox本人认识到"需要改进"并持续修订。在工程界，一位研究者公开承认自己模型的缺点并持续改良，这样的诚实做法相当罕见。使用时必须确认实现的是哪个版本，因为"标准k-ω"这个名字相同但方程在不同时期会有差异。
标准k-ω模型（Wilcox）的数值计算手法
数值实现
标准k-ω的实现中有什么特殊注意事项吗？
ω的壁面边界条件是最关键的。
壁面边界条件
ω的壁面值怎样确定？
在$y^+ \to 0$的极限下：
这是Dirichlet条件，数值实现上有几种方法：
1. 作为壁面第一层单元中心的值设置：$\omega_P = \frac{6\nu}{\beta_1 (\Delta y_1)^2}$
2. 作为壁面Ghost单元的值设置：Wilcox推荐
3. 壁函数切换：$y^+ > 2.5$时$\omega = u_\tau / (\sqrt{\beta^*} \kappa y)$
方法1的话网格越细ω越大，对吧。
对。$y^+ = 1$时$\omega \sim O(10^6)$，$y^+ = 0.1$时$\omega \sim O(10^8)$。这需要注意矩阵的对角占优性。
离散化方案
k-ω推荐的离散化方案是什么？
| 变量 | 推荐方案 | 备注 |
|---|---|---|
| k | 二阶迎风格式 | 带TVD限制函数 |
| ω | 一阶/二阶迎风格式 | ω在壁面附近急剧变化，一阶精度也可接受 |
| 动量 | 二阶迎风及以上 | |
| 交叉扩散项 | 中心差分 | 梯度内积，注意条件分支 |
OpenFOAM中的设置
OpenFOAM的使用方法是什么？
OpenFOAM已实现Wilcox 2006版的kOmega。
```
RAS
{
RASModel kOmega;
turbulence on;
printCoeffs on;
}
```
但实务中kOmegaSST的使用压倒性更多。kOmega主要用于特定基准验证或需要排除SST混合效应的场景。
收敛性改善
收敛困难时有什么技巧？
k-ω与y+=1的相性——壁面解析的意外收益
k-ω模型在壁面附近ω→∞，而壁上ω有解析解存在，因此与y+=1左右的细网格高度相容。得益于这个特性，该模型在直接解析粘性底层、处理转捩流和低Reynolds数区域时表现很强。航空机翼边界层分析长期采用这个模型的历史也来自此。反过来说"选择k-ω就应该力求y+≤1"这一实务指引是从模型的物理特性自然导出的。网格越细k-ω的威力越发挥，这是k-ω用户应该掌握的重要知识。
标准k-ω模型（Wilcox）的实务应用
实务指南
标准k-ω应该在哪些场合使用？有SST的话不是不需要了吗？
大多数情况SST足够了，但标准k-ω在以下情形有应用。
正当使用的情形
| 情形 | 理由 |
|---|---|
| 与Wilcox论文的对比验证 | 需要用同一模型再现结果 |
| 壁面主导的边界层流 | 自由流影响较小时，与SST精度相当 |
| 排除SST混合影响 | 研究用途需要分离混合效应 |
| 历史验证资产 | 过往验证数据基于k-ω时 |
SST vs 标准k-ω的比较
具体怎样不同？
| 特性 | 标准k-ω | SST k-ω |
|---|---|---|
| 壁面附近精度 | 高 | 高（相当） |
| 自由流敏感性 | 中等（2006版改善） | 低（k-ε混合） |
| 逆压力梯度剥离 | 倾向低估 | 改善（SST限制） |
| 滞止点问题 | 存在 | 用生产限制器缓解 |
| DES/DDES化 | 可行但稀有 | 标准方法 |
入口边界条件
入口ω值怎样设置？
一般转换公式：
或从乱流粘性比：
| 参数 | 外部流 | 内部流 |
|---|---|---|
| 乱流强度$I$ | 0.1-1% | 5-10% |
| 长度尺度$l_t$ | $0.01 L_{ref}$ | $0.07 D_h$ |
| 粘性比$\mu_t/\mu$ | 1-10 | 10-100 |
考虑自由流敏感性，应对ω入口值进行敏感性分析验证。尤其1988版结果变化很大。
看来实务中用SST更安全。
对。标准k-ω的存在意义主要在于深化对SST的理解。
k-ω在宇宙飞行器热防护中的应用
返回舱和超音速导弹等航天器的气动加热分析中，k-ω模型是标准选择。这是因为k-ω在冲击波后的高温边界层和壁面热流通量预测中表现相对较好。JAXA、NASA的技术报告中也能看到用k-ω进行超音速、极超音速壁面热流评估的事例。"普通配管分析用的模型竟也在航天器应用"，这展示了k-ω守备范围的广度。
标准k-ω模型（Wilcox）的软件对比
商业工具的支持
各求解器对标准k-ω的支持如何？
| 求解器 | 实现版本 | 备注 |
|---|---|---|
| Ansys Fluent | Wilcox 1998版相当 | 有低Reynolds修正、剪切流修正选项 |
| Ansys CFX | Wilcox版 | 强烈推荐SST。k-ω主要为向后兼容 |
| STAR-CCM+ | Wilcox 2006版 | |
| OpenFOAM | Wilcox 2006版 | kOmega类 |
Fluent中的注意事项
Fluent的标准k-ω有什么选项？
Fluent的k-omega (Standard)有：
- Low-Reynolds Corrections：转捩区的减衰函数
- Shear Flow Corrections：剪切流的扩散项修正
- Transitional Flows：自由剪切层修正
这些默认关闭，需要时手动启用。但Fluent手册也推荐用SST。
实现版本的差异
1988版和2006版差异很大吗？
| 特性 | 1988版 | 2006版 |
|---|---|---|
| 交叉扩散项 | 无 | $\sigma_d/\omega (\nabla k \cdot \nabla \omega)$（有条件） |
| 自由流敏感性 | 大 | 改善 |
| β | 常数 | 由$f_\beta$变动（涡伸效应） |
| Realizability | 无 | 加应力限制 |
2006版的交叉扩散项在形式上类似SST的交叉扩散，但不用混合函数，这是区别。
说到底，选标准k-ω的理由相当有限。
对。除了学术对比研究、与过往验证结果的一致性、教学用途，用SST才是实务上的正确选择。
"Wilcox 1988 vs 2006"——求解器版本迷宫
同样是"标准k-ω"，但ANSYS Fluent、STAR-CCM+、OpenFOAM实现的Wilcox年代版本不同。如OpenFOAM的kOmega是基于Wilcox 2006版，而Fluent旧版接近1988年版本。项目间切换求解器时"设置相同却结果不同"容易发生，接手他人的分析工作时必须查看求解器发行说明来了解模型方程差异。
标准k-ω模型（Wilcox）的前沿研究
前沿话题
标准k-ω的研究发展有吗？
Wilcox本人在2006年后仍继续改进。
Wilcox 2006版的定位
2006版改善了多少？
Wilcox在著作"Turbulence Modeling for CFD" (2006, 第3版)进行了大幅改进。主要改善：
1. 交叉扩散项：大幅降低自由流敏感性
2. 应力限制：类似$\mu_t = \min\left(\frac{\rho k}{\omega}, \frac{\rho a_1 k}{\Omega F_2}\right)$的限制
3. 涡伸补正：β含涡度张量不变量
这些改进使与SST的性能差距缩小，但SST的实绩和普及程度的优势难以追赶。
ω方程的理论意义
为什么用ω而不用ε？
历史上Kolmogorov (1942)最初提出ω（"湍流频率"）。$\omega = \varepsilon/(C_\mu k)$是乱流时间尺度$\tau = 1/\omega$的倒数，在壁面上取有限值（ε在壁面奇异）。
物理意义：
- $k$：乱流能量的"大小"
- $\omega$：乱流能量的"消散速率"（乱流旋转频率）
- $\varepsilon$：乱流能量的"消散量"
ω在壁面的行为比ε更合理，但自由流中ε更稳定。SST就是取两者之长。
可压缩性扩展
k-ω能用于超音速流吗？
Wilcox提出了可压缩性修正。高马赫数的乱流压缩性效应（dilat dissipation）：
$M_t = \sqrt{2k}/a$是乱流马赫数，$a$是音速。$M_t > M_{t0} \approx 0.25$时压缩性效应启动。超音速混合层扩展率的抑制等有效。
了解k-ω的历史和理论基础对理解SST帮助很大。
对。要正确理解SST各项的含义，需要k-ω和k-ε两者的知识。
k-ω的自由流敏感性——烦人的"无穷远处的影响"
标准k-ω最著名的弱点是"自由流敏感性"。仅改变入口边界处ω（比消散率）值，就会导致远处边界层解产生变化。这物理上不合理，相当于"无穷远静止流体的乱流特性影响壁面附近"，违反物理直觉。Wilcox在1998和2006年改版时试图修正这问题，但未完全解决。SST模型的诞生动机之一就是避免这一问题，所以说k-ω的历史弱点反而产生了SST这部杰作。
标准k-ω模型（Wilcox）的故障排除
故障排除
标准k-ω特有的问题有什么？
与SST共通的问题也多，但标准k-ω有其特有问题。
1. 自由流敏感性（1988版）
改变入口ω值，结果变化很大。
原因：1988版k-ω对自由流ω边界值敏感。ω_inlet改变10倍，壁面摩擦系数会变20%以上。
对策：
- 切换到2006版（交叉扩散项改善）
- 切换到SST（根本解决）
- 必须用1988版时，进行ω敏感性分析评估结果不确定性
2. 无转捩全乱流假定
低Reynolds数流中层流区域没被计算出来。
原因：k-ω假定全乱流。即使用Low-Re修正也无法很好预测自然转捩。
对策：
- 用γ-Reθ转捩模型（Transition SST）
- 手动指定层流区域（乱流粘性固定为零）
- 层流-乱流分别建网
3. 后流域的过度扩散
圆柱后流的Karman涡快速消失。
原因：定常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-ω用户间很常见。
关联话题
价值
更详细
报告