汪 威,喻九阳,杨 侠
(武汉工程大学机电工程学院,湖北 武汉 430074)
脉冲流是指按照一定规律如正(余)弦周期变化流动的流体.对于脉冲流速度场分布规律,大量学者都进行了研究,并得到公认的结论.20世纪中前期,Richardson通过对管内往复流动脉冲流的试验发现速度分布存在速度环效应,即最大速度在靠近管壁处而不是在管中心[1];Uchida对此结论进行理论分析,证明速度环效应[2].
但脉冲流是否能强化换热,各专家学者的研究结论都不一致.1990年,Mackly实验证实脉冲流对带内肋片圆管强化效果显著,但对普通圆管直通道没有强化换热的效果[3],俞接成对这两种情况进行数值模拟分析,得到一致的结论[4].
同时,波节管作为一种新型高效换热元件,换热效果是传统直管的1.8~2.5倍[5],并在工业中广泛应用,如若脉冲流在波节管中作用具有强化效果,则换热效率较传统的工业换热器将会有显著的提高[6-7].为此,通过流体数值计算软件fluent对波节管脉冲流的对流换热进行模拟分析,探讨脉冲流是否提高波节管换热效率,并分析其机理.
在数值计算中,采用了二维轴对称方法建立模型,取波节管长度为L=1 820 mm,基管直径(波节管胀扩前的圆管直径)D2=19 mm,波峰处直径D1=25 mm,波长S=20 mm,其中圆弧管段S1=13 mm,直管段S2=7 mm,如图1所示.此外,按照实际模型的结构,在建立波节管几何模型时,在管的进出口端各增加了两个直管段,直径为25 mm,长度为40 mm;在管长方向,除了进出口直管段外,还包含了87个波距的长度.为了便于计算和建模,管壁厚度为零,温度恒定为293 K.流体工作介质为水,入口温度为353 K.
图1 波节管示意图
数值模拟采用结构四边形网格对所建立的几何模型计算区域进行离散,考虑边界层对传热和流动的影响,将壁面附近网格加密,整个模型共划分14.6万个网格单元.计算时选用离散格式,压力与速度的耦合计算采用SIMPLEC方法,压力采用STANDARD,忽略重力的作用,对流项采用二阶迎风格式.流体入口定义为速度输入,并通过UDF输入正弦脉冲流如式(1),出口定义为压力出口,压力为零.
vz=vf[1+Asin(2πnt)]
(1)
图2表示当v=0.2 m/s时,n=10,50 rad/s时波节管流体出口截面平均压力分布.从2图中可以得到,在相同的速度波动条件下,n越大,压力波动也越大.
图2 出口中心处压力波动图
图3 、图4分别是v=0.2 m/s,n=20 r/s第44个波节段的稳态和脉冲流条件下的流函数分布图和速度矢量图.图3中a~h分别表示一个周期下不同相位时的流函数图g为稳态情况下流函数图;图4中a~f表示脉冲流条件下一个周期内不同相位时流体的速度矢量图,i为稳态情况下速度矢量图.通过图形可以表明:稳态流动条件下,此波段内流体流态是固定且无漩涡;而脉冲流情况下,在一个周期内,流体流态不断变化,基本过程如下:首先流体同向流动,速度场无涡旋,当脉冲相位超过180°,即达到负周期时,脉冲流流体在波节靠近流体入口(背速)处产生漩涡,并随着速度的减小逐渐长大,当速度增加时,漩涡又逐渐变小,直至消失.
图3 流函数分布
图4 速度矢量图
图5是v=0.1 m/s,n=20 r/s,x=910处即第44个波节段的稳态和脉冲流条件下的温度分布图.从图5中可以看出,由于脉冲流作用,时均径向温度分布在壁面附近比稳态时的斜率大,降低温度边界层的厚度,从而强化传热.结合图3和图4可得到波节管脉冲流强化传热的基本机理是:由于脉冲流改变了流体的速度分布,产生涡旋,流场紊流增大,从而影响流体径向温度分布.
图5 径向温度分布图
图6给出了波节管在不同速度,不同频率脉冲流作用调节下,波节管换热强化系数变化.从图中可以看出:脉冲流能显著提高波节管的换热效率,在相同的脉冲频率下,换热强化系数随着速度的增加而增加,例如n=20 r/s,速度v=1.0时的R值比速度v=0.1时提高了60.9%;在相同速度条件下,换热强化系数先随着频率的增大而增大,当频率超过一定范围后,R随着频率的增大而减小,以v=1.0为例,n=30 r/s强化传热系数R分别比n=10 r/s,50 r/s的提高了20.5%和8.6%.
图6 强化系数分布图
a.通过fluent数值计算,得到与稳定流动相比脉冲流能显著提高波节管换热效率.
b.脉冲流使得流体在波节管波节处湍流增强并产生漩涡,大大增加流体的扰动,减小边界层厚度,强化波节管换热效果.
c.在强化换热过程中,当脉冲相位在0~180°变化时流场稳定变化,不产生漩涡;当脉冲相位大于180°,流体速度小于平衡值时,波节处流体流场剧烈变化,产生漩涡,并在随着速度的减小,漩涡不断长大,当速度增大时,漩涡逐渐变小,直至消失.
d.换热强化系数在一定速度范围内,随速度的增加而增大;在一定频率范围内,R先随着频率的增大而增大,然后,随着频率的增大而减小.
参考文献:
[1]Richardson E G,Tyler E.The transverse velocity gradient near the mounths of pips in which an atlernating or continuous flow is established [J]. Proceedings of the Physical Society,1929,42:1-15.
[2]Uchida S.The pulsating visous flow superposed on the steady laminar motion of incompressible fluid in a circular pipe [J].ZAMP,1956,7:403-422.
[3]Mamayev V V,NOSOV V S, Syromyatnikov N I.Investigation of heat transfer in pulsed flow of air in pips[J].Heat Transfer-soviet Research,1976,8(3):111-116.
[4]俞接成.脉冲流动和壁面振动传热研究[D]:清华大学,2007.
[5]徐建民,王晓清.波节管管内流动和传热的数值模拟[J].石油化工设备,2008,37(1):4-7.
[6]卢霞,喻九阳.真空电阻凸焊温度场有限元分析[J].武汉工程大学学报,2009,31(5):72-75.
[7]王仕仙,徐建生,卢霞.基于ANSYS的滑动摩擦热结构耦合分析[J]:武汉工程大学学报,2009,31(5):67-70.