APP下载

不同壁面取向下超疏水平面直轨道上的气泡滑移1)

2021-05-30叶煜航凃程旭包福兵汪钰锟杨森森

力学学报 2021年4期
关键词:气液方位角壁面

叶煜航 凃程旭,†,,2) 包福兵 汪钰锟 杨森森

∗(中国计量大学,杭州 310018)

†(浙江利欧泵业有限公司,浙江温岭 317500)

∗∗(浙江大学,杭州 310027)

引言

气液两相系统在自然界和工业过程中广泛存在,包括细胞孵化[1]、矿物浮选[2]、电解[3-6]等.例如在电解过程中基底表面的气泡堆积会导致能耗提高[4-6],因此需要降低气泡在壁面的停留时间.相反,矿物浮选过程中则希望气泡与矿物浆液充分接触使其表面附着更多的待选矿物[2].气泡运动速度的大小影响着化学反应器中的反应效率[7].因此,气泡在壁面附近的运动模式及其操控对上述过程至关重要.

近来,具有特定浸润性几何分布的表面被用于气泡的定向输运,其相关研究获得高度关注,工程应用前景广阔.截至目前,在水环境中主要通过Laplace压差或者浮力两种方式实现气泡自发的定向输运[8].近年来的许多研究中,科研人员通过具有壁面润湿性的几何梯度构建Laplace 压差来实现气泡的输运,具体壁面形式为超疏水铜锥[9-12]、梯形超疏水轨道[13-15]、楔形线性阵列超疏水轨道[16]、超疏水非平行双轨道[8].除了超疏水双轨道外,气泡总是从轨道窄的一端自发移动到宽的一端.在一定几何梯度下的超疏水表面能够实现气泡的抗浮力输运[10,13,17].然而,基于Laplace 压差法定向输运气泡仅限于短距离(∼O(10 mm))低速(∼O(10 mm/s))的定向运输[9,14],未能解决气泡长距离高速输运问题.基于浮力法则可实现气泡的长距离定向连续输运,超疏水轨道可以是直线[18-19]、S 形[20]、螺旋形[21].当壁面的倾角较小时(θ <30◦),气泡的滑移速度较小.随着壁面的倾斜度增大,浮力的作用显著提升,气泡在轨道上的滑移速度增大.当壁面倾角较大(θ >30◦)时,尽管存在局部浮力效应[22],但壁面的超疏水性仍促使气泡在水中滑动速度比气泡自由上升的速度快[18],在不同疏水性的倾斜下表面,毫米级气泡(直径1.79∼3.06 mm)的最大滑移速度发生在60◦θ70◦[18].对于垂直疏水表面(120◦CA125◦),Jeong 和Park[23]的研究发现,在很小的初始距离下气泡可以在这些表面上滑动而不逃逸,但气泡不能在渗透壁或弱亲水性表面黏附和滑动,这与最近的数值结果[24]相矛盾.

上述气泡滑移现象的物理本质是气泡与周围流体的相互作用伴随三相接触线的钉扎和解钉[25-27].接触线在固体表面上的钉扎主要是由于受固体表面粗糙度或表面化学不均匀性的影响[25].对于自由上升气泡,气泡的形状决定了气泡滑移的终了速度,而气泡形状由浮力、形状阻力、表面张力等力的平衡所确定[24].相对于自由上升气泡,贴壁气泡需额外考虑壁面剪切力、接触线钉扎等因素.值得注意的是,基于浮力的简单模型(FB=ρgV)不适用于倾斜壁面的气泡动力学的精确分析[22].目前,倾斜超疏水壁面的倾角对气泡滑移运动的相关研究较少[18],超疏水壁面取向对气泡滑移的影响尚不清楚,而壁面取向是诸多工程结构的重要参数,如精馏塔翅片.为此,本文采用高速阴影成像系统对不同壁面取向及倾角下气泡在超疏水轨道上的运动特性进行了实验研究(雷诺数Re=500 ∼700,韦伯数We=7 ∼13),以期为该气泡操控技术在核反应堆、矿物浮选、池沸腾等领域的应用提供科学依据和创新思路.

