基于气泡动力学分段调控浸润性强化核态沸腾

1、PAGE 16 -基于气泡动力学分段调控浸润性强化核态沸腾引言核态沸腾是一种高效传热方式，在相变换热设备中获得广泛应用，具有传热温差小、传热系数高的优点1-2。在航空航天、燃料电池、微纳电子芯片等产业高速发展的背景下，利用表面复杂微结构3-4与浸润性以强化核态沸腾换热，改善沸腾传热性能5-7，是近年来热门的强化表面换热的研究方向。表面浸润性能够显著影响沸腾表面换热性能8-10。Yim等11制备了表面接触角分别为81和15的普通表面和二氧化钛涂层表面，证明更亲水的二氧化钛涂层表面成核点数量增加，气泡体积减小脱离频率提高，表面传热能力提高64.1%。姜洪鹏等12通过高速相机观察到疏水表面气泡容易汇

2、聚，生长脱离周期较长。Sarode等13发现超疏水性限制板能够吸引并拉离加热器表面的气泡，防止蒸汽聚集在加热器表面，液体通道供应的周围液体可显著降低表面过热与抑制温度波动，传热系数（heattransfercoefficient,HTC）最大可提高82%。浸润性改性可改变气泡动力学特性，从单气泡的生长传热特性研究入手可加深对沸腾换热强化的理解。潘丰等14认为低热通量时的孤单气泡微层立气泡的微层蒸发在沸腾换热过程中起主导作用。Wang等15认为无改性表面上的气泡生长周期为2040ms，而改性表面上的时间小于20ms，传热性能更优。Sarker等16研究发现在低润湿的粗糙表面上，气泡的膨胀和收缩效

3、应造成基底附近的液体传播扰动，导致等待时间延长。在明确了表面亲疏水性分别对沸腾热力学和气泡动力学的影响后，混合浸润性换热表面的强化机理也逐渐清晰。Liang等17通过构建核池沸腾装置，研究了基于超亲水表面的疏水性点状或条纹状混合润湿性表面对核池沸腾传热和临界热通量（criticalheatflux,CHF）的影响，发现润湿性的增加会强化表面的再润湿效应改善表面CHF，低热通量时成核增强，疏水区域传热系数高于亲水区域。Kong等18认为混合润湿表面提供的成核位点密度增大，气泡经过多次聚结后滑向光滑区域，在亲水性区域气泡脱离加速。Shen等19在抛光铜表面沉积具有疏水性质的圆形斑点，制备的101和

4、接近0交替的疏水性和亲水性混合表面，可通过促进气泡成核和扰动效应使核沸腾的起始温度下降35%；过热度达到20K时，与纯铜表面相比HTC最大增强超过300%。Hsu等20发现亲疏水性区域交错次数与面积比值影响铜圆筒壁面过热度和表面传热性能，交错次数越多、亲水区域比值越大使表面成核越容易。除搭建核态沸腾实验台外，由于仿真软件对表面浸润性调控具有明显的便利性，混合浸润性的模拟研究得到了大量关注21-23，主要方法包括格子玻尔兹曼方法、CFD-VOF方法和分子动力学模拟等。Ma等24采用格子玻尔兹曼模型进行池沸腾传热数值模拟，综合利用疏水表面在低过热度以及亲水表面在高过热度时的换热优势，发现柱顶103

5、、底部53的浸润性表面沸腾换热性能最好。Lee等25则探索使用交叉纵横比作为衡量疏水点的形状的变量，研究了亲水性表面上的不同形状疏水点对换热的影响。Chen等26和张龙艳等27采用分子动力学模拟方法，认为光滑表面上更大面积的强亲水区域以及强亲水纳米结构有利于成核，原因是纳米尺度下浸润性增强使固液界面间热阻效应减弱，热量传递效率提高，越容易核化。综上，近年来利用浸润性的调控强化核态沸腾的研究，已经阐明了亲疏水性的影响以及混合浸润性表面的换热优势。但以上方法主要集中在对换热表面的静态修饰，对于不同阶段所展现的换热特点利用不足，因此可以从单个气泡生长脱离阶段的不同时刻入手，有目的性地进行表面浸润性调

6、控。根据气泡动力学需求，亲水性可促进气泡脱离，疏水性可促进气泡生长。根据本课题组对相同微结构表面沸腾气泡的脱离阶段进行亲水性调节研究结果28，本文进而针对生长阶段，尤其是t0.152ms阶段进行疏水性调节，以获得生长阶段疏水性条件的强化效果，为强化核态沸腾换热方法提出新的思路。1模型及方法1.1计算域基本设置计算域底部为厚50m的加热硅板，硅板底面热通量设置为58kW/m2；硅板上方为微柱区域，长宽高分别为50、50和150m，间隔为100m，以阵列形式排布；以平面y=0作为对称轴，流体区域高度设置为1500m，计算域底面面积为1350m675m，满足沸腾计算要求。首先进行纯导热计算，计算域侧

