螺旋扁管折流杆换热器传热强化数值模拟研究

2020-03-26王斯民孙利娟简冠平张早校

高校化学工程学报 2020年1期

王斯民, 孙利娟, 宋晨, 肖娟, 简冠平, 张早校, 文键

(1.西安交通大学化学工程与技术学院, 陕西西安 710049;2.西安交通大学能源与动力学院, 陕西西安 710049)

1 前言

折流杆换热器(rod-baffle heat exchangers，RBHEs)由美国菲利普石油公司于1970年首次提出并制造，用折流栅代替折流板，有防止管束振动，降低壳程压降，减少流动死区和漏流损失等优点[1-4]。很多学者对传统RBHEs(换热管为圆管且四边形布置，折流杆为圆杆)进行了研究。Philiips公司GENTRY等[5-8]和中国郑州大学热能研究中心、董其伍等[9]和马小明等[10]做了大量实验，得到了传统RBHEs壳程努塞尔数Nu和摩擦系数f的经验关联式。在此基础上，研究集中在改变折流杆结构，换热管结构及其布管方式。实验方面：陈文昕[11]得出了三角形布管RBHEs经验关联式，与现有的正方形布管经验关联式比较，结果表明三角形布管方式更优。严良文等[12-13]研究了波形折流杆、三角形布管的 RBHEs，其总传热系数高于折流板换热器。孔松涛等[14]证明了层流时45°菱形折流杆的RBHEs综合传热系数比传统RBHEs高20%。WANG等[15]发现雷诺数Re在2 000 ～ 4 800时，螺纹管RBHEs的传热系数比光管大38%，压降相当。模拟方面：整体模型为主，多孔介质模型、周期性单元流道模型、周期性全截面模型、和“分段模拟，整体综合”为辅。YANG等[16]对比了模拟与实验结果发现了多孔介质模型和整体模型能够同时准确计算传热和压降。刘佳驹等[17-18]对比发现了深槽螺旋波纹管RBHEs综合性能优于光管RBHEs。WANG等[19]证明波形折流杆，三角形布管的RBHEs壳程的综合换热性能比传统RBHEs大10%。YU等[20]提出了六角形夹持和圆弧夹持式 RBHEs，结果表明六角形夹持结构的刚度大，适用于大、重管束，传热效果好，但整体性能指标低于圆弧夹持结构。YOU等[21]使用周期单元流道模型，对比发现带凹槽折流杆的RBHEs的交错排列形式的传热系数h比非交错排列形式提高了41.9%，压降增大一倍，传热系数与压降比(h⋅ΔP-1)更大。YUAN等[22]证明了RBHEs与折流板换热器的传热相当，但压降大幅度降低。WANG等[23]证明双壳程的h⋅ΔP-1比单壳程RBHEs提高了约8.9%。无论实验还是模拟，对RBHEs进行换热管结构改进的文献很少，借鉴现有的纵向流管壳式换热器传热管形式，基于流动传热强化机理，本文提出用螺旋扁管代替圆管安装在RBHEs中。对比了螺旋扁管RBHEs和圆管RBHEs壳程的流动、传热和综合性能。

2 折流杆换热器的整体模拟

2.1 几何模型及网格划分

按照GB151-2014标准，设计并建立了螺旋扁管RBHEs和圆管RBHEs壳程的几何模型，如图1所示，坐标原点在入口正中心。2种模型换热管尺寸不一样，其余结构参数都相同。圆管的直径是16 mm。螺旋扁管的横截面是椭圆，长轴为16 mm，短轴为8 mm，螺距为266.7 mm。管间距为22 mm；壳体内径是147 mm；壳体和换热管长度均为1 000 mm。折流栅间距为200 mm，双排管间布管，沿着流动方向被定义为第1、2、3和4个折流栅，由折流环和折流杆组成。折流环的横截面为圆形，直径为7.3 mm；折流杆的直径为6 mm。螺旋扁管的长轴等于圆管直径，便于将圆管抽出后，螺旋扁管直接插入安装。

图1 折流杆换热器的几何模型Fig.1 Geometry model of RBHEs

先使用ANSYS Meshing对几何模型进行了非结构化网格划分，导入Fluent中转化为多面体网格，大大降低网格数量并提高网格质量。对圆管RBHEs和螺旋扁管RBHEs分别进行了网格独立性验证，结果如图2所示。圆管RBHEs模型在网格数量达到1 213 899后，换热系数h变化低于2.56%，螺旋扁管RBHEs模型在网格数量高于4 928 456时，换热系数h和变化小于1.84%。圆管和螺旋扁管RBHEs最终取定的网格数量分别为1 737 729和5 383 816。

图2 网格独立性验证Fig.2 Validation of mesh independence

2.2 数值方法及模型验证