1 实验方法

1.1 超疏水直线轨道(SHBT)的制作

在洁净的玻璃基底表面贴上预留有后续喷涂区域的掩膜,将纳米SiO2溶液(Glaco,Soft99 Co.,Japan,SiO2初始粒径为30 nm,溶剂为异丙醇)喷涂在玻璃基底上,多次喷涂烘干后在玻璃基底表面形成SiO2纳米涂层,由此得到超疏水轨道.如图1(a)所示,玻璃基底表面附着有一层致密的纳米二氧化硅涂层,表面粗糙度约为100 nm[28].SHBT 表面存在大量的微孔结构,夹带着气体形成气穴和薄气膜,水滴沉积其上表现出超疏水性,而在水环境中呈现出超亲气性,因此超疏水表面在水环境中对气泡具有极强的附着力.详细的超疏水轨道制作方法参见Tu 等[29]最近的研究结果.利用上述方法得到的超疏水轨道其表观接触角如图1(b),约为162◦,滚动角<2◦.

图1 超疏水轨道实物图Fig.1 Physical picture of superhydrophobic trajectory

1.2 实验装置

本研究采用的气泡生成装置和高速阴影成像系统如图2 所示.由注射泵(Harvard Pump 11 Elite,Hollision,MA) 推动注射器往不锈钢喷嘴(喷嘴内径为0.26 mm)注气,从而在水箱底部产生气泡.注射器与喷嘴通过聚氨酯圆管(PU 管)连接,注射泵注射流量设定为0.1 mL/min.实验采用纯净水,水箱为亚克力透明水箱,且其中一个侧面为倾斜面(倾角α(图2)).

图2 实验装置图Fig.2 Experimental setup

气泡从喷嘴处自由脱离后具有高度轴对称性,根据此刻的气泡投影轮廓计算得到气泡的等效直径(Deq) 是2.4 mm、体积(V) 为7.3 µL.高速阴影成像系统由两台高速CCD 相机(FASTCAM Mini UX,Photron)、微距镜头(AT-X Pro,Nikon)和LED 阵列组成,从两个正交视角同步记录气泡的上升,高速相机帧率为2000,相应分辨率为1280×1024 像素.一台相机视场平行于SHBT 中心线且在竖直面内(相机I),另一台相机视场与水箱倾斜面平行(相机II),结合LED 阵列和匀光板构成双视角的高速阴影成像系统,捕捉气泡在超疏水轨道上的滑移.

具有超疏水轨道的玻璃基底通过固定架与转台固定,g为重力加速度(图3),图中y轴沿轨道中心轴线方向(斜向上为正),x轴在水平面内且垂直于y轴(当基底平面为竖直平面时,x轴恰好与基底的法线方向重合),z轴与xy平面垂直,m轴为xz平面与基底平面的交线,玻璃基底平面法线n位于xz平面内,与x轴的夹角为β,通过转动转台使超疏水轨道绕其中心轴线转动,改变轨道表面取向,即改变β(方位角,−90◦β90◦),间隔15◦取一个工况.β>0◦气泡位于轨道上侧,β<0◦气泡位于轨道下侧.图4(a)∼图4(c)为相机I 在3 个重要方位角(−90◦,0◦,90◦)所拍摄的典型气泡图像.图4(a) β=−90◦,气泡位于轨道下方;图4(b) β=0◦,气泡位于竖直面内;图4(c)β=90◦,气泡位置位于轨道上方.转台固定于轨道支撑,使得图3 中轨道与水平面夹角为α 角(轨道倾角,α=45◦,60◦,65◦,70◦,75◦).

图3 倾斜超疏水轨道支撑装置示意图Fig.3 Schematic diagram of inclined superhydrophobic trajectory fixation

