电场和加热器特性对饱和池沸腾传热影响的介观数值方法研究*

胡剑 张森 娄钦

1) (国能浙江宁海发电有限公司,宁波 315600)

2) (上海理工大学能源与动力工程学院,上海 200093)

采用耦合电场模型的相变格子Boltzmann (LB)方法研究了饱和池沸腾传热性能,重点分析了均匀电场作用下加热器表面润湿性以及加热器长度对沸腾过程中气泡生成、合并、断裂等动力学行为的影响以及气泡的动力学行为对池沸腾传热性能的影响.结果表明,电场的作用能否强化沸腾传热与加热器的长度以及润湿性有直接关系.对于亲水表面,当加热器长度 ≤6.25 时,由于加热器尺寸较小,沸腾过程中加热器表面产生的气泡相互作用力弱,此情况下电场的存在使得气泡体积减小,沸腾被抑制.当加热器长度 6.25＜≤9.375 时,均匀电场均能提高临界热流密度(critical heat flux,CHF),且在此加热器长度范围内,CHF 提高的百分比随着电场强度的增大而增大.这是因为 6.25＜≤9.375 时,更长的加热器为气泡的生成提供了充分的空间,气泡之间的相互作用力较强,均匀电场作用下的气泡间距增大,气泡数量增加,且CHF 提高百分比逐渐增大;当＞9.375 时,再润湿阻力随着加热器长度的增大而增大,导致沸腾过程中产生的蒸气在电场力作用下容易被紧贴于加热表面,增加了固体与流体之间的换热热阻,并在气泡根部形成不利于气泡向中间移动的涡,减缓了加热表面热流体与两侧较冷流体的热质交换,CHF 提高的百分比随着加热器长度的增大逐渐减小.对于疏水表面,随着长度的增大,CHF 提高百分比同样为先增大后减小,然而其阈值增大.

1 引言

沸腾是一种典型的气液相变现象,被广泛应用于微纳机电强化换热、高热流密度电子设备冷却、蒸气发电、医学消毒和药物制取、化工杂质去除、食品加工等领域[1−3].如何增强沸腾传热性能以及提高临界热流密度(critical heat flux,CHF)是沸腾研究的核心目的[4−6].从1961 年Clubb[7]提出施加外部电场可以增强流体的传热后,研究者们展开了大量的研究工作以了解电流体动力学(electrohydrodynamic,EHD)增强沸腾传热的机理.20 世纪80 年代以来,EHD 强化沸腾传热的实验研究受到了广泛关注.

在EHD 强化沸腾实验中发现电场的放电源和加热表面特性的不同对沸腾换热具有很大影响.Madadnia 和Koosha[8]利用电极丝施加电场,在直径为1 mm的加热丝上研究了电场作用下孤立气泡的离核直径、成核速率(频率)和成核点密度等特性.发现当典型热流密度为35.4 kW/m2并保持恒定时,电极电压为6.0,6.6 和7.5 kV 时气泡的尺寸和成核频率减小,而电极电压在8 kV 及以上时未观察到成核现象.在5 mm×5 mm 矩形加热器表面,Gao 等[9]研究了沸腾过程中EHD 对R113制冷剂成核及其后续生长动态的影响.他们观察到随着外加电场强度的增大,气泡离开直径(体积)和离开频率减小,而气泡生长时间和等待时间延长.他们还观察到对于单个R113 气泡的情况,电场的存在可增强沸腾传热.实验中除了研究EHD对沸腾过程的影响,还有学者通过在电场中放置单个气泡,在等温条件下研究EHD 对气泡动力学行为的影响.Dong 等[10]调查了直流电场对单个惰性气泡生长、变形和脱离等行为的影响.发现在电场作用下,气泡被电应力的水平分量强烈压缩,并被垂直分量拉长.Zu 和Yan[11]也在EHD 对气泡形变影响的实验研究中观察到相同的现象.

