4.1 湍流建模方法论

4.1 湍流建模方法论

上一章把网格关签完字，本节进入模型关的第一站：在打开湍流模型下拉列表之前，先想明白模型化这件事本身——我们在拿什么换什么。这一节提供的是判断框架，下一节才进具体模型的工具箱。顺序不能反：没有框架的模型选择，只是背口诀。

湍流为什么"算不动"

把纳维-斯托克斯方程直接解到底（直接数值模拟，DNS）在原理上完全可行，坏在成本。湍流能量从大涡逐级传向小涡，最小涡的尺寸与雷诺数的三四分之一次方成反比，三维网格量因此随雷诺数的三次方增长。取工程常见的管道雷诺数十万量级，DNS 所需网格量远超现有任何超算的承受范围。事实上界的残酷在于：这堆小涡对平均流场的贡献其实温和——它们干的活主要是把动量从主流搬向壁面。于是问题转化为：既然算不动全部，能不能只算值得算的部分，把剩下的用模型补上？这就是湍流建模的全部动机。

建模层级：钱与物理的三档交易

RANS（雷诺平均）。把瞬时量分解为平均量加脉动量，脉动对平均量的影响打包成雷诺应力项，再用湍流黏度之类的假设封堵。计算量最小，稳态可解，是工程仿真的绝对主力。代价是丢失全部瞬态细节，且对分离流动、强曲率、旋流的预测依赖模型经验——模型在这些场景经常出错，错得还挺自信。

LES（大涡模拟）。直接解析大尺度涡（几何与流动的主导结构），只对小尺度涡用亚格子模型。物理保真度高，天然给出非定常信息，是声学、燃烧细节、尾流动力学的正选。代价是必须解析到惯性子区，近壁网格要求苛刻，计算量比 RANS 高出几个数量级，且必须跑瞬态。

混合 RANS-LES（DES 家族）。主流区用 LES 保真，近壁或附着流区退化为 RANS 省钱，用一场"分区承包"平衡预算。Fluent 里的 DES、DDES、IDDES、SBES 都属此族，是高保真需求与工程预算之间的现实折中。

三档交易的坐标系是两个轴：你要的物理保真度，和你付得起的计算预算。选型的第一性原理是把问题按轴定位，而不是按"听说哪个模型好用"定位。

图 4-1：湍流建模层级与定位

湍流黏度假设与它的裂缝

RANS 家族的绝大多数成员把雷诺应力简化为湍流黏度乘平均应变率，即 Boussinesq 假设。这个假设把七个分量的雷诺应力问题压缩成求一个标量场（湍流黏度）的问题，代价是默认湍流各向同性。麻烦恰恰出在流动里最工程重要的地方：强旋流、曲率、分离区，湍流天然各向异性，Boussinesq 假设在这里系统性失灵。两个出路：一是换雷诺应力模型（RSM）老老实实解雷诺应力输运方程，各向异性是自带的，代价是多解六七个方程、收敛更娇气；二是继续用两方程模型但打开修正开关（曲率修正、旋流修正），用经验补假设。判断顺序建议：先问流动是否强各向异性，是则 RSM 或修正开关，否则两方程模型足矣。

一份选型方法论案例：两次驳回的经过

背景。某离心风机改造项目，目标是预测整机效率与噪声源强度分布。年轻工程师第一版用标准 k-epsilon 稳态算完，效率偏差百分之六，噪声源云图一片均匀，评审被驳回。

操作。按方法论逐层定位。第一问：需要非定常信息吗？噪声源强度需要，RANS 稳态天然给不了——升级为混合模型是必然。第二问：流动强各向异性吗？叶轮通道带旋转与曲率，是。第三问：预算允许 LES 全域吗？不允许，全环几何网格量超预算。结论定为：近壁与附着区用 SST，尾迹与分离区用 LES 的混合模型（SBES），瞬态推进，采样足够转数做统计。

结果。混合模型下效率预测偏差收窄到百分之一点五，噪声源分布呈现出清晰的叶片通过频率特征，与风洞声学测量的频谱形态一致，方案获得评审通过。

解读。两次驳回的本质都是"模型档位与问题需求错配"：第一版用最低档模型回答需要高保真信息的问题。方法论的价值不是背出哪个模型好，而是三连问逼出需求档位：要什么物理信息、流动是否各向异性、预算几何——三个答案交集就是模型档位。

变式。若项目只要压降与流量，第一版稳态 RANS 本来就够，问题出在需求升级而模型没跟着升级；反过来，若只是初步筛型，上 LES 属于预算错配，钱该花在更多设计方案的对比上。

方法论问答与档位对照

"听说 SST 万金油"对不对。对一半：SST 的近壁处理与分离预测确实是两方程里的全能选手，但强旋流各向异性场景它同样受 Boussinesq 假设拖累——万金油的意思是"默认候选"，不是"免疫误用"。

LES 需要什么网格条件。经验门槛是边界层内展向与流向的分辨率指标（网格单元数在边界层厚度内达到数十量级），比 RANS 的十五层苛刻一个量级；达不到就老实用混合模型，装样子的 LES 比诚实的 RANS 更误导。

稳态还是瞬态怎么定。看关注量：平均性能量（效率、压降）稳态可行；瞬态结构（载荷脉动、声源、火焰位置波动）必须瞬态。判断不了的折中方案是稳态先算、提取关键区域做瞬态快照复核。

三档建模的付出与回报对照：

方法论的最后一块拼图是"模型升级的触发器"。什么信号说明当前档位不够用？三个可观察指标：其一，工程量对模型选择敏感（换一档模型结论翻转），说明物理超出当前模型的能力圈；其二，非定常特征对结论有实质影响（平均量在瞬态解与稳态解间系统性漂移）；其三，实验对表的偏差呈现结构化模式（不是均匀散布而是随工况规律变化）。三者任一出现，就该启动 4.1 的三连问重新定位档位，而不是继续在当前档位里调参——参数救不了档位错配，这是湍流建模的第一戒律。

本节验收点

• 项目需求清单里写明了需要平均量还是非定常量，模型档位据此选定。
• 各向异性判断有依据：旋流、曲率、分离三类特征逐项核对。
• RANS 选型时修正开关（曲率、旋流）状态已记录。
• 模型升级或降级的理由写进了建模说明，而非凭手感。