《武汉工程大学学报》 2014年03期
58-62
出版日期：2014-03-31
ISSN:1674-2869
CN:42-1779/TQ
内置分段式弹簧换热管内流场的数值模拟
0引言管内插入螺旋弹簧作为一种无源技术已被一些学者进行了研究，并证明了该技术能够很好的实现强化传热，但是弹簧结构对管内流体压降与阻力的影响尚不清楚，该研究尚不多见.本文决定对螺旋弹簧进行分段式处理后，通过数值模拟的方法来系统研究分段式弹簧强化传热特性.已有很多学者对内置螺旋弹簧换热管强化传热特性进行数值模拟［13］，多数采用FLUENT软件对流体流动传热的速度场、压力场、温度场进行数值模拟，但是对分段式弹簧换热管数值模拟方法的研究并不多.数值模拟更能给研究人员带来直观的分析数据，且运用起来方便，高效，深受学者的青睐.运用FLUENT软件对不同内置分段式弹簧换热管的速度场进行数值模拟，有利于观察流体流动状态和研究内置分段式弹簧换热管管内流体的传热特性，得到不同长度的分段式弹簧对流动传热性能的影响规律.1建立模型和边界条件1.1计算模型拟采用4种不同尺寸的内置分段式弹簧换热管与光管流动与传热进行三维数值模拟，该模拟选取长度500 mm的管长作为模拟区域，换热管直径为20 mm，弹簧丝径为1 mm，弹簧节距为3 mm，弹簧圈径为16 mm.分段式弹簧的尺寸如表1所示.表1各段弹簧长度Table 1The length of each segment of the springmm序号两端分别内置弹簧长度150210031504200由于计算模型中弹簧呈螺旋结构，并且横截面很小，划分网格有一定的困难，经过反复的尝试，最终选择适应性较好的Tet/Hybird（四面体）类型，该类型网格包含四面体网格单元，但是在合适的位置也包含六面体、锥体，划分方式选用TGrid（主要采用四面体单元）方式，为提高结果的精确度，对弹簧壁面的网格加密.划分网格的局部模型如图1所示. 图1划分网格局部模型Fig.1Mesh generation of local model1.2边界条件流体入口定义为速度入口边界，速度从0.3 m/s增加到1.1 m/s，入口流体的温度为333 K，换热管的壁面定义成无滑移的静止的刚性壁面，壁面温度恒定为293 K，流体出口定义为压力出口边界.采用标准κε湍流模型，压力与速度的耦合采用SIMPLE（基本）方式，对流项采用一阶迎风离散格式，定义收敛条件为残差的绝对值小于1×10-6.2模拟结果2.1流动特性在内置分段式弹簧换热管内，弹簧呈螺旋结构，对流体起到一个引导作用，使得流体随着产生螺旋运动；弹簧在流体的冲击作用下会发生振动，流体受到弹簧的振动作用，轴向速度，径向速度以及弹簧周围流体的切向速度会有所增大，因此弹簧就会产生离心力，对近壁区域流体产生扰动，减薄边界层，从而改善传热效果；弹簧的螺旋结构增加了流体的湍流强度，充分扰动边界层，进一步提高传热效率［46］. 2.24种不同分段弹簧长度下在x=30 mm截面处的径向速度场当流体入口速度为0.3 m/s时，内置不同尺寸分段式弹簧换热管管内x=30 mm截面处的径向速度分布如图2至图6所示.第3期徐建民，等：内置分段式弹簧换热管内流场的数值模拟武汉工程大学学报第36卷 图2光管速度分布图Fig.2Velocity profile of light pipe 图31#管速度分布图Fig.3Velocity profile of 1# pipe 图42#管速度分布图Fig.4Velocity profile of 2# pipe 图53#管速度分布图Fig.5Velocity profile of 3# pipe 图64#管速度分布图Fig.6Velocity profile of 4# pipe从图2至图6可以看出，换热管内置分段式弹簧可以明显提高流体的径向速度，并在换热管壁面区域加强流体扰动，破坏边界层，边界层在弹簧的振动作用下产生分离，使得边界层变薄，从而达到强化传热的效果.同时插入分段弹簧的长度对径向速度影响不大，因此综合流体流动阻力考虑，为了减小弹簧对流体的阻碍，采用长度较小的分段式弹簧（换热管的两端插入弹簧，中间湍流区域没有弹簧）可以更好地提高强化传热效果.2.3弹簧丝径的影响入口速度为0.3 m/s时，丝径分别取1.0 mm、1.4 mm、1.8 mm，内置不同丝径的分段式弹簧换热管管内x=30 mm截面处的径向速度分布如图7至图9所示. 图7丝径为1.0 mm时x=30 mm截面处速度云图Fig.7Speed of cloud images(Wire diameter is 1.0 mm) 图8丝径为1.4 mm时x=30 mm截面处速度云图Fig.8Speed of cloud images(Wire diameter is 1.4 mm) 图9丝径为1.8 mm时x=30 mm截面处速度云图Fig.9Speed of cloud images(Wire diameter is 1.8 mm) 从图7至图9可以看出随着丝径的增大，弹簧的引流作用得到增强，边界层的速度增大，说明起到了扰动作用，并且加强了弹簧周围流体的运动，弹簧又使得流体产生了明显的螺旋运动，进而迫使边界层分离，边界层的厚度随之减小，从而增强强化传热效果，但同时增大了阻力压降.2.4弹簧圈径的影响入口速度为0.3 m/s时，圈径分别取8 mm，12 mm，16 mm，内置不同圈径的分段式弹簧换热管管内x=30 mm截面处的径向速度分布如图10至图12所示. 图10圈径为8 mm时x=30 mm截面处速度云图Fig.10Speed of cloud images(circle diameter is 8 mm) 图11圈径为12 mm时x=30 mm截面处速度云图Fig.11Speed of cloud images(circle diameter is 12 mm) 图12圈径为16 mm时x=30 mm截面处速度云图Fig.12Speed of cloud images(circle diameter is 16 mm)从图10至图12可以看出随着弹簧圈径的增大，流体的速度也会增大，因为圈径大的弹簧离壁面区域近，流体的冲击使弹簧发生振动后，边界层的流体会产生扰动，加上流体在弹簧作用下产生的螺旋运动，使得边界层的流体运动更剧烈，可以较快的从层流阶段过渡到紊流阶段，减薄了边界层的层流底层，降低壁面热阻，进而提高了强化传热效率.3计算结果与数据分析图13可以反映出每段弹簧的长度对强化传热有一定的影响，并与光管进行了比较，内置分段式弹簧对强化传热起到了促进作用.随着雷诺数Re的增加，努赛尔数Nu随之增大，雷诺数Re相同时，随着插入弹簧的长度增加，努赛尔数Nu相应的增加，从而增大强化传热的效果. 图13Nu数随Re的变化曲线Fig.13Nusselt number and Reynolds number variation graph 图14可以反映出，内置分段式弹簧换热管相对于光管，阻力系数会增大，但是随着Re数的增大，阻力系数逐渐减小，且1#管阻力系数最小，说明在两端插入的弹簧长度不宜过长，因为弹簧会增加流体流动阻力，在换热管完全紊流区，可以达到良好的传热效果，因此，采用分段式弹簧可以有效减小阻力，使流体呈现扰动、稳定流、扰动的状态，进一步提高传热效率. 图14阻力系数随Re的变化曲线Fig.14Resistance coefficient and Reynolds number variation graph 4结语a.在内置分段式弹簧换热管的入口与出口端由于存在弹簧的振动，使得流体边界层破坏，加强了流体的扰动，使得两端管内侧边界层变薄，从而提高强化传热的效果.b.分段式弹簧的存在，可以明显起到强化传热的作用.在相同的Re数下，与光管相比，Nu数提高了2~4倍，但是阻力系数也会相应的增加，总体来说还是可以提高强化传热效率.c.内置分段式弹簧强化传热技术适用于低流速下的传热技术应用，并且两端弹簧尽量选择短些，有利于减小压降，兼并弹簧的振动作用，从而实现强化传热.致谢衷心感谢武汉工程大学机电工程学院对本研究工作的支持与帮助.