4. 材料属性与温度相关参数设置

在传热与流体联合仿真中，材料属性往往不是一组固定常数。金属的导热系数会随温度升高而下降，液体的黏度会随温度升高而急剧降低，气体的密度则随温度升高而减小。如果整场都用常温下的单一数值，仿真结果在高温区或大温差工况下会明显偏离实际。本节要解决的核心问题是：如何把随温度变化的导热系数 k(T)、比热容 Cp(T)、密度 ρ(T) 和动力黏度 μ(T) 正确写入材料节点，并把固体域与流体域的材料分配到正确的几何实体上，同时让边界条件与材料属性保持一致。

一、温度相关属性的三种定义方式

COMSOL 中定义温度相关材料属性，本质上是在“材料属性值”这一栏里填入一个关于温度 T 的表达式或函数。常用方式有三种：

• 解析表达式：直接在属性输入框中写 T 的多项式，例如导热系数 k = 45 - 0.02*(T-293.15)，单位 W/(m·K)。适合温度范围窄、有经验公式的情形。
• 插值函数：在“全局定义—函数—插值”中建立一张温度—属性数据表，再在材料属性中调用该函数名，如 k_solid(T)。适合有实验数据点或材料手册表格的情形。
• 分段/阶跃函数：用于相变附近属性突变，例如用阶跃函数让比热容在熔点附近出现一个尖峰，等效处理潜热。

需要特别注意的是温度变量的写法。COMSOL 中温度默认单位为 K，若材料手册给的是摄氏度，必须写成 T-273.15 或 (T-273.15[1/K])，否则量纲会报错。建议在定义函数时就把自变量单位统一为 K，避免后续反复换算。

二、固体域与流体域的材料分配

联合仿真模型通常同时包含固体域（如金属壁、芯片、管壁）和流体域（如空气、水、油）。材料分配的原则是：每个几何域只能被一个材料节点覆盖，但同一材料可以覆盖多个域。实际操作中，在“材料”节点下新建两个材料，分别命名为“固体—铝”和“流体—水”，然后在各自的“选择”栏中指定对应的域。

分配完成后要检查两件事：一是所有域是否都被材料覆盖（未覆盖的域会报“材料未定义”错误）；二是流体域的材料是否同时包含传热所需的 k、Cp、ρ 和流动所需的 μ。如果只做传热不做流动，μ 可以留空；但一旦耦合了层流或湍流接口，μ 缺失会直接导致求解失败。

三、参数关系与典型示例数据

下面给出一组教学示例参数（属于合理假设，非实测数据），用于说明温度相关属性的量级和变化趋势。以水在 293.15 K 到 353.15 K 范围内为例：

可以看到，黏度从 1.002 mPa·s 降到 0.355 mPa·s，降幅接近 65%，而导热系数只变化约 12%。这意味着在自然对流或强制对流问题中，黏度的温度依赖性对流速分布和压降的影响远大于导热系数。若把黏度当作常数，高温区的流量会被严重低估。

四、实操：从材料库到自定义函数

以水为例，COMSOL 内置材料库中已有“Water, liquid”材料，其属性默认多为常数。要改为温度相关，操作路径是：展开该材料的“基本”属性，找到导热系数 k，把输入框中的常数改为用户定义的插值函数 k_w(T)。插值函数的建立步骤为：在“全局定义”下右键“函数—插值”，数据源选择“文件”或“表格”，输入温度列和属性列，单位分别设为 K 和 W/(m·K)，函数名设为 k_w。

对于黏度，若采用指数形式，可写成 μ = 0.001*exp(-0.02*(T-293.15))，单位 Pa·s。这个表达式在 293.15 K 时给出 0.001 Pa·s，在 353.15 K 时约为 0.0003 Pa·s，与上表趋势一致。写完后点击“绘制”按钮，可以直接看到函数曲线，确认没有量纲错误或数量级异常。

五、结果判读与一致性检查

设置完成后，建议先做一次稳态验证：给流体域一个已知入口温度和速度，观察出口温度与压降。若出口温度明显低于手算的换热估算值，优先检查比热容和密度是否被误设为常数；若压降远小于预期，检查黏度是否仍为常温值。另一个有效手段是在“结果”中绘制温度场和速度场，观察高温区黏度降低是否引起流速局部增大——这是温度相关黏度生效的直观证据。

最后要保证边界条件与材料属性在温度范围上一致。例如入口温度设为 353.15 K，而材料插值函数只定义到 333.15 K，求解器会在超出范围处取端点值，结果虽能算出但物理意义存疑。养成在“材料—属性—函数”中检查定义域的习惯，是保证联合仿真可信度的基本要求。