扭曲椭圆管换热器壳程强化传热的数值研究
陈杭生1,2,3, 朱冬生1,2,, 陈二雄1,2, 莫逊1,2, 谭连元4
1.中国科学院广州能源研究所,广州 510640
2.中国科学院可再生能源重点实验室,广州 510640
3.中国科学院大学,北京 100049
4.佛山市步鹿节能科技有限公司,广东 佛山 528000
† 通信作者:朱冬生,E-mail:zhuds@ms.giec.ac.cn

作者简介: 陈杭生(1990-),男,硕士研究生,主要从事余热回收利用、强化传热及节能环保研究。朱冬生(1964-),男,博士,教授,博士生导师,主要从事余热回收利用、强化传热及节能环保研究。

摘要

基于扭曲椭圆管的换热器是一种新型的新风系统换热器,针对扭曲椭圆管及其应用特点,设计了两种不同结构参数的新风系统换热器。应用FLUENT软件,在夏季工况下对两种不同结构参数的新风系统换热器壳程进行模拟分析,并通过与实验数据的对比,验证计算模型的可靠性。结果显示在相同体积流量下,随着壳程开孔面积的增大,对流换热系数 h不断减小,压降Δ p不断减小,综合性能系数 h p1/3变化不明显;随着螺距的减小,对流换热系数 h不断增大,压降Δ p不断增大,综合性能系数 h p1/3也不断增大;流场分析显示,扭曲椭圆管换热器壳程流道内,呈现出明显沿着扭曲椭圆管壁面的螺旋流,使得空气在流道内充分扰动,增强换热效果。

关键词: 扭曲椭圆管; 壳程; 强化传热; 数值模拟
中图分类号:TK124 文献标识码:A 文章编号:2095-560X(2022)06-0519-07
Numerical Study on Shell-Side Heat Transfer Enhancement of Twisted Oval Tube Heat Exchanger
CHEN Hang-sheng1,2,3, ZHU Dong-sheng1,2,, CHEN Er-xiong1,2, MO Xun1,2, TAN Lian-yuan4
1. Guangzhou Institute of Energy Conversion, Chinese Academy of Sciences, Guangzhou 510640, China
2. CAS Key Laboratory of Renewable Energy, Guangzhou 510640, China
3. University of Chinese Academy of Sciences, Beijing 100049, China
4. Foshan Bulu Energy Saving Technology Co., Ltd., Foshan 528000, Guangdong, China
Abstract

The heat exchanger based on a twisted oval tube is a new type of fresh air system heat exchanger. According to the characteristics of twisted oval tubes and their application, two kinds of fresh air system heat exchangers with different structural parameters were designed. The shell-side of two kinds of fresh air system heat exchangers with different structural parameters was simulated and analyzed under summer working conditions by FLUENT software, and the reliability of the calculation model was verified by comparing it with the experimental data. The results showed that under the same volume flow rate, with the increase of the shell-side opening area, the convective heat transfer coefficient h decreased, the pressure drop Δ p decreased, while the change of comprehensive performance coefficient h p1/3 was not obvious. With the decrease of pitch, the convective heat transfer coefficient h increased, and the pressure drop Δ p increased, but h p1/3 also increased. Flow field analysis showed that the spiral flow along the wall surface of the twisted oval tube was obvious in the shell-side flow channel of the twisted oval tube heat exchanger, the air in the flow channel was fully disturbed, and the heat transfer effect was enhanced.

Key words: twisted oval tube; shell-side; heat transfer enhancement; numerical simulation
0 引言

换热器是工程技术中广泛采用的冷热流体交换热量的设备。换热器可以按不同的方式分类, 按换热器操作过程可分为间壁式、混合式和蓄热式三大类。其中, 管壳式换热器是间壁式换热器的一种主要形式[1]。管壳式换热器是诞生至今应用最为广泛的一种换热器, 从发展状况看, 针对管壳式的结构改进包括管程和壳程两个部分。而对于管程结构的改进主要集中在管内插入物与传热管的设计优化[2]

