欢迎访问《新能源进展》官方网站!今天是

扭曲椭圆管内强化传热及熵产特性

  • 冯俊奇 ,
  • 田果 ,
  • 王梦轩 ,
  • 谈莹莹 ,
  • 王占伟 ,
  • 王林 ,
  • 李修真 ,
展开
  • 河南科技大学 建筑能源与热科学技术研究所,河南 洛阳 471023
† 通信作者:李修真,E-mail:

作者简介:冯俊奇(1995-),男,硕士研究生,主要从事对流传热强化及工业节能研究。李修真(1990-),男,博士,讲师,主要从事传热强化新技术与工业节能研究。

收稿日期: 2023-12-05

  修回日期: 2023-12-29

  网络出版日期: 2024-08-30

基金资助

河南省高校科技创新团队项目(22IRTSTHN006)

河南省高校科技创新人才项目(22HASTIT025)

Heat Transfer Enhancement and Entropy Generation Characteristics in Twisted Elliptical Tubes

  • Junqi FENG ,
  • Guo TIAN ,
  • Mengxuan WANG ,
  • Yingying TAN ,
  • Zhanwei WANG ,
  • Lin WANG ,
  • Xiuzhen LI ,
Expand
  • Institute of Building Energy and Thermal Science and Technology, Henan University of Science and Technology, Luoyang 471023, Henan, China

Received date: 2023-12-05

  Revised date: 2023-12-29

  Online published: 2024-08-30

Copyright

版权所有 © 《新能源进展》编辑部

摘要

为探究扭曲椭圆管内强化传热及熵产特性,基于实验验证,在雷诺数(Re)为710 ~ 4 790范围内,采用数值模拟方法着重分析了包含扭曲比(P/D)和长短轴比(A/B)在内的结构参数及Re对扭曲椭圆换热管内传热和流动阻力性能的影响规律,并依据熵产理论分析扭曲椭圆管内的熵产特性,以期揭示扭曲椭圆管内流体对流传热强化机理。结果表明:扭曲椭圆管内对流传热性能、流动阻力及总熵产随Re的增大而分别增大、减小及增大,随P/D的减小或A/B的增大均增大,且总熵产主要由传热过程不可逆损失引起;当A/B = 2.16和P/D = 11.90时,传热系数(h)和压降(Δp)均达到最大值;此外,当P/D = 26.33和A/B = 1.34时,总熵产最小。

本文引用格式

冯俊奇 , 田果 , 王梦轩 , 谈莹莹 , 王占伟 , 王林 , 李修真 . 扭曲椭圆管内强化传热及熵产特性[J]. 新能源进展, 2024 , 12(4) : 473 -481 . DOI: 10.3969/j.issn.2095-560X.2024.04.013

Abstract

To investigate the heat transfer enhancement and entropy generation characteristics of twisted elliptical tubes, based on the experimental verification, the effects of the structural parameters, including twist ratio (P/D), ratio of major to minor axis (A/B), and the flow parameter Reynolds number (Re) on the heat transfer and flow resistance performance in twisted elliptical tubes were analyzed by numerical simulation method in the range of Re from 710 to 4 790. According to entropy generation theory, the entropy generation characteristics of the twisted elliptical tubes were discussed to reveal the mechanism of fluid convection heat transfer enhancement inside twisted elliptical tubes. The results showed that the convective heat transfer performance, flow resistance, and total entropy generation inside the twisted elliptical tubes were increased, decreased, and increased with the rise of Re, respectively; all increased with a decrease of P/D or an increase of A/B, and irreversible losses in the heat transfer process mainly caused the total entropy generation. When A/B = 2.16 and P/D = 11.90, the heat transfer coefficient (h) and pressure drop (Δp) reached their maximum value. Furthermore, the total entropy generation reached the minimum at P/D = 26.33 and A/B = 1.34.

0 引言

