2. RANS方程封闭问题与涡黏假设
在上一节中，我们已经明确：直接数值模拟（DNS）对工程问题几乎不可承受，因此工程上普遍采用雷诺平均（RANS）方法。RANS 的核心思想是把瞬时速度分解为平均量与脉动量：
u_i = U_i + u'_i
p = P + p'其中大写字母表示时间平均量，带撇号表示脉动量。把该分解代入不可压缩 Navier–Stokes 方程并做时间平均，就得到雷诺平均方程：
∂U_i/∂x_i = 0
ρ(∂U_i/∂t + U_j ∂U_i/∂x_j)
= -∂P/∂x_i + ∂/∂x_j [ μ(∂U_i/∂x_j + ∂U_j/∂x_i) ] - ρ ∂(u'_i u'_j)/∂x_j与原始 N–S 方程相比，多出了一项 -ρ ∂(u'_i u'_j)/∂x_j，其中 u'_i u'_j 称为雷诺应力张量。它由速度脉动的二阶统计量构成，代表湍流脉动对平均流动的动量输运。问题在于：未知量从原来的 4 个（三个速度分量加压力）增加到 10 个（4 个平均量 + 6 个独立雷诺应力分量），而方程数仍然只有 4 个。这就是 RANS 的封闭问题——方程组不封闭，必须引入额外的模型假设，把雷诺应力与已知的平均量联系起来。
Boussinesq 涡黏假设
最广泛使用的封闭方式来自 Boussinesq（1877）的类比思想：既然分子黏性通过分子随机运动产生剪切应力，那么湍流脉动是否也能用类似形式表达？于是他提出，雷诺应力可以写成平均应变率的线性函数：
-ρ u'_i u'_j = μ_t (∂U_i/∂x_j + ∂U_j/∂x_i) - (2/3) ρ k δ_ij其中 μ_t 是湍流黏度（涡黏度），k = ½ u'_i u'_i 是湍动能，δ_ij 是克罗内克符号。这个式子就是 Boussinesq 假设。它把 6 个未知的雷诺应力分量，压缩为两个标量未知量：μ_t 和 k。封闭问题由此从“求 6 个应力”转化为“求 2 个湍流量”，难度大幅降低。
需要特别强调：μ_t 不是流体的物理属性，而是流动状态的函数。同一流体在不同流速、不同几何、不同湍流强度下，μ_t 可以相差几个数量级。它反映的是湍流脉动对动量的输运能力，而不是分子间的黏性作用。这一点是理解 Fluent 中所有涡黏模型（k-ε、k-ω、SST 等）的基石。
湍动能 k 与耗散率 ε 的角色
Boussinesq 假设把问题归结为求 μ_t，但 μ_t 本身仍未知。工程上最常用的做法是量纲分析：μ_t 应正比于湍流速度尺度与长度尺度的乘积。在 k-ε 模型中，速度尺度取 √k，长度尺度取 k^(3/2)/ε，于是：
μ_t = ρ C_μ k² / ε, C_μ ≈ 0.09这里 k 是湍动能，代表湍流脉动的能量强度；ε 是湍动能耗散率，代表脉动能量转化为热能的速率。二者各自满足一个输运方程，与平均流方程联立求解。因此 k-ε 模型是“两方程模型”：平均流方程 + k 方程 + ε 方程，共 6 个方程，对应 6 个未知量（U、V、W、P、k、ε），方程组封闭。
在 Fluent 中，k 和 ε 不是直接设定的物性，而是通过边界条件或初始化给定的场量。入口处通常给定湍流强度 I 与水力直径 D_h，Fluent 据此估算：
k = 1.5 (U_ref · I)²
ε = C_μ^(3/4) · k^(3/2) / (0.07 · D_h)例如，入口风速 10 m/s、湍流强度 5%、水力直径 0.2 m，则 k ≈ 1.5 × (10×0.05)² = 0.375 m²/s²，ε ≈ 0.09^(0.75) × 0.375^1.5 / (0.07×0.2) ≈ 1.63 m²/s³。这些值会直接影响近壁区和剪切层的湍流黏度分布，进而影响分离点、再附长度和压降预测。
在 Fluent 中如何观察与判断
求解完成后，Fluent 的“Turbulence”面板会输出 k、ε、μ_t 以及湍流黏度比 μ_t/μ 的场分布。工程上常用 μ_t/μ 判断湍流模型是否合理：在充分发展的管流或边界层中，该比值通常在 10～1000 量级；若在主流区出现 μ_t/μ < 1，说明湍流模型几乎未激活，可能入口湍流参数给得过低，或网格分辨率不足以捕捉剪切层。反之，若近壁区 μ_t/μ 异常大，则可能是壁面处理与 y+ 不匹配，需要检查第一层网格高度。
另一个关键观察量是 k 的分布。在剪切层、分离剪切层和尾迹区，k 应显著高于主流区；若 k 场过于平坦，说明产生项被低估，往往与速度梯度解析不足或模型常数不适用有关。对于 k-ε 模型，标准常数 C_μ=0.09, C_1ε=1.44, C_2ε=1.92, σ_k=1.0, σ_ε=1.3 是默认值，但在强旋流、强曲率或逆压梯度明显的流动中，这些常数可能需要校准，这正是后续“工程校准”环节要处理的问题。
μ_t 概括湍流对平均流的整体影响。k-ε、k-ω、SST 等模型之间的差异，本质上就是 μ_t 的构造方式不同：k-ε 用 k²/ε，k-ω 用 k/ω，SST 则在近壁用 k-ω、远场用 k-ε 的混合形式。
μ_t 当作流体物性输入。在 Fluent 中，μ_t 是求解器根据 k、ε（或 ω）场计算出来的结果，用户只能通过边界条件、初始化和模型常数间接影响它。若发现 μ_t 分布异常，应优先检查入口湍流参数、壁面 y+ 和网格分辨率，而不是试图直接“设定” μ_t。
总结来说，RANS 的封闭问题迫使我们引入模型假设；Boussinesq 假设通过涡黏度 μ_t 把 6 个雷诺应力未知量压缩为 2 个湍流量 k 和 ε；Fluent 中的两方程模型正是围绕 k、ε（或 ω）的输运方程展开求解。理解这条从“不封闭”到“封闭”的逻辑链，是后续选择 k-ε、k-ω 还是 SST 的判断基础。