壳程工作流体为水，假设为不可压缩流体湍流流动，忽略重力和水物性的变化。入口速度边界条件，为0.1、0.2、0.3、0.4、0.5和0.6 m⋅s-1，水的进口温度为300 K。压力出口边界条件。管内为低温沸点工质冷凝，设定管壁温度为350 K。其余壁面均为无滑移绝热壁面。采用标准k-ε湍流模型和标准的壁面函数模型。压力和速度耦合采用SIMPLE算法，动量和能量方程均采用二阶迎风格式，湍动能和湍流耗散率方程均采用一阶迎风格式。残差值除能量方程需小于1×10-6，其余项小于1×10-3认为收敛。标准kε模型的控制方程如下：

连续性方程：

动量方程：

能量方程：

湍动能k方程：

湍动能耗散ε方程：

将折流杆换热器壳程和管程同时模拟，与文献[24]壳程为冷水，管程为热水，折流栅间距为120 mm的工况所得实验数据进行对比。管程和壳程入口边界条件为质量流量和温度，其余设置与上述相同。模拟和实验所得到换热器总换热系数Kh的偏差如图3所示。Kh的最大偏差为5.63%，平均偏差为 4.06%。偏差的主要原因有模型的简化；除换热管的其余壁面绝热、水物性不变的假设；忽略重力的影响等。

图3 实验与模拟结果对比Fig.3 Comparison of experimental and simulation results

3 结果和讨论

3.1 速度分布

图4和图5分别为入口速度0.1和0.6 m⋅s-1时，折流杆换热器X=0，YZ平面速度云图。取速度范围的最小值和最大值速度云图分析，可以分析低速入口条件和高速入口条件对流场的影响。

图4 入口速度vinlet为0.1 m⋅s-1的速度云图Fig.4 Velocity contour at inlet velocity vinlet of 0.1 m⋅s-1

从图4和5中可以看出，圆管RBHEs中的高速区和低速区相对较多，速度大小差别更大，速度分布更加不均匀。在高速入口条件下，2种模型壳程流速不均匀度的差别缩小。从流经每个折流栅后的流场可以看出，圆管RBHEs中折流栅对流体的作用更大。图中流经第1、3和4折流栅时，由于流通截面积的减小，流速会增大，在圆管RBHEs中流速增大更明显。在流经第2个折流栅时，因为所截平面刚好穿过第2个折流栅中间折流杆的中心线，受到折流杆的阻拦，速度降低并波动，在圆管RBHEs中这种波动更明显。

图5 入口速度vinlet为0.6 m⋅s-1的速度云图Fig.5 Velocity contour at inlet velocity vinlet of 0.6 m⋅s-1

图6和7是入口速度为0.1、0.6 m⋅s-1时，2种折流杆换热器Z= -500、-550和 -600 mm，XY平面速度云图。从Z= -500和-550 mm云图中可以看出，螺旋扁管RBHEs的管间流速更大，速度分布更加均匀。由Z= -600 mm的云图可以看出，圆管RBHEs流经折流栅时的速度更大。从图6和7中比较可以看出，高入口速度条件下，2种RBHEs速度不均匀度的差距在减小。

图6 入口速度vinlet为0.1 m⋅s-1的速度云图Fig.6 Velocity contour at inlet velocity vinlet of 0.1 m⋅s-1

图7 入口速度vinlet为0.6 m⋅s-1的速度云图Fig.7 Velocity contour at inlet velocity vinlet of 0.6 m⋅s-1

综上，螺旋扁管RBHEs的管间速度更大，分布更加均匀，因为螺旋扁管的扭转曲面，周期性扰动流体，改变流体的流速和流向，加强了流体的纵向混合。流经折流栅部分，圆管RBHEs的流速更大，速度波动也更大，有3方面的原因：1是折流栅结构相同，圆管入口流通面积小于螺旋扁管RBHEs，折流栅的存在使得流通面积减小相同，但圆管RBHEs的减小率更大。2是圆管RBHEs换热管与折流杆的垂直相切处的流通缝隙更小。3是圆管RBHEs的圆管对流体的扰动没有螺旋扁管强，流体流经折流栅引起的速度波动更明显。

图8和图9是入口速度为0.1、0.6 m⋅s-1时，2种折流杆换热器Z= -500 mm，XY平面湍动能云图。由图可见，螺旋扁管RBHEs换热管椭圆截面两端湍动能很高；圆管RBHEs的湍动能均匀度更好。螺旋扁管RBHEs湍动能大的区域(绿色和红色部分)更加多而且分散，在换热管贴壁处和管间均有分布。因为螺旋扁管扰动流体，形成了复杂的周期性流动和混合，湍动能增大。而圆管RBHEs在入口速度为0.1 m⋅s-1时，湍动能大的区域仅分布在换热管周围，在入口速度为0.6 m⋅s-1时，在壳体贴壁处也存在，而且贴壁处的湍动能增强。随着入口速度的增大，两种换热器湍动能均匀度的差别在缩小。