图4 α=70◦时3 个重要方位角的典型气泡图像(相机I 测得)Fig.4 Typical bubble images captured using Camera I in three important azimuths at α=70◦

超疏水轨道底端在竖直方向上距喷嘴顶端约为4 mm.浮力作用下,气泡从喷嘴口完全脱离后,在壁面与气泡之间的流体开始以一定的速度流动,此时周围流体产生的正压力驱动气泡向壁面迁移[23,30],气泡与壁面之间的流体变薄并最终发生破裂,气泡受到壁面的气膜黏性效应影响附着于轨道表面.实验过程中每个工况重复5 ∼10 次,终了速度误差范围在10%以内.

2 实验方法

2.1 气泡稳定性

针对气泡在不同方位角的SHBT 上滑移的稳定性,本文从气泡形状和气泡上升速度两方面进行分析.

超疏水轨道表面存在一层薄气膜,长期置于水中会造成气膜不稳定导致轨道失效[18,31],此外轨道放入水中的浸没速度对轨道表面气膜的形成至关重要[32],而在本文的实验过程中,连续有气泡在SHBT上滑移,其稳定的终了速度表明气膜的厚度是较为稳定的.由于轨道表面铠甲状气膜的存在[28]和三相接触线在轨道表面的钉扎,使得气泡稳定附着于壁面.当气泡所受浮力克服壁面对气泡的钉扎,气泡将沿着壁面开始滑移[25].SHBT 特殊的表面微观结构(CAH(∆θ <2◦)) 可以较好维持其水下壁面气膜的稳定,极小的亲气角(等价较大的疏水角) 表明其轨道的表面能极低,使得气泡在浮力的作用下能够轻松克服表面缺陷能垒而滑移[33],而滑移的过程实际是三相接触线的滑移−钉扎状态(stick-slip 状态).

如图5,当V=7.3µL 的气泡在W=2.0 mm 的超疏水轨道上运动,气泡整体呈现为具有多个峰脊的半子弹型,Re=500 ∼700.图5 为α=70◦,β=90◦,0◦,−90◦,不同时刻下气泡的阴影成像轮廓(图下侧轮廓线分别为4 个时刻(t1=0 ms,t2=12 ms,t3=24 ms,t4=36 ms)气泡轮廓的叠加,黄色线为壁面).由于气泡表面毛细波的存在,气泡在运动过程中呈现出多脊型[34],根据气泡的气液界面是否存在显著波动可分为波动型和稳定型.波动型如图5(a)和图5(b)所示,气泡在运动过程中,气液界面存在剧烈波动,图中气泡轮廓线重合度较低,而稳定型在气泡运动过程中,形状基本保持稳定,不同时刻下气泡轮廓线高度一致(图5(c)).当轨道倾角较低时(α65◦),随着β从±90◦逐渐接近0◦,气液界面由稳定状态向非稳定状态转变,失稳的临界点均在∼30◦|β|45◦.随着α 的增大,失稳临界点向|β| ∼90◦逐渐靠近,气泡在SHBT 上方(β >0◦)较之气泡在SHBT 下方(β <0◦)更易产生波动,这可能与气液界面在m轴方向受到的静压大小有关.当β=±90◦时,气液界面所受的静压在xz平面上的分量方向与壁面法线方向重合,m轴方向受到的静压分量为0,气液界面较为稳定.随着β 从±90◦逐渐接近0◦,xz平面上的静压分量的方向偏离壁面法线方向,m轴方向受到的静压分量增加,气液界面受力对称性破坏,导致气液界面逐渐失稳.当气泡在轨道正上方滑移时(β=90◦),气泡在轨道法向上的高度(类似于粗糙元特征尺度)大于气泡在轨道正下方(β=−90◦)滑移时的高度,使得后者的气泡形状稳定性高于前者(表1).