7、壁面设置为绝热边界条件，上表面为压力出口，计算域几何尺寸与边界条件如图1所示；成核点距离微柱基底10m，基于之前的实验结果29在过热度达4K时放置种子气泡，半径为7.5m；将三个侧壁面改为压力出口条件后开始核态沸腾计算。计算域内的网格均是六面体网格，为提高计算精度与稳定性，采用如图2所示的分区域网格细化方案。015m高度范围内的基底区域和气泡区域进行二阶加密后最小网格尺寸为2.5m；除气泡区域外，15150m范围内网格大小均为10m，随着高度增加利用边界层网格将高度增加为12、14.4和17.28m，以节省计算资源与时间30。图1图1计算域几何尺寸与边界条件Fig.1Thegeometryan

8、dboundaryconditionsforthecomputationaldomain图2图2分区域网格细化方案28Fig.2Meshschemeswithazonalgridrefinementtechniqueindifferentzone281.2数学模型及求解方法本文基于计算流体力学方法（computationalfluiddynamics，CFD），利用软件Fluent15.0，模拟计算微柱结构表面上单个气泡生长脱离过程，采用流体体积多相流模型（volumeoffluidmodel，VOF）31实现液态水和水蒸气两相的求解以及气液交界面的精确捕捉，基本参数设置见表1。气泡气液界面处

9、的蒸发占总蒸发量的很大一部分，采用能量跳跃模型（energyjumpcondition）32计算气液界面处相变。气泡的微层蒸发在生长阶段初期快速蒸发33，本文采用简化的微层建模思想34，假设微层区域液体不发生流动，且微液层顶部温度保持为饱和温度Tsat不变，计算得到微层高度变化与气相变化率。气液两相间的表面张力采用连续表面张力模型（continuumsurfaceforce,CSF）。将气泡直径的变化与微层蒸发量占比分别与文献对比，误差均在可接受范围内，模型的准确性得到验证，所采用的基本数学模型与具体验证过程可见参考文献30。表1数值模拟涉及工质的物性参数Table1Physicalparam

