SST k-ω模型（Menter）
理论与物理
概述
老师！ SST k-ω模型在CFD中给我的印象是使用最广泛的，为什么它这么受欢迎呢？
Menter (1994) 开发的SST（剪切应力输运）k-ω模型，是一种混合模型，它结合了k-ω模型在近壁区域的优越性，以及k-ε模型在自由流中的稳定性。此外，通过加入考虑湍流剪切应力输运的限制器，其在逆压梯度下的分离预测能力得到了显著改善。
这就像一个一石三鸟的模型呢。
正是如此。可以说它是工业CFD中广泛使用的“默认模型”。
控制方程
请告诉我数学公式。
k方程:
ω方程:
最后一项交叉扩散项（Cross-Diffusion Term）是从k-ε转换时产生的项，由 $F_1$ 混合函数控制。
混合函数:
在壁面附近 $F_1 \to 1$（k-ω行为），在远处 $F_1 \to 0$（k-ε行为）。
涡粘性定义（SST限制器）:
$a_1 = 0.31$，$F_2$ 是第二混合函数。这个限制器近似了Bradshaw假设（$-\overline{u'v'} = a_1 k$），改善了逆压梯度下的分离预测。
$F_1$ 是根据到壁面的距离来切换的呢。
准确地说，它不仅根据壁面距离 $y$，还利用流场的局部量（$k$、$\omega$、$\nu$）进行自适应切换。壁面附近的k-ω区域和远处的k-ε区域被平滑地连接起来。
两个模型的"优点结合"——Menter的构想
Florian Menter于1994年发表的SST模型论文，提出了一种在当时颇具创新性的混合方法：通过混合函数将k-ε的自由流稳定性和k-ω的近壁精度结合起来。这种在翼型周围边界层表现为k-ω、在自由流中接近k-ε行为的切换机制，使其成为航空、汽车、涡轮机械CFD的标准模型。如今它已普及到成为全球CFD求解器的默认设置。
各项的物理意义
- 时间项 $\partial(\rho\phi)/\partial t$：想象一下拧开水龙头的瞬间。最初水流会不稳定地喷溅，过一会儿才变成稳定的水流，对吧？描述这个“变化过程中”的就是时间项。心脏搏动导致血流脉动，发动机阀门每次开闭引起流动变化，这些都是非定常现象。那么定常分析是什么？就是只观察“经过足够时间流动稳定后”的状态——也就是将此项设为零。由于计算成本大幅降低，先用定常分析求解是CFD的基本策略。
- 对流项 $\nabla \cdot (\rho \mathbf{u} \phi)$：把落叶扔进河里会怎样？它会随水流被带到下游，对吧？这就是“对流”——流体运动搬运物质的效果。暖气的热风能到达房间角落，也是因为空气这个“搬运工”通过对流输送热量。这里有趣的是——这项包含“速度×速度”，因此是非线性的。也就是说，流速加快时此项会急剧增强，变得难以控制。这就是湍流的根本原因。常见的误解：“对流和传导差不多”→ 完全不一样！对流是流动搬运，传导是分子传递。效率有天壤之别。
- 扩散项 $\nabla \cdot (\Gamma \nabla \phi)$：有过在咖啡里倒入牛奶后放置不管的经历吗？即使不搅拌，过一会儿也会自然混合。那就是分子扩散。那么下一个问题——蜂蜜和水，哪个更容易流动？当然是水，对吧？因为蜂蜜的粘性（$\mu$）高，所以不易流动。粘性越大扩散项越强，流体的运动就变得“粘稠”。雷诺数小的流动（缓慢、粘稠）中扩散占主导。相反，Re数大的流动中对流占压倒性优势，扩散则成为配角。
- 压力项 $-\nabla p$：推动注射器的活塞，液体就会从针头有力地射出，对吧？为什么呢？因为活塞侧压力高，针头侧压力低——这个压力差成为推动流体的力。水坝放水也是同样原理。天气图中等压线密集的地方会怎样？没错，会刮强风。“有压力差的地方就会产生流动”——这就是纳维-斯托克斯方程压力项的物理意义。这里的误解点：CFD中的“压力”多为表压而非绝对压力。切换到可压缩分析时结果突然出错，原因可能就是混淆了绝对压力/表压。
- 源项 $S_\phi$：受热的空气会上升——为什么呢？因为比周围空气轻（密度低），被浮力推上去了。这个浮力作为源项添加到方程中。此外，燃气灶火焰产生化学反应热、工厂电磁泵对金属熔液施加洛伦兹力……这些都是“从外部向流体注入能量或力”的作用，都用源项表示。忘记源项会怎样？自然对流分析中如果忘记加入浮力，流体就完全不动——就像冬天房间里开了暖气，暖空气却不上浮这种物理上不可能的结果。
假设条件与适用范围
- 连续介质假设：克努森数 Kn < 0.01（分子平均自由程 ≪ 特征长度）时成立
- 牛顿流体假设：剪切应力与应变速率呈线性关系（非牛顿流体需要粘度模型）
- 不可压缩假设（Ma < 0.3时）：将密度视为常数处理。马赫数0.3以上需考虑可压缩性效应
- Boussinesq近似（自然对流）：仅在浮力项中考虑密度变化，其他项使用恒定密度
- 不适用的情形：稀薄气体（Kn > 0.1）、超音速/高超音速流动（需要激波捕捉）、自由表面流动（需要VOF/Level Set等）
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 速度 $u$ | m/s | 入口条件从体积流量换算时，注意截面积单位 |
| 压力 $p$ | Pa | 区分表压与绝对压力。可压缩分析使用绝对压力 |
| 密度 $\rho$ | kg/m³ | 空气: 约1.225 kg/m³@20°C、水: 约998 kg/m³@20°C |
| 粘性系数 $\mu$ | Pa·s | 注意与运动粘性系数 $\nu = \mu/\rho$ [m²/s] 混淆 |
| 雷诺数 $Re$ | 无量纲 | $Re = \rho u L / \mu$。层流/湍流转换的判断指标 |
| CFL数 | 无量纲 | $CFL = u \Delta t / \Delta x$。直接关系到时间步长的稳定性 |
数值解法与实现
数值实现要点
实现SST时需要特别注意哪些点？
有三个重要点。
1. 壁面距离计算
混合函数需要壁面距离 $y$ 对吧。
壁面距离通常在计算开始前通过求解泊松方程获得，或者几何计算到最近壁面单元的距离。对于复杂形状，基于泊松方程的方法（例如OpenFOAM的wallDist）更稳健。
Fluent会自动计算壁面距离。OpenFOAM中通过 fvOptions 的 wallDist 指定。壁面距离的精度通过混合函数影响结果，因此网格在壁面处的正交性很重要。
2. 生成项限制器
$\tilde{P}_k$ 是加了限制器的吗？
是的。在Menter的原论文中，$\tilde{P}_k = \min(P_k,\; 10 \beta^* \rho k \omega)$，限制生成项不超过耗散项的10倍。这是为了防止驻点处的湍流动能过度积累（驻点异常）。
一些求解器提供基于涡量的生成项选项，可以进一步缓解驻点问题。
3. 交叉扩散项处理
$\omega$ 方程中的交叉扩散项在数值上没问题吗？
$\frac{\partial k}{\partial x_j}\frac{\partial \omega}{\partial x_j}$ 可能为正也可能为负。通过 $CD_{k\omega} = \max\left(2\rho\sigma_{\omega 2}\frac{1}{\omega}\frac{\partial k}{\partial x_j}\frac{\partial \omega}{\partial x_j},\; 10^{-10}\right)$ 来防止除零错误。
OpenFOAM中的设置
请告诉我OpenFOAM中的典型设置。
constant/turbulenceProperties:
```
RAS
{
RASModel kOmegaSST;
turbulence on;
printCoeffs on;
}
```
壁面边界条件:
k:kqRWallFunction（若 $y^+ < 1$ 则用fixedValue 1e-10）omega:omegaWallFunctionnut:nutUSpaldingWallFunction（全$y^+$对应）
Fluent中的推荐设置
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 模型 | k-omega SST | 默认开启生成项限制器 |
| 壁面处理 | 若解析则 $y^+ \approx 1$ | 无需增强壁面处理（SST本身支持低Re） |
| URF (k, omega) | 0.7-0.8 | 收敛困难时可降至0.5 |
| 离散化 | Second Order Upwind | k和ω两者 |
混合函数F1的真面目——“此处k-ε，此处k-ω”的边界
SST模型核心的混合函数F1，是根据壁面距离、湍流动能、涡粘性等计算出的无量纲量，在0到1之间平滑变化。F1=1的区域由k-ω主导，F1=0的区域则运行接近k-ε的公式。实现上需要计算壁面距离，OpenFOAM中使用 wallDist 库，但在复杂形状下若此距离计算不准确，可能导致混合出错。当感觉“SST结果不对劲”时，建议先从可视化 wallDist 开始检查。
迎风格式（Upwind）
一阶迎风：数值扩散大但稳定。二阶迎风：精度提高但有振荡风险。高雷诺数流动中必须使用。
中心差分（Central Differencing）
二阶精度，但Pe数 > 2时会产生数值振荡。适用于低雷诺数的扩散主导流动。
TVD格式（MUSCL、QUICK等）
通过限制器函数抑制数值振荡同时保持高精度。对捕捉激波或陡峭梯度有效。
有限体积法 vs 有限元法
FVM：自然满足守恒定律。CFD的主流。FEM：对复杂形状、多物理场有利。SPH等无网格法也在发展中。
CFL条件（库朗数）
显式法：CFL ≤ 1为稳定条件。隐式法：即使CFL > 1也稳定，但影响精度和迭代次数。LES：推荐CFL ≈ 1。物理意义：一个时间步内信息传播不超过一个网格。
残差监控
连续性方程、动量、能量的各项残差下降3~4个数量级可判断为收敛。质量守恒残差尤其重要。
松弛因子
压力：0.2~0.3，速度：0.5~0.7为一般初始值。发散时降低松弛因子。收敛后可提高以加速。
非定常计算的内部迭代
在每个时间步内迭代直至达到定常解。内部迭代次数：5~20次为参考值。若残差在时间步间波动，需重新审视时间步长。
SIMPLE法的比喻
SIMPLE法是一种“交替调整”的方法。先假设求解速度（预测步），然后根据该速度修正压力以满足质量守恒（修正步），再用修正后的压力修正速度——反复进行这种“接传球”以逼近正确答案。类似于两人调整架子水平的作业：一人调整高度，另一人调整平衡，如此反复交替。
迎风格式的比喻
迎风格式是“站在河流中重视上游信息”的方法。站在河里的人看下游无法知道水的来源——这反映了“上游信息决定下游”的物理规律。虽然是一阶精度，但由于能正确捕捉流动方向，稳定性高。
实践指南
实践指南
请告诉我SST模型在实际工作中使用时的最佳实践。
SST通用性很强，但有几个注意事项。
网格要求
壁面的 $y^+$ 应该怎么处理？
对于壁面解析，建议 $y^+ \approx 1$。SST模型本身是低雷诺数模型，因此可以解析粘性底层。如果使用壁面函数，则 $30 < y^+ < 300$ 是典型范围。但请注意，SST的壁面函数是基于Spalding公式的，与标准壁面函数不同。
对于分离流或强逆压梯度流，即使 $y^+$ 稍大（例如~5），SST通常也能给出合理结果，但为了获得最佳精度，$y^+ \approx 1$ 是理想选择。
入口条件
入口的湍流参数怎么设置？
需要指定湍流强度 $I$ 和湍流粘度比 $\mu_t/\mu$，或者直接指定 $k$ 和 $\omega$。
对于内部流动（管道、风洞等）：
- $I = 0.05$（5%）是常见值
- $\mu_t/\mu = 10$ 是合理的初始估计
对于外部流动（汽车、飞机周围）：
- $I = 0.01$（1%）更典型
- $\mu_t/\mu = 1$
如果已知水力直径 $D_h$ 和湍流长度尺度 $l$，可以使用公式：
收敛性
SST的收敛有时很困难，有什么技巧吗？
是的。SST的 $\omega$ 方程在初始阶段可能不稳定。建议：
- 分阶段初始化：先用k-ε或S-A模型获得初始场，再切换到SST。
- 降低松弛因子：将k和ω的URF降至0.5-0.6，稳定后再提高。
- 监视 $\omega$ 的值：确保其不会变得过大（例如 > 1e6），这可能导致数值不稳定。
- 使用伪瞬态（Pseudo-Transient）：对于稳态问题，使用伪瞬态方法可以提高鲁棒性。
后处理
如何验证SST结果的可信度？
检查以下几点：
- $y^+$ 分布：确认壁面附近的 $y^+$ 值符合预期。
- 湍流粘度比 $\mu_t/\mu$：在边界层外应小于1000，通常为10-100。
- 湍流动能 $k$：不应在驻点等处出现非物理的峰值。
- 残差历史：所有方程的残差都应平稳下降至少3个数量级。
- 全局平衡：检查质量、动量、能量的输入输出是否平衡。
SST的“交通信号灯”——何时使用，何时避免
SST模型就像CFD世界的“万能工具”，但并非所有情况都适用。对于强旋转流（如旋风分离器），SST可能过度预测湍流粘度，此时RSM（雷诺应力模型）更合适。对于高度分离的流动（如后台阶），SST表现良好，但DES（分离涡模拟）能提供更精确的瞬态特征。对于自然对流，SST的浮力效应实现需要仔细检查。记住：没有“最好”的湍流模型，只有“最适合”当前问题的模型。
常见错误与排查
| 现象 | 可能原因 | 解决方案 |
|---|---|---|
| 发散（残差爆炸） | 初始条件太差，松弛因子过大 | 降低URF，使用更好的初始场，启用伪瞬态 |
| $\omega$ 值过大（>1e6） | 壁面附近网格过细或边界条件错误 | 检查壁面 $\omega$ 边界条件，确保 $y^+$ 合理 |
| 分离预测不足 | 网格在分离区不够密，$y^+$ 过大 | 细化分离区网格，确保 $y^+ \approx 1$ |
| 驻点处k值过高 | 生成项限制器未生效或设置错误 | 确认求解器中SST的生成项限制器已开启 |
| 结果与实验偏差大 | 入口湍流条件不合理，网格无关性未验证 | 进行网格无关性研究，调整入口湍流参数 |
进阶技巧
- 曲率修正：对于强曲率流动（如涡轮叶片），启用曲率修正可以改善预测。
- 转捩模型：标准SST是全湍流模型。对于需要考虑层流-湍流转捩的问题（如低雷诺数翼型），可以使用SST与转捩模型（如$\gamma$-$\widetilde{Re}_{\theta}$）的耦合。
- DES/LES混合：对于大分离流动，考虑使用SST-DES或SST-LES混合方法以捕捉瞬态特征。
SST模型选择流程图
问题：外部空气动力学（汽车、飞机） → 选择SST（默认设置）。
问题：内部流动（管道、泵） → 选择SST，注意入口湍流条件。
问题：强旋转流（旋风分离器、搅拌槽） → 考虑RSM或SST+曲率修正。
问题：大分离瞬态流（桥梁风振、汽车相关主题
なった
詳しく
報告