图5 SHBT 轨道上,不同时刻气泡形状(α=70◦,t1=0 ms,t2=12 ms,t3=24 ms,t4=36 ms)Fig.5 Bubble shapes at different times in the SHBT trajectory(α=70◦,t1=0 ms,t2=12 ms,t3=24 ms,t4=36 ms)

表1 β=±90◦不同α 下的气泡高度(mm)Table 1 When β=±90◦,the height of bubble under the different α(mm)

通过监测滑移过程中气泡质心的位置随时间的变化,可以确定气泡在超疏水轨道上运动的稳定性.如图6 给出了α=60◦,不同β 下气泡在超疏水轨道上滑移距离随时间的变化,以气泡完整出现在相机视场中作为统计气泡位置的初始时刻.图中纵坐标为气泡沿轨道方向上质心相对于初始时刻位置的位移,横坐标为气泡在轨道上的运动时间.无论是自由上升气泡还是在SHBT 上滑移的气泡,其位移与运动时间(s–t)二者强线性相关(图6),线性函数的斜率(即气泡的运动速度)随方位角的不同而不同,但都大于自由上升气泡.s–t的高水平线性相关说明气泡沿轨道的终了速度波动很小,气泡近似处于受力平衡状态.Jeong 和Park[23]的研究成果也表明气泡在竖直均匀疏水壁面的滑移速度基本恒定.

图6 α=60◦,不同β 下,气泡沿轨道方向滑移距离s 随时间t 变化图Fig.6 α=60◦,the sliding distance of bubble along the trajectory direction varies with time under the different β

2.2 受力分析

气泡在超疏水轨道上运动时受到浮力、形状阻力、轨道表面分裂力(源于接触线在SHBT 表面上的钉扎)及附加质量力的影响.

壁面的存在会导致气泡的局部浮力效应[22],浮力的本质是静水压差,可将静压对气泡的作用力分为竖直作用力FV和水平作用力FH.竖直方向静压力FV是在气液界面A上对静压二重积分,可利用高斯公式,将其转化为气泡的体积分,并减去轨道B面所受竖直方向的静水压力,所以竖直方向静压力FV

方程右边第一项表示等效直径气泡所受竖直方向上的静水压力,第二项表示轨道气固界面(面B) 所受竖直方向的静水压力.式中,ρ 为水的密度,g为重力加速度,为9.81 m/s2,V为气泡等效体积,W为轨道宽度,θf为曲面微元dS法线方向与竖直方向(向上)的夹角,L为气泡与轨道单侧接触线长度.

静水压力对气液界面A在水平方向上的作用力为FH

CD为阻力系数(包括形状阻力和黏滞阻力);Af为迎风面积,实验结果表明气泡滑移时的迎风面积与超疏水轨道的倾角及方位角高度相关.

因为超疏水表面的低表面能,接触线钉扎解钉所需的轨道表面分裂力∆F可采用[25,35]

式中,κ 为欧拉常数,λ 为典型缺陷尺寸,fc为逃脱能垒的临界力,a/λ 为缺陷线密度,kB为玻尔兹曼常数,T为弛豫率,fT为热力,K0为尝试频率,Eb为潜在的能量势垒,ks为气液界面的弹性系数.

对于静止液体中加速的气泡,附加质量力[23-24]如下

其中,Cm为附加质量系数,ρb为气体密度,对于扁球形气泡,Cm取决于气泡的纵横比.

综上,SHBT 上滑移气泡的受力模型如图7 所示,忽略气泡受到的重力,气体与壁面间的黏性阻力(nN量级) 远小于驱动力及形状阻力(约50 µN),从而忽略不计.

故沿轨道方向气泡受力方程为

如前文所述,气泡在轨道上可近似为匀速运动,因此附加质量力可以忽略,此外超疏水轨道滚动角较小,接触角滞后力可忽略不计,接触线耗散弱,轨道表面分裂力量级较小(约0.1 µN)[36].因此式子(7) 可简化为

