3. 液冷系统基础与Fluent建模准备

液冷系统的CFD建模,难点往往不在求解器本身,而在“把真实回路翻译成可计算的几何与边界”。本节围绕液冷回路的构成、冷板与泵阀参数、冷却液物性、流道简化与边界条件设定展开,目标是让读者在打开Fluent之前,已经能明确:我要算的是哪一段、入口给什么、出口给什么、壁面给什么、哪些细节可以简化、简化后误差从哪里来。

一、液冷回路的构成与建模范围界定

一个典型的服务器液冷回路可拆成五个功能段:冷板(cold plate)、快接头与软管、CDU(冷量分配单元,含板式换热器与二次侧泵)、一次侧冷冻水回路、以及膨胀罐与排气阀等辅助件。做Fluent仿真时,通常不会把整条回路都建成三维模型,而是按研究目标切分:

  • 冷板级仿真:只取冷板内部流道+冷板基板+发热芯片热源,入口接软管出口,出口接软管入口,用于评估芯片结温、流阻与流量分配。
  • CDU级仿真:取板式换热器二次侧+泵+管路,用于评估供液温度稳定性与泵功耗。
  • 系统级仿真:用一维系统工具(如Fluent的System Coupling或AMESim)与三维冷板耦合,三维只保留冷板,其余用集中参数替代。

本节以最常见的“单冷板三维仿真”为主线。建模范围一旦确定,后续所有边界条件都围绕这个范围来设定,切忌把整机机柜都建进来却只关心一颗芯片。

二、冷板与泵阀的关键参数

冷板的核心参数包括:流道形式(蛇形、平行微通道、射流冲击、针翅)、流道水力直径、基板厚度与材料、入口/出口接口规格。以某1U服务器CPU冷板为例,教学示例可取:流道为6条平行矩形微通道,单通道宽1.2 mm、高3 mm,基板为铜,厚3 mm,芯片热源面积30 mm×30 mm,发热功率250 W。这些数值属于合理假设的示例,实际项目需以图纸为准。

泵与阀的参数则决定边界条件。泵的关键参数是流量-扬程曲线(Q-H曲线)与工作点;阀的关键参数是流阻系数Kv或局部阻力系数ζ。工程上常把泵简化为“给定质量流量入口”,把阀简化为“局部阻力”或直接忽略。若泵曲线已知,也可用Fluent的Fan/Pump边界或UDF实现,但入门阶段建议先用定流量,收敛后再切换到泵曲线。

判断依据:若仿真目标是“芯片温度是否超标”,用定流量入口即可;若目标是“泵选型是否够用”,则必须用Q-H曲线并让求解器自己找工作点。

三、冷却液物性:别用默认水

液冷系统常用工质为去离子水、乙二醇水溶液(如25% EG)或介电冷却液。Fluent材料库中默认water-liquid的物性在20℃附近,但液冷系统供液温度常在35~45℃,且乙二醇溶液物性与纯水差异明显。教学示例中,25%乙二醇水溶液在40℃时可取:密度约1040 kg/m³、比热约3.75 kJ/(kg·K)、导热系数约0.48 W/(m·K)、粘度约1.2 mPa·s。这些值应作为温度相关函数或至少作为常数输入,并在报告中注明来源。

若忽略温度对粘度的影响,微通道内压降可能被低估10%~20%。因此建议:先用常物性快速收敛,再用分段线性或多项式拟合的变物性重算一次,对比压降与温升差异。

四、流道简化与边界条件设定

真实冷板内部常有圆角、螺纹、密封槽、加工毛刺,这些对宏观流动影响很小,却会让网格量暴增。简化原则是:保留影响流量分配与换热的特征,删除只影响局部应力的特征。具体做法:

  1. 删除小于0.2 mm的圆角与倒角,用直角替代。
  2. 忽略快接头内部螺纹,用等直径圆柱段替代。
  3. 若流道对称,取1/2或1/4模型并设对称边界。
  4. 芯片热源用体热源或面热流密度,不要建真实焊料层。

边界条件设定是本节最关键的实操。以单冷板为例:入口设为mass-flow-inlet,给定质量流量(如0.08 kg/s)与入口温度(如40℃);出口设为pressure-outlet,表压0 Pa;冷板外壁面设为对流换热系数与自由流温度(如h=15 W/(m²·K),T=35℃),或用绝热壁面先看最坏情况;芯片底面设为恒定热流密度q″=250 W/(0.03 m×0.03 m)≈277.8 kW/m²。所有方向约定:入口法向为+x,出口为−x,重力沿−y,全文与图示保持一致。

单冷板液冷仿真:几何简化与边界条件设定 铜基板(3 mm) 6条平行微通道:1.2 mm × 3 mm(简化后直角) 入口 mass-flow-inlet 0.08 kg/s, 40℃ 出口 pressure-outlet 0 Pa 表压 芯片热源 250 W q″ ≈ 277.8 kW/m²(30 mm × 30 mm) 外壁面:h=15 W/(m²·K) T∞=35℃ x y z 25% 乙二醇水溶液 @40℃ ρ=1040 kg/m³ cp=3.75 kJ/(kg·K) k=0.48 W/(m·K) μ=1.2 mPa·s 几何简化原则 删除 <0.2 mm 圆角/螺纹/密封槽 对称面用 symmetry,保留流量分配特征

五、网格与求解设置的衔接

边界条件确定后,网格策略随之明确:微通道内至少布置10~15层边界层网格,第一层高度按y+目标反推;若用k-ε模型,y+可放宽到30~100,但微通道建议用SST k-ω并保证y+≈1。入口段需延长3~5倍水力直径,避免入口效应污染换热结果。教学示例中,水力直径Dh≈1.7 mm,入口延长段取8 mm。

初始化时先用混合初始化,再对入口流量做2~3次局部加密迭代,观察出口质量流量是否与入口守恒(偏差应小于0.5%)。若残差下降到1e-4但出口流量不守恒,优先检查边界条件是否被误设为压力入口或壁面。

常见错误:把芯片热源直接加在流体域上,导致热量没有经过铜基板导热就直接进入冷却液,结温被严重低估。正确做法是热源加在固体域底面,流体与固体之间设coupled壁面。

六、可直接使用的建模准备检查表

检查项 示例值/设定 判断依据
建模范围单冷板+芯片热源目标为芯片结温与流阻
入口边界mass-flow-inlet, 0.08 kg/s, 40℃泵工作点已知
出口边界pressure-outlet, 0 Pa回液接常压CDU
冷却液物性25% EG, 40℃常物性先收敛后变物性复核
热源277.8 kW/m² 底面热流250 W / 0.0009 m²
外壁面h=15 W/(m²·K), T∞=35℃机箱内自然对流估算
网格边界层15层, y+≈1SST k-ω 要求
守恒检查出口/入口质量流量偏差<0.5%边界条件正确性

完成以上准备后,Fluent中的操作顺序是:读入网格→检查最小体积与负体积→设置材料→设置单元区域条件(流体/固体)→设置边界条件→初始化→迭代。每一步的输入都来自本节确定的参数,而不是临时拍脑袋。把这张检查表填满,再开始计算,能避免绝大多数“算完了才发现入口给错”的返工。