搜索
您的当前位置:首页正文

波纹管内的流动与传热强化研究

来源:二三娱乐
维普资讯 http://www.cqvip.com

第33卷第3期 2006正 北京化工大学学报 JOURNAL OF BEIJING UNIVERSITY OF CHEMICAL TECHNOLOGy VoI.33.No.3 2006 波纹管内的流动与传热强化研究 肖金花 钱才富 黄志新 安维峥 (北京化工大学机电工程学院,北京 100029) 摘要:通过数值模拟和实验验证考察了流体在波纹管内的流动与传热情况,研究了不同的流体雷诺数、波形及管 径大小对波纹管强化传热性能的影响,探讨了其强化传热机理。研究发现,在湍流范围内,波纹管的强化传热倍数 达到相同条件下直管的1.78~3.10倍,且最佳强化效果出现在Re=10000附近;文中拟合出的考虑波形影响的波 纹管内传热准则方程,对波纹管的工程设计和应用具有一定的指导意义。 关键词:波纹管;数值模拟;实验验证;强化传热;准则方程 中图分类号:TQ051.5 作为一种新型高效管壳式换热器,波纹管换热 实验所需费用昂贵的问题,另一方面可以取得管内 器已越来越受到关注。它与传统的管壳式换热器结 构基本类似,区别仅在于用横截面周期性变化的波 纹管代替了传统的直管作为换热管。这种换热器具 有传热效率高、不易结垢、热补偿能力强等优点。其 核心换热元件——波纹管,通常分为不锈钢波纹管 流场和温度场的分布细节,探讨传热强化机理。本 文的目的旨在通过对波纹管单管内的流动和传热进 行数值模拟,来研究波纹管的传热强化效果及机理。 1 计算模型 本文首先分析了目前工业上常用的一种波纹 管,其小径D1=19 mm,大径D2=25 mm,波距n= 18 mm,大圆弧半径Rl=15 mm,小圆弧半径R2则 根据几何关系由D1、D2、R和a确定,见图1。 和碳钢波纹管两种。不锈钢波纹管是由薄壁不锈钢 管从管内胀扩加工形成,而碳钢波纹管则由碳钢管 从管外挤压成型,因此波纹管的纵截面可近似认为 由大小圆弧相切组成,如图1所示。 波纹管内的流动可以近似看成是轴对称的管内 不可压缩定常流动,故数值计算的控制方程采用二 f 图1波纹管纵截面图 Fig.1 Longitudinal section of the corrugated tubes 维不可压缩稳态流动的连续性方程、Navier—Stokes 动量方程和能量方程(忽略流体质量力): + P0 (1) 波纹管换热器推出的时间还不是很长,国内外 虽有一些学者对其进行了流动和传热方面的研 (“ + 考):一 + ( + )c2 究u—J,但都不是很深入,关于其强化传热机理的研 究,则更是鲜有报道。对波纹管传热性能的研究大 多以实验为主,随着计算机技术的发展和通用流体 分析软件的出现,人们有可能对流体在波纹管内的 流动和传热进行数值模拟,这一方面可以解决大量 P(“ + 考):一 3y+ ( +等)c3 3T+ 哥)= +矗(寒+熹)c4 “ 湍流模型采用标准矗一£湍流模型[引: a(kui)P = 收稿日期:2005—10—18 P 2 = ( + 一 ㈤ ( + )毒]+ 。 G (6)【6) 基金项目:北京市教育委员会共建项目(XK100100541) 第一作者:女,1979年生。博士生 *通讯联系人 E—mail:qiancf@mail.buct.edu.cn 上式中,忌是湍流动能,£是湍流耗散率, 为 湍流黏性系数, =lD・ ・k 2/£, 、C£l、C£2、O"k和 维普资讯 http://www.cqvip.com

