内置弧形槽扭旋元件圆管中湍流传热特性研究-辽宁化工2022年05期

摘  要：设计出一种新型强化换热元件——弧形槽扭旋元件，并在Re=4 000～10 000范围内对三种弧形槽数量的扭旋元件进行了数值模拟，从平均努塞尔数、温度场分布等方面分析了弧形槽的数量对管道传热特性的影响。模拟结果表明，n=3元件弧形槽的数量和分布位置有利于中心流体与管壁附近流体的混合以此强化传热，在三种元件中具有最佳的换热效果，n=2元件次之，n=1元件的结构不利于槽内流体与外界流体的交汇，换热效果较弱;三种元件引起的管道阻力随弧形槽数量的增加而提升，综合传热性能方面，PEC值均大于1，最大值为n=3元件在Re=4 000时，达到1.33，说明使用弧形槽扭旋元件会使管道传热性能得到不同幅度的提升。

关  键  词：数值模拟;扭旋元件;传热特性;PEC值

中图分类号：TQ053.6 文献标识码： A 文章编号： 1004-0935（2022）05-0629-04

能源的转换和利用是当今社会发展的关键性问题，而热能在各行各业中都占有举足轻重的地位，合理高效的换热设备可以极大促进热能的有效利用[1]，换热管作为换热设备中的重要元件，其传热性能直接决定换热设备的高效与否，因此研究出具有高传热效率C0IcyvKwQ6IZaqgHX2aCvQ==的换热管成为该领域内的热点。

目前强化传热技术主要分为无源技术和有源技术两种[2-3]，本文主要采用管插件无源技术进行强化换热，其作用机理为管插件的合理结构使流体在流经元件时产生径向流动的同时引起二次流[4-7]，此种流动现象会使管道中心流体与壁面附近的高温流体进行充分混合，同时壁面处的边界层也遭到破坏，有利于热量的传递与交换。

由于管插件技术经济高效、易于维修等优点，许多学者对其进行了大量研究。Kumar等[8]在雷诺数为1～25 000的范围内研究了SK型静态混合器的长宽比及单元数量对流体流动和混合的影响，同时对比了元件在不同流态下的混合程度，并通过实验得出了每混合器单位的压降与所用SK元件的数量无关的结论。P. Murugesan等[9]研究了不同扭曲比、深度和宽度比的“V”型双绞扭带对圆管内的传热和摩擦系数的影响，研究发现管内平均努塞尔数和平均摩擦系数随扭曲比、宽度比和深度比的减小而增加。詹三江等[10]在光滑扭带的基础之上研究并分析了开三角形边槽的扭带新结构，实验证明了这种新型结构可以增加换热效果和沿程阻力，并与光滑纽带进行对比，定量分析了阻力系数λ与Nu增加的幅度，并得到了相应的努塞尔数Nu及阻力系数λ的实验关联式。孙达

等[11]对高黏度下中国结式静态混合器的流动特性进行了模拟，发现Re较低时，中国结式静态混合器可以改变流场结构，促进径向混合并减薄边界层，强化高粘流体的混合及传热性能。

本文设计出一种由带有弧形凹槽的波纹板扭旋而成的新型管内扰流元件——弧形槽扭旋元件，探讨了弧形槽扭旋元件强化传热的机理，研究了弧形槽数量的变化对管内传热性能的影响。研究结果可为弧形槽扭旋元件在工程中的应用提供理论基础和技术支持。

1  物理模型及数值模拟

1.1  物理模型

本文研究对象为内含4个扭旋元件的换热管道，相邻元件交错垂直布置，元件直径D，厚度δ，长度H，扭转角为180°，弧形槽的直径d，弧度R =π，弧形槽数量n分别为1、2、3，其中n=2元件结构如图1所示。

图2显示了扭旋元件在管道中的布置情况及三维坐标分布，管长为L，管径D，管道入口段和出口段均保留La。元件及管道的结构参数如表1所示。

1.2  数值模拟方法

1.2.1  模型假设

本文在数值模拟过程中做了如下假设：

1）流体为不可压缩的牛顿流体。

2）流体流动时为稳态。

3）忽略流体物性随温度的变化。

4）忽略自然对流及热辐射。

1.2.2  模拟方法

本文选用的工作流体为水，在雷诺数为4 000～ 10 000的范围内进行研究，选用reliable k-e湍流模型，该模型可以更好地模拟旋转流动和二次流，使用标准壁面函数对近壁区域的计算进行了优化，在求解算法的选择上选用SIMPLEC，方程离散采用二阶迎风格式。

1.2.3  控制方程与边界条件

能量方程：   。   （1）