图8 入口速度vinlet为0.1 m⋅s-1的湍动能Fig.8 Turbulence kinetic energy contour at inlet velocity vinlet of 0.1 m⋅s-1

图9 入口速度vinlet为0.6 m⋅s-1的湍动能Fig.9 Turbulence kinetic energy contour at inlet velocity vinlet of 0.6 m⋅s-1

3.2 温度分布

图10～13是换热器轴线平面(X=0，YZ平面)和径向平面(Z= -500 mm，XY平面)的温度云图。从图10和11中可以看出，流体在螺旋扁管RBHEs中的温度低于圆管RBHEs，因为在螺旋扁管RBHEs中的流体速度整体高于圆管，另外在圆管折流杆换热器中的RBHEs更大，分布更不均匀，这与流场的不均匀度大有关。

图10 入口速度vinlet为0.1 m⋅s-1的温度云图Fig.10 Temperature contour at inlet velocity vinlet of 0.1 m⋅s-1

图11 入口速度vinlet为0.6 m⋅s-1的温度云图Fig.11 Temperature contour at inlet velocity vinlet of 0.6 m⋅s-1

从图12和13可以看出，圆管RBHEs中管间的温度高于螺旋扁管RBHEs，同样是管间流体速度低所致。圆管换热器管间温度高的现象在0.6 m⋅s-1中的现象被减弱，因为在高速入口条件下，圆管RBHEs的管间速度得到增强，螺旋扁管增强管间速度的作用被减弱。螺旋扁管RBHEs的贴近换热管处的温度明显高于相同入口速度条件下圆管RBHEs，因为螺旋扁管附近的湍动能大于圆管附近，能够增强换热，该处流体温度得到提高。

图12 入口速度vinlet为0.1 m⋅s-1的温度云图Fig.12 Temperature contour at inlet velocity vinlet of 0.1 m⋅s-1

图13 入口速度vinlet为0.6 m⋅s-1的温度云图Fig.13 Temperature contour at inlet velocity vinlet of 0.6 m⋅s-1

3.3 流动、传热和综合性能分析

图14为圆管和螺旋扁管RBHEs壳程换热系数h随入口速度变化曲线图。由图可以看出，在低速入口条件下，螺旋扁管RBHEs的壳程传热系数h高于圆管RBHEs，在高速入口条件下，螺旋扁管RBHEs的h低于圆管。在入口速度0.1 m⋅s-1的条件下，螺旋扁管RBHEs的传热系数h比圆管折流杆换热器大42.29%，在入口速度为0.6 m⋅s-1的条件下，比圆管RBHEs小6.63%。因为螺旋扁管RBHEs的换热管对流体形成扰动强于圆管RBHEs，而折流栅部分对流体的扰动弱于圆管RBHEs，低入口速度条件下，加强换热方面更强，高入口速度条件下，削弱换热方面更强。可见螺旋扁管增强流体扰动来增强换热的机理更适合用于流体速度较低时，在流体速度较高时，流体主流区流量增大，螺旋结构影响变小。这与文献[25]研究的单管螺旋扁管代替圆管增强换热的方法更适用于低Re下流动的结论一致。

图14 传热系数h随入口速度vinlet的变化Fig.14 Profiles of heat transfer coefficient h as a function of inlet velocity vinlet

图15 压降ΔP随入口速度vinlet的变化Fig.15 Profiles of pressure drop ΔP as a function of inlet velocity vinlet

不同入口速度条件下RBHEs壳程压降ΔP的变化图如图15所示，由图可知，在相同入口速度条件下，螺旋扁管RBHEs的压降低于圆管RBHEs，减少了约33.60 %。因为折流杆换热器的压降主要由无折流栅时的压降和单独由折流栅引起的压降两部分组成，在螺旋扁管RBHEs中，由折流栅部分引起的压降远低于圆管RBHEs。因为圆管RBHEs中，流经折流栅，流速变化大，后续波动更明显，所造成的流动损失更大。

图16为螺旋扁管RBHEs和圆管RBHEs综合性能评价指标(performance evaluation criteria，PEC)和两者的比值随不同入口速度的变化。PEC为Nu⋅f-1/3，Nu、f分别为换热器壳程的努塞尔数、摩擦系数，是换热器综合性能的评价指标。从图中可以看出，螺旋扁管 RBHEs的综合性能优于圆管RBHEs。随着速度的增大，螺旋扁管与圆管RBHEs的PEC比值从2.52降低到1.65。在0.1 m⋅s-1的入口速度时，螺旋扁管RBHEs的综合性能比圆管RBHEs提高152%，在入口速度为0.6 m⋅s-1时，提高了65%。因为螺旋扁管RBHEs在不同入口速度条件下，虽然不一定能够增大换热系数，但其压降总是低于圆管RBHEs，综合性能更优。

图16 综合性能评价指标PEC及其比值随入口速度vinlet的变化Fig.16 Profiles of comprehensive performance efficiency criteria PEC and PEC ratio as a function of inlet velocity vinlet