毫米级微管径热管的试验与模拟研究

毫米级微管径热管的试验与模拟研究曹甜甜，张志刚（天津城建大学能源与安全工程学院，天津300384）摘要：热管置入式墙体技术是一种新型被动式太阳能利用技术，该技术利用重力热管的传热特性，实现了室内外热量的传输.此技术中内置于墙体的微管径热管的传热能力起着至关重要的作用.为了使微管径热管能最大化利用太阳能，本文采用可视化实验与数值模拟相结合的方法，对其进行了研究.结果表明：蒸发段和冷凝段均相通的毫米级微管径并联集束热管具有良好的均温性；热管充液高度在100~250mm 之间时传热性能最优；在WIHP 所涉及的几何尺寸范围内，8mm 管径热管的综合传热性能比较理想.关键词：毫米级微管径热管；充液高度；数值模拟；气液分布中图分类号：TK172.4文献标志码：A文章编号：2095-719X （2019）01-0015-05Experimental and Simulation Study of Millimeter Micro Diameter Heat PipeCAO Tiantian ，ZHANG Zhigang（School of Energy and Safety Engineering ，TCU ，Tianjin 300384，China ）Abstract ：The technology of the wall implanted with heat pipes is a new passive solar energy utilization technology ,whichrealizes the transmission of indoor and outdoor heat by using the heat transfer characteristics of gravity heat pipe.In this technology ,the heat transfer capability of the micro diameter heat pipe plays a crucial role.In order to maximize the use of solar energy for micro diameter heat pipes ,the combination of visual experiment and numerical simulation is applied to study it in the paper.The results show that the evaporation section and the condensation section that are communicated to the millimeter size parallel micro diameter heat pipe have good temperature uniformity.Heat transfer performance of heat pipe is best between 100~250mm.The comprehensive heat transfer performance of 8mm diameter heat pipe is ideal in the range ofgeometric size involved in WIHP.Key words ：millimeter micro diameter heat pipe ；liquid filling height ；numerical simulation ；gas-liquid distribution 收稿日期：2018-03-12；修订日期：2018-03-27基金项目：天津市自然科学基金（17JCYBJC21400）作者简介：曹甜甜（1992—），女，山西吕梁人，天津城建大学硕士生.通讯作者：张志刚（1961—），男，教授，从事建筑节能与新技术研究.E-mail ：zhangzhigang609@DOI ：10.19479/j.2095-719x.1901015热管是具有较高传热效率的两相热量传输元件，其优势在于利用较小的传热面积和温差就能够达到使用要求.热管已被广泛应用于航空航天、计算机中央处理器CPU 散热、余热回收和冶金工业等诸多工业领域[1-4]，成为解决热量传输等问题的一种有效技术手段.重力热管是利用重力作用使冷凝液回流而进行热量传输的一种热管，它既具有一般热管传热效率高、等温性良好以及环境适应性较强等优点[5]，又由于结构加单、价格低廉被广泛用于热能综合利用和余热回收.中国天津城建大学研究人员2011年首次提出了一种新型被动式太阳能利用技术———热管置入式墙体（the wall implanted with heat pipes ，简称WIHP ），该技术将微管径重力热管与墙体相结合，利用室内外温差进行热量传输，巧妙地解决了外围护结构保温性能与利用天然冷热源的矛盾[6].WIHP 中微管径热管起到传递热量的关键作用，其能否正常工作，直接决定了热管置入式墙体能否正常运行.为了研究微管径热管的传热性能，使之能在WIHP 中最大化利用太阳能，本文采用可视化实验与数值模拟相结合的方法，针对蒸发段和冷凝段均相通的毫米级微管径并联集束热管，分析和探讨了热管内部汽液两相流的流动特征和换热机理.天津城建大学学报Journal o f Tianjin Chengjian University第25卷第1期2019年2月Vol.25No.1Feb.2019天津城建大学学报第25卷第1期1实验介绍和模拟设置1.1实验介绍实验系统主要由加热装置、冷却装置、充液及抽气装置、试验段（热管）和测量装置等组成，如图1所示.加热装置由恒温水箱提供热管蒸发段所需恒温水，冷却装置采用水冷或空冷冷却热管冷凝段，抽气设备采用真空泵，测量装置采集Aglient 34970A 高速数据采集仪采集温度.为进行可视化实验，热管由石英玻璃制作，上部横管端部设充液及抽气口，热管总长500mm ，内径4mm ，壁厚2mm ，根据其形式的对称性，布置了其温度测点.热管形式及测点分布如图2所示，蒸发段与冷凝段长度比为1：1.本实验的目的主要是为了与数值模拟结果作对比，来验证模拟中所建模型的正确性，实验中主要测试热管壁温以及观察热管蒸发段工质流动状态.1.2模拟设置本文建立二维几何模型，采用FLUENT 中VOF 模型进行数值模拟，具体过程分为两部分.（1）GAMBIT 网格划分.为了更好地捕捉边界附近工质的流动和热交换，管壁处采用边界层网格.此外，由于热管内部工质的流动和传热受气液表面张力的影响较大，所以采用非结构化四边形网格[7].局部网格划分如图3所示.（2）FLUENT 设置.FLUENT 设置：求解器形式为Pressure-Base ，选用瞬态模拟，考虑重力，多相流模型选用VOF ，湍流模型为k -着，开启能量方程，添加工质与管壁材料，设置气态为主相，液态为第二项，启用蒸发冷凝模型，设置边界条件，蒸发段为恒壁温，冷凝段为对流换热，为使计算收敛更快，选用PISO 压力-速度耦合算法，格林高斯梯度插值（Green-Gauss CellBased ），设置初始化条件并标记充液高度.2实验研究与数值模拟结果对比分析实验测试和数值模拟中，热管内充液高度（充液高度是指在底部横管充满液体的基础上，立管内所充的液体的高度）为150mm ，压力为0.033MPa ，热管倾角为90°，工质A ，蒸发段加热温度为40℃，冷凝段采用空气冷却，室温为17℃.（1）微管径热管的温度分布.在WIHP 中，重力热管的关键作用是热量的传输，热管工作时的温度分布，体现了整个热管传热性能的好坏，是热管传热性能的重要指标之一.表1为热管各测点的试验测试温度值与模拟温度值，从表中可以看出，试验测试和模拟得到热管上各测点温度的最大差值分别为0.84、0.67K ，二者都能体现热管良好的均温性.由于实测误差和模拟简化造成了测点温度值在数值上有差距，各个测点的试验测试值均较模拟值低2℃左右，模拟误差大约是试验测试值的5.5%左右，模拟与试验值都呈现了相同的规律.图2实验所用热管形式及温度测点分布图1试验系统原理示意图3热管网格划分的局部放大效果恒温水箱过滤器水泵电子称流量计恒温水箱加热水槽冷却水槽水槽热管阀门数据采集仪绝热冷凝段蒸发段98765234116··2019年2月（2）微管径热管的气液分布.为进一步分析热管的传热性能，模拟得到了热管内气液分布规律，并将模拟结果和实验现象进行了对比.在气液分布方面，外界热量通过蒸发段管壁传递给热管内的液态工质，液态工质吸热蒸发，产生饱和蒸汽，饱和蒸汽在管内成蛇形上升，在冷凝段放热冷凝变为液体，液体沿管壁回流至蒸发段，并且回流过程中在管壁上形成一定厚度的液膜.图4为热管试验与模拟现象的对比，从图中可以看出，模拟结果与试验结果基本一致，热管内产生了明显的气液充塞现象.综上分析可知，模拟所得到的热管温度分布、气液分布均合理，且与试验现象和数据基本一致，误差在可接受范围之内，模拟方法合理.3数值模拟分析以下模拟结果均在蒸发段加热温度为40℃，冷凝段空气温度为17℃，管内初始压力为0.033MPa ，工质为A 的条件下得出.3.1不同充液高度条件下微管径热管的传热特性（1）冷凝段温度与充液高度的关系.热管冷凝段的温度分布直接关系到WIHP 传热性能的好坏，图5显示了4mm 热管在不同充液高度条件下热管冷凝段的温度分布情况.其中，温度测点根据热管冷凝段高度等间距设定.从图中可以看出，随着充液高度的增加，冷凝段的温度逐渐升高，尤其是475mm 处测点的温度明显升高，充液高度为0mm 时该测点温度为310.98K ，充液高度为250mm 时为311.72K ，提高了0.74K.在不同充液高度下，随着热管高度的增加，热管壁面温度均逐渐降低.充液高度从0~250mm 所对应的最高测点与最低测点的温差分别为0.51，0.19，0.15，0.09，0.05，0.03K ，可见随着充液高度的增加冷凝段顶端与低端的温度差逐渐降低.当充液高度在0~100mm 之间时，热管冷凝段温度变化较大，且充液高度越小变化越大.当充液高度在100~250mm 之间时，热管冷凝段的温度变化较小，均温性良好.（2）气液分布与充液高度的关系.热管管内是两相流动，当热管传递热量增加时，蒸汽量增加，液膜变厚，因此可通过液膜的厚度分析热管的传热能力.图6为不同充液高度下冷凝段的液膜厚度.图中测点1~10依次由下而上排列，间距为25mm ，其中测点1设置在蒸发段与冷凝段交界处.由于冷凝段顶端横管的存在，立管与横管连接处的液膜厚度会受到横管影响较大，因此只将最高测量点设置在了立管高度为475mm 的位置.从图中可以看出，随着冷凝段高度的降低，液膜厚度逐渐增加，但相同距离液膜厚度增加量逐渐变小，因冷凝液通过重力作用沿管壁回流，当重力大于壁面黏附力时，随着高度的降低冷凝液速度逐渐增大，液膜厚度增加趋势变缓.当充液高度为0mm 时，冷凝段各个测点处液膜均相对较薄.当充液高度为50mm 时，各个测点的液膜厚度均增加0.018mm 左右.当充液高度为100mm 时，液膜厚度进一步增加，但若继续增加充液高度，液膜厚度几乎没有变化.因在稳定状态下，随着充液量的增加，热管蒸发段内液态工质的蒸发量逐渐增加，冷凝段蒸汽的冷凝量也逐渐增加，液膜逐渐变厚，当热管达到传热极限后将无法传递更多热量，因此在冷凝段无法冷凝更多的蒸汽，液膜厚度无法继续增加.充液高度为150mm 和200mm 的液膜厚度曲线有相对较大的波动，因为该充表1热管各测点的试验温度值与模拟温度值测点123456789试验值36.6336.6236.5336.5536.5636.4936.4236.1735.89模拟值38.6038.5838.5738.5638.5838.5638.5438.2737.73（a ）模拟现象（b ）实验现象图4热管模拟与实验现象的对比311.8311.6311.4311.2311.0310.8475300450425400375350325275立管高度/mm0mm100mm 200mm50mm 150mm 250mm图5不同充液高度条件下热管冷凝段的温度分布曹甜甜等：毫米级微管径热管的试验与模拟研究℃17··天津城建大学学报第25卷第1期图7不同充液高度下底部横管及其附近的气液分布液高度下，测点4和6恰好设置在了该充液高度下液膜的波谷和波峰的位置.已有研究结果[8]表明，当充液高度小于干涸临界充液高度时，会发生干涸极限，而当充液高度大于蒸发段长度时，液体在冷凝段易发生过冷现象，甚至会发生周期性爆裂沸腾.因此，为了保持热管稳定有效运行，充液高度应介于干涸临界充液高度和蒸发段高度之间.图7显示了不同充液高度下热管底部横管及其附近的气液分布情况（phase 1代表液体体积分数）.从图中可看出，充液高度为0和50mm 时，并没有液池存在，当充液高度高于100mm 时，随着充液高度的继续增加，底部横管和立管底部的液体体积分数明显增加，液膜更加突出，横管中气泡也较为明显.由此可见，为保证热管连续稳定的运行，充液高度在100~250mm 之间较为合适.0.150.100.050.00475275450425400375350300250立管高度/mm0mm 100mm 200mm 50mm 150mm 250mm图6热管冷凝段液膜与充液高度的关系（a ）0mm（b ）50mm（c ）100mm（d ）150mm（e ）200mm（f ）250mm0.950.900.850.800.750.700.650.600.550.500.450.400.350.300.250.200.150.100.05phase 10.950.900.850.800.750.700.650.600.550.500.450.400.350.300.250.200.150.100.05phase 10.950.900.850.800.750.700.650.600.550.500.450.400.350.300.250.200.150.100.05phase 10.950.900.850.800.750.700.650.600.550.500.450.400.350.300.250.200.150.100.05phase 10.950.900.850.800.750.700.650.600.550.500.450.400.350.300.250.200.150.100.05phase 10.950.900.850.800.750.700.650.600.550.500.450.400.350.300.250.200.150.100.05phase 13.2不同管径条件下微管径热管的传热特性由于受到墙体厚度及内外表面结构的限制，WIHP 所用热管必须是毫米级的微管径热管，而对于毫米级的热管，热管管径的大小会对整个热管的传热性能有较大的影响.为选择合适的管径，本文模拟研究了蒸发段长度比为0.5的情况下，不同管径（内径分别为2，4，6，8mm ）热管的传热效果.（1）冷凝段温度与管径的关系.图8显示了不同管径热管冷凝段温度的分布情况，可以看出，不同管径热管相同高度测点最大温差不超过0.14K ，同一根32518··2019年2月热管冷凝段最大温差不超过0.2K，四种管径的热管都体现了良好的均温性，但就此四种管径的热管而言，随着管径的增加，热管均温性逐渐下降，2mm热管均温性最好，8mm热管均温性最差.随着热管管径的增大，管内脉动现象逐渐变弱，2mm热管具有强烈的脉动现象，而8mm热管中几乎不存在脉动现象.强烈的脉动现象可以使管内蒸汽更易到达热管顶端，同时增加管内的扰动，破坏管壁处的液膜，使得蒸汽可以直接与管壁接触，增强换热效果，因此，脉动现象更有利于热管的均温性.值得注意的是，当管径增加一个数量级，管内的冷凝段温度随管径的变化规律恰好与图8中所显示的相反，即随着管径的增加，冷凝段的均温性越好，冷凝量也增大[9].这是由于对于大管径的热管，管内没有脉动现象，管径越大，蒸汽越容易上升到热管顶部，且冷凝段与外界接触面积增大，更有利于热管与环境换热，壁面温度更易均匀.（2）气液分布与管径的关系.图9为不同管径热管冷凝段液膜厚度随冷凝段高度的变化规律.从图中可以看出，随着热管管径的增加，液膜厚度在整体上逐渐增加.2mm管径热管的液膜厚度随着冷凝段高度的增加波动较大，各测点间液膜厚度相差较大.在300mm处的冷凝段液膜厚度最大为0.221mm，而处于热管冷凝段下部250，275mm处的液膜厚度较薄，分别为0.089，0.075mm，与最大值相差约为1.5mm. 4mm热管冷凝段液膜厚度也存在波动，但波动幅度较小，在275mm处液膜最厚为0.1023mm.6mm热管冷凝段液膜厚度的波动进一步减小，且液膜厚度自上而下逐渐增厚.8mm热管冷凝段液膜厚度基本不存在波动，液膜从上而下几乎呈线性规律逐渐增厚.不同管径热管冷凝段的上部液膜厚度变化趋势基本一致，越靠近热管顶端液膜越薄.总体来说，对于毫米级热管而言，管径越小，冷凝段液膜厚度的波动性越明显，随着热管管径的增大，冷凝段液膜厚度逐渐增大，8mm热管冷凝段底部液膜最厚为0.126mm.虽然在稳定状态下，蒸发段吸收的热量等于冷凝段释放的热量，但是，热管内部狭小的反应空间内，不同时刻的气液分布不同，且没有规律.热管内工质流动的平稳性与壁温的稳定性直接相关，当壁温稳定时，热管内工质的流动平稳性较好，反之，则热管内工质流动会出现波动，导致热管传热性能下降[10].由于随时间的推进，8mm管径的热管内气液分布状态变化较小，液膜波动较小，因此热管温度较为稳定.此外，随着热管内径增大，管内气液充塞现象逐渐减弱，使得横管中堆积的气体减少，液态工质更易回流至横管，立管中液注高度降低，管壁处液膜的面积增加，由于液膜传热系数远大于液池，因此将更有利于热管的传热.4结论本文通过实验验证了采用FLUENT软件中VOF 模型对适用于WIHP的毫米级微管径并联集数热管进行数值模拟的合理性，并分析了充液高度和管径大小对其传热特性的影响，得出如下结论.（1）热管中的不同立管在冷凝段同高度位置的温度基本相同，随着高度的增加，管壁的温度呈降低趋势，整体来看蒸发段和冷凝段均相通的毫米级微管径并联集束热管具有良好的均温性.（2）随着充液高度的增加，冷凝段壁面温度均匀性逐渐提高，毫米级微管径并联集束热管充液高度在100耀250mm之间时能够确保有效稳定运行.（3）WIHP中热管的管径以8mm较为适宜.对于毫米级热管而言，虽然管径的增大会减弱脉动现象，降低热管冷凝段的均温性，但最大温差也仅为0.2K；而增大热管管径的同时，热管温度更为稳定，且液膜图8不同管径热管冷凝段温度的分布情况311.80311.75311.70311.65311.60311.55311.50475300450425400375350325275立管高度/mm2mm4mm6mm8mm0.250.200.150.100.050.00475300450425400375350325275立管高度/mm2mm4mm6mm8mm250图9冷凝段液膜厚度与管径的关系（下转第25页）曹甜甜等：毫米级微管径热管的试验与模拟研究19··2019年2月面积增加、换热增强，因此综合考虑，在WIHP中应尽量选择较大管径的热管.参考文献：[1]FRANCO A，VACCARO M.On the use of heat pipe principle for the exploitation of medium Low temperature geothermal resources[J].Applied Thermal Engineering，2013，59（1）：189-199.[2]王斌斌，仇性启.热管及其换热器在烟气余热回收中的应用[J].工业加热，2006，5（35）：37-41.[3]ZHANG H，ZHUANG J.Research，development and industrial application of heat pipe technology in China[J].Applied Thermal Engineering，2003，75（4）：295-306.[4]El-BAKY M A，Mohamed M M.Heat pipe heat exchanger for heat recovery in air conditioning[J].Applied Thermal Engi-neering，2007，27：795-801.[5]庄骏，张红.热管技术及其工程应用[M].北京：化学工业出版社，2000.[6]ZHANG Z G，SUN Z J，DUAN C X.A new type of passive solar energy utilization technology：the wall implanted with heat pipes[J].Energy and Buildings，2014，84：111-116. [7]王海彦，刘永刚.ANSYS FLUENT15.0流体数值计算方法与实例[M].北京：中国铁道出版，2015.[8]IMURA H，SSASAGUCHI K，KOZAI H.Critical heat flux in a closed two phase thermosiphon[J].International Journal of Heat and Mass Transfer，1983，26（8）：1181-1188.[9]GHORAI S，Nigam K D P.CFD modeling of flow profiles and interfacial phenomena in two-phase flow in pipes[J].Chemical Engineering and Processing Journal，2006，45（1）：55-65.[10]曹丽召.重力热管流动与传热特性的数值模拟[D].北京：中国石油大学，2009.bridge cables at moderate flow turbulence[J].Matec Web of Conferences，2017（107）：1-8.[9]何宇虹.大跨径索桥覆冰形式及气动力特性研究[D].西安：长安大学，2015.[10]刘春城，刘佼.输电线路导线覆冰机理及雨凇覆冰模型[J].高电压技术，2011，37（1）：241-248.[11]张锐.考虑极端气候的季冻区道路结构温度场分布规律研究[D].哈尔滨：哈尔滨工业大学，2012.[12]MAKKONEN L.Modeling of ice accretion on wires[J].Journal of Climate and Applied Meteorology，1984，23（6）：929-939.[13]李新民，朱宽军，刘彬.典型覆冰导线空气动力学特性数值和试验模拟[J].高电压技术，2014，40（2）：427-433. [14]蒋兴良，申强.环境参数对导线覆冰厚度影响的试验分析[J].高电压技术，2010，36（5）：1096-1100.（上接第19页）冯一凡等：斜拉索覆冰形态试验研究25··。