3. 热流体仿真物理基础

服务器机柜热流体仿真的可信度，首先取决于是否选对了控制方程与传热模式。机柜内部同时存在风扇驱动的强制对流、发热芯片与散热器表面的自然对流、固体导热以及壁面之间的辐射换热，若把四类物理过程混为一谈，仿真结果往往会出现“风量对但温度不对”或“局部热点被抹平”的问题。本节从守恒方程出发，说明各传热模式的适用条件，并给出机柜场景中可落地的简化假设。

一、三大守恒方程：仿真的底层约束

热流体仿真的求解对象是空气（或冷板液冷中的冷却液）的速度场、压力场与温度场。质量守恒（连续性方程）要求控制体内密度变化率与流出质量通量之和为零；对机柜内常见的低速不可压流动，可简化为速度散度为零，即进入某控制体的空气体积流量等于流出的体积流量。动量守恒（Navier–Stokes 方程）描述压力梯度、黏性力、体积力与惯性力的平衡，风扇提供的压升、滤网与多孔介质的阻力、浮升力都在此体现。能量守恒则要求控制体内能变化等于对流带入/带出的焓、导热与辐射的净热流，加上内热源（芯片、电源、硬盘）的发热量。

在稳态、不可压、常物性假设下，机柜内空气的守恒方程组可写为：

连续性：  ∇·u = 0
动量：    ρ(u·∇)u = -∇p + μ∇²u + ρgβ(T-T_ref)
能量：    ρc_p(u·∇T) = ∇·(k∇T) + q_vol

其中 ρ 为空气密度，μ 为动力黏度，β 为热膨胀系数，q_vol 为体积热源。注意动量方程中的浮升力项 ρgβ(T−T_ref) 是自然对流的来源；当风扇强制对流占主导时，该项常可忽略，但机柜顶部与盲板附近流速极低，浮升力不可轻易丢弃。

二、共轭传热：固体与流体的耦合界面

共轭传热（Conjugate Heat Transfer, CHT）指固体域导热与流体域对流在同一求解中耦合，界面处满足温度连续与热流连续：T_s = T_f，k_s(∂T/∂n)_s = k_f(∂T/∂n)_f。服务器机柜中，芯片结→封装→散热器基座→鳍片→空气的路径就是典型共轭传热。若只做纯流体仿真而把散热器简化为定热流边界，会高估鳍片效率、低估界面温升。工程上常对散热器做“等效导热+等效对流”的降阶处理，但必须校核其等效热阻与实测或详细模型偏差在可接受范围内。

三、自然对流、强制对流与辐射的适用条件

判断以哪种对流为主，常用瑞利数 Ra 与格拉晓夫数 Gr 与雷诺数 Re 的比值。机柜内空气温升约 10～20 K、特征长度 0.1～1 m 时，自然对流 Ra 约在 10⁶～10⁸，属层流至过渡区；风扇送风速度 1～5 m/s、风道当量直径 0.05～0.2 m 时，Re 约在 3×10³～7×10⁴，属湍流强制对流。一般当 Re²/Gr > 10 时强制对流主导，Re²/Gr < 0.1 时自然对流主导，中间区域需同时保留浮升力项。

辐射的适用条件是存在温差且表面发射率不可忽略。机柜内金属表面发射率通常 0.05～0.3，而 PCB、塑料件、涂漆面板可达 0.8～0.9。当表面温度高于周围空气 15 K 以上、且风速低于 0.5 m/s 时，辐射换热量可占总换热量的 10%～30%，此时应开启辐射模型（如离散坐标法 DOM 或面到面 S2S）。若机柜内风速普遍大于 2 m/s，辐射占比通常低于 5%，可先忽略以节省计算量。

四、机柜场景中的简化假设与校核

完整求解机柜内所有服务器细节在工程上不可行，需做合理简化。常见假设包括：将服务器机箱简化为多孔介质或定阻力模块，用压降–风速曲线（如 Δp = a·v² + b·v）代替真实结构；将风扇简化为风扇曲线或定流量边界；将 PCB 与芯片简化为等效导热块与体积热源；忽略机柜壁面厚度与外部环境，仅保留绝热或定温边界。若机柜内平均风速大于 2 m/s，可先关闭辐射模型；若存在顶部低风速区，则应开启 S2S 辐射并校核该区域温度。

简化是否可接受，需通过网格无关性、能量平衡与关键测点校核。能量平衡要求：所有热源总功率 ≈ 进出风焓差 + 壁面散热。例如总热源 5 kW、风量 600 CFM（约 0.283 m³/s）、空气密度 1.18 kg/m³、比热 1005 J/(kg·K)，则理论温升 ΔT = 5000/(1.18×0.283×1005) ≈ 14.9 K。若仿真出风与进风温差偏离该值超过 10%，说明边界条件或简化假设需要复查。关键测点包括 CPU 散热器基座温度、机柜顶部空气温度、风扇前后压差，三者与实测或经验值的偏差应控制在 5 K 与 10 Pa 以内。

综上，机柜热流体仿真的物理基础可归纳为：以质量、动量、能量守恒为约束，以共轭传热连接固体与流体，以 Re、Gr、Ra 判断对流模式，以表面发射率与温差判断辐射是否可忽略，最后用能量平衡与关键测点校核简化假设的合理性。只有物理模式选对，后续的网格、求解器与后处理才有意义。