除了电场对气泡成核过程的影响,一些学者还研究了电场对整个沸腾曲线以及CHF 的影响.Quan 等[12]利用网状电极在直径约为18 mm 的圆形加热表面对光滑表面和肋表面沸腾过程的电流体动力学进行了实验研究,并探讨了不同壁面过热度情况下EHD 对沸腾换热性能的影响机制.他们发现电对流是低过热区增强换热的主要机制,且电场抑制了气泡成核.在中过热区,电对流效应减弱,电场力阻止气泡从受热表面生长和脱离,导致传热恶化和壁温升高.在高过热区,电场力将气柱解体成更小的气泡,并增大气泡脱离的频率,有助于提高沸腾传热性能.Liu 等[13]研究了10 mm×10 mm 的矩形表面上均匀电场对不同尺寸微针鳍结构表面池沸腾换热的影响,发现蒸气层的破裂有利于CHF 的增强,而场陷阱效应不利于CHF 的增强.因此,电场对CHF 的影响由以上两个因素共同决定.在1 cm×1 cm 的加热表面,Garivalis和Manfredini[14]研究了微重力条件下电场对微结构表面(方柱)和平面表面池沸腾的影响.结果表明,电场和微结构的结合使微重力条件下的CHF 比普通表面提高了114%—144%.

以上在一定的电场强度和加热器尺寸下研究了沸腾换热性能,极大地提高了对该问题的认识.近年来,随着计算机技术和数值方法[15−19]的快速发展,采用模拟方法研究EHD 及其增强沸腾换热受到大量学者的广泛关注.在6 mm×6 mm硅橡胶绝缘疏水表面,Songoro[20]采用有限元法对电场作用下单个水滴的运动进行数值模拟.结果表明,在水平和垂直施加交流电场的情况下,液滴振动的预测共振频率与实验共振频率吻合.采用耦合水平集(level set)和流体体积(volume-of-fluid)的CLSVOF(coupled level set and volume-of-fluid)方法[21],Pandey 等[22]研究了在6.825 mm 加热器表面上电场和过热度对气泡成核方式与传热特性的影响,发现在低过热度下,随着电场强度的增大,气泡释放速率增大,同时气泡间距减小,形成了离散气泡.当过热度较大时,较高的电场强度导致形成大量的蒸气柱.然而,以上方法无法模拟气泡成核,也无法对整个沸腾曲线进行刻画.最近,格子Boltzmann (LB)方法[23]由于可以从微观尺度描述分子之间的相互作用、不需要追踪相界面、自动成核、可以方便描述沸腾的所有阶段等优点被广泛应用于气-液相变的问题研究中[4,5,24−26].采用格子Boltzmann 方法并耦合电场作用力模型,张浏斌等[27]模拟了均匀电场作用下在5 个格子单位的小型加热器表面上,单气泡成核、生长、脱离过程,重点研究了重力加速度和电场强度对气泡动力学的影响,发现了均匀电场有助于低重力加速度下池沸腾气泡的脱离,从而提高换热效率.Yao 等[28]研究了14.4 mm 加热器表面上电场强度对气泡行为的影响,发现了电场可以将气泡变形为扁长形和扁圆形,分别导致气泡加速和减速离开热表面.最近,Li 等[29]用格子Boltzmann 方法研究了电场作用下在柱状结构表面上的沸腾传热性能,发现外加电场对柱状结构表面的池沸腾有正、负两方面的影响,并将润湿性修饰区域应用于柱状结构表面的柱顶,进一步增强了电场作用下柱状结构表面的沸腾换热.

