3. 变比热容条件下的积分热负荷计算
为什么变比热容时不能直接用算术平均
热负荷的严格定义是焓差积分:
Q = ṁ · ∫T₁T₂ cp(T) dT
当 cp 随温度近似线性变化时,算术平均比热容法恰好精确(因为线性函数的积分等于端点值的平均乘以区间长度)。但当 cp(T) 呈明显非线性(如高温烟气、接近临界区的流体),算术平均会产生系统性偏差。本例用一组教学示例数据,分别用分段线性积分、二次多项式积分和算术平均法计算,并量化偏差。
示例工况与比热容数据
设热流体质量流量 ṁ = 2.5 kg/s,进口温度 T₁ = 500 °C,出口温度 T₂ = 200 °C。已知该流体在 200–500 °C 区间的比热容实测点如下(示例假设数据,单位 kJ/(kg·K)):
| T (°C) | 200 | 250 | 300 | 350 | 400 | 450 | 500 |
|---|---|---|---|---|---|---|---|
| cp | 1.080 | 1.105 | 1.135 | 1.170 | 1.212 | 1.260 | 1.315 |
注意 cp 的增量在高温段更大(200→250 增 0.025,450→500 增 0.055),说明曲线向上弯曲,算术平均会低估积分值。
三种方法的计算过程
方法一:分段线性积分(梯形法)。每段用两端 cp 的算术平均乘以温差,累加:
ΔT = 50 K,共 6 段
段1: (1.080+1.105)/2 × 50 = 54.625
段2: (1.105+1.135)/2 × 50 = 56.000
段3: (1.135+1.170)/2 × 50 = 57.625
段4: (1.170+1.212)/2 × 50 = 59.550
段5: (1.212+1.260)/2 × 50 = 61.800
段6: (1.260+1.315)/2 × 50 = 64.375
∫cp dT = 54.625+56.000+57.625+59.550+61.800+64.375 = 353.975 kJ/kg
Q = 2.5 × 353.975 = 884.94 kW
方法二:二次多项式拟合积分。对 7 个点做最小二乘二次拟合(以 T 为变量,单位 °C),得到示例拟合式:
cp(T) ≈ 0.0000008·T² + 0.00025·T + 1.0300
(拟合仅用于演示,实际工程应使用物性数据库或厂家数据)
∫[200→500] cp dT = [0.0000008/3·T³ + 0.00025/2·T² + 1.0300·T] 从200到500
= (0.0000002667·1.25e8 + 0.000125·2.5e5 + 1.03·500)
− (0.0000002667·8e6 + 0.000125·4e4 + 1.03·200)
= (33.33 + 31.25 + 515) − (2.133 + 5.0 + 206)
= 579.58 − 213.13 = 366.45 kJ/kg
Q = 2.5 × 366.45 = 916.13 kW
方法三:算术平均比热容法:
cp,avg = (1.080 + 1.315)/2 = 1.1975 kJ/(kg·K)
Q = ṁ · cp,avg · (T₁ − T₂) = 2.5 × 1.1975 × 300 = 898.13 kW
结果对比与判读
| 方法 | ∫cp dT (kJ/kg) | Q (kW) | 相对分段线性基准偏差 |
|---|---|---|---|
| 分段线性积分(基准) | 353.98 | 884.94 | — |
| 二次多项式积分 | 366.45 | 916.13 | +3.5% |
| 算术平均比热容法 | 359.25 | 898.13 | +1.5% |
判读要点:本例中算术平均法比分段线性基准高约 1.5%,而二次拟合高约 3.5%。两者偏差方向不同,说明偏差大小与符号取决于 cp(T) 的弯曲方向和拟合质量。工程上应以分段线性梯形法为保守基准(段数足够时误差可控),多项式拟合仅作为平滑外推的辅助手段,且必须核对拟合残差。当 cp 变化幅度超过约 10%(本例 1.080→1.315,变化 21.8%)时,算术平均法的偏差不可忽略,应改用分段积分。
常见错误:① 把温度单位混用(°C 与 K 在温差上等价,但在多项式拟合中 T 的绝对值会影响系数,必须与拟合式约定一致);② 分段数过少(如只分 2 段)导致梯形法本身误差偏大,建议每段温差不超过 50–100 K,或使每段 cp 变化小于 2%;③ 直接用算术平均法却不做偏差评估,在高温烟气、近临界流体等场景可能低估或高估热负荷数个百分点,进而影响换热面积选型。
可复用的计算校核清单
- 确认 cp(T) 数据来源与温度单位,记录数据点数量与间距。
- 先算算术平均法结果,作为快速估算与量级校核。
- 用分段梯形法(段温差 ≤ 50–100 K)计算积分基准值。
- 若需平滑曲线,做多项式拟合并检查残差;拟合阶数不宜过高以免振荡。
- 比较各方法偏差,偏差 > 2% 时以分段积分为准并在计算书中注明。
- 最终 Q 与冷流体侧吸热量做能量平衡校核,偏差应在合理热损失范围内。