故气泡运动主要受(气液界面静压产生的)驱动力及形状阻力影响.联立方程式(3)、式(4)和式(8),可得气泡滑移速度与迎风面积的关系

由实验结果反映的气泡迎风面积与倾角和方位角的相关性,结合式(4)和式(9)可知,倾角和方位角对气泡滑移速度的影响机制在于,倾角和方位角的改变导致了气泡迎风面积和有效浮力的变化,而后两者是气泡滑移速度的关键影响因素.

2.3 轨道倾角和方位角对气泡滑移速度影响

图8 给出了不同方位角下,气泡沿轨道方向的滑移速度随着α 增加的变化曲线.对于45◦α75◦,不同方位角下气泡的滑移速度相对于自由上升气泡(0.2 ∼0.3 m/s) 始终维持在一个较高的水平(0.53 m/su0.66 m/s).气泡等效直径Deq=2.4 mm,在轨道上的滑移速度最快可达0.66 m/s,明显大于其他类型超疏水表面的结果[18].Maleprade 等[37]采用尼龙绳控制气泡倾斜运动,速度也仅为0.1 m/s.Wu 等[18]使用的Cu2O 超疏水轨道(θ ≈160◦)和聚四氟乙烯超疏水轨道(PTFE,θ ≈175◦),粉末粒径约为74µm,轨道表面粗糙度为微米级,两种超疏水倾斜表面上的气泡最高速度分别约为0.42 m/s(Cu2O 超疏水轨道)和0.53 m/s(PTFE 超疏水轨道).Glaco Soft 99 溶液静置干燥形成的超疏水轨道,表面粗糙度约100 nm[28],本研究利用阴影成像测量经Glaco Soft 99 喷涂干燥的轨道表面所附着的气膜厚度在微米级,与Landau-Levich 公式计算的理论值[32]很接近.无论是静置还是喷涂得到的水平超疏水表面,气泡在其上铺展时,三相接触线的铺展演化过程都非常相似,最高运动速度都可达1.0 m/s[28-29,38].对表面微观形貌而言,由于Cu2O 轨道和纳米SiO2轨道表面粗糙度差异较大,使得Cu2O 轨道的钉扎作用强于纳米SiO2轨道[25-26,34],而且三相接触线在Cu2O 轨道移动时其在轨道法线方向上需克服的浸润距离(约102µm)[34]远大于纳米级的SiO2轨道,因此气泡在两种轨道上运动存在明显的速度差.相同大小的气泡在超疏水轨道上呈多脊型,气泡的迎风面积远小于单丝上滑移的类球形气泡[37].

图8 不同β 下,气泡滑移速度u 随α 变化图Fig.8 Bubble slip velocity changes with the trajectory inclination angle for different azimuth angles

此外,气泡滑移速度随α 增大的变化规律根据β的不同可归纳为3 种基本模态:

(1) 当β=−90◦,随着α 的增大,气泡受到的驱动力Fdriven增强,气泡沿轨道方向滑移速度增大,在α=70◦时,滑移速度达到最大值u=0.6 m/s.随着倾角继续增大,滑移速度下降.该特性与Wu 等[18]实验结果一致.β=90◦,气泡滑移速度随着α 的增大先增大后减小,其规律与β=−90◦相似,变化趋势更为剧烈.特别当β=±90◦,气泡滑移速度在α >70◦开始下降,图9 的气泡迎风面积数据表明这很可能是由于迎风面积的增大导致的.其余角度下的气泡迎风面积难以提取,仅能获取到β=±90◦,β=±45◦,β=0◦这5 个典型的角度下的气泡迎风面积.

图9 不同β 下气泡迎风面积与α 关系Fig.9 Relation between the frontal area of bubble and α under different β

(2)当β=0◦,气泡位于竖直面内,滑移速度随着倾角的增大而减小.显然也是受到气泡的迎风面积的影响(图9),随着α 的增大气泡迎风面积明显增大,形状阻力对气泡运动的影响占主导.

