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)200250300350400450500
cp1.0801.1051.1351.1701.2121.2601.315

注意 cp 的增量在高温段更大(200→250 增 0.025,450→500 增 0.055),说明曲线向上弯曲,算术平均会低估积分值。

c_p(T) 曲线:分段线性积分 vs 算术平均 温度 T (°C) c_p (kJ/(kg·K)) 200 250 300 350 400 450 500 1.05 1.15 1.25 1.35 c_p,avg=1.1975 分段线性 c_p(T) 积分(真实焓差) 算术平均比热容法

三种方法的计算过程

方法一:分段线性积分(梯形法)。每段用两端 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.98884.94—
二次多项式积分366.45916.13+3.5%
算术平均比热容法359.25898.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 与冷流体侧吸热量做能量平衡校核,偏差应在合理热损失范围内。