在全球大力开展节能减碳的科技运动中,研究者们已意识到提高能源利用效率是解决日益增长的能源需求和有限能源间矛盾的有效途径之一[1]。换热器作为一种应用最为广泛的能量转换设备,其换热性能的优劣直接影响能量传输效率的高低。管壳式换热器仍然是目前应用最为普遍的换热器,其性能的提高很大程度上依赖于换热管内外空间结构和流体流动模式的改善[2]。在此背景下,国内外学者对形式多样的新型强化传热管如螺旋管[3]、椭圆管[4]、交叉缩放椭圆管[5]等管型开展了理论研究和技术研发。而扭曲椭圆管[6]凭借结构简单、生产成本低且传热性能优异等优点,吸引了许多学者的研究兴趣。
研究者们对扭曲椭圆管内流体流动特性和传热性能开展了实验和数值模拟研究。YANG等[7]在雷诺数Re为600 ~ 55 000范围内,以水为工质实验研究了扭曲椭圆管内层流、过渡流和湍流时的传热和流动特性,得到了传热系数和流动阻力的准则关联式。还有学者对扭曲椭圆管外流动传热特性进行了研究。思勤等[8]对扭曲椭圆管换热器壳侧和管侧工作介质分别为水和油的传热性能进行了实验研究,得到了壳侧和管侧努塞尔数和摩擦系数的实验关联式。LIU等[9]和DONG等[10]通过对实验数据的拟合分别得到了在纵向流情况下扭曲椭圆管束对流传热性能和流动阻力特性计算的准则关联式。这些扭曲椭圆管内外的准则关联式对设计换热器具有理论指导意义。
还有一些学者利用场协同原理和扭曲椭圆管内二次流的演变规律,对强化传热性能进行分析以探索对流传热强化机理。马芳芳等[11]利用计算流体力学(computational fluid dynamics, CFD)技术对扭曲椭圆管内介质为水的传热和流动阻力特性进行了数值模拟,得到不同工况下的温度场、压力场及速度场,并基于场协同原理对其温度与速度梯度的协同效果及压力与速度梯度的协同效果进行比较和分析,获得不同工况下传热与流动阻力特性随长短轴比变化的规律,分析了影响扭曲椭圆管内流体传热效果与流动阻力损失的因素及其传热强化的机理。CHANG等[12]发现扭曲椭圆管内二次流强度以及传热和流阻特性均与长短轴比和扭曲比紧密相关,长短轴比越大或扭曲比越小,二次流强度越强且流动阻力也越大,进而传热强化效果就越显著。近些年,众多学者将热力学定律和火积理论用于评价扭曲椭圆管内对流传热性能。CHEN等[13]比较了采用熵产理论和火积耗散理论对传热强化的优化效果,研究表明,在相同的限制条件下,熵产最小的优化导致最佳的热功转换效率,而火积耗散极值原理的优化导致最佳的对流传热效率。EBRAHIMI等[14]采用热力学第二定律探究了微型扭曲椭圆管内强化传热机理,结果表明,较低的能量不可逆损失导致扭曲椭圆管内传热增强。
综上所述,目前针对对流传热优化设计的研究主要是基于热力学第一定律,即从能量的数量上进行传热强化评价,而对流传热中包含温差的热量传递和流体流动摩擦阻力损失,这些过程均为不可逆过程,往往会引起能量在质上的贬值。传热及流动中的不可逆损失所引起熵产已得到重视,并被用于评估热力系统热力学性能和热力设备的结构优化[15,16,17]。如ZHENG等[18]通过熵产分析对肋槽换热管的热工水力性能进行了数值研究,结果表明,肋槽换热管的最佳设计参数组合可以使熵产最小且综合传热性能最优。尽管熵产最小化原理已普遍应用于多种强化换热管的优化设计,但将熵产理论用于分析扭曲椭圆管内对流传热特性的研究还较为少见。为此,本文从热力学第二定律的角度,利用熵产最小化原理,分析不同结构参数的扭曲椭圆管内总熵产的变化规律,以期揭示扭曲椭圆管的传热强化机理。

1 数值计算方法

1.1 物理模型

扭曲椭圆强化传热管结构如图1所示,采用Solidworks软件对其建立模型,其由内径(d)为20 mm的铝质光滑圆管加工而成,对于具有不同长短轴比的扭曲椭圆管所对应的当量直径(D)取值分别为19 mm和17 mm。由中心复合设计法得出扭曲比和长短轴比的设计参数,具体参数见表1。测试段(L2)作为计算区域,长度为1 000 mm,为了保证在测试段入口处为充分发展的湍流并且消除出口回流对测试段的影响,设计了上游段(L1)和下游段(L3),长度分别为100 mm和160 mm。长短轴比(A/B)定义为椭圆横截面长轴和短轴的比值;扭曲节距(P)定义为椭圆形横截面扭转360°所对应管轴方向的长度;扭曲比(P/D)定义为扭曲节距和当量直径的比值。
Fig. 1 Schematic diagram of enhanced heat transfer tube structure

图1 强化传热管结构示意图

Table 1 Values of design parameters

表1 设计参数的取值