以上研究表明,电场对沸腾过程的气泡动力学行为和传热性能的影响较为复杂,其与电极施加方式、加热表面设置、过热度等有密切联系.然而,大部分关于电场增强池沸腾换热的研究都是在特定的加热器长度基础上进行的,很少有研究人员考虑电场对不同长度加热器池沸腾传热的影响.相关研究表明: 加热器长度对气泡动力学、流场分布和沸腾传热性能有显著影响.Rainey 和You[30]和Lee 等[31]在实验研究中发现,随着长度增大,加热器侧面冷流体与加热表面上方热流体的对流阻力增加,从而导致CHF 逐渐降低.Zhang 等[32]采用二维格子Boltzmann 方法模拟不同长度加热器的池沸腾过程也得出类似结论.施加电场时,加热器长度对气泡动力学行为以及沸腾传热性能的影响尚不清晰.因此,本文将研究不同长度加热器下均匀电场对气泡动力学行为、沸腾曲线、以及临界热流密度的影响.由于实验工作中不同尺寸加热器样本与电极设置较为复杂,且基于放电源与电场之间的复杂关系,很难保证施加于不同加热器样本的电场保持一致.而LB 方法在研究气相相变问题时具有不需要追踪界面、成核自发发生、流固耦合易于处理等优点[33].本文采用二维伪势LB 和电场的耦合模型,数值模拟均匀电场作用下不同长度加热器的池沸腾过程,详细研究不同长度加热器下均匀电场对沸腾曲线,气泡动力学行为和CHF 的影响.

2 数值方法

2.1 格子Boltzmann 方法

在LB 模型[34]中,流体的密度分布函数为

式中wi是与离散速度模型相关的权重系数.在(1)式中,力项Fi表达式如下所示:

其中 Δu=Fδt/ρ,F为相互作用力,其包含粒子间相互作用力Fint,外力Fext以及电场力Fe,具体表达式如下:

式中粒子间相互作用力Fint可以写为

其中ψ(x)为有效质量,其与局部密度和压力有关:,其中G是相互作用强度,局部压力p由状态方程得到,本研究中采用Redlich-Kwong (R-K)状态方程,对应的压力形式如下:

其中g为重力加速度,ρave为每一步计算时整个计算域的平均密度.方程(4)中电场力Fe将在下节具体给出.

温度分布函数gi(x,t)其对应的LB 演化方程为[34]

其中U是实际的流体速度.(8)式中的源项ϕ表示为

其中cυ是流体的比热容.流体的宏观量可由密度分布函数和温度分布函数的矩得到:

特别指出,模拟中气液界面处的物性χ(如热扩散系数、介电常数、动力黏度、比热容)由(15)式根据气液份额给出：

二维格子Boltzmann 方法的D2Q9 模型[35,36]对模拟沸腾传热过程具有良好的适用性[5,29].因此,本文采用D2Q9 模型,其离散速度为

其对应的权重系数为w0=4/9;wi=1/9(i=1—4);wi=1/36(i=5—8).

2.2 电场模型

根据电流体动力学理论,电介质在电场中所受的电场力Fe为[37]

其中ρe是自由电荷密度,ρ是流体密度,E是电场强度,ε0为真空介电常数,ε为流体介电常数.等式(17)右边第一项为库仑力,理想流体中的自由电荷为零,因此可以忽略库仑力.第二项为介电泳力,表示电场中介电常数分布不均匀所产生的力.在两相流中,介电泳力主要由气相与液相的介电常数差引起,这决定了气泡在电场中的运动方向.第三项为电致伸缩力,对于不可压缩流体,此项可以忽略.因此流体所受电场力Fe为

对于电场强度E可由下式求解:

电场强度E表示为电势V的梯度,即E=-∇V,于是(19)式可以改写为

为了求解(20)式,本文采用文献[26,38]中的LB方程:

3 物理模型和计算参数

