体相微马达双气泡聚并驱动的界面演化机制.pdf

微纳米马达是建立微流体环境与宏观操控的桥梁,气泡微马达的驱动速度高,这一优势在实际应用中不可替代.管式气泡马达适用于复杂场景但能量转化率低,气泡驱动的 Janus 微球马达效率高但仅适用于气液界面附近.鉴于此,本文提出通过双气泡聚并方式驱动 Janus 微球马达的新体系,调和了高能量转化率与界面受限的矛盾.在实验中,借助高速摄像记录了双气泡聚并驱动微马达的100s 级过程,气泡聚并紧邻微球发生,通过释放的能量驱动微球显著运动,其融合过程是独特的可动曲壁受限下的气液界面演化问题.进一步结合伪势格子 Boltzmann 数值方法探究了气泡聚并驱动的流体动力学机制.研究结果揭示了不同时段气泡聚并的细节,给出了气泡颗粒尺寸比等因素对微球位移、初始动能转换率的影响,确认了双气泡聚并释放表面能的高效驱动机制.关键词：双气泡聚并,微马达,界面能,伪势格子 Boltzmann 方法PACS：47.63.mf,07.10.Cm,65.40.gpDOI:10.7498/aps.72.202306081引言微马达因其体积微小,更易执行精细化的操作,在微机器人操控1、靶向药物运输2和生物传感3等领域具有广泛应用前景.微马达进行自驱动的形式多样4,5,如化学方式驱动、增加外部声场、光场6、电场、磁场驱动,以及两者相结合的混合方式等.由于化学驱动方式的实验形式简洁,驱动效率相对较高,本文采用 Janus 微球生成气泡的化学驱动方法.Janus 微马达可以通过非对称的表面催化反应将周围环境溶液中的化学能转化为机械能,驱动自身运动.相比于其他利用物理场梯度的驱动方式(热泳7、扩散泳8等),气泡驱动具有更强的驱动力,可以实现更高的微马达速度.尽管化学方式自驱动微马达的能量都来自溶质分解所释放的化学能,并最终转化为热能耗散掉,但不同驱动机制的能量转化形式和中间路径决定了微马达的能量转化率.理论上,当研究问题的尺度从宏观尺度下降至微纳尺度,体积力的影响会大幅下降,而表面力及线力会居于主导,这一尺度效应在微马达的设计问题中仍适用,意味着为获得高的能量转化率,微马达应尽量利用储存在界面中的能量.尽管同为气泡驱动微马达,但现有研究中却存在着两类性质截然不同的气泡驱动模式:一是利用气泡与气液界面瞬间融合、界面完全消失所释放出的能量;二是利用完整的气泡反弹推动微马达,下面分别进行说明.H2O2Janus 微球一般采用电子束蒸发镀膜制备,可以大批量获得.将 Janus 微球置于低浓度 溶*陕西省自然科学基础研究计划(批准号:2020JM-479)、中国博士后科学基金(批准号:2022M720148)和国家自然科学基金(批准号:52208126)资助的课题.通信作者.E-mail:通信作者.E-mail:2023中国物理学会ChinesePhysicalSocietyhttp:/物理学报ActaPhys.Sin.Vol.72,No.16(2023)164703164703-1=Eo/Et107102液中,当微球粒径大于气核生成临界粒径(约 20m)时,借助 Pt 表面的催化反应会析出过饱和氧气泡,发生微气泡推进9.王雷磊等1通过实验揭示了如何利用 Janus 微球实现气泡融合驱动微马达.这个过程发生在近气液界面处,随着气泡上浮或由于气液界面的薄膜蒸发,这个气泡会与气液界面融合,释放出界面能量,推动微球运动.Wang 等1012也观察到了类似的气泡推进方式.经计算,这种驱动方式的能量转化率(,其中 Eo为输出机械能,Et为总机械能)约为 量级13.有趣的是,在自然界中也观察到了类似的利用两相间界面能的力学机制.例如,苔藓类植物14利用囊腔中的高压气体、蕨类植物15利用“空化弹弓”发射孢子,它们的能量转化率则要高得多,约为 量级.1010另一类气泡驱动微马达为经典的管式微马达,它采用气泡反弹的驱动方式16.微管结构可通过应变诱导纳米薄膜的自卷曲、基于模板的电沉积或双光子激光光刻制备17,工艺相对复杂,制备的中空微管结构两端通常是几何不对称的18.不同于Janus 微球生成的气泡位于球体外部,管式微马达Pt 催化层在微管内壁,气泡形成于管内,随着气泡不断地生成,后续气泡会排挤前面的气泡从微管的一端释放到溶液中19.这里,管式微马达对内部气泡的保护作用要远大于 Janus 微球对外部气泡的作用,意味着管式马达的工作模式可以不受外部环境的限制,适用于复杂的溶液场合.实验还观察到,在推动管式微马达后,气泡整体仍完整地存在于体相溶液中,形成长长的气泡尾链.而我们知道气泡大部分的能量贮存在界面中,因此管式微马达并未充分利用这部分界面能,这就限制了管式微马达能量转化率,Wang 等20计算了管式微马达的能量转换率约为 量级.可见,一方面现有人工微马达的能量转化率远不如自然界中的微生物,仍有很大的提升空间;另一方面,不同形式微马达的优缺点也非常明显.Janus 微球马达加工简单,可以更好地利用气泡界面消失所释放的能量,能量转化率较高,但气泡融合需在临近气液界面处才能触发;管式微马达气泡生成于微管内部,适用场景不局限于气液自由面,体相中甚至固液界面附近也能有效驱动,但其制备复杂,主要利用了气泡反推的能量,能量转化率偏低.综上所述,是否可以设计一种新的驱动模式,能够兼顾上述两种微马达的优势,使其在气液界面、体相、固液界面处均能高效驱动?鉴于此,基于 Janus 微球体系提出了新颖的双气泡聚并驱动方式,既释放了部分气泡界面能量,保证了较高的能量转化率,又不需要气液界面触发的气泡融合条件,并保有了简洁的微球结构.另外,双气泡聚并驱动的方式还减小了气泡周期性生长、融合给流场带来的扰动,驱动过程更加稳定,定向性更好.这一新的驱动方式为微马达在复杂场景的应用提供了可能.本文借助高速摄像捕获了双气泡聚并驱动微马达的过程,通过观测两个气泡聚并前后的界面演化细节,结合伪势格子 Boltzmann 方法数值模拟,阐明了双气泡聚并驱动微球的流体动力学机制.2双气泡聚并驱动微马达实验 2.1 实验设置H2O22H2O2Pt 2H2O+O2H2O2H2O2H2O2H2O2VH2O2/VH2O实验中,气泡生成是通过 Janus 颗粒(直径在30100m 之间)一侧的 Pt 催化剂分解溶液中的()而产生的.溶解氧过饱和后析出产生氧气泡,如图 1(a)所示.借助倒置显微镜及高速摄像(PhantomV2512)记录实验过程.先滴加 反应溶液在载玻片上形成液膜,再蘸取 Janus 颗粒悬浮于液膜间.考虑到气泡在纯水中维持的时间极短难以拍摄记录,但甘油黏稠且不与 水溶液反应.本实验在 水溶液的基础上添加了甘油配制出了不同浓度的反应试剂,增大溶液黏度,延长聚并发生过程.实验中主要通过控制 的浓度(10%13%,)控制双气泡的生成.两个气泡依托成核点不断生长,界面相互接近,接触之后便会发生气泡聚并,此过程中颗粒被推开.气泡不断地生成、聚并、脱离颗粒,实现了对颗粒的驱动.rneck为了后续更好地对气泡聚并的演化过程进行描述,对双气泡的几何坐标进行定义,如图 1(b)所示.在初始时刻,两气泡接触点的位置为 O 点,长轴方向规定为 Z 方向,两气泡颈部生长的方向规定为 r 方向,生长长度的一半为颈部半径 .H2O2V甘油/VH2O2水溶液首先给出基本的气泡聚并过程的实验观察,如图 2 所示.两个直径分别为 88.0m 和 87.6m 的气泡在甘油 水溶液(30%,)中发生聚并.这里的气泡通过 Janus 微球产生,与物理学报ActaPhys.Sin.Vol.72,No.16(2023)164703164703-2rneck微球脱离后,稳定存在于溶液中,与溶液中另一个气泡相遇后,被记录下来.此实验给出了基础的无临近颗粒影响的双气泡聚并过程.图 2(a)为气泡聚并的初始时刻,从这个时刻开始,开始出现.图 2(b)中气泡聚并颈部均会受到黏性力、惯性力以及表面张力的共同作用,且两侧所受力相等,出现 Z 方向上严格对称的演化过程.当气泡颈部半径与气泡半径的比值为 1 时,可认为颈部生长过程完成.故图 2(c)中气泡聚并颈部的生长在 10s 已基本完成.图 2(c)和图 2(d)为气泡回缩为球形的振荡过程.Rb1Rb2E=4R2bRb4.84 1084.14 108E对气泡聚并前后的表面能变化进行分析.表面能根据两气泡半径 ,以及所处溶液的表面张力系数g 计算得到:,表示气泡半径,其中表面张力系数的测定采用受限液滴表面张力仪.初始两气泡表面积之和为 m2,聚并完成后单气泡表面积变为 m2;双气泡聚并后表面积减小,表面能降低.气泡表面能 的释放量为聚并前的 14.6%.对气泡聚并前后不同尺寸气泡的表面能释放量进行统计,如图 3 所示,表面能释放量随气泡半径呈二次函数关系.这里数据离散的可能原因在于理论推导中按照水的表面张力系数进行计算,而实验中溶液的初始配比及持续的化学反应使其浓度发生变化,造成实际表面张力系数会存在一定偏差.H2O2实验实验考虑气泡聚并过程与临近颗粒的相互作用情况.图4(a)(c)分别为半径24m 和25m的气泡在 水溶液中的聚并演化过程.在 026s 气泡已经完成了颈部的生长;在图 4(d)和图 4(e)过程中颗粒在气泡的界面振荡中被弹开,明场光源控制端过氧化氢水溶液气泡Janus颗粒2H2O2 2H2O+O2Pt电源高迼相机PhantomV2512(a)(b)气相液相neck图1(a)实验装置示意图;(b)气泡聚并形状演化的方向定位Fig.1.(a)Schematicdiagramoftheexperimentalsetup;(b)orientationoftheevolutionofbubbleaggregationandshape.(a)(b)(c)(d)图2气泡聚集与界面演化的时间序列(a)0;(b)4.34s;(c)10.5s;(d)43.39sFig.2.Timeseriesofbubbleaggregationandinterfaceevolution:(a)0;(b)4.34s;(c)10.5s;(d)43.39s.02.04.06.08.000.51.01.52.0=2.89T10-3 PaSs理论值=2.06T10-3 PaSs=1.00T10-3 PaSsD/(10-9 J)b/(10-5 m)D=0.3b2图3气泡聚并后表面能释放量随气泡半径变化图Fig.3.Variation of surface energy release with bubbleradiusafterbubbleaggregation.物理学报ActaPhys.Sin.Vol.72,No.16(2023)164703164703-38.1 1011109出现了明显位移,同时气泡也被弹开;在图 4(f)中聚并后的气泡仍存在溶液中,对体相溶液的扰动较小.经计算,图 4 实验中气泡表面能的释放量约为 J,其中一部分用于推动颗粒,转化为颗粒动能,能量转化率大约在 量级.H2O2V甘油/VH2O2水溶液实验实验进一步考虑在实际相关应用中,会出现气泡在较高黏度溶液中聚并的情况,如人体血液的黏度约是水的 45 倍.本组实验将半径为36m 的 Janus 颗粒置于甘油 水溶液(10%,)中,生成两个半径分别为 47m和 49m 的气泡(图 5).与图 2 相比较,增加了颗粒对气泡一侧界面演化的约束,在 5s 时气泡颈部的生长基本不受颗粒的影响,见图 5(b);而在10s 处,靠近颗粒一侧的颈部半径明显小于另一侧,说明在水动力作用下颈部的生长已经受到了颗粒的影响,见图 5(c).然而,参考黄色标记线,可以发现颗粒并没有显著的移动.净位移主要出现在图 5(d)(f)中,气泡界面在 Z 方向回缩,r 方向上拉伸.界面在惯性力作用下在 r 方向拉伸的同时,也受表面张力的作用,在不断收缩,使气泡趋向于稳定状态.从相对于绿色的标记线位置可以看出,这里气泡的反弹位移不明显,可以忽略.当时间达到 200s 后,推动颗粒的净位移为25.3m,这个时间段内,颗粒的平均速度高达0.12m/s.气泡表面能的释放量为聚并过程未发生时的 11.2%.由于增大溶液黏度后,气泡聚并后与颗粒弹开需要的时间较长.区别于之前的能量转化率,定义“初始动能转换率”专门来评价气泡聚并推动颗粒的能量利用率,即颗粒与气泡未发生脱离之前,气泡聚并赋予颗粒的动能 Ev与气泡聚并前后表面能差值 DEs之比,即i=Ev/E,(1)Ev=mv2/2,(2)E=4(R2b1+R2b2 R2b3),(3)Rb1Rb2Rb3其中,m 为颗粒质量,v 是气泡聚并开始到结束时间段内的颗粒最大驱动速度,分别为聚并前两气泡的尺寸,为气泡聚并后的尺寸.为了与前人对能量转化率的计算相对比,能量转化率的表达式为e=Ev/E,(4)Ev=mv2/2,(5)vEv其中,表示从气泡生长到聚并结束时间段内的颗粒平均速度,是其对应的颗粒动能,Es为气泡聚并前最大表面能.能量转化率是从微马达的运动周期来考虑,包含了催化反应气泡生长阶段.气泡生长的时间持续较长,大约为 40ms,目前关注于(a)50 mm(b)(c)(d)(e)(f)图4Rb1=24m,Rb2=25m,Rp=27m 时,实验拍摄的气泡聚并驱动颗粒运动的时间序列(黄色虚线标示初始时刻颗粒的位置,绿色虚线标示初始时刻聚并气泡中心)(a)0;(b)13.3s;(c)26.6s;(d)39.9s;(e)66.6s;(f)533.0sFig.4.ExperimentallyfilmedtimeseriesofbubblesaggregatinganddrivingparticlemotionatRb1=24m,Rb2=25m,Rp=27m(Yellowdashedlinemarksthepositionoftheparticleattheinitialmoment,andthegreendashedlinemarksthecenteroftheaggregationbubbleattheinitialmoment):(a)0;(b)13.3s;(c)26.6s;(d)39.9s;(e)66.6s;(f)533.0s.(a)50 mm(b)(c)(d)(e)(f)图5Rb1=47m,Rb2=49m,Rp=36m 时,实验记录气泡聚并驱动颗粒运动序列图(a)0;(b)5s;(c)10s;(d)20s;(e)30s;(f)200sFig.5.ExperimentallyrecordedbubbleaggregationanddrivenparticlemotionsequencediagramatRb1=47m,Rb2=49m,Rp=36m:(a)0;(b)5s;(c)10s;(d)20s;(e)30s;(f)200s.物理学报ActaPhys.Sin.Vol.72,No.16(2023)164703164703-4103109气泡聚并过程界面形变带来的影响,颗粒平均动能无法准确表达文章主题.而初始动能转换率是其中气泡聚并界面形变让颗粒产生位移短时间内的转换率.由此计算得出,图 5 所示过程的初始动能转换率约为 量级,能量转化率约为 量级.实验实验考虑气泡与颗粒尺寸差异较大的情况.图 6 中颗粒的半径为 18m,气泡的半径分别为 7.8 和 10.4m,此类实验中颗粒的尺寸比气泡尺寸大得多,颗粒没有发生明显位移,而气泡被弹开.将气泡颗粒作为一个系统,若不考虑附近溶液对系统的能量耗散,这个系统机械能守恒,气泡聚并释放能量,气泡表面能转化为气泡动能,使气泡弹离颗粒.根据动量守恒,由于颗粒质量远大于气泡,故颗粒几乎未被推动,气泡被明显弹开.在图 6(a)(e)过程中,虽然没有推动颗粒,但由于气泡大小不一,较大气泡侧对颗粒的作用力更大,颗粒仍具有朝向左上的作用力,在图 6(f)中,气泡自身不断形变从颗粒脱离、弹开.对比实验 II 可见,当气泡与颗粒的相对尺寸小于某个临界值,便无法实现自驱动.|Rb1/Rb2 1|0.15Rb=(Rb1+Rb2)/2从实验实验 III 可以看出,气泡聚并过程对颗粒的驱动作用与它们的相对尺寸密切相关.首先,对比图 5 与图 6(a)(c)过程,尺寸相对小的颗粒会更早地影响气泡颈部的生长,颗粒与颈部界面相互作用,使颗粒更容易被推动,产生位移;而与气泡尺寸相比,较大的颗粒在颈部生长过程中基本不受影响.其次,在气泡的表面能快速释放的阶段,小气泡释放的能量自然也比较小.因此,存在一个合理的颗粒、气泡相对大小范围,使颗粒产生有效位移,以达到较高的初始动能转换率.鉴于上述对气泡与颗粒相对尺寸影响气泡聚并驱动颗粒运动的实验,对低黏度溶液中的双气泡聚并驱动颗粒运动的位移及其的初始动能转换率与气泡、颗粒尺寸比值的关系进行统计,并只关注两气泡尺寸相近的情况(),为两气泡的平均半径.Rb/RpRb/RpRb/Rp1.4 Rb/Rp 1.7如图 7 所示,当 在 0.71.4 之间时,颗粒被推动的位移随气泡与颗粒尺寸比出现近似线性增大.然而当 =1.7 时,颗粒并没有随之出现大幅位移,其原因可通过数值模拟来解释.据实验统计,当气泡半径在 1080m 范围,气泡与颗粒尺寸比在 0.71.7 之间,气泡聚并演化为球形,界面演化对颗粒的驱动速度在 0.010.40m/s内.同时注意到,在 值为 0.71.4 内,所对应微马达的初始动能转化率也随气泡与颗粒尺寸比呈线性增大.分析其原因,这个尺寸比范围内,颗粒足够大几乎可以完全捕获气泡聚并在颗粒方向上释放的表面能,所以其动能转化率主要由气泡与颗粒的大小来决定.根据动量守恒,小质量气泡速度大,大质量颗粒速度小,所以气泡尺寸小时,颗粒获得的初速度较小,故其初始动能转化率和位移都较小.从图 7 可推测,在 之间存在临界状态,颗粒尺寸的减小导致捕捉到的气泡聚并释放的表面能减小,初始动能转换率开始降低.(a)50 mm(b)(c)(d)(e)(f)图6Rb1=7.8m,Rb2=10.4m,Rp=18.0m 时,实验记录气泡聚并驱动颗粒运动序列图(a)0;(b)5s;(c)10s;(d)20s;(e)30s;(f)200sFig.6.ExperimentallyrecordedbubbleaggregationanddrivenparticlemotionsequencediagramatRb1=7.8m,Rb2=10.4m,Rp=18.0m:(a)0;(b)5s;(c)10s;(d)20s;(e)30s;(f)200s.0.51.01.52.0-2024680246810初始动能转换率颗粒位移/10-2/(10-5 m)b/p图7双气泡聚并驱动颗粒运动的位移、初始动能转换率随气泡与颗粒尺寸比值的变化Fig.7.Variationofdisplacement,initialkineticenergycon-version rate with bubble to particle size ratio for doublebubbleaggregationanddrivingparticlemotion.物理学报ActaPhys.Sin.Vol.72,No.16(2023)164703164703-53数值模拟根据实验,气泡界面的非对称演化推动了临近颗粒.为了探究更多颗粒限制气泡聚并的细节,并且认识到双气泡聚并的界面演化驱动颗粒的运动机理,采用伪势格子 Boltzmann 方法(伪势 LBM)对靠近颗粒的双气泡聚并过程以及驱动颗粒运动的细节进行研究,通过对界面演化过程的分析来揭示双气泡驱动颗粒运动的流动机理.篇幅所限,暂不研究溶液浓度对双气泡生成的影响.3.1 伪势 LBM 模型采用精确差分法作用力模型21和 Carnahan-Starling(C-S)状态方程以获取较大密度比的 LBM伪势模型22.其分布函数形式为23f(x+et,t+t)f(x,t)=1f(x,t)feq(x,t)+F(x,t),(6)ffeqet=(0.5)c2stcsc2s=c2/3其中,为离散粒子的运动方向,为 离散方向上的密度分布函数,为 离散方向上的局部平衡态分布函数,x 为粒子位置,为离散方向上的粒子移动速度,为时间步长,为松弛时间.运动黏度 与松弛时间的关系:,为格子声速,.粒子局部平衡态分布函数为feq=1+e uc2s+(e u)22c4su22c2s,(7)0=4/914=1/958=1/36其中,为平衡态函数对应的权重,与 u 分别为流体的密度和速度.对于 D2Q9 模型 ,.宏观密度 和流场的实际速度可以通过如下方法求得:=f,(8)u=ef.(9)F 为 方向上粒子之间的相互作用力,N 个方向上的总力表示为F(X)=g(X)c2sN=1w(|e|2)(X+e)e.(10)g(X)w(|e|2)式中,为粒子间相互作用强度,为有效密度,是用于计算各向同性相互作用力的权重.通过引入非理想流体状态方程,可以得到有效密度的表达式为=2gc2s(p c2s),(11)式中压力 p 由 C-S 状态方程来求解:p=RT1+b/4+(b/4)2(b/4)3(1 b/4)3 a2,(12)其中,a=0.4963(RTc)2/Pc,b=0.1873RTc/Pc,a=1,b=4R=1Tc和 Pc分别是临界温度和临界压力.一般为方便计算,规定,和 .颗粒的运动方程由牛顿第二定律得到,包括平动方程和转动方程:MdUdt=F(t),Jddt=T(t),(13)其中 M 是颗粒质量,U 是颗粒速度,F(t)是颗粒受力,J 为转动惯量,是角速度,T(t)是颗粒的扭矩.对于固体动边界的处理,采用半步长反弹格式及局部演化方法24:fi(x,t+t)=fi(x,t)2iuw ec2s,(14)uw=U+(x+12et xc),(15)iuwxc其中 i 和 分别表示流入和流出的方向,为边界速度,为型心位置.利用动量交换法获得颗粒边界上的表面力和力矩:F=xBce2efi(x,t)iuw ec2s,T=xBce(x xc)2efi(x,t)iuw ec2s,(16)e其中,B 代表边界周围的格点,表示格子速度,其方向从当前流体格点指向固体格点.速度 U 和角速度 通过解方程获得.对于固体壁面的亲疏水作用采用流体与固体壁之间的黏附力表征25.为了验证模拟方法的可靠性,模拟计算前进行了网格无关性验证以及 Laplace 定律的验证.3.2 模拟结果图 8 显示了计算区域,其四面均采用反弹边界条件,两气泡中心连线的中点位于模拟区域的中物理学报ActaPhys.Sin.Vol.72,No.16(2023)164703164703-6L=2 104x=1.1 106t=1.8 107T/Tc=0.0658l/g=68.8心,两气泡相界面间距为 0.Janus 颗粒的位置与两气泡相切.此模型的边界长度 m,均匀划分为 201201 个网格,单个网格的长度 m,时间步长 s,温度与临界温度比值 ,液相与气相的密度比 .模拟中的设置均为格子单位.Rb=20Rp=20模拟设置两气泡半径 m,颗粒半径 m,如图 9 所示.在图 9(a)中,0 时刻,气泡接触打通了相界面,气泡内的流体开始互通.在奇点处受到黏性力、惯性力和表面张力的共同作用,聚并开始发生,颈部开始生长.图 9(c)中气泡颈部两侧的形态出现差异,颗粒对颈部的进一步演化形成阻碍,此时颗粒没有明显的位移.图 9(d)(f)演化过程中,利用气泡相界面形变将颗粒推出.以上现象与实验中所观察到的情况比较相符,颗粒出现净位移的时间都在气泡颈部生长结束之后.通过模拟结果计算可知,颗粒的净位移为 5.78m,相当于推进了半个颗粒的长度.为了清楚气泡聚并与颗粒的作用机制,监测了气泡聚并过程中流场压强与速度场的变化.由图 9(a)可知,在 t=9.1s 左右,气泡开始出现聚并,颈部两侧出现白色的高压区,在远离颗粒侧颈部的压强略低于近颗粒侧的压强,可见,颗粒制约了这一侧高压向溶液中扩散.同样,颈部生长的界面形变诱导附近流场,对颗粒有向 r 方向的作用力.气泡内部气体流速较高,有向外扩张的趋势.但是在界面表面张力的作用下,流体速度指向界面内.然而可以看到在聚并气泡外侧,有一圈高压区域,这是方法本身带来的虚拟压力.图 9(b)中气泡聚并颈部一侧界面紧贴颗粒,高压区消失.另一侧由于颈部的快速生长,带动附近流体向 r 负方向流动.t=36.6到 s(图 9(c),远离颗粒侧的颈部生长已经完成,周围流场中的流体速度明显降低,由此反映出,气泡聚并颈部的生长速度逐渐降低.Janusbbparticle反弹边界反弹边界反弹边界反弹边界图8反弹边界条件液相流场中双气泡数值模拟的计算域及边界条件Fig.8.Computationaldomainandboundaryconditionsforthenumericalsimulationofdoublebubblesinliquid-phaseflowfieldwithreboundboundaryconditions.(a)(b)(c)(d)(e)(f)/(mSs-1)0.100.090.080.070.060.050.040.030.020.01/Pa0.020-0.02-0.04-0.06-0.08-0.10-0.1240 mm图9不同时刻近颗粒气泡聚并演化过程压强及速度场变化(a)9.1s;(b)18.3s;(c)36.6s;(d)73.1s;(e)310.8s;(f)914.1sFig.9.Variationofpressureandvelocityfieldsduringnear-particlebubbleaggregationandevolutionatdifferentmoments:(a)9.1s;(b)18.3s;(c)36.6s;(d)73.1s;(e)310.8s;(f)914.1s.物理学报ActaPhys.Sin.Vol.72,No.16(2023)164703164703-7此时的颗粒依然未被明显推动.图 9(d)中,Z 方向上的气泡两端开始向内回缩,周围流体向 O 点流动.而由于受到颈部生长惯性的作用以及气泡内部Z 方向流体受到挤压也向 r 的反方向运动,在远离颗粒侧的界面已经与原有的两个气泡的界面融合为一个具有更大曲率半径的界面.在贴近气液界面处流体速度梯度较大,随 Z 方向上的流体,向中心O 流动.t=310.8图 9(d)(e)说明聚并过程中周围流体会持续作用,气泡在 r 方向上持续形变,颗粒被气泡推出.这个时段中颗粒出现明显净位移,也是整个演化过程中对驱动贡献最大的时间段.到 s 左右,Z 方向上流体向中心 O 的流动速度减小,说明Z 方向上的界面形变趋于稳定.根据速度场显示,颗粒周围更多流体向 r 方向运动,说明此时的颗粒正以较大的速度向 r 方向运动.在气泡界面将颗粒推动到最远处后,颗粒也会在流场的影响下,继续向 r 方向移动.气泡界面与颗粒相接触,呈现出水滴状,但在表面张力的作用下,气泡在 Z 方向回缩,与颗粒分离,即图 9(f)所示.界面回弹势必会引起颗粒与气泡之间的流体,向 r 反向流动,这时颗粒也会在这部分流体的影响下,向 r 反向有小距离的回拉.回弹后的气泡已基本演变为一个球体.Rb/Rp=4Rb/Rp 1Rb/RpRb/Rp图 10 为实验与模拟中不同时刻气泡聚并推动颗粒的形态演化变化对比分析,发现数值模拟与实验的演化形态基本一致.颗粒拥有最大位移的时间均在 t=100200s 之间,即气泡聚并后在表面张力作用下 Z 方向开始向中间收缩,在惯性作用下 r 方向拉伸阶段,是驱动颗粒使之产生较大净位移的主要方式.通过数值模拟计算了不同气泡与颗粒尺寸比的初始动能转换效率以及驱动颗粒的位移,如图 11所示.当 时,大气泡驱动小颗粒会提高驱动位移,从而提高驱动速率,但是这个增量不足以抵消颗粒质量本身对动能的影响,造成了气泡尺寸与颗粒尺寸的比值越大,转换率越低.但是当 时(图中绿色虚线左侧部分),上述规律并不成立,小气泡驱动大颗粒,首先释放的能量低,其次,从动量守恒考虑,颗粒质量远大于气泡,小气泡被弹开.模拟结果表明,颗粒的位移随气泡与颗粒尺寸比的增加呈现近似线性增长,在 =1 处的初始动能转换率达到最大,随着 逐渐接近 0,几乎无法驱动颗粒.0 ms0 ms7 ms5 ms50 mm27 ms10 ms62 ms20 ms154 ms30 ms210 ms200 ms(a)(c)(e)(b)(d)(f)图10实验与数值模拟气泡聚并驱动颗粒形态演化对照橙色框内为数值模拟结果(Rb=20m,Rp=10m),蓝色框内为实验结果(Rb1=41m,Rb2=39m,Rp=34m)Fig.10.Comparisonofexperimentalandnumericalsimulationsofbubbleaggregation-drivenparticlemorphologyevolution.Theorangeboxshowsthenumericalsimulationresults(Rb=20m,Rp=10m),andtheblueboxshowstheexperimentalresults(Rb1=41m,Rb2=39m,Rp=34m).物理学报ActaPhys.Sin.Vol.72,No.16(2023)164703164703-84讨论与分析以观察到的实验现象为基础,通过数值模拟方法研究了微米量级的双气泡聚并驱动颗粒的多相流体动力学机制,阐明了双气泡聚并驱动微马达的可行性.实验中观察到,气泡发生聚并之后,没有立即发生破裂,而是在溶液中稳定存在一段时间,这个现象出现的概率随着溶液黏度的增大而增大.没有立即出现气泡溃灭或与界面融合的情况,也说明气泡聚并微马达适用于体相溶液中,且自由界面对气泡聚并过程的影响暂时可以忽略.对比实验与数值模拟结果,获得了定性一致的结果,说明对气泡聚并的力学机制分析的正确性,但其中部分情况存在定量的微小差异.这原因可能在于,由于完整的气泡过程包括化学反应-气泡生长-气泡聚并-气泡脱离的连续过程,处于降低模拟难度与复杂性的考虑,在模拟中仅关注了气泡聚并这一关键环节进行研究,因此在数值模型的建立上进行简化.数值模拟中,没有考虑实际成核点位置与气泡生长的影响,根据实验和实验,从气泡表面能的释放量来看,对比无颗粒时的气泡聚并过程和气泡附着于颗粒生长过程中发生聚并,气泡生长对推动颗粒的运动影响不明显.通过对颗粒周围速度场的观察,本文提出的双气泡聚并的驱动方式,这一现象可类比 Jens 等26提出的自然界中“偶极泳动”中“Pusher”的一类,图 12 中“马达”在后面形成的流场推动颗粒运动,出现“Pusher”的特征流场.这类似自然界中微生物通过鞭毛摆动,使得左右两侧流体向中间汇集,再向前后推动,进而实现泳动.本文提出的体系中气泡界面的聚并演化过程实际上实现了类似鞭毛的非对称运动效果,构造出了类似的流场.40 mm(a)(b)图12(a)“力偶极子”微马达流场26;(b)双气泡聚并微马达速度场Fig.12.(a)Flowfieldofthe“forcedipole”micromotor26;(b)velocity field of the double bubble aggregation andmicromotor.在实验过程中发现,当两气泡尺寸差距较大,能够改变颗粒驱动的方向,在以后的研究中,若能实现对气泡大小的实时量化控制,那么对气泡聚并微马达的方向和速度调控会更加精确.然而,深刻的理解气泡聚并与颗粒之间相互作用的全部细节,需要进一步数值模型研究.目前模拟中采用二维多相流模型,可以在现象上给出流动机理的解释,精确的定量对比需要在三维模拟的基础上进行.5结论本文提出用双气泡聚并方式驱动 Janus 微球马达的新颖驱动方式,通过实验研究与数值模拟对气泡聚并驱动颗粒的力学机理与运动机制进行了研究,得到下列结论.1)气泡聚并导致的界面振荡是驱动微马达运动的主要原因,其能量转化率介于管式微马达反弹驱动与气泡和界面融合驱动之间.2)当气泡颗粒尺寸比在 0.71.4 内,气泡聚并驱动颗粒的位移随气泡与颗粒尺寸比的增大而增大,初始动能转换率基本一致.3)在黏度较小的溶液中,气泡聚并驱动颗粒产生净位移的主要阶段是气泡聚并颈部生长完成后,界面推动以及界面推出后,惯性作用的结果;在黏度相对较大的溶液中,气泡聚并驱动颗粒产生净位移的主要阶段只有气泡聚并颈部生长完成后,界面推动的结果.012345012345/10-5012345/(10-5 m)b/p初始动能转换率颗粒位移图11气泡与颗粒尺寸比与初始动能转换率、颗粒位移之间的关系Fig.11.Relationshipbetweenbubbletoparticlesizeratioandinitialkineticenergyconversionrate,particledisplace-ment.物理学报ActaPhys.Sin.Vol.72,No.16(2023)164703164703-9参考文献 WangLL,ChenL,ZhengX,YuZX,LWC,ShengMJ,WangLN,NiePC,LiHY,GuanDS,CuiHH2022Small1822038721LiMY,WuJR,ZhengYY2020Chin.Sci.Bull.654123(inChinese)李茂垚,吴建荣,郑元义2020科学通报6541232LiM,XiN,WangYC,LiuLQ2021IEEE Trans.Biomed.Eng.681303LiuL,WangDW,RaoW2021Micromachines122804WangW,DuanWT,AhmedS,MalloukTE,SenA2013Nano Today85315FengYW,JiaDL,YueHE,WangJ,SongWP,LiLQ,Zhang A M,Li S,Chang X C,Zhou D K 2023 Small 1922075656BarabanL,StreubelR,MakarovD,HanLY,KarnaushenkoD,SchmidtOG,CunibertiG2013ACS Nano713607PavlickRA,SenguptaS,McFaddenT,ZhangH,SenA2011Angew.Chem.Int.Ed.5093748ZhengX,CuiHH,LiZH2017Chin.Sci.Bull.62167(inChinese)郑旭,崔海航,李战华2017科学通报621679WangSJ,WuN2014Langmuir30347710ManjareM,YangB,ZhaoYP2012Phys.Rev.Lett.10912830511WeiMJ,ChenL,WuT,ZhangHY,CuiHH2017Acta.Phys.Sin.66164702(inChinese)魏梦举,陈力,伍涛,张鸿雁,崔海航2017物理学报6616470212WangLL,ChenL,ZhangJ,DuanJM,WangL,LiZH,ZhengX,CuiHH2018Langmuir341042613VanLeeuwenJL2010Science32939514Noblin X,Rojas N O,Westbrook J,Llorens C,Dumais J2012Science.335132215NaeemS,NaeemF,ManjareM,LiaoF,BolaosQuionesVA,HuangGS,LiY,ZhangJ,SolovevAA,MeiYF2019Appl.Phys.Lett.11403370116ZhangPY,Hu