3. 热仿真基本概念与物理模型
数据中心热仿真的本质，是在计算机中求解机房或机柜内部空气的流动与传热过程。无论使用哪种商业软件，其底层都依赖同一套物理模型：以计算流体力学（CFD）为框架，用质量守恒、动量守恒和能量守恒三个方程描述空气的速度、压力与温度分布，再通过湍流模型处理工程中无法直接解析的湍流脉动，最后用边界条件把真实机房“翻译”成数学问题。理解这四层结构，才能判断仿真结果是否可信，而不是只看云图是否好看。
一、计算流体力学基础：把连续空间离散化
CFD 的核心思想是“离散化”。空气在机房中是连续介质，其状态由速度矢量 u=(u,v,w)、压力 p、温度 T 等场变量描述，满足一组偏微分方程。解析解几乎不存在，因此需要把计算域划分成大量小控制体（网格），在每个控制体上把微分方程近似为代数方程，再迭代求解。
常用的离散方法为有限体积法（FVM）：对每个网格单元，守恒方程写成“流入量 − 流出量 + 源项 = 累积量”的平衡形式。以质量守恒为例，单元内空气密度不变时，进入单元的质量流量必须等于流出的质量流量。网格质量直接决定结果可信度：在机柜进风口、地板送风口、冷热通道交界处，速度与温度梯度大，需要局部加密；而在大空间静压区可适当放粗。工程上常要求关键区域网格尺寸在 2～5 cm 量级，并做网格无关性验证——把网格加密一倍，若关键测点温度变化小于 0.5 ℃，可认为网格已足够。
二、三大守恒方程：仿真的物理骨架
质量守恒（连续性方程）：对不可压缩空气，∇·u=0，即速度场的散度为零。它约束了气流不能凭空产生或消失，是压力场求解的基础。
动量守恒（Navier–Stokes 方程）：描述空气在压力差、黏性力和体积力作用下的加速与减速。机房中主要驱动力是风机压差与浮升力。浮升力项常用 Boussinesq 近似：密度变化只保留在重力项中，ρ−ρ₀≈−ρ₀β(T−T₀)，其中 β 为热膨胀系数，约 1/300 K⁻¹。这一项决定了热空气上浮、冷空气下沉的自然对流行为，是冷通道封闭、热通道封闭效果评估的关键。
能量守恒：描述对流换热与导热。空气经过服务器时吸收热量，温度升高；其关系可简化为 Q=ṁcpΔT，其中 ṁ 为质量流量，cp 为定压比热（约 1005 J/(kg·K)）。例如某机柜风量 0.3 m³/s、空气密度 1.2 kg/m³，则 ṁ≈0.36 kg/s；若服务器发热 5 kW，则温升 ΔT≈5000/(0.36×1005)≈13.8 ℃。这个估算在仿真前就能判断送风温度是否足够低。
三、湍流模型：工程精度的折中
机房气流雷诺数通常在 10⁴～10⁶，属于充分湍流。直接数值模拟（DNS）需要极小网格和极短时间步，工程上不可行。常用做法是雷诺平均（RANS），把瞬时量分解为时均量与脉动量，额外引入湍流黏度 μt 封闭方程。
工程中最常用的是 k–ε 模型：求解湍动能 k 与耗散率 ε，适合大空间、高雷诺数的主流区域，计算稳定、成本低，但对分离流和近壁面精度一般。对机柜进风、送风口射流等强曲率区域，可用 k–ω SST 模型，它在近壁面用 k–ω、远场用 k–ε，兼顾边界层与自由剪切流。若关注温度分层与浮升力主导的流动，可选用 RNG k–ε 或 Realizable k–ε，它们对浮升力和旋流的修正更好。选择原则是：先保证边界条件与网格合理，再考虑湍流模型升级；模型差异通常小于边界条件误差。
四、边界条件类型：把真实机房翻译成数学问题
边界条件是仿真误差的最大来源。常见类型包括：
- 速度入口：给定送风口风速与方向，如地板穿孔砖出风 2.5 m/s、垂直向上。
- 质量流量入口：给定风量，适合已知风机曲线的场景。
- 压力出口：给定回风静压，如吊顶回风 0 Pa 表压。
- 壁面：无滑移壁面，可设绝热或给定热流密度。
- 多孔跳跃：把服务器、穿孔地板简化为多孔介质，用阻力系数 Δp=K·v² 描述压降。
- 热源：服务器按体积热源或面热源给定发热功率，如每 U 空间 200 W。
五、从物理模型到仿真判断
把上述模型组合起来，一次机房热仿真就是：几何建模 → 网格划分 → 选择湍流模型 → 设置边界条件 → 迭代求解 → 后处理判读。判读时重点关注：机柜进风温度是否超过 ASHRAE 推荐上限（如 A1 级 27 ℃）、冷通道是否存在局部热点、回风温度是否与空调负荷匹配。若某机柜进风温度仿真值 32 ℃，而现场实测 30 ℃，差异可能来自服务器多孔阻力系数偏大或送风量偏小，应优先复核边界条件而非直接改湍流模型。
需要强调的是，三大守恒方程与湍流模型属于“物理骨架”，边界条件属于“工程输入”。骨架决定仿真能算什么，输入决定仿真算得准不准。实际项目中，边界条件的误差往往远大于湍流模型之间的差异。因此，在完成一次机房热仿真后，应把机柜进风温度、冷通道温差、空调回风温度与现场实测做对比，若偏差超过 2 ℃，优先检查送风量、服务器阻力系数和热源功率是否设置合理，再考虑更换湍流模型或加密网格。只有物理模型与边界条件同时可信，仿真结果才能用于指导冷通道封闭、送风温度调整或机柜布局优化。