如图1 所示,整个计算区域设置为Lx×Ly=400×1000 格子单位,红色部分为加热台,位于底部中心,高度H=30 格子单位,长度为LH,加热器表面具有不同润湿性,由润湿性边界条件[39]描述,计算区域内红色以外部分为流体区域.另外,在加热台下表面及计算域上边界施加大小为V0的电压,并在计算域上界面接地线,以保证上界面电势为0,加热台上表面的电势为V0.由于加热器上表面与上边界平行,类似于两平行带电平板,会在蓝色阴影区域形成电场强度为E=V0/(Ly-H) 的均匀电场.在计算区域的左右边界,对密度分布函数和温度分布函数采用周期边界格式;在计算区域的上边界,对密度和温度分布函数采用适用于自由边界的对流边界格式,以消除边界处可能给内部流场带来的影响;在下边界对密度分布函数采用半步长反弹格式,对温度分布函数在加热台上边界采用流固耦合边界条件[40],考虑固体导热率以及流体物性差异,为确保流固边界温度和热通量传递的连续性,加热台左右两边及流体侧下边界采用绝热边界条件,在加热台内部下边界设置等温加热方式.对于理想电介质模型,上下边界采用非平衡外推边界,左右边界为周期边界条件[26].初始时刻,流体温度为饱和温度Tsat,密度为Tsat对应的饱和密度ρl,加热器下边界温度为Tb.在数值模拟中,R-K 状态方程中的参数设置为a=2/49,b=2/21,R=1.0,计算得Tc=0.1961,pc=0.1784.取饱和温度Tsat=0.9Tc,则饱和液体和饱和气体的密度分别ρl=5.426,ρv=0.8113.比热容为cpl=cυl=4.0,cpv=cvυ=2.0,动力黏度µl=0.3257,µv=0.0487,热扩散系数αl=0.05,αv=0.06,液体和蒸气的导热比为λl/λv=11.15,流体的潜热hfg=0.624,表面张力系数σ=0.095,重力加速度g=(0,-0.00008).文献[24]在研究加热器物性对CHF 影响时,对流固耦合边界条件加热器物性取(ρcp)s/(ρcp)l=0.4,1.0,1.5,发现加热器物性对CHF 值影响不明显.因此,本文取(ρcp)s/(ρcp)l=1.0,即(ρcp)s=21.704,αs=1.5,则固体/液体和固体/气体的导热系数比为λs/λl=30,λs/λv=334.4.真空介电常数取ε0=2.2360,蒸气与液体的介电常数比为εv/εl=0.4472.需要特别指出的是,以上单位均为格子单位,其与实际物理单位的转换关系如表1 所示.

表1 格子单位与物理单位转换Table 1.The unit conversion from lattice unit to physical unit.

图1 物理问题示意图Fig.1.Diagram of the physical problem.

为了便于与前人工作进行对比,如无特别说明均采用无量纲单位,所采用的特征长度、特征速度、特征速度和特征热流密度如下:

经过无量纲化得到的无量纲长度、无量纲时间、无量纲热流分别定义为L*=L/l0,t*=t/t0,q*=q/q0.另外绘制沸腾曲线需要的时间和空间平均热流如下:

其中ta到tb时间段是沸腾开始后的很长一段时间,q是空间平均热流:

4 模拟结果

4.1 无电场作用时气泡行为的变化对不同长度加热器表面沸腾传热性能的影响

在研究均匀电场对不同加热器长度的沸腾传热影响之前,首先对不同加热器长度下的沸腾传热进行研究.本节模拟了不同长度加热器的沸腾过程,考虑到整个计算域尺寸以及物理模型设置,格子单位长度取LH=50,100,120,150,200,300,对应无量纲长度=3.125,6.25,7.5,9.375,12.5,18.75.为了研究润湿性的影响,本文分别考虑了亲水表面和疏水表面,气泡与加热表面静态接触角分别设置为θ=45°,θ=120°.

图2(a)和图2(b)分别展示了亲水表面和疏水表面不同长度加热器下的沸腾曲线.从图2(a)可以观察到对于亲水表面,当=3.125 时,随着壁面过热度增大,空间时间平均热流密度越来越大,整个沸腾曲线并没有出现过渡沸腾阶段,即不存在临界热流密度.为了揭示这一现象出现的原因,图3给出了=3.125 时不同加热温度下同一时刻的气泡形态.从图中可以发现，当壁面过热度较小,即Tb=0.98Tc时,加热表面每次仅生成一个气泡.当壁面温度升高到Tb=1.00Tc时,尽管壁面过热度的增大导致生成气泡的尺寸增大,然而由于加热器尺度较小,加热表面同样只能生成一个气泡,即加热表面生成的气体始终以孤立的形式存在,并没有发生由于气泡的合并导致的沸腾换热恶化的现象.而当壁面温度升高到Tb=1.02Tc,加热表面形成稳定的蒸气柱,蒸气柱顶端持续形成气泡并继续上升,此沸腾状态为膜态沸腾.由于加热器尺寸较小,在加热表面形成的小尺寸蒸气膜在浮力和重力的共同作用下整体处于上升状态.由以上结果可知,由于加热器尺寸较小,没有为出现多个气泡提供充足的空间,因此没有出现因为多气泡合并和脱离速率降低导致的过渡沸腾阶段,此时随着壁面过热度的增加,沸腾由核态沸腾阶段直接转变为膜态沸腾阶段,空间时间平均热流密度在整个沸腾阶段逐渐升高.