设计参数 取值
Re 710、1 500、2 750、4 000、4 790
A/B 1.34、1.50、1.75、2.00、2.16
D/mm 19、19、17、17、17
P/D 26.33、23.53、19.12、14.70、11.90

1.2 控制方程

在三维稳态数值模拟中,采用Renormalization group(RNG)k-ε湍流模型对流体进行计算,控制方程为质量、动量和能量方程[19],计算式如下。
质量方程:
$\frac{\partial \left( \rho {{u}_{i}} \right)}{\partial {{x}_{i}}}=0$ (1)
动量方程:
$\frac{\partial \left( \rho {{u}_{i}}{{u}_{j}} \right)}{\partial {{x}_{i}}}=-\frac{\partial p}{\partial {{x}_{i}}}+\mu \frac{\partial }{\partial {{x}_{j}}}\left( \frac{\partial {{u}_{i}}}{\partial {{x}_{j}}}+\frac{\partial {{u}_{j}}}{\partial {{x}_{i}}} \right)$ (2)
能量方程:
$\begin{align} & \frac{\partial }{\partial {{x}_{i}}}\left( \rho {{T}_{\text{m}}} \right)+\frac{\partial }{\partial {{x}_{i}}}\left( \rho {{u}_{i}}{{T}_{\text{m}}} \right)=-p\frac{\partial {{u}_{i}}}{\partial {{x}_{i}}}+ \\ & \lambda \frac{\partial }{\partial {{x}_{j}}}\left( \frac{\partial {{u}_{i}}}{\partial {{x}_{j}}}+\frac{\partial {{u}_{j}}}{\partial {{x}_{i}}} \right) \end{align}$ (3)
湍流动能k方程和湍流动能耗散率ε方程[20]的计算式如下。
湍流动能k方程:
$\begin{align} & \frac{\partial \left( \rho {{u}_{i}}k \right)}{\partial {{x}_{i}}}=\frac{\partial }{\partial {{x}_{j}}}\left[ \left( \mu +\frac{{{v}_{\text{m}}}}{{{\sigma }_{k}}} \right)\frac{\partial k}{\partial {{x}_{j}}} \right]+ \\ & {{\eta }_{t}}\frac{\partial {{u}_{i}}}{\partial {{x}_{j}}}\left( \frac{\partial {{u}_{i}}}{\partial {{x}_{j}}}+\frac{\partial {{u}_{j}}}{\partial {{x}_{i}}} \right)-\rho \varepsilon \end{align}$ (4)
湍流动能耗散率ε方程:
$\begin{align} & \rho \frac{\partial \left( {{u}_{j}}\varepsilon \right)}{\partial {{x}_{j}}}=\frac{\partial }{\partial {{x}_{j}}}\left[ \left( \mu +\frac{{{v}_{\text{m}}}}{{{\sigma }_{\varepsilon }}} \right)\frac{\partial \varepsilon }{\partial {{x}_{j}}} \right]+{{c}_{1}}\rho S\varepsilon - \\ & {{c}_{2}}\rho \frac{{{\varepsilon }^{2}}}{k+\sqrt{v\varepsilon }} \end{align}$ (5)
式中:${{c}_{1}}=\max \left[ 0.43,\frac{\eta }{\eta +5} \right]$;$\text{ }\eta \text{=}\frac{Sk}{\varepsilon }$;$\text{ }S=\sqrt{2{{S}_{i,j}}\text{ }{{S}_{i,j}}}$;${{S}_{i,j}}=\frac{1}{2}\left( \frac{\partial {{u}_{i}}}{\partial {{x}_{j}}}+\frac{\partial {{u}_{j}}}{\partial {{x}_{i}}} \right)$;c2 = 1.9;σε = 1.2;σk = 1.0;μ= cμρk2

1.3 边界条件及计算方法

数值模拟计算时边界条件设定如下:入口处给定均匀速度,流体温度设置为303.15 K,出口处为压力出口;测试段管壁面为无滑移壁面且壁温恒定为343.15 K,上游段和下游段壁面均为绝热壁面;工作介质设置为不可压缩流——水(热物理性质见表2),不考虑重力和黏性耗散的影响;选择SIMPLE算法将动量方程和连续性方程分别离散以得到满足连续性条件的速度场和压力场;动量、能量、湍流动能及湍流动能耗散率的求解采用二阶迎风格式;收敛的残差标准设置为1.0 × 10-6
Table 2 Thermophysical properties of working fluid

表2 工质热物理性质

参数 取值
比热容cp/[J/(kg∙K)] 4 174
动力黏度μ/(Pa∙s) 5.49 × 10-4
密度ρ/(kg/m3) 988.1
普朗特数Pr 3.54
热导率λ/[W/(m∙K)] 0.648

