3. 冷却液物性获取与定义

在液冷电子系统的CFD建模中，冷却液物性不是一组可以随手填写的常数，而是决定换热能力、压降水平和温度场分布的关键输入。同一套冷板模型，仅把工质从纯水换成50%乙二醇水溶液，在相同流量下换热系数可能下降15%～25%，压降却明显上升。因此，本节的核心任务是：明确需要哪些物性、这些物性如何随温度变化、从哪里获取或拟合、以及如何正确写入求解器。

一、CFD真正需要的物性集合

对于单相液冷，求解器通常需要以下物性：密度ρ、定压比热容c_p、导热系数k、动力黏度μ（或运动黏度ν）。若涉及自然对流或浮升力影响，还需热膨胀系数β；若涉及湍流普朗特数修正，则需普朗特数Pr。对于两相或相变工质，还必须补充饱和温度、汽化潜热、饱和蒸气压曲线以及表面张力。

需要特别强调的是：物性必须与温度绑定。水的黏度从20℃到60℃会下降约45%，若在模型中固定为20℃的值，高温区流量分配和压降预测都会失真。因此，正确做法是建立“温度—物性”曲线或分段多项式，而不是单一常数。

二、三类典型工质的物性特征

纯水是液冷系统的基准工质。其比热容高、导热系数大、黏度低，但冰点高、电导率敏感。在5～90℃范围内，水的密度约从1000 kg/m³缓降至965 kg/m³，c_p在4.18 kJ/(kg·K)附近小幅变化，k约从0.58升至0.68 W/(m·K)，而μ从约1.0 mPa·s快速降至0.3 mPa·s。水的物性数据成熟，可直接采用IAPWS-IF97标准或NIST数据库。

乙二醇水溶液用于防冻和抑制腐蚀，常见体积浓度为30%、50%、60%。其特点是：冰点显著降低，但黏度大幅上升、比热容和导热系数下降。以50%乙二醇水溶液为例，20℃时μ约为4～5 mPa·s，是纯水的4倍以上；c_p约3.3 kJ/(kg·K)，k约0.40 W/(m·K)。浓度越高，防冻能力越强，但换热性能越差。工程上常用厂商提供的浓度—温度—物性表，再插值成曲线。

介电液（如氟化液、矿物油、硅油）用于直接接触冷却，最大优势是绝缘。其典型特征是：密度大（约1600～1800 kg/m³）、c_p低（约1.0～1.1 kJ/(kg·K)）、k低（约0.06～0.13 W/(m·K)）、μ随温度变化剧烈。介电液的换热能力远弱于水，因此往往需要更高流量或更大换热面积。其物性通常由供应商提供，且不同牌号差异明显，不能互相套用。

三、温度依赖关系的数学处理

把物性表变成求解器可用的函数，常用三种方式。第一种是分段线性插值：将温度区间划分为若干段，每段用线性函数近似。例如水的黏度在20～40℃可写为μ(T)=1.002−0.0245(T−20) mPa·s，在40～60℃另取一段。这种方式简单稳定，适合大多数工程仿真。

第二种是多项式拟合：用二次或三次多项式拟合整段曲线。例如50%乙二醇水溶液的密度可拟合为ρ(T)=1065−0.42T−0.0012T²（示例拟合式，仅用于说明形式）。多项式光滑但外推风险大，必须限制在拟合温区内使用。

第三种是查表+插值：直接读取NIST或厂商表格，在求解器中以表格形式输入。这种方式精度最高，但需要确保温度节点覆盖实际工况范围，且节点间距不宜过大。

四、纳米流体与相变工质的特殊处理

纳米流体是在基液（水、乙二醇）中加入纳米颗粒（Al₂O₃、CuO、TiO₂等）形成悬浮液。其导热系数通常高于基液，但黏度也显著增加，且存在稳定性问题。CFD中处理纳米流体物性时，常用经典模型估算：导热系数用Maxwell模型，黏度用Einstein模型或Batchelors模型，密度和比热容按体积分数加权。需要特别注意的是，纳米流体的物性对颗粒体积分数、粒径和团聚状态非常敏感，文献关联式差异较大，建议以实验数据为准，并在报告中注明所用模型和适用范围。

相变工质（如氟化液在沸腾冷却中）的处理更复杂。除单相物性外，还需定义饱和温度—压力曲线、汽化潜热、表面张力以及沸腾换热系数。CFD中通常采用两相模型（如VOF、Mixture或Eulerian），并配合相变模型（如Lee模型）定义蒸发/冷凝速率。此时物性不再是单值函数，而是随相态和温度跳跃变化，必须分别定义液相和气相物性。

五、实操：把物性写入求解器

以ANSYS Fluent为例，物性定义路径为：Materials → Create/Edit → 选择fluid → 在Density、Cp、Thermal Conductivity、Viscosity下拉框中选择“piecewise-linear”或“polynomial”，然后输入温度—物性数据点。若使用UDF，可用DEFINE_PROPERTY宏定义μ(T)和k(T)。

/* 示例：50%乙二醇水溶液黏度UDF（分段线性，仅示意） */
DEFINE_PROPERTY(mu_eg50, c, t)
{
  real T = C_T(c, t);          /* 温度，单位K */
  real Tc = T - 273.15;        /* 转为℃ */
  real mu;
  if (Tc <= 20.0)
    mu = 0.0050 - 0.00008 * (Tc - 0.0);
  else if (Tc <= 60.0)
    mu = 0.0034 - 0.00005 * (Tc - 20.0);
  else
    mu = 0.0014 - 0.00002 * (Tc - 60.0);
  return mu;                   /* 单位：Pa·s */
}

定义完成后，建议做一次“物性敏感性检查”：在相同边界条件下，分别用常数物性和温度依赖物性各跑一次，对比出口温度、压降和芯片结温。若差异超过5%，说明物性温度依赖不可忽略，必须保留曲线定义。

六、物性数据记录表示例

为保证可追溯性，建议为每个工质建立一份物性记录表，至少包含以下字段：