微圆柱阵微通道内沸腾换热及气泡运动特性研究_左玉清.pdf

第 39 卷第 1 期2023 年 2月上海电力大学学报JournalofShanghaiUniversityofElectricPowerVol39，No 1Feb2023DOI:10 3969/j issn 2096 8299 2023 01 009收稿日期:2022-02-11通信作者简介:杨艳霞(1985)，女，博士，副教授。主要研究方向为流动、传热传质及数值计算方法。E-mail:yangyanxia tyut edu cn。基金项目:国家自然科学基金(5197632);山西省基础研究计划项目(202203021211167)。微圆柱阵微通道内沸腾换热及气泡运动特性研究左玉清，杨艳霞，赵贯甲(太原理工大学 电气与动力工程学院，山西 太原030024)摘要:基于流体体积(VOF)模型对恒定热流密度下的微圆柱阵通道内的流动沸腾换热特性展开数值研究，并利用几何重构法捕捉气液两相界面迁移变化，研究了流道内流速的分布和微圆柱高度对气泡分布、沸腾换热性能及沸腾不稳定性的影响。结果表明:微圆柱的存在增加了气泡核化点密度;气泡运动增加了工质与加热壁面接触的可能性;微圆柱高度较低时，易发生气塞现象，加热面散热不均导致局部异常过热，换热性能下降;微圆柱阵通过阻碍气泡反向流动来降低沸腾不稳定性，但微圆柱高度过高会造成换热系统振荡，影响传热可靠性。关键词:微圆柱阵;微通道;沸腾换热;气液两相流;气泡中图分类号:TK124文献标志码:A文章编号:2096 8299(2023)01 0054 07Study on Boiling Heat Transfer and Bubble Motion inMicro-channel of Micro-rib ArrayZUO Yuqing，YANG Yanxia，ZHAO Guanjia(College of Electrical and Power Engineering，Taiyuan University of Technology，Taiyuan030024，China)Abstract:The heat transfer characteristics of flow boiling in the cylindrical micro-fin channel withconstant heat flux at the bottom were numerically studied，which was based on the VOF flow modelcoupled with boiling heat and mass transfer model To capture the movement of vapor-liquid interfaceby the geo-reconstruction method，this article studied the partial velocity variation of two-phase inter-face Simultaneously，it analyzed the effect of rib height on the morphology of bubbles，flow boiling heattransfer characteristics and boiling instability in the micro-channel The results show that the presence ofmicro-columns increases the density of bubble nucleation point The bubbles movement increases thepossibility of contact between working medium and heating wall When the height of rib is low，the cap-illary force generated by the secondary channel is weak and air bubbles are not easy to escape This canlead to locally abnormal overheating and cause the heat transfer to deteriorate eversing flow of thebubbles is blocked by micro-fin array，and the boiling instability is reduced However，over-hight micro-fin can lead to oscillate heat transfer system，thus affecting the reliability of heat transferKey words:micro-rib structure;micro-channel;boiling heat transfer;vapor-liquid two-phaseflow;air bubbles左玉清，等:微肋阵微通道内沸腾换热及气泡运动特性研究微通道散热器具有可靠运行、微型化和高能效的优势1，在农业、航空航天、太阳能等领域得到广泛应用。随着技术的不断发展，急剧增加的热负荷使设备的安全性和寿命受到挑战，因此迫切需要高效的散热手段。TUCKEMAN D B 等人2 在 30 多年前提出的微通道沸腾换热可以有效解决该问题，即液体通过吸收大量热量生成气泡从而有效增强换热性能。目前，已有很多学者在微通道沸腾两相流领域展开了研究。文献 3研究发现，相比常规尺寸管，过冷液体在微通道中流动沸腾的传热特性和冷却性能有较大提高。文献 4-5 在微通道流动沸腾实验中观察到气泡的多种流型，同时发现流型转变会影响沸腾换热性能和传热不稳定性。段芮等人6 通过实验观察到对于翅片较密的沸腾管，排液能力较差会导致换热系数不佳。LIU H等人7 利用高速摄像捕捉气泡运动轨迹，发现较高热通量会导致气泡滑动速度更高且气泡间更易合并。KISHNAMUTHY S 等人8 通过微针肋通道的沸腾换热实验发现，与光滑平面相比，微针肋表面的流动沸腾换热显著增强，且两相传热系数与质量流量、热流密度均无关。HUA J Y 等人9研究了微通道沸腾的传热特性与气泡动力学特性的内在联系，从气泡出现的数量和顺序等角度分析三角形针翅具有较高传热系数的原因。LEE W 等人10通过沸腾模拟计算发现，锥形空腔表面微结构加速了气泡的生长和脱离进程，有效提高了沸腾换热性能。LIN Y H 等人11模拟研究了平行微通道 2 个单气泡的生长，发现增加质量通量能够抑制气泡流动逆转，降低沸腾不稳定性。PAJAPATI Y K 等人12对分段翅片微通道沸腾过程进行了数值研究，结果表明，二次通道的存在使气泡的生长速度更快，换热系数更高。陈宏霞等人13采用流体体积(Volume of Fluid，VOF)模型模拟研究了微肋表面单个气泡的动力特性，表明微肋结构对于气泡的生长和脱离具有明显的促进作用。KIM J S 等人14研究了水平管沸腾生成的气泡形核位置分布和气泡行为，采用新型热分配模型更加真实地反映了气泡滑动现象。LI X L 等人15模拟研究了圆柱微通道内液滴的两相流传热，发现随着滑移长度的增加，液滴两相流中的传热效果也会增强。综上所述，对于微通道沸腾换热的研究主要基于整体宏观方面，而由于微结构的复杂性，目前对局部微观的沸腾换热研究尚少，因此本文对圆形微圆柱阵通道内的流动沸腾过程进行数值模拟，通过捕捉气液相界面的运动行为，研究分析局部气泡运动特性，分析微圆柱高度对微通道内气泡分布、沸腾换热特性及沸腾不稳定性的影响，以期为散热器的结构优化提供思路。1物理模型1 1几何模型微通道的几何模型如图 1 所示。水平微通道内交叉排布7 排10 列微圆柱阵。微通道基底板和微圆柱的材质均为铜，沸腾工质为水。流体流入微通道，通过沸腾换热将加热面热量带走。图 1 中，k为基底板厚度，mm;l 为微圆柱高度，mm;L、W、H分别为微通道的长度、宽度、高度，mm;d 为微圆柱直径，mm;SL为微圆柱横间距，mm;ST为微圆柱纵间距，mm。图 1几何模型55上海电力大学学报2022 年本文取微通道的 L W H 为 16 0 mm 4 2 mm 0 3 mm，h1为 0 05 mm，d 为 0 6 mm，SL=ST=1 mm，4 组 微 圆 柱 高 度 l 分 别 为0 05 mm、0 10 mm、0 15 mm、0 20 mm。1 2控制方程本文做出以下假设:水和水蒸气为不可压缩流体;水和水蒸气在饱和压力及温度下的物性参数为常数;通道尺寸过小，不考虑重力影响。VOF 模型通过求解一组动量方程，并跟踪整个计算域中气液两相的体积分数来对两相进行建模，能够更加精确地追踪气液相界面。VOF 模型的建立基于以下条件:2 种或多种流体(或相)不互溶，则液相和气相体积分数的总和应当保持恒定，即l+v=1(1)式中:l、v 液相、气相的体积分数。基于 VOF 模型的液相和气相的连续性方程为lt+(vl)=Sl(2)vt+(vv)=Sv(3)式中:t 时间，s;梯度算符;v 气液两相的平均速度，m/s;S 液相转移到气相的质量源项，kg;l、v 液体密度和气体密度，kg/m3。动量方程为(v)t+(vv)=p+(v+vT)+g+Fvol(4)式中:混合相的密度，kg/m3;p 气液相界面的局部压力，Pa;混合相动力黏度，Pas;g 重力加速度;Fvol 体积力，N。气泡与微通道尺寸处于同一量级，因此气泡所受表面张力16 对于沸腾流动是非常重要的影响因素。将表面张力作为体积力加入动量方程中，其形式为Fvol=ll0 5(l+v)(5)式中:表面张力系数，N/m;l 两相界面的液相曲率。由于温度变化幅度较小，表面张力采用常数，数值取 0 059。能量方程为(E)t+v(E+p)=(T)+Q(6)式中:E 气液两相流体的内能，J/kg;气液两相的导热系数，W/(mK);T 气液两相共用的温度，K;Q 通过界面的体积能量源项，J/kg。沸腾换热需要考虑质量和能量在气液两相界面处的传递。质量传递取决于局部温度和饱和温度之间的温差，如果局部温度高于饱和温度，就会发生沸腾，液相质量转移到气相，否则会发生凝结。因此，质量源项17 为S=rlll(Tl Tsat)Tsat，Tl Tsatrvvv(Tl Tsat)Tsat，Tl Tsat(7)式中:rl、rv 液相和气相的质量传递强度因子，取 1 s1;Tl、Tsat 液体的温度和饱和温度，K。相应的能量源项17-20 为Q=hlvS(8)式中:hlv 液相传递到气相的潜热，J/kg。模拟计算中，质量源项和能量源项通过用户自定义函数与 VOF 模型进行耦合。1 3网格划分及网格无关性验证使用 ICEM 软件对微圆柱阵微通道进行结构化网格划分。为提高计算效率，对每一个微圆柱都采用 O-block 方式生成，网格质量检测值达到0 6 以上，以保证计算精度。选取 t=40 ms 时的沸腾换热系数进行网格无关性验证，数据如表 1所示。根据表 1 数据计算可得，当网格数大于55 万时，偏差在 2%以内，因此本文采用的网格数量为 55 万。表 1网格无关性验证网格数/104沸腾换热系数/W(m2K)1网格数/104沸腾换热系数/W(m2K)11315 0505525 5132520 0527926 0013922 46565左玉清，等:微肋阵微通道内沸腾换热及气泡运动特性研究1 4边界条件微通道内沸腾换热的边界条件设置如下:入口边界，速度 v=04 m/s，温度 T1=370 K;出口边界，压力出口处的回流温度 T2=373 15 K;流体区边界，上边界及四周边界均为绝热边界;固体区边界，加热底面为恒定热流密度 0=300 kW/m2;基底板侧壁面均为绝热壁面。求解方法采用 PISO 算法和几何重构法，模拟计算的时间步长取 1 0 107s。2数值模拟结果与分析2 1模型验证为了验证模型的可行性，对三维矩形平直微通道内的流动沸腾换热进行模拟，与文献 21 中参数设置一致，计算 3 组入口流速下的流动沸腾换热系数，并与文献 21 中的模拟值和实验值进行比较，结果如表 2 所示。表 2本文计算结果及与文献 21 的数据对比入口速度/(ms1)沸腾换热系数/W(m2K)1计算值模拟值实验值误差/%模拟误差 实验误差0 1087 35884 92778 1462 91180 1586 46285 18978 6131 5990 2088 20585 83179 2752 7113其中，误差的计算公式为:误差=(模拟结果 文献结果)/文献结果。由表 2可以得出，本文计算结果与文献 21模拟值的平均误差为2 4%，与实验值平均误差为 11 0%，吻合程度良好，证明本文所采用的 VOF 模型及传热传质模型是合理可靠的。2 2模拟结果与分析2 2 1气泡动力特性分析l=0 1 mm、t=40 ms 时的微通道内气泡分布云图如图 2 所示，蓝色部分为气泡。由图 2 可以看出，通道入口段气泡较少，这是因为入口处流体不断冲刷加热壁面使气泡脱离;小气泡主要集中在微圆柱周围，这是由于微圆柱的存在提供了大量的气化核心，有助于气泡的生成;通道出口段多为大气泡，这是因为从加热面上图 2微通道内气泡分布云图脱离的气泡随着流动方向(沿 x 正方向)汇集在通道尾端，小气泡聚并成气泡块，堵塞通道。l=0 05 mm、t=39 5 ms 时，微通道内流体的速度矢量图如图 3 所示。由图 3 可以看出，3 13 mm处布置有微圆柱阵，流体阻力作用增大，微圆柱周围流体速度骤降至 0 2 m/s 以下，整体流速偏低;在 x=13 mm 处流体进入平直通道，速度提高至 0 45 m/s 以上。图 3(b)和图 3(c)是图 3(a)内方框区域 A 的放大图，可以看出气液两相界面处流体速度远大于周围，微圆柱间气泡发生聚并现象。t=39 5 ms 时 2 个气泡分别生长在相邻的微圆柱上，t=40 ms 时气泡发生合并，且速度沿着合并的方向加大。这是由于随着气泡的生长、聚并，其体积不断膨胀，推动流体流动。同时，由于微通道尺寸过小，气泡与微通道处于相同量级，对流体的扰动更加剧烈。这个现象说明，气泡运动是影响通道内部流速变化的主要因素之一。不同微圆柱高度下微通道内的沸腾现象图片如图 4 所示。4 组微通道皆出现了许多贯通整个通道高度方向的塞状气泡。图 4(a)是流体区上壁面截面处的气泡分布，可以看出，l=0 15 mm比 l=0 10 mm 的微通道内塞状气泡明显减少，l=0 05 mm 的微通道内甚至出现了拉伸大气泡。由图 4(b)可以看出，微圆柱高度的增高推后了塞状气泡出现的位置，且微圆柱高度越高，大气泡数量越少。这是因为较高的微圆柱之间形成明显的次级通道，次级通道内产生的毛细力加快了流体的汇入，气泡跃离频率增加，相邻气泡不易合并，通道可以及时排出气泡，不易造成气塞现象。75上海电力大学学报2022 年图 3通道内气液两相流体的速度场图 4不同微圆柱高度微通道内的沸腾现象85左玉清，等:微肋阵微通道内沸腾换热及气泡运动特性研究2 2 2微圆柱高度对沸腾换热性能的影响t=40 ms 时努塞尔数(Nu)沿流动方向的变化曲线如图 5 所示。图 5加热壁面 Nu 值沿流动方向变化曲线图5 中，0 3 mm 处为平直通道，Nu 值较小;3 13 mm 布置有微圆柱阵，Nu 值明显增大，并呈现不同程度的波动;13 16 mm 进入平直通道，Nu值大幅减小。结合图2 可知，由于微圆柱的存在导致生成气泡增多，次级通道产生的毛细力促进流体汇入，加快了气泡脱离，Nu 值增大，换热能力增强。同时，l=005 mm 条件下，x=3 mm 比 x=2 mm 处的 Nu 值增幅较小，仅增大了 10%;而 l=0 20 mm时的 Nu 值则增大了 58%;在微圆柱位置处，l=020 mm 时的Nu 值整体高于其他通道。这是因为微圆柱高度越高，次级通道产生的毛细力作用越强，加快了周围流体的汇入，推动了气泡脱离，提高了液体与加热面接触的可能性，增强了换热效率。在不同微圆柱高度下，沸腾换热系数随时间的变化曲线如图 6 所示。图 6沸腾换热系数随时间的变化曲线由图 6 可以看出，在 t=30 ms 之后流体完全流过通道，沸腾换热系数大致趋于稳定，随着微圆柱高度的增高，沸腾换热系数逐渐增大。沸腾换热系数 h 的计算公式为h=qw-iTw-i Tl(9)式中:qw-i 通过固液界面的面积加权平均热通量，kW/m2;Tw-i 固液界面的面积加权平均温度，K;Tl 流体的质量加权温度，K。平均含气率随微圆柱高度的变化曲线如图 7所示。其中，平均含气率取微通道在 40 50 ms间含气率的平均值。随着微圆柱高度的增高，通道含气率减小。结合图 3 可知，由于沸腾过程中气液两相界面的运动增强了流体的扰动，各处工质实现了与加热面的接触传热，充分吸收热量，提高了换热效率，沸腾换热系数增大，从而对加热壁面起到显著的降温效果;微圆柱高度降低，平均含气率增大，结合图 4 可知，塞状气泡增多，堵塞微通道，使壁面出现干涸现象，造成局部异常过热，以过热高温点为中心向周围壁面扩散热量，进而影响整个微通道壁面的冷却，沸腾换热系数减小，导致传热效果减弱。图 7平均含气率随微圆柱高度的变化曲线2 2 3沸腾不稳定性分析微通道内沸腾产生的气泡易发生反向流动，会对整个微通道的换热系统造成振荡。流动沸腾不稳定会直接导致提前达到临界热流密度，致使传热效果减弱，而微圆柱阵的存在可以有效抑制气泡的反向流动，提高沸腾的稳定性 22-23。一般将通道进出口的气液两相压降作为衡量流动是否稳定的标准。不同微圆柱高度下，进出口压降随时间的变化95上海电力大学学报2022 年趋势如图8 所示。在 t=30 ms 之后流体完全流入微通道，压降保持在一定范围内的振荡。微圆柱高度进出口压降 p 的波动范围分别为158 181 kPa(l=0 05 mm)，1 72 2 07 kPa(l=010 mm)、227 266 kPa(l=0 15 mm)、3 14 355 kPa(l=020 mm)。这表明随着微圆柱高度的增高，压降也相应增大，流动沸腾不稳定性升高，换热系统可靠性降低，影响了换热器的安全运行。图 8进出口流动压降随时间的变化曲线3结论本文对含有微圆柱阵的微通道内流动沸腾的气泡分布规律及流道内速度变化进行了研究，并分析了微圆柱高度对沸腾换热特性和不稳定性的影响。主要结论如下。(1)气液两相界面流体速度大于周围，次级通道的存在促进了流体的汇入，加速了气泡脱离，较少发生气塞现象。(2)随着微圆柱高度的增高，平均含气率减小，换热系数增大，换热性能增强。(3)微圆柱越高，进出口流动压降波幅越大，不稳定性提高，换热系统振荡明显，可靠性下降。参考文献:1过增元 当前国际传热界的热点 微电子器件的冷却J 中国科学基金，1988，2(7):20-252TUCKEMAN D B，PEASE F W Corrections to high per-formance heat sinking for vlsiJ IEEE Electron Device Let-ters，1981，2(8):213-2133PENG X F，WANG B X Cooling augmentation with micro-channeled structuresJ Journal of Thermal Science，1993，2(2):98-1104FAYYADH E M，MAHMOUD M M Flow boiling heattransfer of 134a in multi-micro-channelsJ InternationalJournal of Heat and Mass Transfer，2017(3):422-4365HUH C，KIM J，KIM M H Flow pattern transition in a stabil-ity during flow boiling in a single micro-channel J Interna-tional Journal of Heat and Mass Transfer，2007，50(5):1049-10606段芮，马虎根，刘君，等 表面增强型强化传热管换热特性实验研究J 上海电力学院学报，2006，22(1):84-887LIU H，CHEN D，QIN J，et al Visual experimental study onbubble behavior from nucleation site in vertical rectangularchannelJ Applied Thermal Engineering，2021，201:20-31(下转第 90 页)06上海电力大学学报2022 年Conference on Genetic and Evolutionary Computing IEEEComputer Society，Shenzhen，China，2010:142-1457史欢，黄晓胜，曾文君，等 磁控电抗器的动态无功补偿装置J 中国电力，2011，44(1):32-368曹积欣，王冠军，李建民 风电并网电力系统无功补偿动态性能研究 J 电力电容器与无功补偿，2012，33(3):16-229王琳基，王苏潭，韩群勇，等 谐波状况下电容柜损坏原因分析及采取措施 J 电气应用，2008(8):76-78 10 陈浩 电容补偿柜在谐波状况下优化运行若干问题的探讨D 福州;福建农林大学，2006 11 林因 基于 ADPSS 仿真平台的 220 kV 变电站静止无功补偿装置的建模及动态特性研究 J 上海电力学院学报，2018，34(2):173-179 12 朱洪波，宋颖巍，刘岩，等 风电场无功补偿计算方法与容量配置的研究J 现代电力，2011，28(6):73-76 13 JAISWAL S P，SHIVASTAVA V，PALWALIA D K Im-pact of semiconductor devices on voltage stability of distribu-tion system C/Proceedings of the National Conference onMaterials and Devices(NCMD)Greater Noida，Inidia，2018:581-589 14 杨狄赛，童向亚 基于单片机的发电机功率因数测量系统设计 J 赤峰学院学报(自然科学版)，2015，31(2):40-42 15 刘懿，王大志 无功功率补偿技术 J 科技创新导报，2008(34):58 16 纵榜峰 基于集成运算放大器的电压比较器实验综述报告J 科学技术创新，2020(4):19-20(责任编辑谢冉)(上接第 60 页)8KISHNAMUTHY S，PELES Y Cross flow boiling in mi-cro-pin fin heat sinks C/Asme 2007 International Mechani-cal Engineering Congress and Exposition，2007，10:11-159HUA J Y，HOU W X，LI G，et al The experimental study onthe boiling bubble behaviors and heat transfer performance ofthe deionized water across various pin fins J Heat and MassTransfer，2021，57(10):1-14 10 LEE W，SON G，JEONG J J Numerical analysis of bubblegrowth and departure from a micro-cavityJ NumericalHeat Transfer Part B Fundamentals，2010，58(5):323-342 11 LIN Y H，LUO Y，LI W，et al Numerical study of flow re-versal during bubble growth and confinement of flow boilingin micro-channelsJ International Journal of Heat and MassTransfer，2021，177:121-135 12 PAJAPATI Y K，PATHAK M，KALEEM K M A compara-tive study of flow boiling heat transfer in three different con-figurations of micro-channelsJ International Journal ofHeat and Mass Transfer，2015，85(2):711-722 13 陈宏霞，孙源，肖红洋，等 微柱结构表面核态沸腾单气泡的数值模拟 J 化工进展，2019，69(12):4979-4989 14 KIM J S，CHO H K Advanced boiling heat transfer model fora horizontal tube with numerical analysis of bubble behav-ioursJ International Journal of Heat and Mass Transfer，2021，175:121-135 15 LI X L，HE L Q，QIAN P，et al Heat transfer enhancement ofdroplet two-phase flow in cylindrical micro-channelJ Applied Thermal Engineering，2021，11:64-74 16 BACKBILL J U，KOTHE D B，ZEMACH C A continuummethod for modeling surface tension J Journal of Computa-tional Physics，1992，100(2):335-354 17 DING S T，LUO B，LI G A volume of fluid based methodfor vapor-liquid phase change simulation with numerical os-cillation suppressionJ International Journal of Heat andMass Transfer，2017(3):348-359 18 YANG Z，PENG X，YE P Numerical and experimental investi-gation of two phase flow during boiling in a coiled tube J In-ternational Journal of Heat and Mass Transfer，2008，51(5/6):1003-1016 19 WEI J，PAN L，CHEN D，et al Numerical simulation of bubblebehaviors in subcooled flow boiling under swing motion J Nuclear Engineering and Design，2011(8):2898-2908 20 ZHOU Z，SHI J，CHEN H H，et al Two-phase flow overflooded micro-pillar structures with engineered wettabilitypattern J International Journal of Heat and Mass Transfer，2014，71:593-605 21 乔纬韬 两相流微通道换热特性的仿真与实验研究D成都:电子科技大学，2016 22 QU W L，MUDAWA I Measurement and correlation ofcritical heat flux in two-phase micro-channel heat sinksJInternational Journal of Heat and Mass Transfer，2004，47(10):2045-2059 23 KOSA A，KUO C J，PELES Y Suppression of boiling flowoscillations in parallel micro-channels by inlet restrictors JJournal of Heat Transfer，2006，128(3):251-260(责任编辑白林雪)09