图2 (a)亲水表面和(b)疏水表面不同长度加热器的沸腾曲线Fig.2.Boiling curves of heaters with different lengths: (a) Hydrophilic surfaces;(b) hydrophobic surfaces.

图3 =3.125,t*=66.99 时刻,不同加热温度下的气泡形态(a)Tb=0.98Tc;(b) Tb=1.00Tc;(c) Tb=1.02TcFig.3.The bubble morphology at =3.125,t*=66.99 moments with different heating temperatures: (a)Tb=0.98Tc;(b)Tb=1.00Tc;(c) Tb=1.02Tc.

图4 t*=55.83,Tb=0.98Tc 条件下,不同长度加热器下沸腾的气泡形态和流场 (a) =6.25;(b)=18.75Fig.4.Bubble morphology and flow field for boiling with different length heaters under t*=55.83,Tb=0.98Tc conditions: (a) =6.25;(b) =18.75.

当表面疏水时不同加热器长度对沸腾曲线的影响规律(如图2(b)所示)与亲水表面时类似,在加热器长度较小时,即=6.25 时沸腾曲线同样未出现过渡沸腾阶段.随着加热器长度的增大,核态沸腾阶段和膜态沸腾阶段热流密度均减小,CHF同样由逐渐降低到维持相对稳定.然而当表面疏水时,CHF 由逐渐降低到维持相对稳定的阈值增加到为=12.5.

4.2 均匀电场作用下气泡行为的变化对亲水表面沸腾传热的影响

本节模拟了在均匀电场作用下亲水表面不同长度加热器的池沸腾传热过程.加热器长度和亲水表面润湿性与4.1 节相同,即=3.125,6.25,7.5,9.375,12.5,18.75,θ=45°.均匀电场强度E取0.05442,0.10884,0.16327,这里无电场作用时的结果(E=0)用于对比分析电场作用对沸腾过程的影响.

图5 为不同长度加热器在不同电场强度下的沸腾曲线,其中ONB (onset of nucleate boiling)表示核态沸腾起始点.从图5(a)和图5(b)可以发现，当=3.125,6.25 时,均匀电场对核态沸腾阶段和膜态沸腾阶段均表现为抑制作用,随着电场强度的增大,核态沸腾和膜态沸腾的热流密度有所下降.为了揭示导致这一现象的原因,图6 为=6.25,电场强度E=0,E=0.16327 作用下,核态沸腾(图6(a))和膜态沸腾(图6(b))气泡形态对比.从图6(a)可以发现,核态沸腾时电场强度E=0.16327 作用下的气泡尺寸比E=0 的气泡尺寸小.沸腾过程主要通过饱和液体相变吸热过程进行热量传递,气泡尺寸减小意味着由饱和液体相变为气体的总质量减小,即潜热换热减小,从而导致核态沸腾阶段热流密度降低.从图6(b)可以发现，在膜态沸腾阶段,蒸气在加热表面铺展成膜.无电场作用时(E=0)气体由两边向加热器中心聚集并在浮力的作用下向上运动,在加热表面中心位置形成向上的蒸气柱,铺展在加热表面蒸气膜的气液界面比较平滑.当存在电场作用时(E=0.16327),同样在加热表面形成蒸气膜,但在电场力的作用下,铺展在加热表面蒸气膜的气液界面发生形变,蒸气膜更加均匀地铺展在加热表面,增大了蒸气膜的整体厚度,导致气膜施加在整个加热表面的热阻增大,膜态沸腾的热流密度降低.可以得出结论:当加热器尺寸较小时,核态沸腾阶段均匀电场减小气泡尺寸,降低了沸腾换热过程中的潜热换热,从而降低热流密度;膜态沸腾阶段均匀电场使加热表面蒸气膜更加均匀地铺展,增大了热阻从而降低膜态沸腾阶段热流密度.