1.4 参数定义

平均对流传热系数、当量直径、雷诺数、努塞尔数、流动摩擦系数和压降公式列于表3
Table 3 Parameter definition formula for twisted elliptical tubes

表3 扭曲椭圆管参数定义公式

变量名称 公式 序号
对流传热系数 $h={q}/{\left( {{T}_{\text{w}}}-{{T}_{\text{m}}} \right)}\;$ (6)
当量直径 $D={4{A}'}/{U}\;$ (7)
雷诺数 $Re={\rho {{v}_{\text{m}}}D}/{\text{ }\!\!\mu\!\!\text{ }}\;$ (8)
努塞尔数 $Nu={hD}/{\lambda }\;$ (9)
摩擦系数 $f={2D\Delta p}/{\left( {{L}_{2}}\rho v_{\text{m}}^{2} \right)}\;$ (10)
压降 $\Delta p={{p}_{\text{in}}}-{{p}_{\text{out}}}$ (11)

1.5 湍流模型验证

采用Fluent软件对三维计算区域进行数值模拟。为确保上述所建立模型及其数值计算方法的准确性与可靠性,选择合适的湍流模型对扭曲椭圆管内传热和流动特性进行模拟分析。在本文确定的工况下,基于RNG k-ε、Realized k-ε、Standard k-ω及Shear stress transport(SST)k-ω四种湍流模型进行数值计算,将数值模拟结果与实验结果[21]对比,如图2所示。由图可知,当采用RNG k-ε湍流模型时,Nuf的数值计算结果与实验结果最为接近,且与实验结果的平均偏差分别为8.5%和2.4%,在工程所允许的误差范围内,这说明采用RNG k-ε湍流模型对所建立模型进行数值模拟是准确且可靠的。
Fig. 2 Comparison of numerical simulation and experimental results

图2 数值模拟和实验结果对比

1.6 网格生成和独立性验证

使用计算流体力学前处理软件ICEM CFD对扭曲椭圆管模型进行网格划分。考虑模型结构的不规则性以及管内流动的湍动特性,对扭曲椭圆管进行网格划分时,管内近壁面处设置第一层网格高度的y+ 值接近1,具体的加密网格划分如图3所示。为避免网格数量对数值计算结果的影响及确定其结果的准确性,有必要对模型结构进行网格无关性验证,此处以管型(P/D = 19.12、A/B = 1.75)为例进行网格无关性分析。在Re为2 750时,选取四套不同尺寸的网格系统(粗、中细、细和更细)进行计算模拟,网格数分别为767 088、773 932、880 102、1 084 666。Nuf随网格数的变化如图4所示。由图可得,细(880 102)和更细(1 084 666)的网格系统计算出的Nuf相对偏差为1.1%和4.1%。在保证节约计算时间及计算结果准确条件下,综合考虑计算效率和求解精度,选用网格数为880 102的网格划分方法进行数值模拟计算。
Fig. 3 Grid division

图3 网格划分

Fig. 4 Variation of Nu and f with the number of grids

图4 Nuf随网格数的变化

2 结果与讨论

2.1 流场和温度场分布

在分析设计参数对管内传热和流动特性的影响之前,首先要探究管内流动曲线、速度分布和温度分布。利用Tecplot软件对扭曲椭圆管内流场和温度场进行可视化后处理,以反映流动和传热特性。
Re为4 000时,不同节距的扭曲椭圆管(A/B = 2.00)内局部流动曲线和速度分布如图5所示,在管中心位置流体沿管轴线方向流动且速度随P的减小而增大。由于存在壁面剪切应力,管内中心处的速度大于近壁面处的速度,流线由管中心直线形状逐渐向管壁面变化为螺旋扭曲形状,随P减小,管壁附近流体流动曲线螺旋化程度更显著。流体在扭曲流道引导下,受管壁面扰动的影响,在管壁附近形成螺旋流,且在离心力的作用下,流动方向发生周期性改变。此外,可以看出螺旋流随P的减小更加剧烈,主要是由于P越小,在一定长度的扭曲椭圆管上的节距数量越多,流体湍动程度增强,对壁面附近的边界层扰动作用增强,强化了传热效果。
Fig. 5 Local flow curve and velocity distribution in twisted elliptical tube

图5 扭曲椭圆管内局部流动曲线和速度分布

Re为4 790时,扭曲椭圆管(A/B = 2.16,P/D = 11.90)轴向不同横截面的温度分布如图6所示。由图可知近壁面区域的温度较高,且随流体从入口端沿轴线朝向出口端流动时,管中心区域的流体温度显著升高,说明扭曲椭圆管内较大径向速度形成的螺旋流改变了管内流场结构,使温度场分布发生了变化,促进了中心区域和近壁面区域之间流体的混合,流体从近壁面向中心区域传递更多的热量以实现管内强化传热。在中下游横截面上的温度分布更均匀,主要归因于近壁面的温度梯度更大,增大了流体湍流强度进而增强了流体的扰动。此外,由于在长轴端形成的涡流方向相反,横截面上的温度分布形状大致呈现“S”形。
Fig. 6 Temperature distribution of twisted elliptical tube on different axial cross sections

图6 扭曲椭圆管在轴向不同横截面上的温度分布

2.2 结构参数对传热和流动特性的影响

在传热计算中,h和ΔP通常表示传热性能和流动阻力损失。为揭示扭曲椭圆管的长短轴比对传热性能和流动特性的影响,当P/D = 19.12时,考察不同长短轴比的扭曲椭圆管内h和ΔPRe的变化,如图7所示。不同长短轴比的扭曲椭圆管内h和ΔP都随Re的增大而增大,说明传热性能增强必然使压降损失增大。对于给定的Re,h和ΔP均随长短轴比的增大而增加,表明较大A和较小B能够使切向和径向速度快速交替变化并生成更多的小尺度涡流,促进核心区域和近壁面区域之间流体的混合以提高传热效率。但由于速度在切向和径向方向快速变化,流体沿着流动方向形成较大的压力梯度,且在近壁面区域集中混合,导致更大的流动阻力损失。
Fig. 7 Variation of h and ΔP with Re in twisted elliptical tubes with different A/B when P/D = 19.12

图7 P/D = 19.12时,不同A/B的扭曲椭圆管内h和ΔPRe变化

A/B = 1.75时,不同扭曲比的扭曲椭圆管内流体h和ΔPRe的变化如图8所示。根据数值模拟结果,不同扭曲比的扭曲椭圆管内hRe的增加和扭曲比的减小而增加。类似地,对流传热性能的增强往往引起流动阻力的增大,主要归因于管壁的扭曲程度越大,近壁面流体与壁面冲刷碰撞的强度越大,削弱了流体沿流动方向的动能,却增大了局部流动的涡流强度进而强化传热性能。上述所得结果与KIM等[22]的研究结论一致。
Fig. 8 Variation of h and ΔP with Re in twisted elliptical tubes with different P/D when A/B = 1.75

图8 A/B = 1.75时,不同P/D的扭曲椭圆管内h和ΔPRe的变化

2.3 熵产分析

任何被动强化传热技术均会增大流动功耗,不可避免地增加熵产。为提高强化传热效率,并获得传热过程中引起的不可逆与换热管传热性能之间的关系,利用BEJAN[15]提出的熵产最小化方法分析扭曲椭圆管内由传热和流动过程引起的不可逆损失以估算复杂传热和流动过程的能量耗散,该方法是一种从热力学第二定律角度对特定热力学系统优化的方法。在熵产最小化的情况下对流动通道或热力学系统进行优化设计,由于流体在水平扭曲椭圆管内流动,重力对流动的影响较小,且速度大小几乎不变,势能和动能变化忽略不计。沿着流动方向长度为dz的扭曲椭圆管微元体为热力学系统,热力学第一定律和第二定律表示如下。
$\dot{m}\text{d}h\text{=}{q}'\text{d}z$ (12)
$\text{d}{{{S}'}_{\text{gen}}}=\dot{m}\text{d}s-{{q}'\text{d}z}/{\left( T+\Delta T \right)}\;$ (13)
由热力学关系式得:
${\text{d}h}/{\text{d}z}\;=T\left( {\text{d}s}/{\text{d}z}\; \right)+\left( {1}/{\rho }\; \right)\left( {\text{d}P}/{\text{d}z}\; \right)$ (14)
$\Delta {{{S}'}_{\text{gen}}}={\text{d}{{{{S}'}}_{\text{gen}}}}/{\text{d}z}\;$ (15)
$\Delta {{{S}'}_{\text{gen}}}={{q}'\Delta T}/{{{T}^{2}}\left( 1+\tau \right)}\;-\left( {{\dot{m}}}/{\rho T}\; \right)\left( {\text{d}P}/{\text{d}z}\; \right)$ (16)
压降表示为:
$\text{d}P=-\left( {2f{{G}^{2}}}/{\rho D}\; \right)\text{d}z$ (17)
质量流速:
$G={{\dot{m}}}/{{{A}'}}\;$ (18)
无量纲温差$\tau \text{=}{\left( {{T}_{\text{w}}}-{{T}_{\text{in}}} \right)}/{{{T}_{\text{in}}}}\;$忽略不计,q′和ΔTSt形式表示为:
$St={{{q}'}}/{U\Delta T{{c}_{p}}G}\;$ (19)
代入等式(16)化简得:
$\Delta {{{S}'}_{\text{gen}}}={{{{{q}'}}^{2}}D}/{\left( 4{{T}^{2}}\dot{m}{{c}_{p}}St \right)}\;+{2{{{\dot{m}}}^{3}}f}/{\left( {{\rho }^{2}}TD{{{{A}'}}^{2}} \right)}\;$ (20)
通常,等式(20)是由传热和摩擦效应影响的熵产组合,由式(21)和式(22)表示如下。
${{\left( \Delta {{{{S}'}}_{\text{gen}}} \right)}_{\Delta T}}={{{{{q}'}}^{2}}D}/{\left( 4{{T}^{2}}\dot{m}{{c}_{p}}St \right)}\;$ (21)
${{\left( \Delta {{{{S}'}}_{\text{gen}}} \right)}_{\Delta P}}={2{{{\dot{m}}}^{3}}f}/{\left( {{\rho }^{2}}TD{{{{A}'}}^{2}} \right)}\;$ (22)
S′gen)ΔT和β(ΔS′gen)ΔP分别表示由传热和流动摩擦力引起的不可逆损失。PAOLETTI等[23]对由传热引起的不可逆损失相对于总熵产的大小提出了不可逆性分布参数Be,一个无量纲参数定义如下:
$Be\text{ =}{\text{ }{{\left( \Delta {{{{S}'}}_{\text{gen}}} \right)}_{\Delta T}}}/{\Delta {{{{S}'}}_{\text{gen}}}}\;$ (23)
图9(a、b)为A/B = 1.75和P/D = 11.90时,扭曲椭圆管内总熵产随Re的变化情况,对于给定的Re,扭曲椭圆管内总熵产随P/D的减小或A/B的增大而增大;对于相同的A/BP/D,总熵产随Re的增大均增大,说明总熵产变化受结构参数和流动参数共同作用的影响。当P/D = 19.12时,不同A/B的扭曲椭圆管内 (ΔS′gen)ΔT和β(ΔS′gen)ΔPRe的变化如图10所示,图中可见,(ΔS′gen)ΔT和β(ΔS′gen)ΔPRe的增加分别减小和增大,且β(ΔS′gen)ΔT均大于 (ΔS′gen)ΔP,说明熵产受传热过程的影响更大,这是由于Re增加使流体传热性能增强,流场中的温度分布更加均匀,整体流动温度降低,从而降低了由传热引起的不可逆损失。与此同时,Re增加会导致更大的流体黏度和速度梯度,增大了扭曲椭圆管内流动摩擦引起的不可逆损失。
Fig. 9 Variation of ΔSgen in a twisted elliptical tubes with different Re: (a) P/D (A/B = 1.75); (b) A/B (P/D = 11.90)

图9 扭曲椭圆管内ΔS′genRe的变化:(a)P/DA/B = 1.75);(b)A/BP/D = 11.90)

Fig. 10 Variation of (ΔSgen)ΔT and (ΔSgen)ΔP with Re in a twisted elliptical tubes with different A/B when P/D = 19.12

图10 P/D = 19.12时,不同A/B的扭曲椭圆管内β(ΔS′gen)ΔT和β(ΔS′gen)ΔPRe的变化

A/B = 1.50和P/D = 19.12时,扭曲椭圆管内BeRe的变化情况如图11所示。图中BeRe的增加而减小,表明流动摩擦对总熵产的影响逐渐增加,而传热对总熵产的影响降低。Be均接近1,说明传热对总熵产的影响远大于流动摩擦对总熵产的影响,A/B = 1.34和P/D = 26.33对应的Be最小。结果表明,Be随长短轴比的增大或扭曲比的减小而增大。说明流体流动时由传热引起的不可逆损失远大于由流动摩擦引起的不可逆损失[24]
Fig. 11 Variation of Be in a twisted elliptical tubes with different Re: (a) P/D (A/B = 1.50); (b) A/B (P/D = 19.12)

图11 扭曲椭圆管内BeRe的变化:(a)P/DA/B = 1.50);(b)A/BP/D = 19.12)

综上,根据熵产最小化原理对数值模拟结果进行分析以揭示扭曲椭圆管内传热强化机理。结果表明,扭曲椭圆管内提高传热性能的同时增大了由传热和流动过程引起的不可逆损失。

3 结论

基于实验验证,采用数值计算方法分析结构和流动参数对扭曲椭圆管内对流传热性能及流动阻力特性的影响,进而采用熵产理论揭示扭曲椭圆管内强化对流传热机理,在Re为710 ~ 4 790范围内,具体结论如下:
(1)不同结构参数的扭曲椭圆管内h和ΔPRe增大均增大;对于给定的Re,h和ΔPA/B的增大或P/D的减小而增大,而且在A/B = 2.16和P/D = 11.90,h和ΔP均达到最大值。
(2)不同结构参数的扭曲椭圆管内总熵产均随Re的增大而增大;对于一定的Re,总熵产随A/B的增大或P/D的减小而增大;(ΔS′gen)ΔT和(ΔS′gen)ΔPRe增大分别增大和减小。
(3)当P/D = 26.33和A/B = 1.34时,总熵产达到最小。当P/D = 11.90和A/B = 2.16时,Be均接近1,说明由传热引起的不可逆损失远大于由流动摩擦引起的不可逆损失。
符号表:
A 半长轴,mm
A′ 扭曲椭圆管横截面面积,m2
B 半短轴,mm
Be 不可逆性分布参数
D 当量直径,mm
d 圆管内径,mm
dS′gen 熵产率,W/K
ε 湍流动能耗散率,%
s 比熵产率,W/(kg∙K)
τ 无量纲温差,%
cp 比热容,J/(kg∙K)
ρ 密度,kg/m3
Pr 普朗特数
λ 热导率,W/(m∙K)
μ动力黏度,Pa∙s
f 流动摩擦系数
G 质量流速,kg/(m2∙s)
h 对流传热系数,W/(m2∙K)
k 湍流动能,m2/s2
L 管长,mm
L1 上游段长度,mm
L2 测试段长度,mm
L3 下游段长度,mm
$\dot{m}$ 质量流量,kg/s
Nu 努塞尔数
P 扭曲节距,mm
q 壁面与流体间的平均热流量,W/m2
q′ 换热管单位长度的热流密度,W
Re 雷诺数
St 斯坦顿数
T 温度,K
Tw 壁面温度,K
Tm 流体平均温度,K
U 横截面润湿周长,m
$\nu $ 运动黏度,m2/s
vm 管中流体的平均速度,m/s
Δp 测试段进口和出口的压差,Pa
ΔT 壁面和流体温度差,K
ΔS′gen 总熵,W/K
S′gen)ΔT 由传热产生的熵,W/K
S′gen)ΔP 由流体摩擦力产生的熵,W/K
下角标:
in 入口
out 出口
m 平均
w 壁面
[1]
FENG C, WANG M. Analysis of energy efficiency and energy savings potential in china’s provincial industrial sectors[J]. Journal of cleaner production, 2017, 164: 1531-1541. DOI: 10.1016/j.jclepro.2017.07.081.

[2]
ZHANG J F, LI B, HUANG W J, et al. Experimental performance comparison of shell-side heat transfer for shell-and-tube heat exchangers with middle-overlapped helical baffles and segmental baffles[J]. Chemical engineering science, 2009, 64(8): 1643-1653. DOI: 10.1016/j.ces.2008.12.018.

[3]
MISURATI K A, QUAN Y K, GONG W, et al. Contrastive study of flow and heat transfer characteristics in a helically coiled tube under uniform heating and one-side heating[J]. Applied thermal engineering, 2017, 114: 77-84. DOI: 10.1016/j.applthermaleng.2016.11.168.

[4]
HAN H, HE Y L, LI Y S, et al. A numerical study on compact enhanced fin-and-tube heat exchangers with oval and circular tube configurations[J]. International journal of heat and mass transfer, 2013, 65: 686-695. DOI: 10.1016/j.ijheatmasstransfer.2013.06.049.

[5]
MENG J A, LIANG X G, CHEN Z J, et al. Experimental study on convective heat transfer in alternating elliptical axis tubes[J]. Experimental thermal and fluid science, 2005, 29(4): 457-465. DOI: 10.1016/j.expthermflusci.2004.04.006.

[6]
CHENG J L, QIAN Z Q, WANG Q. Analysis of heat transfer and flow resistance of twisted oval tube in low Reynolds number flow[J]. International journal of heat and mass transfer, 2017, 109: 761-777. DOI: 10.1016/j.ijheatmasstransfer.2017.02.061.

[7]
YANG S, ZHANG L, XU H. Experimental study on convective heat transfer and flow resistance characteristics of water flow in twisted elliptical tubes[J]. Applied thermal engineering, 2011, 31(14-15): 2981-2991. DOI: 10.1016/j.applthermaleng.2011.05.030.

[8]
思勤, 夏清, 梁龙虎, 等. 螺旋扁管换热器传热与阻力性能[J]. 化工学报, 1995, 46(5): 601-608.

[9]
LIU S J, YIN Y D, TU A M, et al. Experimental investigation on shell-side performance of a novel shell and tube oil cooler with twisted oval tubes[J]. International journal of thermal sciences, 2020, 152(4): 106290. DOI: 10.1016/j.ijthermalsci.2020.106290.

[10]
DONG X Y, JIN X D, LI P Y, et al. Experimental research on heat transfer and flow resistance properties in spiral twisted tube heat exchanger[J]. Applied thermal engineering, 2020, 176: 115397. DOI: 10.1016/j.applthermaleng.2020.115397.

[11]
马芳芳, 云和明, 郭训虎, 等. 扭曲椭圆管强化换热与流动减阻特性的场协同分析[J]. 高校化学工程学报, 2017, 31(6): 1301-1310. DOI: 10.3969/j.issn.1003-9015.2017.00.045.

[12]
CHANG K S, CHOI J S, Kim J S. Laminar fluid flow in a twisted elliptic tube[J]. KSME journal, 1988, 2(1): 44-51. DOI: 10.1007/BF02944075.

[13]
CHEN Q, LIANG X G, GUO Z Y, et al. Entransy theory for the optimization of heat transfer-a review and update[J]. International journal of heat and mass transfer, 2013, 63(1): 65-81. DOI: 10.1016/j.ijheatmasstransfer.2013.03.019.

[14]
EBRAHIMI A, ROOHI E. Numerical study of flow patterns and heat transfer in mini twisted oval tubes[J]. International journal of modern physics C, 2015, 26(12): 1550140. DOI: 10.48550/arXiv.1811.02451.

[15]
BEJAN A. Entropy generation minimization[M]. New York, 1996, 79(3): 1191-1218.

[16]
KO T H. Thermodynamic Analysis of optimal mass flow rate for fully developed laminar forced convection in a helical coiled tube based on minimal entropy generation principle[J]. Energy conversion and management, 2006, 47(18-19): 3094-3104. DOI: 10.1016/j.enconman.2006.03.006.

[17]
TANDIROGLU A. Effect of flow geometry parameters on transient entropy generation for turbulent flow in circular tube with baffle inserts[J]. Energy conversion and management, 2007, 48(3): 898-906. DOI: 10.1016/j.enconman.2006.08.022.

[18]
ZHENG N B, LIU P, FENG S, et al. Numerical investigations of the thermal hydraulic performance in a rib-grooved heat exchanger tube based on entropy generation analysis[J]. Applied thermal engineering, 2016, 99: 1071-1085. DOI: 10.1016/j.applthermaleng.2016.02.008.

[19]
DARBANDI M, ABDOLLAHPOUR M S. An innovative cost-effective numerical-analytical procedure to derive heat correlations for shell and tube heat exchangers equipped with smooth and dimple tubes[J]. Numerical heat transfer, 2024 (early access). DOI: 10.1080/10407782.2023.2301521.

[20]
GU H D, CHEN Y P, WU J F, et al. Performance investigation on twisted elliptical tube heat exchangers with coupling-vortex square tube layout[J]. International journal of heat and mass transfer, 2020, (151): 119473. DOI: 10.1016/j.ijheatmasstransfer.2020.119473.

[21]
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.

[22]
KIM H R, KIM S, KIM M, et al. Numerical study of fluid flow and convective heat transfer characteristics in a twisted elliptic tube[J]. Journal of mechanical science and technology, 2016, 30: 719-732. DOI: 10.1007/s12206-016-0127-4.

[23]
PAOLETTI S, RISPOLI F, SCIUBBA E. Calculation of exergetic losses in compact heat exchanger passages[J]. American society of mechanical engineers, advanced energy systems, 1989, 10(2): 21-29.

[24]
GUO J F, CHENG L, XU M T. Multi-objective optimization of heat exchanger design by entropy generation minimization[J]. Journal of heat transfer, 2010, 132(8): 081801. DOI: 10.1115/1.4001317.

文章导航

/