1. 引言
热交换器工业企业中占有很重要的地位 [1]。由于管壳式冷却器结构简单、应用普遍、生产成本低、适用性强、处理量大,因而多应用于热交换换热设备中 [2]。而对于管壳式冷却器壳侧两相传热特性及压降的研究,可以完善冷却器传热过程分析以及流阻压降特性,对于提高冷却器换热效率以及在节约能源、提供换热效率方面具有重要意义 [3] [4]。
近年国内外许多学者对冷却器两相流进行了研究,黄国权等 [5] 针对壳侧流体介质为高温饱和蒸汽情况下,对不同介质入口速度和温度情况下冷却器壳侧传热系数、湍流动能和压降变化的模拟结果进行了比较,从而得出流体介质入口速度和温度对冷却器壳侧传热特点和传热性能的影响规律。文献 [6] [7] [8] [9] 建立了管壳式水冷冷凝器的换热模型基于该模型编制了管壳式冷凝器的仿真程序,并对冷凝器的性能进行模拟,对仿真结果进行了分析。Yang [10] 研究了垂直螺旋折流板管壳式冷却器壳体内丙烷的传热特性,得到了垂直螺旋板管壳式冷却器壳侧的流型和传热系数,并与实验数据进行了对比分析。
冷却器在工作过程中某些特定温度下会发生相变,从而导致冷却器的换热特性及压降分析十分复杂 [4] [11]。本文给出了某压缩机出口管壳式冷却器温度及压降理论计算的具体表达式,考虑了冷却器壳程热流体在传热过程中发生的相变影响,应用MATLAB软件编程计算壳程出口温度及压降值并与已有实验结果进行了对比验证,进一步研究了各参数对其结果的影响。
2. 几何模型及其数学描述
2.1. 几何模型
某供抽气机组压缩机出口使用管壳式冷却器,该冷却器外结构如图1所示。冷却器总长为5.840 m,管程入口及出口在冷却器端部圆面分布,壳程入口及出口分布于冷却器上部。冷却器由两个支架支撑固定。冷却器管束部分排列结构示意图如图2所示,具体为横排32根、竖排42根、41根交替排列。其中,管外径为18 mm,长度为5000 mm,数量1328根,管束排列为等边三角形排列形式。为改善换热效果,壳程使用折流板改变介质流向,且同时起到支承管束的作用,该冷却器内部有三个折流板,其形状为单弓形。冷却器内部局部结构三维示意图如图3所示,其相应的几何参数如表1所示。
![](//html.hanspub.org/file/2-2650235x9_hanspub.png?20210914082547100)
Figure 1. Outline structure diagram of the shell-and-tube cooler
图1. 管壳式冷却器外形结构图
![](//html.hanspub.org/file/2-2650235x11_hanspub.png?20210914082547100)
Figure 3. 3D schematic diagram of internal local structure
图3. 内部局部结构三维示意图
![](Images/Table_Tmp.jpg)
Table 1. Structural geometric parameters of the shell-and-tube cooler
表1. 管壳式冷却器结构几何参数
2.2. 数学模型
2.2.1. 传热过程
传热方程可由下述方程给出:
(1)
(2)
其中,QH为换热量,W;k为整个传热面上的平均传热系数,W∙(m2∙℃)−1;A为传热面积,m2;Δtm为两种流体间的平均温差℃。c1和c2分别为热流体、冷流体的定压比热容,kg∙(kg∙℃);qm1和qm2分别为热流体及冷流体的质量流速,kg∙s;
和
为热流体进出口温度,
和
为冷流体进出口温度,单位均为℃。
冷却器性能通过有效度及热容量比测量,有效度εH由。方程给出:
(3)
其中,qH为冷却器的实际换热量,qHmax为冷却器的最大可能换热量,单位为W。热容量比C*由下式给出:
(4)
式中,Cmax为较大换热量;Cmin为较小换热量;Ch为热流体热容量;Cc为冷流体热容量,单位均为W/℃。
若冷却器的总传热面积为F1,总传热系数为K1;壳程进口温度tkin,壳程出口温度tkout,壳程流量Mk,比热cpk,热容量Ck = Mkcpk;管程流体进口温度为tgin,管程出口温度为tgout,流量Mg,比热cpg,热容量Cg = Mgcpg。已知传热面积、传热系数及热容量,可求出传热单元数NTU:
(5)
其中
(6)
壳程出口温度表达式为:
(7)
由于壳程流体存在部分冷凝现象,水平管外冷凝时的表面换热系数为:
(8)
式中,nm为管束上冷凝时传热系数的修正系数;Bm为管壁上冷凝液膜的组合物性参数;rs为潜热,J/kg;Δt0为壳侧流体温度与管外污垢层外表面温度间的对数平均温差,℃;d0为管外径,m;ρ为冷凝液膜的密度,kg/m3;λ为冷凝液膜的导热系数,W/(m∙K)−1;vm为冷凝液膜的运动粘,m2/s。
发生冷凝现象时,通常纯净的冷凝液体到壁面之间的换热量为
(9)
式子中,hcon为冷凝时的换热系数,Tsat为冷凝液体在给定压力时所对应的饱和温度,Tw是壁面温度。
由于管程流体未发生冷凝现象,首先计算得到管程流体出口温度,从而得到管程流体的换热量,通过管程换热量及管外冷凝时的表面换热系数,并对其出口温度计算公式进行修正,从而得出适用于该冷却器的壳程出口温度数学公式:
(10)
2.2.2. 压降分析
管壳式冷却器壳程压降分为三部分,即中心区(交叉流区)压降;缺口压降;壳侧的入口和出口处的压降。中心区(交叉流区)压降可表示为:
(11)
其中,Nr,cc为流过截面的有效管排数;Gc为交叉流区质量流量;fid为单位管道的理想范宁摩擦因子;gc是牛顿第二定律中的比例常数;μw是管壁处液体的粘性,μm是平均温度下液体的粘性;ρs为流体密度。
缺口压降与雷诺数有关:
(12)
其中,pt为管距,Nr,cw为缺口内交叉流过管排的有效管数,Lb为折流板间距,Dh,w为缺口处水力直径。
壳侧的入口和出口区压降可表示为:
(13)
冷却器的壳程压降表示为:
(14)
其中Nb为折流板数;ξb为旁流压降修正因子;ξl为折流板与壳侧及管与折流板侧的泄露修正因子;ξs为进出口折流板截面不均匀流板间距产生的压降修正因子。
2.3. 理论计算具体表达式
基于上述冷却器温度及其压降理论公式,应用所研究冷却器的几何参数及相应的物料参数,可计算得到冷却器壳程出口的温度与压力,并与已有该冷却器的实际工况测量数据进行了对比验证。
在实际工况中,根据冷却器的进出口温度实测数据,发现温度存在较大差异,存在冷凝现象。因此,计算中由入口边界条件得到管程出口温度及总换热量,并通过实际工况中的参数推导出管外冷凝时的表面传热系数α,如下式
(15)
冷却器的壳程出口温度tkout:
(16)
上式中,Pkin为壳程入口压力;tkin为壳程入口温度,tgin为管程入口温度;Ck为壳程热流体比热容,εH为冷却器有效度。Ck可表示为:
(17)
其中Wy为压缩机的功率,Ak为关于壳程入口温度及压力的多项式,表示为:
(18)
εH可表示为:
(19)
其中传热系数k及热容量比C*可分别表示为:
(20)
(21)
Pgin为管程入口压力。
壳程出口压降通过推导计算获得,带入相关参数,化简后的出口压降ΔPs为:
(22)
其中At为壳程入口温度的函数,Ak2为壳程入口温度及压力的函数:
(23)
(24)
3. 计算结果与分析
基于以上管壳式冷却器的数学模型,应用MATLAB软件建立了管壳式冷却器进出口温度及压降计算仿真程序,其流程图如图4所示。
3.1. 计算结果与实验结果的对比分析
给定冷却器结构参数,流体物性参数,其管程及壳程进口温度和压力采用已有实测数据,可计算得到壳程出口温度及压降,并与已有实验结果进行对比分析。图5为实验测得的壳程进口温度随时间变化曲线,利用该温度作为冷却器入口的初始条件,图6给出了考虑冷凝和不考虑冷凝壳程出口温度理论计算结果与已有实测值的对比曲线。可以看出壳程出口温度变化趋势与入口温度变化趋势一致,均随时间增加逐渐增大。考虑冷凝现象的计算结果较不考虑冷凝的结果更接近与实测值。由于冷却器壳程空气流动情况比较复杂,存在热对流等现象,本文计算中仅热传导,两者仍存在微小差异。
图7给出了该冷却器壳程进口压力实测曲线,图8给出了壳程出口压力计算结果与已有实测值的对比曲线,可以看出两者具有很好的一致性。图9给出了壳程压降随时间变化曲线,可以看出压降的变化与压力有关,压降随进口压力的增大而增大。壳程进出口温度及压力的计算值和实测结果为33,848个数据点,表2列出了其部分时间段的具体数值。
![](//html.hanspub.org/file/2-2650235x41_hanspub.png?20210914082547100)
Figure 5. Shell side inlet temperature curve
图5. 壳程进口温度变化曲线
![](//html.hanspub.org/file/2-2650235x42_hanspub.png?20210914082547100)
Figure 6. Shell side outlet temperature curve
图6. 壳程出口温度变化曲线
![](//html.hanspub.org/file/2-2650235x43_hanspub.png?20210914082547100)
Figure 7. Shell side inlet pressure curve
图7. 壳程进口压力变化曲线
![](//html.hanspub.org/file/2-2650235x44_hanspub.png?20210914082547100)
Figure 8. Shell side outlet pressure curve
图8. 壳程出口压力变化曲线
![](//html.hanspub.org/file/2-2650235x45_hanspub.png?20210914082547100)
Figure 9. Shell side pressure drop curve
图9. 壳程压降变化曲线
![](Images/Table_Tmp.jpg)
Table 2. Partial calculation results and measured data of the shell-and-tube cooler
表2. 冷却器部分计算结果与实测数据
3.2. 各参数对冷却器性能的影响分析
3.2.1. 各参数对壳程出口温度的影响
由于冷却器性能参数对壳程出口温度及压降影响较大,为定量分析该影响作用,利用控制变量法分析壳程出口温度及压降随各参数的变化规律。图10给出了冷却器壳程出口温度随管程介质密度的变化曲线,可以看出其壳程出口温度随冷却器管程介质密度的增加而增加,呈非线性变化。图11给出了换热管外径对壳程出口温度的影响曲线,可以看出冷却器换热管的结构参数换热管外径同样影响着冷却器的换热特性和换热系数,冷却器壳程出口温度随换热管的外径增大而减小。
![](//html.hanspub.org/file/2-2650235x46_hanspub.png?20210914082547100)
Figure 10. Effect of tube medium density on shell side outlet temperature
图10. 管程介质密度对壳程出口温度的影响
![](//html.hanspub.org/file/2-2650235x47_hanspub.png?20210914082547100)
Figure 11. Effect of outer diameter of heat exchanger tube on outlet temperature of shell side
图11. 换热管外径对壳程出口温度的影响
3.2.2. 各参数对壳程压降的影响
冷却器中各重要的部件包括换热管、折流板的结构参数以及其布置方式均对冷却器流阻压降有显著的影响。图12和图13分别给出了换热管横向及纵向间距对壳程流阻压降的影响,可以看出随着换热管间距的增大,壳程的流阻压降逐渐减小,横向间距较纵向间距对其影响更为显著。图14和图15分别给出了折流板间距以及其缺口高度对壳程流阻压降的影响,增大折流板间距或者折流板的缺口高度均会显著降低壳程流体出口压降。
![](//html.hanspub.org/file/2-2650235x48_hanspub.png?20210914082547100)
Figure 12. Effect of transverse tube spacing on pressure drop
图12. 横向管间距对压降的影响
![](//html.hanspub.org/file/2-2650235x49_hanspub.png?20210914082547100)
Figure 13. Effect of longitudinal tube spacing on pressure drop
图13. 纵向管间距对压降的影响
![](//html.hanspub.org/file/2-2650235x50_hanspub.png?20210914082547100)
Figure 14. Effect of baffle spacing on pressure drop
图14. 折流板间距对压降的影响
图16给出了冷却器壳程流体流量对壳程流阻压降的影响,随着壳程流体进口流量的增大,壳程压降逐渐增大,呈非线性变化。
4. 结论
本文研究了管壳式冷却器两相流动传热及压降特性,给出了适用于该冷却器壳程出口温度及压降的计算表达式,并考虑了其壳程热流体的相变影响,数值结果与已有实测数据进行了对比验证,进一步分析了结构参数对其壳程出口温度及压降的影响,得到下述结论。
![](//html.hanspub.org/file/2-2650235x51_hanspub.png?20210914082547100)
Figure 15. Influence of baffle notch height on pressure drop
图15. 折流板缺口高度对压降的影响
![](//html.hanspub.org/file/2-2650235x52_hanspub.png?20210914082547100)
Figure 16. Effect of fluid flow on pressure drop
图16. 流体流量对压降的影响
冷却器壳程出口温度随管程介质密度的增加而增加,同时随换热管的外径增大而减小。
换热管和折流板的间距及其缺口高度对壳程的流阻压降影响显著,壳程压降随换热管横向及纵向管间距、折流板间距及缺口高度的增大而减小,随壳程流体流量的增大而增大,呈非线性变化。
由于实际工况中壳程热流体流动及换热情况比较复杂,为考虑换热器与外界的热交换及内部实际两相流动情况,在之后的分析中应考虑到包括热对流的实际换热模型及其它两相流动因素的影响,进一步提出更精确的壳程出口温度及压降理论模型。