扭曲椭圆管是一种高效强化传热的换热管, 经钢带捆扎后的扭曲椭圆管管束在最大变径凸点处相互接触, 可形成自支撑结构, 减少振动[3, 4]。刘世杰等[5]利用FLUENT软件, 研究了扭曲椭圆管的不同结构参数对管内传热与压降性能的影响, 结果显示, 由于受到离心力的作用, 流体在扭曲椭圆管内产生了垂直于主流方向上的二次流, 促进了边界层与主流区的径向混合, 强化了传热。郝洪亮等[6]对采用扭曲椭圆管和H型翅片管的换热器进行了试验测试, 结果表明, 扭曲椭圆管换热器壳侧传热系数明显优于H型翅片管换热器, 且空气流速越大, 优势越明显, 但同时压降也会增大。TAN等[7]对扭曲椭圆管换热器整体性能的研究表明, 扭曲椭圆管换热器在管侧低流量和壳侧高流量情况下均具有较好的换热性能。LI等[8]采用扭曲椭圆管作内管, 可使材料的综合传热性能提高24.0% ~ 39.0%。内扭曲椭圆管的左右扭转方向对壳程传热的强化有相同的影响。DZYUBENKO[9]采用数值计算的方法对扭曲椭圆管束间流动和换热情况进行了研究, 得出了三扭曲椭圆管管束中速度场和温度场的计算方法。陈文静[10]通过数值模拟的方法研究扭曲管扭转方向对管壳式换热器壳侧换热性能的影响, 存在最佳的短长轴比和扭矩能够使得换热器的换热性能更优。

基于扭曲椭圆管的管壳式换热器是一种紧凑式换热器, 由其组成带热回收功能的新风系统, 国内外尚未有报道。本文应用FLUENT软件对换热器壳程强化传热特性展开数值研究, 模拟分析不同螺距、不同壳程开孔面积对换热器壳程对流换热系数、压降的影响, 并分析其温度流场和速度流场, 研究结果可为扭曲椭圆管管壳式换热器的工程应用提供一定的参考。

1 换热器计算模型与边界条件
1.1 物理模型

图1为扭曲椭圆管结构示意图, 扭曲椭圆管由圆管经过压扁扭转而成。图2为扭曲椭圆管换热器结构示意图, 为了简化模型, 本文物理模型由9根扭曲椭圆管组成, 分为3排, 每排3根。图3为扭曲椭圆管换热器截面示意图。

图1 扭曲椭圆管结构示意图Fig. 1 Twisted oval tube structure diagram

图2 换热器结构示意图Fig. 2 Heat exchanger structure diagram

图3 简化后模型截面示意图Fig. 3 Structure diagram of the simplified model section

利用三维专业软件Inventor进行本文所研究的物理建模, 采用圆管的直径D为10 mm, 扭曲椭圆管长轴为12 mm, 短轴为7.77 mm, 总长度为700 mm, 螺距P分别为80 mm、100 mm、120 mm, 壳程开孔宽度d为37.8 mm不变, 壳程开孔长度L分别为110 mm、130 mm、150 mm时, 开孔面积(Ax = L × d)分别为4 158 mm2、4 914 mm2、5 670 mm2, 表1为模型管型尺寸及壳程开孔面积。

表1 模型管型尺寸及壳程开孔面积 Table 1 Tube size and shell-side opening area of the model

采用FLUENT软件进行数值分析, 根据本研究模拟的流动和传热特性, 及各种湍流模型的适用性, 压力和速度耦合采用SIMPLE算法, 湍流模型采用基于扩散系数和湍流动能的Standardk-ε 模型。

为方便模拟及相关数值计算, 进行以下假设[11]:①管程和壳程的风量相等; ②换热器的传热过程处于稳态; ③相关热物性参数在模拟过程中恒定不变; ④流体为不可压缩的理想流体; ⑤热交换界面密封严密, 不存在串风及漏风现象。

1.2 边界条件

研究对象边界条件设置为:壳程介质为空气, 入口设置为体积流量入口, 恒温T = 308.15 K, 出口设置为自由压力出口; 管程介质为空气, 入口设置为体积流量入口, 恒温T = 300.15 K, 出口设置为自由压力出口; 壳体设置为绝热边界。

1.3 网格无关性验证

模型采用ANSYS 18.0中Fluent模块自带的Mesh功能进行处理, 为保证近壁面处计算精度以及气液界面的准确捕捉, 对近壁面液膜流动区域进行加密, 以保证计算结果的精度, 图4为扭曲椭圆管换热器结构化网格示意图。

图4 结构化网格示意图Fig. 4 Structured grid diagram

新风系统温度效率是一个很重要的参数, 其表征了新风系统换热效果的好坏。