图5 亲水表面不同长度加热器在不同电场强度下的沸腾曲线 (a) =3.125;(b) =6.25;(c) =7.5;(d) =9.375;(e) =12.5;(f) =18.75Fig.5.Boiling curves of heaters of different lengths on hydrophilic surfaces under different electric field strengths: (a) =3.12;(b) =6.25;(c) =7.5;(d) =9.375;(e) =12.5;(f) =18.75.

图6 =6.25 时,在均匀电场强度 E=0,0.16327 作用下 (a) Tb=0.97Tc,核态沸腾状态;(b) Tb=1.02Tc,膜态沸腾状态的气泡形态Fig.6.The bubble morphology of (a) Tb=0.97Tc,nucleated boiling state;(b) Tb=1.02Tc,film boiling state under the action of uniform electric field strength E=0,0.16327,=6.25.

图7 不同长度加热器在电场强度 E=0,0.10884 下的气泡形态对比 (a) =7.5;(b) =9.375;(c) =12.5;(d) =18.75Fig.7.Comparison of bubble morphology of heaters with different lengths at electric field strength E=0,0.10884: (a) =7.5;(b) =9.375;(c) =12.5;(d) =18.75.

为了进一步研究在均匀电场作用下,加热器长度与CHF 的关系,图8 展示了不同加热器长度下均匀电场对CHF 提高的百分比.首先,在相同加热器长度下,随着电场强度的提高,CHF 提高的百分比越来越大.此现象与图5 和图7 中规律一致,即电场强度越大,加热表面气泡越不易合并,会使沸腾在更高过热度下达到CHF,从而使CHF 越来越大.另外,当≤9.375 时,在同一电场强度下,随着加热器长度的增大,CHF 提高的百分比逐渐增大;当 9.375＜≤18.75,在同一电场强度下,随着加热器长度的增大,CHF 提高的百分比逐渐减小.为了解释这两种现象,下面分别选取符合这两种条件的工况进行研究.

图8 不同电场强度下亲水表面加热器长度与临界热流密度的关系Fig.8.Relationship between hydrophilic surface heater length and critical heat flow density at different electric field strengths.

图9 展示了在电场强度E=0.10884,加热器长度分别为=7.5 和=9.375 时,CHF 对应加热温度(Tb=0.995Tc)下空间平均热流密度随时间的变化,为保证沸腾过程处于充分发展阶段,图中展示了t*=20 时刻以后的数据.可以发现当=7.5 时,q*基本在0.10 和0.20 之间波动,而=9.375 的q*整体处于0.15 和0.20 之间,其沸腾传热性能提高.为了进一步分析引起该现象的原因,图10(a),(b)分别展示了在t*=69.32、t*=71.46 时刻两种加热器长度下的气泡形态演变、流场分布和气泡所受电场力的方向以及位置(图中白色箭头).从图中结果可以发现,当=7.5 时(如图10(a)所示)t*=69.32 时刻加热表面有3 个气泡生成,这3 个气泡在t*=71.46 时刻合并为一个大气泡.通过分析以上两个时刻的流场图和气泡的受力可以发现,虽然均匀电场具有横向压缩气泡阻止气泡在加热表面合并的作用,但由于加热台长度较小导致气泡间距很小,并且根据流场的方向得知气泡倾向于向加热表面中心部分聚集,最终仍然会出现气泡合并现象.气泡在加热表面合并会导致其施加在加热表面的换热热阻增大,使热流密度陡然降低.当=9.375 时(如图10(b)所示),t*=69.32时刻加热表面同样生成3 个气泡,由于此时加热台长度增大,生成的3 个气泡之间的间距较大,在电场力作用下这3 个气泡在t*=71.46 时刻并未发生合并,并且在加热器两端生成又生成两个新的气泡,此时沸腾换热热阻较小,并且由于新气泡的产生,沸腾换热潜热增加,从而使q*维持在较高水平.由此可以得出结论: 当 6.25＜≤9.375,更长的加热器为加热表面气泡的生成提供了充分的空间,均匀电场作用下的气泡间距增大、气泡数量增加.加热表面气泡不易合并,降低换热热阻,同时增加了气化核心,提高了潜热换热.