模拟所采用的边界条件如下：

壁面热边界条件为：

其中， =343 K， =293 K

扭旋元件表面采用无滑移速度边界条件。

1.3  网格无关性验证

模型采用ICEM进行网格划分，为准确模拟近壁面处的流动换热状况，对靠近壁面的网格进行了加密处理，分别对三种不同的弧形槽扭旋元件进行了网格无关性检验。Nu与网格数量的关系如图4所示，当n=1，n=2，n=3三种元件的网格数分别达到1 795 731、1 836 308和1 968 024时网格数量对模拟结果的影响较小，综合考虑计算精度和计算机资源，本文选取上述三种网格数量进行模拟分析。

2  结果分析与讨论

2.1  平均努塞尔数与温度场分析

为了分析弧形槽扭旋元件中凹槽数量对管道传热性能的影响，本章对凹槽数量分别为1个、2个、3个的弧形槽扭旋元件做了数值模拟，计算结果如图3所示，可以看到在相同工况下，三种元件的传热性能相对于光管都有不同幅度的提升，n=3元件传热性能最佳，n=2元件次之，两者整体差距不大，随Re的升高其差值有逐渐增大的趋势，传热效果最差的是n=1元件，较n=2和n=3元件分别低了12.2%～15.6%和16.5%～20.2%，其主要原因是凹槽数量的变化对流场和温度场产生不同的扰动作用。

三种不同数量弧形槽元件的截面温度云图如图4所示，可以看到，三种元件都有不同的低温流体聚集区，其中n=1元件流体温度最低点主要位于弧形槽内，此外，弧形槽凸侧临近壁面的区域也有两处面积较小的低温区，从该元件弧形槽的分布位置就可以解释此种现象，其只有单个弧形槽且位于管道中心线处，当流体进入管道后在该元件的作用下产生径向旋转流动，但中心处流经单弧形槽的流体合速度最低，因此其产生的二次流作用较弱且不易与管壁处的流体交汇混合，这种现象导致了弧形槽中流体温度明显低于其他位置。

n=3元件虽有着上下对称的结构，但其温度分布并未像n=2元件一样呈现对称的现象，其原因是流体流经中间槽时会受其影响产生分离流动，导致整个流场对称性失衡，而n=2元件中间为平板结构，并不会对流体的流向产生影响。此外，n=3元件于其他两种元件不同的是，其上下两侧弧形槽中并没有出现流体低温区，上侧弧形槽中反而温度高于周边流体，这种现象是由流体的旋向及与管道壁面的距离决定的，上侧弧形槽距离管道壁面最近且槽中的流体流向管壁处，这种现象使得弧形槽内的流体与壁面高温区频繁交换混合，因此槽内流体温度较高

2.2  阻力特性分析

三种不同数量弧形槽扭旋元件的阻力系数如图5所示。

可以看到，阻力系数最大的为n=3元件，其所引起的阻力为光管的2.63～2.95倍， n=2元件次之，但两者整体差距不大。弧形槽的存在会引起动能损耗导致阻力增大，但结果表明n=2元件和n=3元件相差并不大，考虑这是由于n=3元件两侧的弧形槽分布位置所导致，该元件上下两个弧形槽紧靠壁面，而壁面处流体的速度相对较大，会将更容易带出槽内的流体使得槽内流体与外界流体汇合时阻力减小，因此虽然其有3个弧形槽，但整体阻力并未出现剧增的现象;n=1元件由于单个弧形槽所引起的动能损耗较小，所以其阻力也低于其他两种元件。

2.3  综合传热性能研究

PEC值评价方法可以用来表征换热效果与流动阻力的综合影响[12]，表达式[13]为

。  （2）

图6显示了三种不同弧形槽数量的元件的PEC值，可以看到，n=2元件的综合传热性能明显优于n=1元件，而n=3元件则呈现较大波动，其PEC值在1.28～1.33之间，且雷诺数处于4 000～6 300时其综合传热性能表现较好，但随着雷诺数增大PEC值有较大的下降幅度，Re=6 300～10 000时已经低于n=2元件，而n=1元件在所测Re范围内PEC值低于其他两种元件。由此可以看出，在合适的工况下选择最佳性能的元件是非常有必要的。

3  结 论

本文研究了不同弧形槽数量的扭旋元件传热性能与阻力特性等，并与光管进行对比以此体现弧形槽扭旋元件的强化换热效果，分析了管内流场及温度分布情况，揭示了不同弧形槽数量的元件对流动换热影响的机理，研究结论如下：

1）装有弧形槽扭旋元件的管道较光管具有更佳的换热效果，同时阻力也相应增大，三种弧形槽扭旋元件的PEC均大于1，说明其综合传热性能优于光管。