(3)当β=±45◦时,Re∼O(102),高We∼O(10)数的近壁气泡运动过程中,惯性力起主导作用,气泡开始出现形状不稳定[39],且气泡界面压力在垂直于运动方向平面(xz平面) 内的分力偏离壁面法线方向,破坏了气液界面对称性,气液界面产生波动.气泡的迎风面积Af随着轨道倾角的变化存在一定的波动(图9),与气泡的滑移速度呈负相关.

图8(b)表明±75◦β±90◦,气泡滑移速度随α 的变化规律与模态I 相似,0◦β±15◦,滑移速度随着倾角的变化与模态II 相似.±30◦β±65◦,气泡滑移速度随着α 的变化与模态III 相似.

有趣的是,不同β 气泡滑移速度的变化曲线,都在α=70◦处开始交汇.通过对比β=±90◦,β=±45◦,β=0◦下,气泡迎风面积随α 的变化规律(图9),结果表明在α65◦时,不同方位角下气泡的迎风面积存在较大的差异;当α70◦时,不同方位角下气泡迎风面积基本重合(图9;不同方位角下气泡迎风面积的比值参见附表2).因此,在α70◦时,β 对气泡滑移速度影响较小,此时不同β 下的气泡滑移速度非常接近.

气泡沿着超疏水轨道的滑移速度随着轨道方位角的增加呈现出2 种不同的变化趋势.不同倾角下,滑移速度随方位角的变化特性如图10 所示.

图10 不同轨道倾角下,气泡沿轨道滑移速度u 随β 变化图Fig.10 Sliding velocity of bubble along the trajectory varies with the azimuth angle for different trajectory inclination angles

第一种u-β 分布:当轨道倾角45◦α65◦时,气泡滑移速度随着β 从−90◦逐渐增大至90◦,其值先增大后减小,近似为关于β=0◦的单峰分布.α=45◦,β=0◦时,气泡的滑移速度最大(约0.66 m/s).较小倾角下(α<70◦),β 直接决定了气泡滑移时的基本形态,不同的方位角下气泡的形态不同(图11,黄线为壁面,正视图(相机I)).气泡的迎风面积受β 的影响(图9),随着β=±90◦趋近于β=0◦,气泡的迎风面积逐渐减小.

图11 α=45◦,不同方位角下气泡形态图Fig.11 When α=45◦,the shape of bubble at different azimuths

第二种u-β 分布:当α70◦时,随着方位角变化,气泡的滑移速度围绕某一定值小幅波动.图12(黄线为壁面,正视图(相机I))中所示不同方位角下,气泡运动形状基本一致,迎风面积Af的值较为相近(图9),对方位角的改变不敏感.方位角的变化对大倾角(α70◦) 下气泡滑移速度的影响较小,且当α=70◦时不同方位角下的气泡滑移速度最为相近.

图12 α=70◦,不同方位角下气泡形态图Fig.12 When α=70◦,the shape of bubble at different azimuths

气泡沿轨道方向上的驱动力Fdriven是气泡滑移的主要动力,也是出现上述两种u-β 分布的关键因素之一.当α 恒定时,式(3)中气泡驱动力Fdriven中第一项ρgVsin α 为一定值,第二项ρgWL2f(α,β,θ1,θ2)/2随着β 的改变而改变,即随着β 从−90◦到90◦,第二项的值先减小后增大(拐点为β=0◦).因此,驱动力先增大,在β=0◦时达到最大值(ρgVsin α),然后减小.当β=±90◦,局部浮力效应[22]最明显,气泡驱动力最小.