图9 均匀电场强度 E=0.10884 作用下,Tb=0.995Tc沸腾过程的空间平均热流密度随时间的变化Fig.9.Spatially averaged heat flow density with time for Tb=0.995Tc boiling process under uniform electric field strength E=0.10884.

图10 均匀电场强度 E=0.10884 作用下,Tb=0.995Tc沸腾状态的气泡形态演变、当前时刻气泡所受电场力和当前时刻流场分布 (a) =7.5;(b) =9.375Fig.10.Evolution of bubble morphology,electric field force on the bubble at the current moment and flow field distribution at the current moment under the action of uniform electric field strength E=0.10884,Tb=0.995Tc boiling state: (a) =7.5;(b) =9.375.

图11 均匀电场强度 E=0.16327 作用下,Tb=0.995Tc沸腾过程的空间平均热流密度随时间的变化Fig.11.Spatially averaged heat flow density with time for Tb=0.995Tc boiling process under uniform electric field strength E=0.16327.

图12 均匀电场强度 E=0.16327 作 用下,Tb=0.995Tc沸腾状态的气泡形态演变、当前时刻气泡所受电场力和当前时刻流场分布 (a) =12.5;(b) =18.75Fig.12.Evolution of bubble morphology,electric field force on the bubble at the current moment and flow field distribution at the current moment under the action of uniform electric field strength E=0.16327,Tb=0.995Tc boiling state: (a) =12.5;(b) =18.75.

根据前面的分析和讨论,对于亲水表面,均匀电场对不同长度加热器池沸腾的影响不尽相同.当≤6.25 时,由于加热表面产生气泡较少,均匀电场对池沸腾的影响表现为减小气泡尺寸从而抑制沸腾.当 7.5≤≤18.75 时,均匀电场总是能够提高CHF,且CHF 随着电场强度的增大而增大.另一方面,相同均匀电场强度对不同长度加热器池沸腾的CHF 提高的百分比也有较大的差异.≤9.375 时,长度的增大为电场力作用下的加热表面气泡更大的间距提供了空间条件,气泡不宜合并,均匀电场对CHF 提高的百分比逐渐增大.＞9.375 时,过大的加热器长度具有更长的水平流线区域,均匀电场作用下此区域蒸气容易贴附于加热表面,沸腾换热热阻增大,导致CHF 提高的百分比降低.

4.3 均匀电场作用下气泡行为的变化对疏水表面沸腾传热的影响

本节研究疏水表面均匀电场作用下,加热器长度对池沸腾传热性能的影响.为了与亲水表面的结果做对比,这里加热器长度和均匀电场强度与上节亲水表面情况的数据相同,而疏水表面润湿性取θ=120°.

图13 为疏水表面不同长度加热器在不同电场强度下的沸腾曲线.观察发现当=3.125 时均匀电场对整个沸腾阶段均有抑制作用,当≥6.25时,均匀电场对较高过热度下的沸腾有促进作用,并且其此时CHF 随电场强度的增大而有所均有提高,这与亲水表面得到的结果类似.不同的是当=6.25 时亲水表面的CHF 被抑制(图5(b)),而疏水表面的CHF 被促进.为了分析引起这一现象的原因,图14 展示了=6.25时,加热温度Tb=0.99Tc的疏水表面气泡动力学行为.从图中结果可以发现当电场强度为零时,由于疏水表面相变速率较大[33],t*=80.39 时刻形成了一个占据大部分加热表面的气泡,蒸气施加在加热表面的热阻较大,导致热流密度减小.当E=0.16327,在相同时刻加热表面的气泡在电场力作用下发生横向挤压形变,并在两侧生成两个小气泡,成核点的增加进一步使沸腾换热性能增强,因此CHF 被促进.

