浅析污水换热器的特点和设计计算方法

浅析污水换热器的特点和设计计算方法

①受条件限制污水换热器只能以小温差换热。例如在北方严寒地区,冬季污水进口温度只有10℃左右,而受热泵机组蒸发温度的限制,中介水温度不能太低,因而污水最大温降只能是3~4℃,两级换热每级的温差小于3℃。
②由于温差小,换热系数小,污水的降温需要较长的流程,通常为30 m左右。这就需要对壳管式换热器的内部结构作全新的设计,包括管程的安排、壳空间的结构设计等等,以保证全程逆流换热和换热器内部各处均有理想的换热温差。
③换热器内污水黏度大[2],换热面污染严重[3],导致污水侧的阻力大于清水侧,而换热系数小于清水侧。
④由于存在堵塞与污染的危险,无法像针对清水那样采取波纹管、内肋片、内插物等增强换热的设施。
  

污水-水换热器的设计计算方法
  
 计算参数分类
    一般而言,换热器设计应追求使用较少的耗材、具有较小的体积、消耗较少的水泵用电来获得预期的换热效果。虽然上述要求相互矛盾,但人们对一般的换热器早已有了成熟的设计经验。而污水换热器要达到上述要求就绝非易事,必须精确计算、精打细算,才能避免耗材过多、体积过大或泵耗过高,才能在各项指标中找到一个最佳至少较佳的平衡点。
    与耗材、体积、泵耗这三大指标相关的参数不仅有结构参数,诸如管径、管长、管间距、管程数等等,也有运行参数,诸如流速、流量、压降、各节点进出口温度、换热量等等。为了给设计计算一个清晰的思路,对这些参数作如下分类。
    1)已定或可定参数
    ①换热量Q。依据需求而定。
    ②污水流量V与温降(或温升)Δt。这两个参数是在设计污水源热泵系统时确定的。依据公式Q =ρcVΔt(式中c为比热容)可知,确定其中一个量,另一个量可求。
    ③传热平均温差Δtm依据污水初始温度、中介水温度、蒸发温度等在系统设计中确定,一般为3~4℃。
    ④管径d。依据经验与定型管材的规格确定,一般选取外径dw为20 mm,内径di为16 mm。
    ⑤管间距。考虑壳空间的流通面积,一般壳管式换热器管子的排列间距由s/dw=1.5(式中s为管心距)确定,换热管按正三角形排列。在这种布管方式下,管空间与壳空间之比为0.4∶0.6。
    2)计算取值参数
    ①特征温度。由于污水和中介水温度变化不大,其特征温度均可按进口温度选取。
    ②污水黏度。由实验测得。课题组选用旋转圆筒黏度计测试了不同温度下的污水黏度,测试结果表明换热器设计中,可按污水黏度为清水黏度的2.5倍左右计算。
    ③壳侧流速。一般换热器设计中,壳侧的推荐流速为管侧的1/2[4]。然而由于污水换热器中装有纵向折流板,在同样的流速下,压降必然高于一般换热器,因而此处推荐流速u2=u1/3(式中u2为壳侧流速,u1为管内流速)。
    ④污垢热阻。工程经验表明,考虑污垢热阻时,换热器总传热热阻为两侧对流换热热阻的1.2倍。
    ⑤实验验证糙粒高度Ks取0.001 m时按照希弗林松公式计算所得阻力系数与实测值相近,最大误差在10%以内,表1给出了计算值与实测值的比较。
   
    3)待定参数
    除了上述参数,换热器设计还要涉及许多重要参数,诸如传热系数、换热面积等等。仔细分析这些参数,发现它们并不是各自独立的,而是互相关联的,再经分析知道它们当中基本的参数只有3个,即管根数N、每根管的管长L(含折返)、管内流速u1。其余参数均可由这3个基本参数计算出来。
    ①中间计算参数
    传热系数K取决于管侧表面传热系数h1、壳侧表面传热系数h2,而h1, h2分别取决于管内流速u1和壳侧流速u2。
    ②效果参数
    管侧阻力Δp关系到功耗,由沿程阻力系数f,L与u1确定;换热面积A关系到投资,由L,u1与N确定;换热器管程数n和换热器组数需综合考虑N,L及设备允许占地空间后确定。
    2.2　设计所用方程及公式
    设计涉及的所有参数与基本参数之间的关系可由以下公式计算。
   
    3)换热计算
    ①污水侧表面传热系数h1一般来说,管内流速u1在0.8~1.3 m/s之间,对应的最小雷诺数Remin=3 920>2 300,最大雷诺数Remax=6 370<104,由此可知管内流动属于过渡流,选用的计算公式为：
    （4）
    式中　λ,Nu, Re, Pr分别为导热系数、努塞尔数雷诺数和普朗特数。
 
   中介水侧表面传热系数h2
    有挡板的壳侧流道复杂,难以采用简单的公式进行准确的计算。特别是对污水换热器,壳侧往往需要增设纵、横较多的挡板。关于壳侧表面传热系数较为准确的计算,笔者将另文详述。此处暂粗略认为壳侧的换热为外掠圆管对流换热,并按经验选取系数φ=0.7,取壳侧最小流通截面处流速u2为管内流速的1/3,此时Remin>1 000,Remax<2×105,选用的计算式为：