2）弧形槽的数量和分布位置对换热效果的影响较大，n=3元件较其他两种有着较高的平均努塞尔数，较n=1和n=2元件最高增加了4%和15.9%。

3）装有弧形8d83ec67d92aa5ec3f19b7233b89c415fecab118394f87422c9165a895433f50槽扭旋元件管道的阻力系数随凹数量的增加而提升，n=1元件所引起的阻力较小，而n=2元件和n=3元件阻力系数提升明显且两者差距较小。

4）在Re=4 000～6 300时n=3元件的PEC值最高，达到1.33;Re=6 300～10 000时则n=2元件表现较好，n=1元件的综合传热性能在所测Re范围内均不如其他两种元件。

参考文献：

[1]LIN Y. Performance investigation on an air source heat pump system with latent heat thermal energy storage[J]. Energy 2021 239（PA）.

[2]CHANG S W YANG T L LIOU J S. Heat transfer and Pressure Drop in Tube with Broken Twisted Tape insert[J]. Therm Fluid 2007 32： 489-501.

[3]HANG W JAN Y J LIOU J S. Turbulent Heat Transfer and Pressure Drop in Tube Fitted with Serrated Twisted-tape[J]. Int. J. Therm 2007 46： 506-5181.

[4]]张晓屿 刘志春 刘伟. 内置螺旋片的强化传热管的数值模拟研究[J]. 工程热物理学报 2011 32 （003）： 448-450.

[5]郭剑 杨昆 刘伟. 圆管插入十字形扭带强化传热数值模拟[J]. 工程热物理学报 2009 30 （7）： 1216-1218.

[6]AGARWAL S K  RAO M R . Heat transfer augmentation for the flow of a viscous liquid in circular tubes using twisted tape inserts[J]. International Journal of Heat & Mass Transfer 1996 39 （17）： 3547-3557.

[7]齐辉. 静态混合器的流体力学性能和传热性能研究[J]. 化工装备技术 2006 27 （003）： 66-67.

[8]KUMAR V SHIRKE V NIGAM K. Performance of Kenics static mixer over a wide range of Reynolds number[J]. Chemical Engineering Journal 2008 139 （2）： 284-295.

[9]MURUGESAN P MAYILSAMY K SURESH S et al. Heat transfer and pressure drop characteristics in a circular tube fitted with and without V-cut twisted tape insert[J]. International communications in heat and mass transfer 2011 （3）： 329-334.

[10]詹三江 林清宇 孙瑞娟 等. 换热管内置开边槽扭带的阻力与传热特性研究[J]. 化工科技 2012 20 （1）： 6-9.

[11]孙达 仇枭逸 奚桢浩 等. 中国结式静态混合器的流体力学特性[J]. 石油化工 2020 49 （10）： 40-46.

[12]GUO Z Y ZHU H Y LIANG X G. Entransy-a Physical Quantity Describing Heat Transfer Ability [J]. International Journal of Heat and Mass Transfer 2007 50 （13）：2545-2556.

[13]张傲宇，丁桂彬，王宗勇.同心双扭旋元件对管内层流传热特性的影响[J].当代化工，2021，50（01）：229-232.

Study on Turbulent Heat Transfer Characteristics in Circular

Tube With Built-in Arc Groove Torsional Element

ZHANG Chun-mei LIU Bin LI Yu-xuan LI Quan YANG Peng-da LI Shuai

（School of Energy and Power Engineering Shenyang University of Chemical Technology Shenyang Liaoning 110142 China）

Abstract： A new type of enhanced heat exchange element-arc groove torsional element was designed and the torsional elements with three kinds of arc grooves were numerically simulated in the range of Re = 4 000～10 000. The influence of the number of arc grooves on the heat transfer characteristics of the pipeline was analyzed from the aspects of average Nusselt number and temperature field distribution. The simulation results showed that the number and distribution position of n = 3 element arc-shaped grooves were conducive to the mixing of the central fluid and the fluid near the pipe wall so as to strengthen the heat transfer. Among the three elements n = 2 element had the best heat transfer effect followed by n = 2 element. The structure of n = 1 element was not conducive to the intersection of the fluid in the groove and the external fluid and the heat transfer effect was weak; The pipe resistance caused by the three elements increased with the increase of the number of arc grooves. In terms of comprehensive heat transfer performance the PEC value was greater than 1 the maximum value was 1.33 when n  was 3 and Re was 4 000 indicating that the use of arc groove torsional elements improved the heat transfer performance of the pipe by different extents.

Key words： Numerical simulation; Torsion element; Heat transfer characteristics; PEC value