温度交换效率计算公式:

${{\eta }_{\text{wd}}}\text{=}\frac{{{T}_{\text{xj}}}-{{T}_{\text{xc}}}}{{{T}_{\text{xj}}}-{{T}_{\text{pj}}}}\times 100\text{ }\!\!%\!\!\text{ }$ (1)

式中:${{\eta }_{\text{wd}}}$为温度交换效率, %; Txj为新风进风干球温度, K; Txc为新风出风干球温度, K; Tpj为排风进风干球温度, K。

当网格数量从6.5 × 106增加到7.5 × 106时, ${{\eta }_{\text{wd}}}$从55.21%变为54.96%, 仅减小了0.25%, 为了兼顾计算精度和效率, 将该模型的网格数量确定为6.5 × 106

图5 温度交换效率${{\eta }_{\text{wd}}}$随网格数量的变化Fig. 5 Temperature exchange efficiency${{\eta }_{\text{wd}}}$changes with the number of grid

2 模拟结果准确性验证

为了验证本文计算模型的准确性, 搭建了新风系统测试平台, 对尺寸结构长轴A = 12 mm、短轴B = 7.774 mm、螺距P = 120 mm的新风系统, 新风量与回风量相同, 选取两室法来进行${{\eta }_{\text{wd}}}$的测试[12], 图6为实验装置示意图, 图7为实验测试平台。在夏季标准工况下新风系统稳定运行30 min后, 采用热电阻温度采集模块进行连续30 min测试, 分别测试新风出口温度Txc、排风进口温度Tpj、新风进口温度Txj, 每隔1 min测试1组数据, 取平均值进行计算${{\eta }_{\text{wd}}}$。实验结果与模拟结果的对比如图8所示, 最大误差为8.5%, 尚在工程应用允许的误差范围内。分析误差产生的主要原因是温度探头设置位置在整机的进出口, 与新风系统换热器内部的进出口位置有一段距离, 新风在到达换热器之前会有部分的换热, 导致实验值比模拟值效果更好, 造成温度误差。

图6 实验装置示意图Fig. 6 Diagram of the experimental setup

图7 实验测试平台Fig. 7 The experimental test platform

图8 温度交换效率${{\eta }_{\text{wd}}}$实验值与模拟值对比Fig. 8 Comparison of experimental and simulated temperature exchange efficiency ${{\eta }_{\text{wd}}}$

3 结果与分析
3.1 壳程开孔面积Ax对传热与阻力的影响

图9所示为新风系统换热器壳程换热系数h随壳程开孔面积变化的性能曲线, 从图中可以看出, 随着体积流量的增大, 不同壳程开孔面积的h均增大, 且以近似线性的规律变化。在相同体积流量下, 随着壳程开孔面积的增大, h不断减小。随着体积流量的增加, h也随着增大, 近似于线性增长。在体积流量为18 m3/h时, 壳程开孔面积4 158 mm2h比壳程开孔面积5 670 mm2h大5.19%。

图9 壳程开孔面积Ax对换热系数h的影响Fig. 9 Effect of shell-side opening area Ax on heat transfer coefficient h

图10显示, 壳程压降随着体积流量的增大而不断增大, 且其增长的幅度逐渐增大。由于雷诺数的增大造成入口速度的提高, 此外壳程平均流速增加, 引起压降增大。由图也可以看出, 在相同体积流量下, 壳程开孔面积越小, 阻力越大。在体积流量为27 m3/h时, 壳程开孔面积为4 158 mm2的压降比壳程开孔面积5 670 mm2的压降大18.73%。

图10 壳程开孔面积Ax对压降的影响Fig. 10 Effect of shell-side opening area Ax on pressure drop

在实际工程应用中, 需要同时兼顾压降和换热系数的大小, 本文采用壳程对流换热系数h与压降Δ p1/3的比值作为新风系统换热器的综合性能系数。图11为三种不同壳程开孔面积的新风系统换热器综合性能系数变化曲线图。在同一壳程开孔面积下, 综合性能系数随着体积流量的增大而增大, 但壳程开孔面积为4 914 mm2和5 670 mm2时, 综合性能系数相差不大。

图11 综合性能系数与体积流量之间的关系Fig. 11 Relationship between comprehensive performance coefficient and volumetric flow

