4. 湍流模型选择与壁面处理

机柜级热仿真中，空气流过服务器、交换机、电源模块时，雷诺数通常落在 10³～10⁵ 量级，属于充分发展的湍流。湍流会显著增强动量与热量掺混，直接决定机柜进出风温差、局部热点位置和风扇工作点。因此，湍流模型不是“选一个默认项”即可，它决定了近壁面速度梯度、边界层分离、尾流掺混和换热系数 h 的预测精度。本节围绕三类常用模型——标准/Realizable k-ε、k-ω SST、雷诺应力模型（RSM）——说明它们在机柜尺度与芯片尺度上的适用边界，并给出壁面函数与网格分辨率的匹配方法。

一、三类模型的物理能力差异

k-ε 模型求解湍动能 k 与耗散率 ε，采用涡黏假设，把雷诺应力表示为平均应变率的线性函数。它的优势是鲁棒、收敛快、对自由剪切流（如机柜顶部排风、机房回风）预测稳定；缺点是近壁面处理依赖壁面函数，对逆压梯度、分离流和强各向异性流动预测偏差较大。机柜内服务器进风口、盲板边缘、线缆束后方常出现分离与回流，k-ε 容易高估掺混、低估局部高温。

k-ω SST 模型在近壁区使用 k-ω 形式，在远场切换为 k-ε 形式，并通过剪切应力输运限制涡黏系数。它对逆压梯度边界层、分离点和再附长度更敏感，适合机柜内挡板、服务器面板开孔、风扇入口等存在局部分离的区域。代价是近壁网格要求更严，y⁺ 通常希望控制在 1 附近，否则壁面函数与模型假设不匹配。

雷诺应力模型直接求解六个雷诺应力分量，放弃各向同性涡黏假设，能刻画二次流、强旋流和浮力驱动流动。机柜内热空气上升、冷通道下沉形成的自然对流叠加强迫对流，以及风扇尾旋，都是 RSM 的潜在优势场景。但 RSM 方程多、刚性大，收敛慢，对网格质量和初场更敏感，通常用于局部关键区域或作为校核手段，而非整柜全尺寸常规计算。

二、壁面函数与网格分辨率的匹配

壁面处理的核心是 y⁺ 与壁面函数的一致性。y⁺ 定义为 y⁺ = ρ u_τ y / μ，其中 u_τ 为摩擦速度，y 为第一层网格中心到壁面的距离。机柜内空气物性近似取 ρ=1.18 kg/m³、μ=1.85×10⁻⁵ Pa·s。若目标 y⁺=30，壁面剪切速度 u_τ 约 0.5 m/s，则首层网格高度 y ≈ y⁺ μ /(ρ u_τ) ≈ 30×1.85×10⁻⁵/(1.18×0.5) ≈ 0.94 mm。若目标 y⁺=1，则首层高度约 0.03 mm，边界层内至少需要 10～15 层，总厚度覆盖到 y⁺≈100～200。

标准壁面函数假设第一层网格落在对数律层（y⁺≈30～300），此时壁面换热通过经验关系计算，对 y⁺ 不敏感但无法解析黏性底层。增强壁面处理（Enhanced Wall Treatment）或 SST 的低 Re 形式则要求 y⁺≈1，直接积分到壁面，能更准确给出壁面温度梯度和换热系数。实际工程中常见错误是：用 k-ε + 标准壁面函数，却把首层网格做到 0.1 mm，导致 y⁺≈3，落在缓冲层，壁面函数与网格双重不匹配，换热系数可能偏差 20% 以上。

三、机柜尺度与芯片尺度的差异化选择

机柜尺度关注的是流量分配、进出风温度和整体热分布。此时计算域大、网格量受限，优先选择 Realizable k-ε 配合标准壁面函数，首层网格按 y⁺≈50 设计，可在数小时内完成整柜稳态计算。若发现某台服务器进风温度预测偏高、局部回流明显，可对该区域做局部加密并切换为 k-ω SST，采用分区建模或嵌套网格，而不是整柜全部改用 SST。

芯片尺度关注的是结温、壳温和散热器热阻。散热器翅片间距常为 1～3 mm，翅片间流速低、边界层厚，必须解析边界层。推荐 k-ω SST + y⁺≤1，边界层 10 层以上，增长比 1.2。对于封装内部、焊球和基板，特征尺寸进入毫米以下，空气流速极低，可能已接近层流，此时应改用层流模型或低 Re 修正，强行套用湍流模型反而引入虚假掺混。

四、一个可复现的校核流程

以某 42U 机柜为例，单柜功率 8 kW，风量 3000 m³/h，进风温度 25 ℃。第一步，用 Realizable k-ε + 标准壁面函数做整柜初算，首层高度 1 mm，目标 y⁺≈50，得到各服务器进风温度与风量分配。第二步，提取进风温度最高的服务器，将其内部风道单独建模，改用 k-ω SST，首层高度 0.02 mm，边界层 12 层，目标 y⁺≈1，重新计算该服务器内 CPU 散热器换热。第三步，对比两种模型在散热器表面换热系数 h 上的差异：若差异小于 10%，说明整柜 k-ε 结果可用于工程判断；若差异大于 20%，则需在整柜模型中对该区域引入修正或局部替换模型。

示例参数（教学假设，非实测）：
机柜：42U，8 kW，风量 3000 m³/h，进风 25 ℃
整柜模型：Realizable k-ε，标准壁面函数，首层 1.0 mm，y⁺≈50
局部模型：k-ω SST，低 y⁺ 壁面处理，首层 0.02 mm，12 层，y⁺≈1
对比量：散热器表面平均 h、CPU 壳温、服务器出风温度
判断阈值：h 差异 <10% 可接受；>20% 需局部修正或模型替换

五、常见误区与处理建议

• 误区一：整柜和芯片都用默认 k-ε。结果往往是整柜流量分配尚可，但芯片结温偏低，因为 k-ε 高估了翅片间掺混和换热。
• 误区二：盲目追求 y⁺=1 却未加密边界层。首层很薄但层数不足，边界层外直接跳到粗网格，梯度重建误差反而更大。
• 误区三：忽略浮力。机柜内温差 10～20 ℃ 时，热空气上升会改变局部流场，k-ε 需开启浮力修正，SST 和 RSM 对浮力更敏感，应同时检查格拉晓夫数 Gr 与 Re 的比值。
• 建议：先做网格无关性验证，至少三套网格（粗、中、细），观察 h 和关键温度随 y⁺ 的变化；再固定网格，对比 k-ε 与 SST 的差异。两者都稳定后，才把结果用于风扇选型或机柜布局决策。

总结来说，机柜尺度优先保证流量与温度场整体趋势，k-ε 是效率最高的起点；芯片尺度必须解析边界层，k-ω SST 配合 y⁺≤1 是更可靠的选择；RSM 适合强旋流、强浮力或作为局部校核。壁面函数与网格分辨率必须成对设计，y⁺ 落在缓冲层是换热预测偏差的主要来源之一。把模型选择、y⁺ 目标和网格层数写进仿真规范，比事后反复调参更能保证结果可复现。