10、etersemployedinthenumericalsimulations材料密度/(kg/m3)比热容/(J/(kgK)热导率/(W/(mK)黏度/(kg/(ms)气液界面表面张力系数/(N/m)单晶硅2330766148液态水958.45664215.50.67220.0002820260.06164水蒸气0.59762079.80.02461.2210-30.061642讨论与分析2.1气泡动力学规律基于目前研究结果28，接触角为48表面的沸腾气泡在t=0.152ms时刻，具有最大蒸发功率；而改变表面浸润性，不仅此气泡最大蒸发功率的时刻不同，气泡生长阶段到脱离阶段的转折时刻不同、气泡整

11、个周期也发生相应变化。为表征不同浸润性对生长阶段沸腾速率的影响，本文统一选取同一时刻t=0.152ms时，利用Fluent软件直接调整固体壁面本征接触角改变表面浸润性，对比其沸腾规律的区别。图3展示了接触角为48表面在t=0.152ms时调整为接触角为20时气泡形态随时间的演变。通过观察可见t=0.152ms调控后，气泡的生长趋势主要由横向生长为主转变为纵向生长为主，气泡呈现缩脖子的蘑菇状；亲水性可促进液体沿微结构表面的吸入，缩小气泡与壁面的接触面积，促进气泡向上生长及脱离。图4展示了不同浸润性表面气泡体积随时间变化规律。t=0.152ms时，接触角分别为48和60的表面上气泡基本相同为0.0

12、65mm3，90和110表面的气泡体积分别为0.077mm3和0.087mm3，相同时间内气泡沸腾速率分别增大18.5%和33.8%。在48-20的表面上，气泡于t=1.120ms脱离，最终时刻体积为0.40mm3。60-20表面的气泡于t=1.13ms脱离，体积为0.43mm3，较前者增大7.5%；90-20表面的气泡于t=1.170ms脱离，体积为0.57mm3，增大42.5%。对于表面接触角为110的表面，尽管在相同计算域中气泡未实现完全脱离，但其在t0.152ms时间内已然展示出了疏水性的提高显著增大了气泡体积，加速了气泡生长速度；至t=0.495ms，气泡体积为0.47mm3，与原4

13、8、原60以及原90表面上的气泡体积相比分别增大74.1%、62.1%和34.3%。以上均表明生长阶段调控表面疏水性可明显增大气泡体积。图3图348-20表面上气泡形态随时间t的变化Fig.3Variationofbubbleshapeon48-20withtimet图4图4气泡体积Vbubble随时间t的变化Fig.4VariationofbubblevolumeVbubblewithtimet利用VOF模型计算过程可时时追踪气液界面，进而利用曲面积分的数值方法直接导出气泡与换热表面间的接触面积，图5(a)展示的是气泡与换热表面基底的接触面积Ab-s随时间的变化规律。根据之前研究35，气泡的

14、单个周期可分为生长阶段和脱离阶段，气泡在基底的接触面积减小后不再回升的时刻，为开始进入脱离阶段的临界时刻，即生长阶段和脱离阶段的分界时刻。基于以上分析，由图3(a)可知在48-20的表面上气泡于t=0.45ms进入脱离阶段，60-20表面上的气泡在t=0.42ms进入脱离阶段，提前6.7%；疏水性较强的原90和原110表面上曲线规律性不如原48和60表面明显，但仍可认为气泡分别于t=0.37ms和t=0.27ms进入脱离阶段,分别提前17.7%和40.0%。生长阶段疏水性增强了气泡在微结构间的铺展，t=0.15ms时48表面上底面接触面积为0.10mm2，60、90和110表面上底面接触面积分

15、别为0.11、0.18和0.23mm2，分别增大了10.0%、80.0%和130.0%。由于微结构的存在能够改善表面沸腾换热性能，因此与微结构间有更大接触面积意味着具有更加高效的微层和气液界面蒸发功率。图5图5气泡与壁面接触面积随时间t的变化Fig.5Variationofcontactareabetweenbubbleandsurfacewithtimet图5(b)展示的是气泡与微柱侧壁面接触面积Ab-p随时间的变化规律。气泡在生长阶段不断包覆微柱并向周围迅速扩张，接触面积曲线近似直线增加。表面疏水性增强使气泡的扩张更加显著，110表面上的曲线斜率是48和60曲线斜率的近两倍。对于48表面，

16、气泡与侧壁面接触的面积在t=0.152ms之前迅速增大，变为20之后短暂下降。原因如2.2节所述气泡与微柱表面间形成的薄液膜，由周围饱和液体沿界面内流产生，减小了气泡与侧壁面的接触面积；由于在t=0.152ms时气泡均处于强烈的生长阶段，气泡继续包覆微柱使接触面积重新增加至t=0.495ms，之后随气泡向上逐渐加快脱离微柱逐渐减小为0。在t=0.15ms时，90和110表面上的底面接触面积相比48表面分别增大了0.8倍和1.3倍，与侧壁面的接触面积则分别增大0.35倍和0.9倍，说明气泡在微结构底面的扩展受浸润性影响更明显，生长阶段疏水性导致更大的底面铺展面积。2.2浸润性引起的界面曲率变化及

17、受力分析图6显示了在t=0.15ms时，48、60、90和110表面y=0截面上气泡温度云图与形态对比。通过对比气泡在微结构间部分的占比大小，可发现疏水表面使气泡更易横向铺展，是由微柱间隙的黏附力增大导致。图7显示了在t=0.50ms时，原48、60和90表面以及t=0.495ms时110表面y=0截面上气泡的温度云图与形态对比，可见原疏水性越强的表面上气泡包裹的微柱越多，其处于热边界层区域内的部分越多，更有利于热量的交换。图6图6温度分布云图(t=0.15ms，y=0截面)Fig.6Temperaturedistributionwithy=0att=0.15ms图7图7温度分布云图(t=0.

18、5ms，y=0截面)Fig.7Temperaturedistributionwithy=0att=0.5ms为分析亲疏水性分段调控对微柱结构间气泡界面受力的影响，将气泡受力以微柱柱顶高度为界分为方向竖直向下的Fsa与竖直向上的Fsb，如图7(a)所示。简化认为Fsa和Fsb可分别代表壁面吸附力和浮力的物理作用，气泡所受合力可表示为：Fs=Fsa+Fsb(1)根据Young-Laplace公式，气泡弯液面所受的表面张力Fs与曲率半径R的关系为：Fs=2RFs=2R(2)R=1R=1(3)式中,为表面张力系数，N/m；为界面曲率。对比观察图7(a)、(b)中红色圆圈部分，原60表面上的气泡与原48

19、表面相比，位于两个柱顶之间的气液界面下沉，曲率增大曲率半径R减小，沿着气泡下半界面方向向上的Fsb增大，壁面吸附力Fsa因亲水性而减弱，导致合力Fs增大，促使气泡向上运动加快脱离。亲水性调控使气泡周围饱和流体沿界面回流形成新薄液膜区域，对比观察图7(a)与图7(b)、(c)，可发现该区域取代了部分气泡与微柱壁面直接接触形成的干烧区，使相变换热量增加、蒸发速率提高，黑色圆圈所标识的微柱间隙内出现了部分低温区域。通过观察图7(d)，薄液膜区域在原110表面微柱侧壁上并不明显，但相似低温区域在右侧微柱上半部分的黑色圆圈内出现，说明110-20的表面上仍存在一定程度的薄液膜效应。2.3成核点过热度的演

20、变规律经以上分析，分段调控浸润性对气泡整个周期内的动力学和热力学过程具有重要影响。疏水性使气泡在微结构间横向扩展以加强换热，亲水性使界面曲率提高加快脱离，微液膜效应增强使得基底壁面温度迅速降低，表面过热度明显下降。图8为成核点处过热度随时间的变化曲线，壁面过热度在加载种子气泡后的短时间内波动，随后因气泡体积增大吸热能力增强而不断降低；当气泡进入脱离阶段之后，随着气泡与基底接触面积的减小与蒸发功率的减弱，过热度重新回升。原48、60和90表面上成核点过热度分别从t=0.330ms、0.500ms和0.620ms开始回升，至气泡完全脱离微结构为止。原48和原60表面上的回升速度分别为2.0K/ms

21、和1.5K/ms，90表面上则为2.5K/ms，说明气泡的脱离速度更快。60表面上成核点处温度始终高于48表面，但根据气泡体积与微结构的接触面积变化以及下文所述的蒸发功率规律，仍可认为其整体换热情况强于48表面。图8图8成核点过热度T随时间t的变化Fig.8VariationofsuperheatTatthenucleationsitewithtimet2.4气液界面和微层蒸发功率图9展示了气泡的微层蒸发功率qml和气液界面蒸发功率qlv随时间的变化规律。气泡微层的蒸发和气液界面的蒸发是沸腾气泡传热的两个主要来源。通过对比分析蒸发速率的变化，可研究核态沸腾传热受浸润性分段调控的影响。图9图9气

22、泡蒸发功率随时间t的变化Fig.9Variationofevaporationratewithtimet微层部分的蒸发对沸腾气泡传热具有重要贡献。由图9(a)微层蒸发功率曲线可以看出，经过在气泡生长初期的快速增长达到最大值后，在生长后期和脱离阶段逐渐减小，最后在气泡完全脱离时减小到零。根据2.1节中分析，生长阶段疏水性增强了气泡在微结构间的铺展，与基底的接触面积增大，微层蒸发功率相应增强。在t=0.150ms时，48表面微层蒸发功率为0.30W，60、90和110表面上微层蒸发功率分别为0.32、0.49和0.66W，分别增大6.7%、63.3%和120%，与底面接触面积的增大规律相符合。气液

23、界面处相变蒸发量受微结构整体过热度影响很大。如图9(b)曲线所示，在生长阶段随着气泡界面在微结构间的扩展，气液界面换热增强，之后气泡主体界面逐渐向上远离微结构，周围液体过热度逐渐减小，相变过热度及气液界面蒸发功率减小。在气泡完全脱离微柱之前，气液界面处的换热速度已逐渐减小至较低水平。对于48表面，接触角在t=0.152ms变为20后，如图7中所示的薄液膜导致了气液界面蒸发功率的提高，t=0.200ms时原60、原90和原110表面上气液界面蒸发功率比原48表面上分别增大2.0%、27.3%和61.5%。2.5单周期内换热性能的强化为了评估分段调控浸润性所产生的强化换热效果，定义核态沸腾表面上的

24、总瞬时蒸发功率为：q(t)=Qml+Qlvt=t0qmldt+t0qlvdttqt=Qml+Qlvt=0tqmldt+0tqlvdtt(4)一个气泡完整的周期包括生长、脱离和等待阶段。研究结果29表明气泡在微柱表面上的等待时间很短，一个气泡的周期时间t0可近似等于生长和脱离阶段的时间，气泡的离开频率f0可认为近似等于1/t0。因此表面上平均蒸发功率可以表示为：q=q(t0)t0f0=q(t0)q=qt0t0f0=qt0(5)由图10可知，气泡在48-20的换热表面上的平均蒸发功率为0.36W，60-20和90-20表面上则为0.39和0.48W，分别提高8.3%和33.3%。图10图10气泡平均蒸发功率随时间t的变化Fig.10Variationofaverageevaporationrateqqwithtimet值得注意的是，生长阶段和脱离阶段内的热流均因疏水性调控而明显提高，在48-20的表面上生长阶段内的平均热流为0.57W，60-20和90-20的表