前文提到, 温度交换效率表征了新风系统换热效果的好坏, 不同壳程开孔面积下, 换热器温度交换效率与体积流量的关系如图12所示。温度交换效率随着体积流量的增大而减小。在相同的体积流量下, 壳程开孔面积增大, 换热器的温度交换效率呈现减小趋势。

图12 温度交换效率${{\eta }_{\text{wd}}}$与体积流量之间的关系Fig. 12 Relationship between shell-side temperature exchange efficiency ${{\eta }_{\text{wd}}}$ and volumetric flow

3.2 螺距P对传热与阻力的影响

图13、图14展示了新风系统换热器扭曲椭圆管Ph和Δ p的影响。由图可知, 壳程h和Δ p均随着P的增大而减小。显然P的增大不利于壳程h的提高, 减小P能对提高h起到积极的作用。从图中还可以看出, Δ p随着体积流量的增加而增大。在相等体积流量前提下, P数值越小时, 单位长度的扭曲椭圆管螺转圈数越多, Δ p受到体积流量的影响就越大。

图13 P对换热系数的影响Fig. 13 Effect of pitch Pon heat transfer coefficient

图14 P对压降的影响Fig. 14 Effect of pitch P on pressure drop

图15为三种不同P的新风系统换热器综合性能系数变化曲线。从图中可以发现, 在相同体积流量下, P= 80 mm时, 综合性能系数最大。在体积流量为18 m3/h时, 当P从120 mm变化到80 mm时, 综合性能系数增加4.4%。

图15 综合性能系数与体积流量之间的关系Fig. 15 Relationship between comprehensive performance coefficient and volumetric flow

不同P下, 换热器温度交换效率与体积流量的关系如图16所示。可以得出, 不同的P下, 新风系统的换热器温度交换效率随体积流量的增大而减小, 并且下降的趋势一致。在相同的体积流量下, P增大, 换热器的温度交换效率呈现减小趋势。

图16 温度交换效率${{\eta }_{\text{wd}}}$与体积流量之间的关系Fig. 16 Relationship between shell-side temperature exchange efficiency ${{\eta }_{\text{wd}}}$ and volumetric flow

4 流场显示与分析

图17展现了扭曲椭圆管换热器和普通圆管换热器壳程流道内的流线分布, 截面位置依次为z = 150 mm、z = 230 mm、z = 310 mm、z = 390 mm、z = 470 mm、z = 550 mm。

图17 换热器壳程流线分布特点Fig. 17 Streamline distribution on shell-side of heat transfer

在普通圆管换热器壳程流道内, 空气均为纵向的平行流, 流动状态相对稳定, 没有垂直于主流方向上的扰流和流体间相互干扰的现象。而在扭曲椭圆管换热器壳程流道内, 空气呈现出了明显的沿着扭曲椭圆管壁面的螺旋流, 这种螺旋流的存在使得空气在扭曲椭圆管换热器壳程流道内的流动状态与普通圆管换热器壳程流道内的流动状态大为不同, 温度分布也因此而改变。螺旋流不仅造成了空气在垂直于主流方向上剧烈的混合扰动, 而且形成了对扭曲椭圆管壁面的冲刷作用, 扰动了热边界层, 增强了换热效果。

图18是在同一风量18 m3/h下, 不同壳程开孔面积的换热器出口处(截面z = 150 mm)速度云图对比情况。扭曲椭圆管换热器的横截面有明显的速度边界层, 在截面大的区域流速大, 呈现出多个“ 椭圆” 形状分布, 符合流体向阻力小的方向流动。随着壳程开孔面积的增大, 出口处的流体高速度区域范围变得越小。

图18 不同壳程开孔面积Ax在出口处的速度云图Fig. 18 Velocity cloud of different shell-side opening area Axat the outlet

图19是在同一风量18 m3/h下, 不同P的换热器在出口处(截面z = 150 mm)温度云图对比情况。由于P的不同, 同一截面处, 扭曲椭圆管的方向也有所不同, 但在长轴凸点附近的速度较小, 对流换热系数小, 因此相应的温度边界层一致显示为厚度大。当P减小时, 单位长度的扭曲椭圆管扭转圈数越多, 壳侧的流体和主流方向的流体也会经过多次的翻转造成充分混合, 增强换热效果。

图19 不同螺距P在出口处的温度云图Fig. 19 Temperature cloud of different pitch P at the outlet

5 结论