气泡在小倾角(45◦α65◦)轨道上滑移,由于α 较小,第一项(ρgVsin α) 占比较小,所以驱动力中第二项的值对驱动力的大小影响较大,故气泡的u-β分布在小倾角下表现为近似于驱动力关于β=0◦的(Fdriven-β)单峰正弦分布.对于大倾角(70◦α75◦)轨道,第二项随α 的增大单调递减(f(α,β,θ1,θ2)的值变化见附件表3),同时α 的增大使得第一项(ρgVsin α)的值增加,因此第二项的占比进一步减小,对驱动力大小的影响较小.当α=90◦(壁面竖直),驱动力的第二项为0,即方位角的变化对气泡滑移时的驱动力没有影响.另一方面,当α70◦,不同β 下气泡滑移时的迎风面积基本一致(见图9),表明大倾角下Af与β 无关.综述两方面,大倾角下气泡在轨道上滑移的驱动力和迎风面积在不同的β 下非常接近,这很可能是大倾角下气泡滑移速度随β 的增加基本保持不变的原因(第二种u-β 分布).气泡在超疏水轨道上的运动主要受浮力及形状阻力的影响,而气泡的形状阻力可以通过阻力系数(CD)进行量化.对于自由上升的气泡(Re∼O(102)),其CD≈0.25[40],而本文中超疏水轨道上运动的多脊型气泡阻力系数CD≈0.1 ∼0.3,变化范围较大,且不同方位角下的阻力系数与倾角存在明显的正相关(图13,红色虚线为自由上升气泡的CD).

图13 不同β 下,多脊型气泡CD 随α 增大的变化Fig.13 Relation of bubble CD with α under the different β

3 结论

本文采用具有正交视场的高速阴影成像系统,研究了毫米气泡(Deq=2.4 mm,Re=500 ∼ 700,We=7 ∼13) 在不同轨道表面取向下(方位角β),倾斜超疏水轨道(W=2 mm)上的滑移特性,重点分析了方位角(−90◦β90◦)和倾角(45◦α75◦)对气泡形态和滑移速度的影响.结果表明,气泡黏附到超疏水轨道上后会经历短暂的加速,而后迅速达到近似匀速的三相接触线钉扎−解钉运动状态.由于表面毛细波的存在气泡滑移时形状为多脊型,根据气液界面波动程度的不同滑移气泡可分为波动型和稳定型.稳定型气泡只出现在较小倾角且较大方位角的超疏水轨道上(α <70◦,45◦|β|).气液界面上的静压在m轴方向上的分量大小对气泡滑移过程中的形状稳定性有决定性的影响,当该分量为零(β=±90◦),气泡形状最稳定.当气泡在轨道正上方滑移时(β=90◦),气泡在轨道法向上的高度大于气泡在轨道正下方(β=−90◦)滑移时的高度,使得后者的气泡形状稳定性高于前者.

方位角对气泡滑移速度的影响存在2 种趋势:

(1) 45◦α65◦,气泡滑移速度近似为关于β=0◦的单峰分布(β=0◦时,气泡滑移速度最大),此时沿轨道方向的驱动力起主导作用;

(2) 70◦α75◦,气泡滑移速度对方位角的改变不敏感,此时气泡的迎风面积是影响滑移速度的主要因素.

综上,壁面取向(方位角β)及轨道倾角(α)通过改变气泡沿轨道方向的驱动力和气泡滑移时的迎风面积影响气泡在超疏水表面上的滑移速度和气液界面稳定性.这一新的认识有望为该气泡操控技术在核反应堆、矿物浮选、池沸腾等领域的应用提供科学依据和创新思路.

猜你喜欢

气液方位角壁面
二维有限长度柔性壁面上T-S波演化的数值研究
高温壁面润湿性对气层稳定性及其壁面滑移性能的分子动力学研究
运载火箭气液组合连接器动态自动对接技术
壁面滑移对聚合物微挤出成型流变特性的影响研究
微重力下两相控温型储液器内气液界面仿真分析
气液固多相流对法兰接缝处的腐蚀行为研究
基于停车场ETC天线设备的定位算法实现
壁面喷射当量比对支板凹腔耦合燃烧的影响
基于新型C4D的小管道气液两相流流型辨识方法
无处不在的方位角