3. 传热学基础回顾

管路系统中的传热问题,几乎不会以单一传热方式孤立出现。以一根输送高温流体的裸管为例:管壁内部靠导热把热量从内壁传到外壁;管外空气或水靠对流把热量带走;同时外壁还向周围环境辐射散热。三种方式同时发生、彼此串联或并联,最终决定管外壁温度、散热量以及流体沿程温降。仿真建模时若只保留对流而忽略辐射,在高温或真空/低对流工况下会明显低估散热,导致出口温度预测偏高。

三种传热方式在管路中的物理位置

先明确各方式发生在哪个环节,这是搭建热阻网络的前提。

  • 导热:发生在管壁、保温层、结垢层等固体内部。圆筒壁稳态导热热阻为 R = ln(r₂/r₁)/(2πλL),其中 λ 为材料导热系数,L 为管长。管壁薄、金属 λ 大(碳钢约 45~50 W/(m·K)),热阻通常很小;保温层 λ 小(岩棉约 0.04 W/(m·K)),是主要热阻来源。
  • 对流:发生在流体与固体壁面的交界处,分管内强制对流和管外自然/强制对流。热阻为 R = 1/(hA),h 为表面传热系数。管内水流速 1~2 m/s 时 h 可达 2000~5000 W/(m²·K),而管外静止空气自然对流 h 仅 5~25 W/(m²·K),两者可差两个数量级。
  • 辐射:发生在管外壁与周围物体(墙体、地面、天空)之间,热阻为 R = 1/(h_r A),辐射换热系数 h_r = εσ(T_s²+T_sur²)(T_s+T_sur),ε 为发射率,σ = 5.67×10⁻⁸ W/(m²·K⁴)。温度越高,辐射占比越大。

耦合关系:串联为主,局部并联

从管内热流体到管外环境,热量依次穿过“管内对流 → 内壁结垢导热 → 管壁导热 → 保温层导热 → 外壁结垢导热 → 管外对流”,这是一条串联热路。而管外壁向环境的散热,对流与辐射是并联的——同一壁温下两条路径同时向外散热,总热阻为两者并联:1/R_out = 1/R_conv + 1/R_rad。这正是耦合关系的关键:辐射不是附加项,而是与对流并列的分流支路。

下面这张热阻网络图把上述结构完整表达出来,节点为温度点,支路为热阻,串联段与并联段一目了然。

管路传热热阻网络(串联 + 管外并联) T_f1 热流体 T_w1 内壁 T_w2 外壁 T_s 保温外表面 T_f2 环境 R_conv,in R_wall R_insul R_conv,out R_rad 并联 热流 Q 方向(由热流体指向环境) 总热阻 R_total = R_conv,in + R_wall + R_insul + [ 1/(1/R_conv,out + 1/R_rad) ] 总传热系数 K = 1 / (R_total · A_ref),A_ref 取外表面或内表面需统一 注:结垢热阻 R_foul 可分别串入内壁、外壁对应位置

总传热系数 K 与热阻网络的对应

工程上更常用总传热系数 K 描述整体传热能力,其定义与热阻网络完全等价:K = 1/(R_total·A_ref)。以圆筒壁为例,若以管外表面 A_o 为基准,则

1/(K_o·A_o) = 1/(h_i·A_i) + ln(r_o/r_i)/(2πλL) + 1/(h_o·A_o)

其中管外 h_o 为对流与辐射的合成:
h_o = h_conv,out + h_rad
(因两者并联,换热系数直接相加)

注意 A_i 与 A_o 不相等,这是圆筒壁与平壁的关键差别。仿真中若把 K 直接按平壁处理,薄壁大管径时误差可达百分之几到十几。

教学示例:保温管散热估算

设一根碳钢管,内径 d_i = 0.10 m,外径 d_o = 0.114 m,外包岩棉保温层厚 50 mm(外径 d_ins = 0.214 m),管长 L = 10 m。管内热水 T_f1 = 90 ℃,h_i = 3000 W/(m²·K);管外环境 T_f2 = 20 ℃,自然对流 h_conv,out = 10 W/(m²·K),外表面发射率 ε = 0.9。碳钢 λ = 50 W/(m·K),岩棉 λ = 0.04 W/(m·K)。

先算各串联热阻(以每米管长计,单位 K/W):

R_conv,in = 1/(h_i·π·d_i) = 1/(3000×π×0.10) ≈ 0.00106
R_wall    = ln(0.114/0.10)/(2π×50) ≈ 0.00042
R_insul   = ln(0.214/0.114)/(2π×0.04) ≈ 2.505

管外表面温度先假设 T_s ≈ 40 ℃(313 K),环境 T_sur = 293 K:
h_rad = 0.9×5.67e-8×(313²+293²)×(313+293) ≈ 5.9 W/(m²·K)
h_o = 10 + 5.9 = 15.9 W/(m²·K)
R_out = 1/(h_o·π·d_ins) = 1/(15.9×π×0.214) ≈ 0.0936

总热阻 R_total ≈ 0.00106 + 0.00042 + 2.505 + 0.0936 ≈ 2.600 K/W(每米)。每米散热量 q' = (90−20)/2.600 ≈ 26.9 W/m,10 m 管段约 269 W。可见保温层热阻占总热阻 96% 以上,是控制环节;管内对流和管壁热阻几乎可忽略。若去掉保温层,R_insul 消失,R_out 按裸管 d_o = 0.114 m 重算约 0.176 K/W,总热阻降到约 0.18 K/W,散热量增大约 14 倍——这就是保温层在热阻网络中的决定性作用。

辐射项 h_rad 依赖表面温度,而表面温度又由热平衡决定,因此含辐射的管路传热是弱非线性问题。仿真中通常先假设 T_s 迭代 2~3 次即可收敛;若 T_s 与环境温差小于 20 ℃,辐射占比通常低于 15%,可先忽略再校核。

仿真建模中的三个易错点

  1. 面积基准不统一:K 必须绑定 A_ref。管内 h_i 用 A_i,管外 h_o 用 A_o,混用会导致 K 偏差。建议统一折算到外表面。
  2. 忽略辐射或重复计入:管外 h_o 应等于 h_conv + h_rad,而不是取两者较大值或相乘。低温小温差工况可忽略辐射,但高温烟气管道必须计入。
  3. 结垢热阻漏项:运行一段时间后内壁结垢 R_foul 可达 0.0002~0.001 m²·K/W,虽数值小,但与管内对流热阻同量级,长期仿真不可忽略。

把上述热阻网络直接映射到仿真工具中,就是一组串联热阻加一个并联节点。理解了这个结构,后续设置边界条件、判断哪一侧是控制热阻、解释出口温度对保温厚度或流速的敏感度,都会变得有据可依。