1、1编号热 3017高 频 焊 螺 旋 翅 片 热 管 换 热 器 特 性 的 研 究徐 通 明 牟 楷南 京 化 工 学 院中 国 工 程 热 物 理 学 会第一届热管会议传 热 传 质 学 会一九八三年八月2符 号 说 明A 换热器表面积 米 2Ah 换热器蒸发段翅片管管外表面积 米 2Ac 换热器冷凝段光管管外表面积 米 2 Cp 定压比热 千卡/ 公斤CDo 光管外径 米Di 热管内径 米d v 热管内蒸汽通道直径 米Dev 容积当量直径 米 NFV 流动净自由容积 米 3F 流体阻力系数Gm 通过最小截面积的最大质量流量 公斤/米 2时Gc 重力换算系数gc=1.27108 公斤时 2
2、h 管外换热系数 千卡/米 2时Chh 以 Ah 为基准的管外换热系数 同上hc 以 A c 为基准的管外换热系数 同上hie 以热管内蒸发段面积 为基准的管内蒸发沸腾换热系数 同上J 热的机械功当量 公斤米/千卡Lg 气流流过的长度 米Nu 努赛尔数3P 锅炉内饱和蒸汽压公斤/厘 米 2p v 热管内蒸汽压差 同上Q 热流量 千卡/时Q 热流密度 千卡/米 2时R 热阻 C/W 或 千卡/W 或米 2 C/千卡Reh 传热计算雷诺数Rep 阻力计算雷诺数R 汽化潜热 千卡/ 公斤热管间距 米SL 热管间距 米T 温度 Ctw 壁面温度 Ctf 气流温度 Cthi 换热器热流体进口温 度 Ct
3、ci 换热器冷流体进口温度 Cthcp 换热器热流体、出口平均温度 CTv 热管内蒸汽温度 Ct m 对数平均温度 CU 总传热系数 千卡/ 米2时CW 气流速度 米/秒 管壁导热系数 千卡/米2时CSTSL4 l 液体导热系数 同上 气体粘度 公斤/米秒 w 壁温下气体粘度 同上 l 液体密度 公斤/米 3 U管内蒸汽密度 同上 f翅片效率 换热器效率5一、引 言热管是一种高效传热元件,由热管组成的换热器,由于其传热效率高。阻力损失小,结构紧凑,运行周期长等特点,常作废热回收装置,在国外已获得广泛的应用。近几年来,国内对热管和热管换热器研究和实用计算公式正在蓬勃兴起。除了需要对热管的单管性能
4、进行试验以外,还必需对热管换热器整体性能进行研究。尤其在为提高传热率而广泛采用高频焊接螺旋翅片管作为管壳后,对这种换热器的整体性能研究更为迫切。本文是对气流横掠高频焊接螺旋翅片管管束的换热规律和阻力特性进行了实验研究。通过对四十四个工况下的总传热系数的测定,并分离出管外换热系数以及阻力试验,得出相应的计算式。用实验得出的准则方程式和计算式,设计式,设计了两台热管换热器,对此两台热管换热器的性能进行了现场实测,实测数据与设计数据合得较好。二 、 传 热 特 性 的 研 究1实验模型分析:热虹吸管具有传热率高,结构简单,工作可靠和价格便宜等优点,目前地面应用的热管换热器很多采用热虹管作为传热元件。
5、热虹吸管的传热模型如图 1 所示(1) 。热虹吸管的传热过程是由七个传 热环节所组成,即 蒸发段管外换热热阻hR16图 1 热虹吸管传热模型蒸发段管壁导热热阻蒸发段管内蒸发、沸腾热阻。它是由两部份组成(1) 过量充液量形成的4.023.1.04.25.02.7.3.6.3. eaevpep qRgch 膜LCLaLe iodenR02ieR137液池沸腾换热热阻,沸腾换热系数 hp通过量纲分析,可用下式计算 1 (1)(2) 回流冷凝液层流液膜蒸发换热热阻,它可以看成是冷凝换热热的逆过程,液膜传热系数 h f可用下式计算(2) (2) 蒸汽由蒸发段流向冷凝段的轴向流动热阻,这热阻与其他各项热阻
6、比较,数量级小得多3,因此可以忽略不计。冷凝段凝结换热热阴,可根据 Nusselt 冷凝理论求出平均传热系数(3)JTRv4 314312/efef RgcihR15314312/ecec Rgh8其中(4)冷凝段管壁导热 热阻冷凝段管外换热热阻上述七个传热环节构成了从高温热源到低温冷源的全部传热过程,这七个传热环节是一个完整的串联热路如图 2 所示,由此可得总热阻的计算式(5)若忽略蒸汽压差流动过程的热阻 R4,则(5)式可以写成R=R1+R2+R3+ R5+ R6+R7 (6)其中 R1与 R7是管外热阻,R 2,R 3,R 5与 R6是管内热阻图 2 热阻和热路2.实验数学关系式的建立今
7、以蒸发段管外总表面积为基准,则总传热率的方程式为(7)这个方程也是总传热系数 U 的定义方程。总热阻与各热阻之间的关系可写成/.4LcqRecideng006chR17127iRimhtUAQ chchcihhhfh ARARAUR 11653200 9(8) 根据式(8) ,等式左边总传热系数 U 可以根据实验测得的 Q、tm 和已知的 Ah计算出来。等式右边第二项中的 R2与 R6, 由于碳钢的导热系数随温度变化甚微,所以当热管尺寸一定时,R 2与 R6可以认为是常数。等式右边的第三项与第四项中,R 3可以用式(1)与式(2)进行计算,R 5可以用式(3)与式(4)进行计算,单管试验表明用
8、式(1)一式(4)计算误差较小。等式右边的第六项中 hc可以采用大容积沸腾换热系数的经验公式(4)进行计算,当饱和蒸汽压在 1-40 公斤/厘米 2时其中 p 与 q 均可由实验测得。等式右边第一项中的 hh是我们所要寻求的量,从传热学知道,气流横掠管束时的换热准则方程为(10)(10)式亦可写成式 (11) 将(11)代入(8)式,经整理后可得其 中(12)(13)15.07.3hc31rmehupARN3100.rmehhrhpARd chcihhichinrfmeh ARAAdlCpdCUAR 1215303101 10(12)两边取对数( 14)如果设则(14)式可写 成(15)式(1
9、5)即为实验数学关系式,y 与 x 成线性关系,通过实验测得四十四对 x,y,采用线性回归的方法,就可以求得 A 与 m,从而求得高频焊接螺旋翅片管的管外换热的准则方程式。3.实验装置与方法试验用的热管,光管外径为32.5 毫米,翅片高为 14.75 毫米,翅片间距为 6。33 毫米,蒸发段长度为 500 毫米,冷凝段长度为 30毫米,冷凝段为光管,热管作成无吸液芯的碳钢-水热虹吸管。试验用热管换热器由 14 根热管组成,按等边三角形排列,间距sL与 sT 为 65 毫米,共排成 9 排。由燃油炉产生高温烟气经配风口配入适当的空气后横掠热管管束的加热段,再经引风机排空。冷水用注水泵连续地向热管
10、管束的冷凝段供水,吸热后变成蒸汽经汽水分离及节流膨胀后进入板式冷凝器冷凝成水。 3101rfuehun pdCUlRmlAl baypdCUlRxmbAlarfnehu310111烟气进、出口温度由镍铬-镍硅热电偶与数字电压表示测定。烟气流量是用笛形管和倾斜式微压计测量。蒸汽的温度和压力用水银温度计和压力表测量,通过节流膨胀测定的干度,并计量蒸发率。从公式(11)中可以看出,管外换热系数 hh是几何参数、气流速度和物性参数的函数。在试验过程中几何参数是不变的,而物性参数 、p 和 在烟气成份变化很小的情况下主要取决于烟气的温度。所以可以认为(16)根据以上分析,试验时维持某一个烟气流速和烟气平
11、均温度作为一个试验工况。实验时烟气温度分别选用240C、280 C、320 C、350C、和 380C,烟气流速是在Reh=6000-14000 范围内选取 6-9 个测点。4.试验结果通过 44 个试验工况下测得的数据经过计算得出 y-x 的关系曲线fht,12如图 4 所示。从图中看出有较好的线性关系。用线性回归的方法得出的气流横掠高频焊接螺旋翅片管的管外换热 系数的准则方程为 (17)本准则方程的定性尺寸为光管外径 d o,定性温度为烟气平均温度 t f。适用范围 tf=240380C、Re n=600014000。三阻力特性的研究设计热管换热器时,气流通过换热器的压力损失将导致能量损失
12、,直接影响换热器的经济效益。根据相似原理管内流动的准则方程式(6)为(18)或 (19) 由在水平管内流动压差 与管长成正比,因此(20)最后可以写成管内流动沿程阻力的统一计算式为(21)对于气流通过换热器内翅片管束时,压力损失可以认为与气流密度成正比与 容积当量直径成反比。S、L、Jameson(7)提出的经验公式为3168.017.rehuRN211dlRpEeeu2dlRFpe2dlfp13(22)并推 荐(23)(23)式是特定情况下得出的,对于高频焊执着螺旋翅片管管束叉排的情况下,f 与雷诺数 Rep的关系还需用实验确定。而在北紊流区与 Rep 的关系一般具有以下形式为了用实验方法来
13、找出系数 B 与指数 n,将(24)式代入(22) ,经整理后得(25)两边取对数后得(26)今设 则(26)式成为y=a+bx (27)同样用所选的 44 个试验工况下得出实验数据计算出 44 组 x1,y。用线性回归的方法得出阻力系数 f 与雷诺数 Rep 的关系式(28)度验所得的 r 与 Rep 的关系表示在图 5 上。其中雷诺数 Rep用下式计算(29) 25.0 6.04.014.02Re38pf SLTDevvgcLGfwmnpfRe6.04.04.02Re SLTDevLgevGcpBwmn 6.04.04.02e SLTevLgevclpll wmnn 6.04.04.02R
14、e SLTDevLgevGclypxblawmnn 206.Re37.1pf DevGmRe14特 尺寸为 Dev,定性温度为气流温度 tf,适用范围为 Rep=2400-6000, t f=240C-380C四讨论1.实验精度分析气流横掠翅片管管束管外换热的准则方程式 Nu=0.1370Re+ 0.6338prt的相关系数 r=0.912,说明线性关系较好,剩余标准差 S=0.0364。说明按3s 取分布范围,置信区间为 y=a+bx3s 即 99.7%数值落在此区间内。气流横掠翅片管管束的阻力系数关系式 =1.307Rep-0.2206 它的相关系数 r=0.90,说明线性程度仍较好,剩余
15、标准差 s=0.0139,取3s 为分布范围则置信区间为 y1=a1+b1x13s 说明有 99.7%的数值落在此区间内。2.应用范围本文得出的两个关系式,上面已指出了气流温度范围与气流 Re的范围,除此之外,还受几何条件的限制,由于我组在设计每小时回收 10 万到 60 万千卡热量的热管换热器,普遍采用上述几何尺寸的热管与等边三角形的叉排形式,所以度验结果在这种几何条件下是完全适用的。本文还没有作管径、顺排、管间距、翅片高与翅片间距等因素的影响。所以对于与本实验几何条件不同的热管换热器,原则上不能应用本文的结果,但是我们从实际换热器的测试中发现,对于几何条件不同(顺排除外)的热管换热器本文所
16、得的关系式也具有良好的参考价值。3.与实测结果对比用本实验所得的关系式分别设计了两台热管换热器,一台是江苏省如东县化肥厂的热管余热锅炉,一台是常州市第四印染厂的热15管空气预热器,其设计值与实测值对比如下表常州市第四例子染厂热管空气预热器如东化肥厂热管余热锅炉设计值 QUAP(热)AP(冷)24825 千卡/时15.20 千卡/米 2时C17.9mmH2O12.5mmH2O233100 千卡/时(按排烟温度 220C30 千卡/米 2时C33mmH2O实测值 QUAP(热)AP(冷)24082 千卡/时14.69 千卡/米 2时C15.1mmH2O8.3mmH2O163700 千卡/时(在烟温
17、为 208C31.82 千卡/米2时C40 mmH2O以上两台热管换热器分别于 1982 年 7 月与 9 月通过了技术鉴定。4.与其他文献对比在资料(8)中对翅片管管外换热系数推荐的准则方程是,而本实验所得的准则方程是,系数略低与前者,指数略高于前者。在资料(9)中推荐了翅片管束阻力系数的计算式为 ,而资料(7)推荐的计算式为 而本实验所得的计算式为 与前者较为接近。而资料的对比曲线分别表示在图 6 与图 7 上。这种又罗接近又有差异的情况,本文认为是于几何条件,试验范围,试验方法与测度精度等因素的差异所造成的。5.实验方法比较25.0Re3.1f 25.0Re38.f206.Re3.1f1
18、6气流横掠管束的换热系数的确定,由于准确测量壁温比较困难,一般均采用从总传热系数中将管外换热系数分离出来的方法。目前采用的有威尔逊法,修正威尔逊法,等雷诺数法及逐次逼近法。最近还有用非稳态法测定管外换热系数。以上这些方法有的在试验时控制因素较多,造成实验条件苛刻。有的数学处理比较繁锁,非稳态一般要采用微处理机进行控制和数据处理。本文所采用的方法,实验程序比较简单,实验条件容易实现,数据处理也不复杂。6.结论本文所得出的气流横掠高频焊接螺旋翅片管管外换热的准则方程在气流温度从 240。 -380。范围内,Re 在 6000-14000 范围内,几何条件与试验用管束相似的情况下可以适用,以满足工业
19、用热管换热器设计的精度要求。对于一般的叉排的热管换热器的设计具有参考价值。同样所得出的阻力系数的方程式 f=1.307Re-0.2206在气流温度自 240-380 和雷诺数 Re 在 2400-6000 范围内均可适用。1M.shiraishi “Investigation of Heat Transfer Characteristics of a Two-phase Closed Thermosyphon” 4 th International Heat Pipe Conference 1981.2.S.W.Chi “Heat Pipe Theory and Practice”McGraw
20、-Hill 19763.P.D.D.unn, D.A.Reay “Heat pipes”4.同济大学“传热学”建筑工业出版社 19805.J.P.Holman “Heat Transfer” McGraw-Hill Book Company 3168.0PrRe137.Na1719766.清华大学热能工程系“流体力学机械工业出版社 19807.S.L.Janeson,ASME trans 67.1945 p633-6418.Robert.H.Derry”Chemical Engineers Handbook”(Fifth Edition)19739.A.Rgunter,ASME trans 67.1945 p633-660