通过FLUENT软件研究扭曲椭圆管新风系统换热器壳侧的强化传热特性, 分析对流换热系数、压降在不同壳程开孔面积、不同螺距下受到不同的影响, 得出以下结论:

(1)在本文研究范围内, 相同体积流量下, 随着壳程开孔面积的增大, 对流换热系数h不断减小, 压降Δ p不断减小, 但壳程开孔面积从4 914 mm2变化到5 670 mm2时, 综合性能系数h/Δ p1/3变化不明显。

(2)相同体积流量下, 随着螺距的减小, 对流换热系数h不断增大, 压降Δ p不断增大。螺距从120 mm变化到80 mm时, 对流换热系数h从36.9增加到40.1, 综合性能系数h/Δ p1/3增加4.4%。

(3)由扭曲椭圆管换热器组成的新风系统的温度交换效率随着壳程开孔面积的增大而减小, 随着螺距的减小而增大。

(4)流场云图显示, 由于扭曲椭圆管结构的特殊性, 壳程流道内沿着扭曲椭圆管壁面的螺旋流对传热有明显的强化作用, 长轴凸点的存在对流场的分布有所影响, 因此扭曲椭圆管的长短轴对其强化传热的影响可作为下一步的研究方向。

参考文献
[1] 杨世铭, 陶文铨. 传热学[M]. 4版. 北京: 高等教育出版社, 2006. [本文引用:1]
[2] 杜文静, 赵浚哲, 张立新, . 换热器结构发展综述及展望[J]. 山东大学学报(工学版), 2021, 51(5): 76-83. DOI: 10.6040/j.issn.1672-3961.0.2020.423. [本文引用:1]
[3] TAN X H, ZHU D S, ZHOU G Y, et al. Experimental and numerical study of convective heat transfer and fluid flow in twisted oval tubes[J]. International journal of heat and mass transfer, 2012, 55(17/18): 4701-4710. DOI: 10.1016/j.ijheatmasstransfer.2012.04.030. [本文引用:1]
[4] 莫逊, 朱冬生, 张洁娜. 扭曲椭圆管在MVR系统降膜蒸发器上的应用研究[J]. 化学工程, 2016, 44(9): 24-28. DOI: 10.3969/j.issn.1005-9954.2016.09.005. [本文引用:1]
[5] 刘世杰, 朱冬生, 张洁娜, . 三维变形管MVR蒸发器的强化传热数值研究[J]. 高校化学工程学报, 2018, 32(6): 1307-1313. DOI: 10.3969/j.issn.1003-9015.2018.06.009. [本文引用:1]
[6] 郝洪亮, 林城迪, 朱鸿飞, . 三维变形管用于燃气轮机进气加热器上的传热特性研究[J]. 节能与环保, 2019(9): 44-45. DOI: 10.3969/j.issn.1009-539X.2019.09.019. [本文引用:1]
[7] TAN X H, ZHU D S, ZHOU G Y, et al. Heat transfer and pressure drop performance of twisted oval tube heat exchanger[J]. Applied thermal engineering, 2013, 50(1): 374-383. DOI: 10.1016/j.applthermaleng.2012.06.037. [本文引用:1]
[8] LI X Z, WANG L, FENG R, et al. Study on shell side heat transport enhancement of double tube heat exchangers by twisted oval tubes[J]. International communications in heat and mass transfer, 2021, 124: 105273. DOI: 10.1016/j.icheatmasstransfer.2021.105273. [本文引用:1]
[9] DZYUBENKO B V. Influence of flow twisting on convective heat transfer in banks of twisted tubes[J]. Heat transfer research, 2005, 36(6): 449-460. DOI: 10.1615/HeatTransRes.v36.i6.20. [本文引用:1]
[10] 陈文静. 扭曲椭圆管管壳换热器壳程强化换热研究[D]. 兰州: 兰州交通大学, 2020. DOI: 10.27205/d.cnki.gltec.2020.001007. [本文引用:1]
[11] 张涛, 朱晓军, 彭飞, . 近壁面处理对湍流数值计算的影响分析[J]. 海军工程大学学报, 2013, 25(6): 104-108. DOI: 10.7495/j.issn.1009-3486.2013.06.019. [本文引用:1]
[12] 周玥. 两室式新风换气机实验平台建设及换气机初步测试分析[D]. 哈尔滨: 哈尔滨工业大学, 2008. [本文引用:1]