图13 疏水表面不同长度加热器在不同电场强度下的沸腾曲线 (a) =3.125;(b) =6.25;(c) =7.5;(d) =9.375;(e) =12.5;(f) =18.75Fig.13.Boiling curves of different lengths of heaters on hydrophobic surfaces under different electric field strengths: (a) =3.125;(b) =6.25;(c) =7.5;(d) =9.375;(e) =12.5,(f) =18.75.

图14 =80.39 时刻,Tb=0.99Tc时 (a) E=0;(b) E=0.16327 的疏水表面对应的气泡动力学行为Fig.14.The bubble dynamics behavior of the hydrophobic surface: (a) E=0,(b) E=0.16327 at the =80.39,Tb=0.99Tc.

图15 给出了不同电场强度下疏水表面加热器长度与CHF 的关系.可以发现在同一加热器长度下,CHF 随着电场强度的增大而增大.相同电场强度下,当＜12.5,CHF 提高百分比随着长度的增大逐渐增大,而当≥12.5,随着长度的增大逐渐减小.对比图8 发现,亲水表面和疏水表面关于加热器长度与CHF 的关系具有相同的趋势,不同的是CHF 提高百分比的最大值所对应的加热器长度不同,如E=0.10884 时,亲水表面的阈值为=9.375,疏水表面的阈值为=12.5.另外,均匀电场对疏水表面CHF 的提高程度远高于亲水表面,比如当E=0.16327 时,所有列举的加热器长度,亲水表面最大提高25.30%,疏水表面最大提高73.07%. 这是因为一方面疏水表面的CHF 本身就小于亲水表面;另一方面,疏水表面的气泡基底比亲水表面大,气泡容易在加热表面合并,而均匀电场提高CHF 的机理为通过施加于气液界面的电场力阻止气泡合并,使沸腾在更高过热度下达到CHF.这两个原因导致均匀电场对疏水表面的CHF 提高程度大于亲水表面.

图15 不同电场强度下疏水表面加热器长度与临界热流密度(CHF)的关系Fig.15.Relationship between hydrophobic surface heater length and critical heat flow density (CHF) at different electric field strengths.

5 结论

采用格子Boltzmann 直接模拟数值方法研究了均匀电场作用下不同加热器表面润湿性以及加热器长度的池沸腾传热过程.重点研究了气泡行为变化、电场强度、加热器长度以及临界热流密度之间的关系,并通过气泡动力学行为、气泡所受电场力和流场分布对其内在机理进行综合分析,得出以下结论.

1) 无电场作用时,在沸腾过程中生成的气泡之间的相互作用随着加热器长度的增大先增大再不变,对应的临界热流密度也随着加热器长度的增大先降低再趋于稳定,且亲水表面临界热流密度保持稳定时对应的加热器长度值小于疏水表面的情况.

2) 当存在均匀电场时,不论是亲水表面还是疏水表面,在核态沸腾阶段初期,电场的存在会小幅度降低热流密度.然而在核态沸腾阶段的中后期,电场的存在会增大热流密度,且疏水表面电场作用增大的热流密度比例大于亲水表面的情况.

3) 随着加热器长度的增大,均匀电场对亲水表面和疏水表面CHF 提高百分比均先增大后减小.电场强度E=0.05442,0.10884,0.16327 时,亲水表面CHF 分别在=12.5,9.375,9.375 达到增大的最大百分比,而疏水表面在=9.375,12.5,12.5 达到增大的最大百分比.

本文的研究在微观尺度揭示了电场力和加热器特性对池沸腾换热性能的影响机理,所得到的结论可为微电子器械加热器的设计提供理论基础.