第3期 肖金花等:波纹管内的流动与传热强化研究 ・69・ 为常数,其值分别是:c =0.09,C l=1.44, 2= 1.92,O"k=1.0, =1.3。 波纹管的流动和传热边界条件是:入口给定流 速,并假定入口横截面上各点流速均相同;出口相对 压力为零;壁面无滑移,60"C恒温。波纹管内的流动 介质采用常温下的水,近似处理为不可压缩的常物 性流体,且忽略其重力。 为进行对比,本文还对一长度与波纹管长度相 同,且内径等于波纹管小径的直管进行了同样条件 下的数值计算。 2 数值方法 本文采用有限单元法对上述控制方程组进行离 散化,应用流线二阶迎风来离散对流项,应用SIM— PLE算法进行压力修正,代数方程采用超松弛法进 行迭代求解。计算采用的网格划分如图2所示。其 中近壁处网格加密。 图2计算网格示意图 Fig.2 Computing mesh diagram 3 计算结果及讨论 根据流体物性,以波纹管小径为特征尺寸得出 流体流动的特征参数: RP: ;Pr: ;N“: H K K 3.1流动特性 图3和图4为入口雷诺数Re=56614时,计算 所得的充分发展段波纹管横截面速度及湍流动能分 布与直管的比较,所取横截面如图2所示,图中y= 0为管中心线。由图3可知,波纹管靠近管中心区 域的横截面流速比直管有较大提高,且在波纹管的 波峰处,截面流速出现了负值,这意味着在波峰处产 生了回流。图4则说明波纹管横截面的湍流动能要 远远大于直管,尤其在管中心区域附近,这意味着波 纹管内流体的湍动程度相对于直管来说要剧烈得 多,因此其应该具有很好的传热效果。 图5所示为Re=56614时波纹管波峰截面处 流速方向角的分布及其与直管的对比。方向角的定 义为arctan(v/U)。由图中可以看出。波纹管横截面 上靠近波峰处流动方向角为负值,这也意味着波峰 内存在负向流速,图6所示的流线图也说明了这一 点。图5同时还表明,由于截面负向速度区域的存 在,波纹管相对于直管也存在较大的流动阻力。 横截面速度大小/m.s 图3横截面速度分布比较(Re=56614) Fig.3 Comparison of cross—sectional velocity magnitudes at Re=56614 图4横截面湍流动能分布比较(Re:56614) Fig.4 Compariosn of cross-sectional turbulence kinetic energies at Re=56614 图5波纹管波峰截面上速度方向分布 Fig.5 Velcoity direction distribution in section 3 of the two tubes 3.2传热特性 为研究波纹管的传热性能,本文在进行数值计 维普资讯 http://www.cqvip.com

・70・ 北京化工大学学报 2006正 量纲Nu数则随Re的增大而减小。这说明,在Re =10000附近是波纹管强化传热效果最佳的一段区 域,此时波纹管的平均传热系数达到了直管的2.8 倍。而Re=10000左右正是直管中的流动从过渡 流进入湍流的临界区域。实验研究[ 一]表明,流体 图6 Re=56614时波纹管充分发展段流线图 Fig.6 Streamline distibutiron in the corrugated tube 在波纹通道中从层流过渡到湍流时的临界雷诺数要 比直管低得多,这就意味着,当直管中的流动刚进入 湍流状态时,波纹管中的流动已处于中等或高度湍 at fully developed regions when Re=56614 算的同时,还进行了实验验证。实验用的波纹管由 外径20mm的薄壁不锈钢管从管内胀扩加工形成, 其公称规格是25 mm×0.8 mm×2500 mm,实测数 据是波峰直径24.5 mm,波谷直径20.2 mm,波距 18.2 mm。实验工质采用常温下的水。实验时在管 外壁包上带有电阻丝的保温层,以施加恒定热流边 界,在管入口处改变流量,以测量不同流量下的进出 口温度和压降等数据。 图7所示为波纹管与直管的努塞尔数N“随雷 窆 图7波纹管与直管努塞尔数对比 Fig.7 Comparison of Nu values for the two tubes 诺数Re的变化关系。由图可知,波纹管的N“数比 直管的大很多,强化传热效果比较显著。 图8则显示了波纹管与直管的N“数比值 图8波纹管内Nu /Nu。随融的变化关系 Fig.8 Relationship between Nu c/Nu s and Re (N“ /N“ ,称之为无量纲N“数)随RP的变化关 系。由图8可知,当4718<Re<10000时,无量纲 N“数随RP的增大而增大,而当RP>10000后。无 流状态,故此时波纹管相对于直管的强化传热效果 最佳。随着管内湍流进一步发展,直管和波纹管内 都出现了大量而剧烈的旋涡,这些旋涡不断冲刷管 壁,使层流边界层减薄,传热热阻减小,故此时边界 层对直管和波纹管传热性能的影响已相差不大,因 此波纹管相对于直管的强化传热倍数也逐渐减小, 但最小也达到直管的1.9倍左右。 图7和图8还表明,数值计算结果与实验结果 非常吻合。 3.3结构参数对传热强化的影响 波纹管的纵截面波形是由大小相切的圆弧组 成,因此其大小圆弧半径不同,管内的流动和传热情 况也不尽相同。此外,波纹管的基管直径也会影响 成型后的波纹管内的流动与传热,因此本文研究了 工程中常用的三种不同管径的波纹管。 图9所示为波纹管的基直管径一定、而大小圆 弧半径改变时,波纹管的无量纲N“数随Re的变化 图9相同管径不同波纹半径时 Nu /Nu。随Re的变化关系 Fig.9 Relationship between Nu c/Nu s and Re at the same pipe diameter and different wave radii 关系。由图9可知,此时波纹管的强化传热趋势一 致:且最佳强化传热区仍然是在Re=10000左右, 这意味着波纹管基管直径一定时,其最佳强化传热 区的位置与波形无关。此外,图9还表明,波纹管的 大小圆弧半径越小,即截面变化越剧烈,其强化传热 效果也越高,最大强化倍数可达到3.1倍左右。 维普资讯 http://www.cqvip.com

第3期 肖金花等:波纹管内的流动与传热强化研究 图10所示是波纹管基管直径不同但大小圆弧 半径相同时,波纹管内无量纲N“随Re的变化关 系。由图10可知,不同基管直径对波纹管内流动和 善 窆 图10不同管径时N“ /Nu 随Re的变化关系 Fig.10 Relationship between Nu c/Nu s and Re at different pipe diameters 传热的影响不大;波纹管的总体强化传热趋势一 致,都存在一个最佳强化区,但最佳强化区出现的位 置稍有不同;随波纹管基管直径的增大,波纹管的最 佳强化传热区稍有后移,且最大强化倍数也略有下 降。这些结论可为波纹换热管的工程应用提供一定 的依据。 3.4波纹管管内传热准则方程的拟合 根据上文所做的研究和分析,本文拟合出了基 于传统直管传热准则方程的波纹管管内的传热准则 方程,并给出方程的形式为 Nu =f・Nu (7) 上式中,N“ 和N“ 分别为波纹管和传统直管 的努塞尔数,其中Nu 可由经典的Dittus.Boelter关 联式[ ]表示。厂为多项式系数,包含了Re和波纹 管大小圆弧半径Rl和R2的影响,故,可写成: f=,(Re)f(Rl,R2) (8) 式(8)中,没有包含波纹管管径的影响,这是因 为由图10可知,管径的影响不大,因此可从图10中 近似取中间一根曲线来分段拟合厂(Re): 5000<Re<11000时,,(RP)=Al+AzRe(9) 11000<Re<75000时,厂(Re)=A 3+A 4Re (10) ,(Re)确定后,根据图9可拟合出,(Rl,R2), 即 ,(Rl,RE)=A5+A6 Rl+A7 R2 (11) 式(9)到(11)中,Al到A7为带单位的拟合系 数,其值和结构参数的单位相关,当结构参数单位为 mm时,拟合系数的值分别为:Al=2.02646,A2= 1.11827×10一 ,A 3=2.91323,A4=一3.50207× 10一,A5=2.2065,A6=一0.01216,A7=~ 0.00411。 公式(7)的应用范围是:Dl=19~32 mm,D2= 25--42 mm,4718<Re<95350,5 mm<Rl,RE<15 mm。在此范围内,数值计算结果与拟合公式(7)的 计算结果最大误差小于8%。和一些实验拟合关联 式不同的是,这里提出的准则方程考虑了波纹管的 不同管径和波形尺寸的影响,因此工程应用性更广。 4 结论 (1)横截面上湍流动能大于相同条件下的直管, 波峰内存在回流; (2)波纹管能显著强化管内传热,其强化倍数达 到相同条件下直管的1.78~3.10倍,且最佳强化传 热区在Re=10000附近; (3)波纹管的管径对其强化效果影响不大,在管 径一定时减小波纹半径能增大其强化传热倍数; (4)拟合出了波纹管管内传热准则方程,这对于 波纹管的工程设计和应用具有一定的参考价值。 符号说明 U, ——直角坐标系下x和Y方向的流体速度分量,m・sI1 T——流体温度.℃ ——流体压力,Pa ——流体动力粘度,Pa・S p——流体密度,kg・m-3 C ——流体定压比热,kJ・kgI1・℃ ——流体导热系数,w・m_2・℃I1 a——对流传热系数。W・m ・℃I1 De——波纹管当量直径 Re——Rey110lds数Pr——Prandtl数 参考文献 [1]俞惠敏,蔡业彬.几种换热管强化传热性能实验分析与 比较[J].流体机械,2003,31(6):7—10. [2] 商福民,毕庆生,张志正,等.波纹管与光管流体传热 性能对比实验研究[J].长春工程学院学报,2003,4 (3):33—35. [3】RussG,Beer H.Heat transfer and flow field in a pipe with sinusoidal wavy surface I:Numerical investigation [J].International Journal ofHeat and Mas Transfer。 1997。40:1061—1070. [4]Russ G,Beer H.Heat transfer and flow field in a pipe with sinusoidal wavy surface II:Experimental investiga— 维普资讯 http://www.cqvip.com

・72・ 北京化工大学学报 Japan,1985,18:550—555. 2006矩 tion[J].International Journal of Heat and Mass Trans— fer。1997,40:1071—1081. [7] Nishimura T,Kajimoto Y,Kawamura Y.Mas transfer characteristics in a channel with a wavy wall[J].Journal of Chemical Engineering of Japan,1985,19:142—144. [5]Nishimura T,Ohori Y,Kawamura Y.Flow characteris— tics in a channel with symmetric wavy wall for steady lfow[J].Journal of Chemical Engineering of Japan, 1984,17:466—471. [8]徐文熙,徐文灿.粘性流体力学[M].北京:北京理工 大学出版社。1989. [6]Nishimura T,Ohori Y,Kajimoto,et a1.Mass transfer characteristics in a channel with symmetric wavy wall for [9]柴诚敬,张国亮.化工流体流动与传热[M].北京:化 学工业出版社,2000. steady flow[】].Journal of Chemical Engineering of Study of Fluid Flow and Heat Transfer in Corrugated Tubes XIAO Jin—hua QIAN Cai—fu HUANG Zhi—xin AN Wei—zheng (College of Mechanical and Elcterical Engineering,Beijing University of Chemical Technology,Belling 100029,China) Abstract:Fluid flow and heat transfer in a corrugated tube have been numerically and experimentally studied. The influence of different Reynolds Number,wave shape and pipe diameter has been investigated.It was found that in the turbulent region。the average convection film coefficient in a corrugated tube is 1.78~3.10 times that in a straight tube under the same conditions,and the highest heat transfer occurs when Re is about 10000. To facilitate engineering design and application of corrugated tube heat exchangers,a formula for the calculation of the average convection film coefficient of a fluid in a corrugated tube is also given. Key words:corrugated pipe;numerical simulation;experimental verification;heat transfer intensification;for— mula for heat transfer 我刊开通“网上投稿系统" 我刊“网上投稿系统”已开通,校内外作者均可通过该系统在线投稿,并实时查询稿件状态。 作者投稿请先登陆我校期刊中心网站(WWW.journa1.buct.edu.cn),然后点击“作者在线投稿”栏目,提 供作者信息,进行注册。接着按照页面提示,分别提交稿件各相关信息,逐步完成投稿步骤。最后,对提交的 稿件进行确认,完成最终投稿。 该投稿系统还具有在线审稿、远程编辑、在线期刊等功能,不仅方便了校内外作者投稿及专家审稿,缩短 了稿件送审周期,也为编辑部稿件自动管理、即时发布预印本提供了保障,极大地提升了我刊现代化管理水 平。 欢迎校内外读者登陆我刊网站,欢迎作者使用我刊网上投稿系统后提出建议,以便编辑部改进系统,使 之更好地为读者、作者服务。 (北京化工大学学报(自然科学版)>编辑部 

因篇幅问题不能全部显示,请点此查看更多更全内容

Top