气流作用下同轴带电射流的不稳定性研究1)

李帅兵 杨 睿 罗喜胜 司 廷

(中国科学技术大学近代力学系，合肥230027)

气流作用下同轴带电射流的不稳定性研究1)

李帅兵 杨 睿 罗喜胜 司 廷2)

(中国科学技术大学近代力学系，合肥230027)

通过对气体驱动同轴电流动聚焦的实验模型进行简化，开展了电场力和惯性力共同作用下同轴带电射流的不稳定性理论研究.在流动为无黏、不可压缩、无旋的假设下，建立了三层流体带电射流物理模型并得到了扰动在时间域内发展演化的解析形式色散关系，利用正则模方法求解色散方程发现了流动的不稳定模态，进而分析了主要控制参数对不稳定模态的影响.结果表明，在参考状态下轴对称模态的最不稳定增长率最大，因此轴对称扰动控制整个流场.外层气流速度越高，气体惯性力越大，射流的界面越容易失稳.内外层液--液同轴射流之间的速度差越大，射流越不稳定.表面张力对射流不稳定性起到促进作用.轴向电场对射流不稳定性具有双重影响：当加载电场强度较小时，射流不稳定性被抑制;当施加电压大于某一临界值时，轴向电场会促进射流失稳.临界电压的大小与界面上自由电荷密度和射流表面扰动发展关系密切.这些结果与已有的实验现象吻合，能够对实验的过程控制提供理论指导.

同轴射流,带电射流,电流动聚焦,射流不稳定性,过程控制

引言

电流体力学(electrohydrodynamics,EHD)是流体力学和电动力学之间的交叉学科.随着纳米技术、微机电系统和生物技术的迅猛发展，EHD研究在纳米材料、微流动、生物医学等领域显示了巨大的应用价值.电雾化(electrospray，ES)是EHD的一个重要分支，可描述为利用单轴或复合结构的毛细管，在电场力作用下促使流体界面拉伸成微细射流[13].在ES中，带电液体射流周围的气体是静止的.不同于ES，近年来另一种有效的微细液体射流产生方法被提出，称之为流动聚焦( fl ow focusing，FF)，也就是从毛细管流出的液体由高速运动的气体驱动，经小孔聚焦后形成稳定的锥射流结构[35].可以看出，在FF中液体射流周围的气体是高速运动的.FF和ES形成的液体射流将破碎成微纳米量级的液滴、颗粒或胶囊，具有相似的射流不稳定性机理.Ga˜n´an-Calvo等[6]和司廷等[7]将FF与ES相结合，提出了电流动聚焦(electro- fl ow focusing，EFF)技术，能够使雾化更稳定，制备的产物尺寸更小，因而具有重要的应用前景.

单轴结构的ES,FF和EFF相对简单，目前实验和不稳定性理论研究都已取得一定进展[5,89].Si等[1011]采用时间线性不稳定性方法分析了各个参数对不稳定模的影响，采用时空线性不稳定性方法分析了参数对鞍点的位置以及雷诺数--韦伯数(Re-We)平面上绝对和对流不稳定性的转换曲线的影响，均取得了一定成果.Li等[12]采用时间不稳定性方法进行了轴向电场作用下有黏及无黏单轴射流不稳定性研究.理论和实验结果吻合，一定程度上证明了不稳定性理论分析方法的可靠性.

同轴结构的ES、FF和EFF变得复杂，不仅多了一层流体，还增加了一层界面[1315].目前相应的实验研究已经开展，但是同轴带电射流不稳定性的理论和实验研究仍十分有限[5,1623].在ES方面，根据加载电场方向的不同，Li等分别就径向[18]、轴向[19]电场作用下同轴射流无黏稳定性进行了理论分析，取得了一定成功.根据射流流体性质的不同，Li等[19]将轴向电场同轴射流分成4种情况，并针对每种情况进行了射流的不稳定性分析.本文是这一工作的延续，选取内层为导体外层为电介质的同轴射流在轴向电场作用下的情况，进一步研究惯性力和电场力共同作用下的同轴带电射流不稳定性，为气体驱动同轴EFF实验提供理论指导.

本文采用的不稳定性研究方法为模态展开法(即正则模法).在稳定性分析过程中，当考虑黏性对流体影响时，控制方程存在非线性项，使得求解过程变得较复杂[24].李广滨等[25]在单轴EFF研究中，采用无黏假设研究了表面张力、外层气体流速等参数对射流的影响，取得了一定研究成果，对实验结果进行了定性预测.本文采用类似的基本假设和无黏情况下的时间不稳定性分析方法，将小扰动展开成Fourier级数的形式，推导并求解出解析形式的色散方程，从而分析扰动在时间域中的发展，开展同轴结构的EFF理论研究.

1 理论模型

1.1 模型简化及基本假设

完全按照实验情况建立物理模型需要考虑的因素较多且十分复杂，因此在开展理论研究之前需要对物理模型进行简化.在前期工作中，Li等[1819]建立了同轴ES中的同轴带电射流模型，Li等[12]建立了单轴EFF中的气体作用下的带电射流模型，取得一定进展.本文在以上两个模型基础上，开展同轴射流在外部气流裹挟下向下游运动的理论建模和不稳定性分析.由于复合液体射流流出小孔后在一段距离内稳定运动，直径几乎保持不变，可以将其近似为圆柱形内部液体，环形中间层液体以及外部高速气流，如图1所示.

与Li等[12]的假设类似，本模型给出下列假设或近似并对符号含义进行说明.同轴射流由半径为R1的圆柱射流(内层流体)和半径为R2的环形射流(中间层流体)组成；同轴射流外部为压强差推动的气体(外层流体).本文假设流动无旋、无黏、不可压，忽略流体间温度差异，假设三种流体等温且不互溶，界面无质量传递.忽略磁场对射流的作用，假设内层流体为导体，电导率为σ1，电容率为ε1；中间层流体为电介质，电容率为ε2.内层流体以及中间层流体之间界面为“界面1”，中间层流体以及外层气体射流之间界面为“界面2”.内层、中间层以及外层流体密度分别为ρ1,ρ2,ρ3，速度为u1,u2,u3.流体无扰动情况下的速度型为均匀速度型；射流处在场强为E0的理想轴向电场下.

图1 气流作用下液--液同轴带电射流的简化模型，其中1,2,3分别表示内层、中间层和外层流体Fig.1 Simpli fi ed model of electri fi ed coaxial liquid-liquid jet in a co fl owing gas stream.1,2 and 3 stand for the inner,middle and outer fl uids

下面从量纲分析的角度解释忽略重力、体电荷以及流体黏性的合理性.在流体力学中无量纲参数Fr数 (惯性力/重力)以及 Bo数 (重力／表面张力)与重力相关.取一组同轴射流实验[26]的参数作为参考状态：速度 U0=2m/s；长度 L0=30µm；密度ρ0=1000kg/m3；黏度 µ0=10−3kg/(s·m)；气体压力差∆pg=5kPa，液体流量速度Q0=50mL/h，表面张力系数γ0=7.17×10−2N/m,重力加速度g=9.8m/s2.可得

可见在同轴射流中，重力相对于惯性力以及表面张力均为可忽略的小量.雷诺数量级为102，当液体流量速度或气体压力差继续增大时，雷诺数量级可以增大到103.因此，本文不考虑重力以及黏性影响是合理的.

水的电导率为 σw=1.4×10−4S/m，介电常数εwε0≈ 7.1×10−10C2/(N·m2).可得电松弛时间以及对流运动特征时间为

可见电松弛时间相对于对流运动特征时间很小.又因内层流体为导体，中间层流体为电介质，由电磁场知识可知，在轴向电场作用下，电荷有充足时间松弛到内层和中间层流体间的界面上并保持平衡状态.因此，可以认为内层、中间层及外层流体中不存在体电荷，电荷存在于导体和电介质之间的界面(即内界面)上.而电荷的这种分布一方面会在最内层流体中提供径向电场，另一方面会和轴向电场相互作用产生切向电应力.在不考虑黏性剪切力情况下，该切向电应力将孤立存在[16].为此本文假设模型中自由电荷密度足够小，以致可以忽略其对射流稳定性的影响.

1.2 控制方程及边界条件

由于假设流体为不可压流体，射流满足连续性方程为：▽·ui=0.其中i=1,2,3分别表示内层、中间层以及最外层流体.因为流动无旋，故可引入速度势函数ϕi

忽略重力，电场力，磁场力以及黏性作用的动量方程为

忽略磁场和体电荷的Maxwell方程为

本文假设射流处于定常状态，故只需对模型提出边值条件.

(1)射流对称轴上的扰动量有限以及无穷远处的扰动为零其中为速度势扰动量，为电势扰动量.

(2)交界面上的运动学边界条件

假设扰动后界面i(i=1,2表示界面1,2)离开平衡位置距离为ηi，则因分界面为一个物质面，可得

展开上式，去掉二阶及二阶以上的小量得到线性化后的运动学边界条件

(3)交界面动力学边界条件

界面1上法向应力有压力、电场力以及表面张力，三力满足平衡关系

同样在界面2上有

其中n为界面上法向单位矢量，T为水动力学张量Th与电应力张量Te之和T=Th+Te

其中ε为介电常数，δ为Kronecker函数.

将Th以及Te代入平衡方程(8)和式(9)得

(4)电场边界条件

Ei(i=1,2,3)分别表示内层、中间层以及外层流体中的电场强度，εi表示介质的绝对介电常数，Di=εiEi表示电位移矢量，qs为界面1上的自由电荷密度，σi为物质的电导率，Ji=σiEi表示电流密度.忽略磁场作用，由法拉第感应定律得

上式也是电场分量在切向方向连续的条件.

由高斯定律得

由界面上电流密度连续条件可得

1.3 基本流场和电场

假设流动无黏，采用 Kelvin-Helmholtz速度型,Ui(r,θ,z)=(0,0,Uz).

基本电场为轴向均匀电场：Ei=E0ez，i=1,2,3.

2 色散关系

2.1 线性稳定性分析

当界面受到任意小扰动时，它们将离开平衡位置发生微小的位移，新位置可以表示为：ri=Ri+ηi(i=1,2表示界面1和界面2).

内层、中间层以及外层流体中，流体速度、压力、电势等物理量可以分解成基本量加扰动量的形式式中 uir,uiθ,uiz,pi,Vi,ϕi分别表示流体 i(i=1,2,3分别表示内层、中间层和外层流体)中的径向速度分量、周向速度分量、轴向速度分量、压力、电势以及速度势.为相应物理量的扰动量；Uiz,Pi0,Vi0,Φi,Qis分别表示无扰动情况下的轴向速度分量、压力、电势、速度势、自由电荷密度.

本文采用正则模法处理稳定性问题.假设内外界面受到的是任意无限小的扰动.这些扰动相互独立，每一个扰动都可以分解成Fourier级数，写成扰动量乘以 exp[ωt+i(kz+nθ)]的形式

其中ℜ=ωt+i(kz+nθ).扰动量是r的函数.

其中

将ui代入动量方程(2)得

由式(4)可引入势函数Vi：Ei=−▽Vi(i=1,2,3)又由式(3)可得

由式(12)～式(14)得

2.2 色散方程及无量纲化

选择中间层射流半径R2，密度ρ2，无扰动速度U2，表面张力系数γ2以及空气介电常数ε3作为特征尺度.设定无量纲参数：α=R2k为无量纲波数；β= ωR2/U2为无量纲复频率；ρ12= ρ1/ρ2为内层流体与中间层流体密度比；ρ32=ρ3/ρ2为外层流体与中间层流体密度比；U12=U1z/U2z为内层流体与中间层流体的速度比；U32=U3z/U2z为外层流体与中间层流体的速度比；a=R1/R2为无扰动时界面1与界面2的位置比；Γ=γ1/γ2为界面1与界面2的张力系数比；ε13=ε1/ε3为内层流体与外层流体的介电常数比；ε23=ε2/ε3为中间层流体与外层流体的介电常数比；We= ρ2U22R2/γ2为 Weber数；Eu= ε3E20/ρ2U22为电Euler数.

将 p1,p2,p3,E1,E2,E3代入动力学边界条件(10)和(11)，经过运算，最终得到无量纲化的色散方程

同时可得扰动后界面偏离平衡位置的距离比

其中φ1−φ2表示界面运动的相位差，φ1−φ2=0表示界面同向运动，φ1−φ2=π表示界面反向运动.

当E0=U3z=0时色散方程转化为

此为无黏不带电同轴射流色散关系，与李广滨[26]所得结果相同.

3 理论结果

3.1 无量纲参数的参考值

无量纲色散关系式(22)是关于无量纲复频率β的四次代数方程，给定一组(n,α)可得两对四个复特征值 β.由 ℜ =ωt+i(kz+nθ)及 β= ωR2/U2可知，β的实部βr为正值对应于射流的不稳定状态.

L´opez-Herrera等[27]以水作为导体以及葵花籽油作为电介质开展了同轴ES实验.相应试剂的物理参数见表1.葵花籽油与空气之间界面表面张力系数γ2为 32.8mN/m，葵花籽油与水之间界面张力系数γ1为16.6mN/m.实验装置包括内外直径为200µm和 400µm的内部针管以及内外直径为 420µm和800µm的外部针管.这样内层流体为半径约100µm的圆柱射流，中间层流体为半径约10µm的环形射流.作为中间层流体的葵花籽油流速U2≈1m/s，又因R2≈110µm.加在同轴射流上的电压约为1kV，电极距离约为1cm[19]，可得E0≈105V/m.空气的绝对介电常数为 εair=8.85× 10−12C2/(N·m2). 根据实验参数，本文拟设定无量纲参数的参考值近似为ρ12=1.19,ρ32=0.001,ε13=80,ε23=3.4,Γ =0.51,a=0.9,U12=1.25,U32=2,We=3,Eu=10−4.

表1 实验试剂的物理特性Table 1 Physical property of the tested material

3.2 周向波数的影响

n代表周向波数，n=0表示轴对称不稳定性，n≥1代表第n类非轴对称不稳定性.由图2可知在参考状态下只存在轴对称不稳定性.为方便，本文主要围绕轴对称不稳定模态开展相关的研究.

图2 周向波数n对扰动增长率βr的影响Fig.2 E ff ect of the azimuthal wave number n on the temporal growth rate βrof perturbations

3.3 电欧拉数的影响

电场对射流稳定性影响是通过电欧拉数Eu来体现的，对色散方程赋不同的Eu，在给定(n,α)情况下得到不同的β.通常会有一至两个β实部为正，对应同轴带电射流的不稳定状态.色散关系的曲线存在一个最大扰动轴向波数及其对应的最大扰动增长率：(αmax,βrmax).同时存在若干区间 (α1,α2)，使得对于任意α∈(α1,α2)，对应的扰动增长率实部βr＞0.本文采用李芳[16]的模态分析法将不稳定波数区间划分为不同模态进行稳定性分析.将实部为正的无量纲复频率β代入关系式(23)中发现κ2并不严格地等于0和π.如图3所示：相位差约为0的模态称为类反对称模，相位差约为π的部分称为类对称模，而相位差随着轴向波数不同在0和π之间连续变化的对应过渡模.

图3 (a)类反对称模,(b)类对称模,(c)过渡模Fig.3(a)parasinuous mode,(b)paravaricose mode,(c)transitional mode

当Eu≈10−4甚至更小时，存在类反对称模、类对称模以及过渡模 3种模态[16]，所得结果与文献[19]所得结果相近.此时类反对称模对应的扰动增长率显著大于另两种模态的扰动增长率起主导作用.图4(a)显示的是Eu＜10−3时的类反对称模.由图4(a)可知电欧拉数Eu从0增长到10−3时，扰动增长率βr变化很小，且扰动增长率βr随着Eu增大而减小.可见界面电荷密度以及加载轴向电场很小时，轴向电场虽然能够抑制射流不稳定性但其作用较小.司廷等[7]在其单轴射流实验中研究了电场对射流稳定性的影响，发现当加载电压V0从0kV提高到0.5kV时射流破碎长度增大，表明加载电场较小时电场强度的增大有助于抑制射流不稳定性.本文结果与实验结论吻合.

当Eu在10−3和10−2之间时只出现类反对称模.此时不同 Eu下 βr与 α关系如图 4(b)所示.随着Eu的增大，最大扰动增长率βrmax减小，即轴向电场的增大对扰动不稳定性仍有抑制作用.值得关注的是，图4(b)中当Eu→0.1时，Eu的变化对最大扰动增长率的影响迅速减小，即∆βrmax/∆Eu逐渐减小，所以Eu→ 0.1时电场对不稳定性的抑制作用减小.而当Eu＞0.1时(图4(c))，最大扰动增长率βrmax显著变大，且对应的最不稳定波数αmax变大，不稳定模由长波不稳定性向短波不稳定性转换.

图4 电欧拉数Eu对扰动增长率βr的影响Fig.4 E ff ect of the electric Euler number Eu on the temporal growth rate βr

由上述分析可知在Eu=0.1附近存在临界点Eu=Ecri，当Eu≪ Ecri时3种不稳定模态共存，类反对称模起到主导作用；当0≪Eu＜Ecri时只存在类反对称模；在这两种情况中轴向电场对射流的不稳定性起到抑制作用.当Eu＞Ecri时轴向电场能够促进射流的不稳定性.司廷等[7]在其实验中发现，加载电场小于0.5kV时能够抑制射流不稳定性，但加载电场大于0.5kV时有助于促进射流不稳定性.在本文所列参数下，临界点Ecri在0.1附近远大于实验中Eu的临界值.原因在于：一方面理论中忽略了界面上的自由电荷；另一方面当电压较大时非轴对称运动剧烈，破碎的液体会变形、扭曲并进一步破碎，这些情况均会导致临界点Ecri减小.

3.4 外层气流的影响

外层流体的流速对射流不稳定性的影响是通过改变U32来实现，如图5所示.当U32＜10时最外层流体对射流不稳定性的影响较小；U32→0时，扰动增长曲线(α,βr)以U32=0时的扰动增长曲线为界，并不断趋近于此边界.当U32逐渐变大后，最大扰动增长率βrmax增大且射流不稳定模态由长波不稳定性向短波不稳定性演化.可见外层流体低速流出时，对射流不稳定性影响较小；当其以较高流速裹挟同轴射流流出小孔时，更易使射流破碎.这类似于带电单轴射流情况中最外层气流流速对射流不稳定性的影响[25].

图5 外层流体与中间层流体流速比U32对增长率βr的影响Fig.5 E ff ect of the velocity ratio U32between the outer layer fl uid and middle layer fl uid on the growth rate βr

比较图4(c)及图5发现，不同Eu下射流扰动规律与不同U32下射流扰动规律相似：长波区域(α,βr)呈线性关系；在短波区域Eu及U32对射流的不稳定性影响显著.Eu及U32对带电同轴射流的影响主要体现在短波区域.相较而言射流不稳定性更依赖于Eu的变化，但由于实验条件限制加载较高Eu不切实际，故通过提高外层气流速度更方便控制射流不稳定性.实验中，往往通过改变气体压力差∆pg来控制流速u3.定常状态下忽略重力对射流的影响，由伯努利方程得：.司廷等[11]在研究单轴流动聚焦时发现，气体压力差∆pg对液体射流的形态影响较大，随着∆pg的增大，在频闪灯下，实验图像表明扰动波长λ不断减小，振幅不断变大.李广滨[26]在研究带电单轴射流以及不带电同轴射流的过程中发现，随着∆pg增大，射流内外径减小，液体射流稳定性越差，扰动由长波不稳定性向短波不稳定性过渡.本文结果与上述结论相符.

3.5 同轴射流速度比的影响

同轴射流内层流体与中间层流体之间的速度比由U12表示随着U12变化同轴射流速度比对稳定性的影响可分为3个阶段.

内层流体流速相对于中间层流体流速很小时，即当U12较小时，如图6(a)所示，扰动增长曲线有两种模态：在长波区存在类对称模(左峰)，短波区存在过渡模(右峰)，此时最大扰动增长率βrmax往往出现在短波区域，过渡模起到主导作用.βrmax随着内层射流流速增大而减小，类对称模逐渐起到主导作用.

内层流体流速与中间层流体流速比进一步增大，当U12～1时只出现类反对称模.如图6(b)，U12变化并未对最大扰动增长率和不稳定波长区间产生较大影响，其影响主要在长波区.可见同轴射流内外层流体流速相差较小时，也即同轴射流内外层流体压力差较小时，压力差的变化对射流稳定性影响较小.

图6 同轴射流间速度比U12对增长率βr的影响Fig.6 E ff ect of the velocity ratio U12between the coaxial jet on the growth rate βr

图6 同轴射流间速度比U12对增长率βr的影响(续)Fig.6 E ff ect of the velocity ratio U12between the coaxial jet on the growth rate βr(continued)

当内层流体流速相对于中间层流体流速较大时(U12＞1)，在长波区出现类对称模，短波区出现过渡模.如图6(c)所示，随着流速U12的增大，过渡模的影响愈来愈显著，最大扰动增长率βrmax变大，射流不稳定性由长波不稳定性向短波不稳定性转化.

可见内外层同轴射流流速相差较大时，存在两个不稳定模态，随着流速差减小，双模中过渡模的作用衰减，类对称模不稳定作用凸显.当内层以及中间层同轴射流流速差较小时仅存在类对称模.随着流速差增大，射流不稳定由长波不稳定性向短波不稳定性过度，最大扰动增长率βrmax增大，射流不稳定越明显.

3.6 表面张力的影响

当同轴射流实验针头内外径以及实验材料确定时，两层界面之间We1与We的关系只依赖于射流流速U12的大小.上文已研究了同轴射流速度比对射流不稳定性的影响，在此我们仅讨论界面2上表面张力的影响，即考虑不同We下增长率βr的变化.

图7 韦伯数We对增长率βr的影响Fig.7 E ff ect of the Weber number We on the growth rate βr

如图7(a)所示，当We≤1时存在类反对称模以及类对称模，两种模的不稳定波数区间相近，但前者起到主导作用.当We＞1时，只存在类反对称模，此时不同We数下(α,βr)关系如图7(b)所示.随着We不断增大，表面张力作用逐渐减弱，扰动增长率βr逐渐减小，射流轴对称不稳定性被抑制.说明表面张力对时间不稳定性起到促进作用，这与单轴情况[25]波数α＜1时的结果相同.

4 结论

通过对气体驱动同轴电流动聚焦的实验模型进行简化，构建了三层流体带电射流物理模型.利用正则模方法研究射流不稳定性，得到解析形式的色散方程.数值求解色散方程，分析了电欧拉数Eu，外层气流流速U32，同轴射流速度比U12，韦伯数We等控制参数对射流稳定性的影响.得到如下结果：

(1)轴向电场对射流不稳定性具有双重影响.电场强度存在临界值Ecri，当电场小于临界值Ecri时，轴向电场能抑制射流的不稳定性；当电场大于临界值Ecri时，轴向电场会大幅增强射流不稳定性.临界值Ecri的大小与界面自由电荷密度以及射流表面扰动发展关系有关；

(2)气体压力差∆pg增大，外层气体射流流速提高，惯性力变大，能够促进短波区域射流不稳定性，但对长波区域稳定性影响有限.此结论与实验现象及单轴情况下相应结果吻合；

(3)同轴射流流速差越大，不稳定性越显著，过渡模主导作用越强.流速差由大向小过渡时类对称模的不稳定性逐渐占据主导地位.当最内层以及中间层同轴射流流速相差很小时，只存在类反对称模，此时流速差变化对射流不稳定性几乎无影响；

(4)表面张力对射流不稳定性起到促进作用，与单轴情况所得结论相同.

下一步我们将尝试利用稳定性分析方法，研究黏性情况下各物理参数特别是Re对射流不稳定性的影响，并探讨实验中观察到的各种现象的物理机制.

1 Ferna´ndez de la Mora J.The fl uid dynamics of Taylor cones.Annu Rev Fluid Mech,2007,39:217-243

2 Melcher JR,Taylor GI.Electrohydrodynamics:a review of the role of interfacial shear stresses.J Fluid Mech,1969,1:111-146

3 Barrero A,Loscertales IG.Micro-and nanoparticles via capillary fl ows.Annu Rev Fluid Mech,2007,39:89-106

4 Gan˜a´n-Calvo AM.Generation of steady liquid microthreads and micron-sized monodisperse sprays in gas streams.Phys Rev Lett,1998,80:285-288

5 司廷,李广滨,尹协振.流动聚焦及射流不稳定性.力学进展,2017,47:201706(Si Ting,Li Guangbin,Yin Xiezhen.Flow focusing and jet instability.Advances in Mechanics,2017,47:201706(in Chinese))

6 Gan˜a´n-Calvo AM,Lo´pez-Herrera JM,Riesco-Chueca P.The combination of electrospray and fl owfocusing.JFluid Mech,2006,566:421-445

7 司廷,李广滨,田瑞军等.电场作用下流动聚焦的实验研究.力学学报,2011,43(6):1030-1036(Si Ting,Li Guangbin,Tian Ruijun,et al.Experimental study of the fl ow focusing under an electric fi eld.Chinese J Theo Appl Mech,2011,43(6):1030-1036(in Chinese))

8 司廷.流动聚焦的实验和理论研究.[博士论文].合肥：中国科学技术大学,2009(Si Ting.Experimental and theoretical investigation on fl ow focusing.[PhD Thesis].Hefei:University of Science and Technology of China,2009(in Chinese))

9 Gan˜a´n-Calvo AM,Montanero JM,Mart´ın-Banderas L,et al.Buildingfunctionalmaterialsforhealthcareandpharmacyfrommicro fl uidic principles and fl ow focusing.Adv Drug Deliv Rev,2013,65:1447-1469

10 Si T,Li F,Yin XY,et al.Modes in fl ow focusing and instability of coaxial liquid-gas jets.J Fluid Mech,2009,629:1-23

11 Si T,Li F,Yin XY,et al.Spatial instability of co fl owing liquid-gas jets in capillary fl ow focusing.Phys Fluids,2010,22:112105

12 Li GB,Luo XS,Si T,et al.Temporal instability of co fl owing liquidgas jets under an electric fi eld.Phys Fluids,2014,26(5):054101

13 Si T,Feng HX,Luo XS,et al.Formation of steady compound conejet modes and multilayered droplets in a tri-axial capillary fl ow focusing device.Micro fl uid.Nano fl uid,2015,18:967-977

14 Si T,Li GB,Wu Q,et al. Optical droplet vaporization of nanoparticle-loaded stimuli-responsive microbubbles.Appl Phys Lett,2016,108:111109

15 Si T,Yin CS,Gao P,et al.Steady cone-jet mode in compound fl uidic electro- fl ow focusing for fabricating multicompartment microcapsules.Appl Phys Lett,2016,108:021601

16 李芳.同轴带电射流的稳定性研究.[博士论文].合肥：中国科学技术大学，2007(Li Fang.Stability investigation on coaxial electrifi ed jets.[PhD Thesis].Hefei:University of Science and Technology of China,2007(in Chinese))

17 Gan˜a´n-Calvo AM,Gonza´lez-Prieto R,Riesco-Chueca P,et al.Focusing capillary jets close to the continuum limit.Nat Phys,2007,3(10):737-742

18 Li F,Yin XY,Yin XZ.Instability analysis of a coaxial jet under a radial electric fi eld in the nonequipotential case.Phys Fluids,2006,18:037101

19 Li F,Yin XY,Yin XZ.Linear instability of a co fl owing jet under an axial electric fi eld.Phys Rev E,2006,74:036304

20 Loscertales IG,Barrero A,Guerrero I,et al.Micro/nano encapsulation via electri fi ed coaxial liquid jets.Science,2002,295:1695-1698

21 Utada AS,Lorenceau E,Link DR,et al.Monodisperse double emulsions generated from a microcapillary device.Science,2005,308:537-541

22 Mart´ın-Banderas L,Flores-Mosquera M,Riesco-Chueca P,et al.Flow focusing:a versatile technology to produce size-controlled and speci fi c morphology microparticles.Small,2005,7:688-692

23 Herrada MA,Montanero JM,Ferrera C,et al.Analysis of the dripping-jetting transition in compound capillary jets.J Fluid Mech,2010,649:523-536

24 Li F,Yin XY,Yin XZ.Instability of a viscous co fl owing jet in a radial electric fi eld.J Fluid Mech,2008,596:285-311

25 李广滨,司廷,尹协振.电场作用下无黏聚焦射流的时间不稳定性研究.力学学报,2012,44(5):876-883(Li Guangbin,Si Ting,Yin Xiezhen.Temporal instability study of an inviscid focused jet under an electric fi eld.Chinese J Theo Appl Mech,2012,44(5):876-883(in Chinese))

26 李广滨.复合流动聚焦的实验和理论研究.[博士论文].合肥:中国科学技术大学,2016(Li Guangbin.Experimental and theoretical investigation on co- fl ow focusing.[PhD Thesis].Hefei:University of Science and Technology of China,2016(in Chinese))

27 L´opez-Herrera JM,Barrero A,L´opez A,et al.Coaxial jets generated from electri fi ed Taylor cones.Scaling laws.J Aerosol Sci,2003,34:535-552

INSTABILITY STUDY OF AN ELECTRIFIED COAXIAL JET IN A COFLOWING GAS STREAM1)

Li Shuaibing Yang Rui Luo Xisheng Si Ting2)
(Department of Modern Mechanics，University of Science and Technology of China，Hefei 230027，China)

Instability study of electri fi ed coaxial jet coupling the electric and inertial forces is performed based on the simpli fi ed experimental model of gas-driven coaxial electro- fl ow focusing.Under the assumption that the fl uids are inviscid,incompressible and irrotational,a triple-layer electri fi ed fl uid jet model is established and an analytical dispersion relation in the temporal regime is obtained.The dispersion equation is solved by the normal mode method,the unstable modesofthe fl owarecalculatedand thee ff ectsof mainly controllable parameterson the unstable modes areanalyzed.The results indicate that the axisymmetric mode dominates the complete fl ow as the maximum growth rate of the axisymmetric modeisthelargestamongallunstablemodes.Asthevelocityofoutergasstreamincreases,theinertialforcecande fi nitely promote the jet instability.The jet will become more unstable as the velocity di ff erence between the inner and outer liquid jets increases.The surface tension also promotes the jet instability.The axial electric fi eld has two-fold in fl uence on the axisymmetric jet instabilities.There is a critical value for the axial electric voltage which is related to the free electric charge density at the interface and the perturbation propagations on the jet surfaces.The applied axial electric fi eld can suppress the jet instability when its intensity is smaller than the critical value;otherwise,the applied axial electric fi eld can promote jet instability.These results are in good agreement with the existing experimental results and can provide guidance on the process control of experiments.

coaxial jet,electri fi ed jet,electro- fl ow focusing,jet instability,process control

O361.5

A

10.6052/0459-1879-17-082

2017–03–12收稿，2017–05–26 录用,2017–06–05 网络版发表.

1)国家自然科学基金(11472270,81327803,11621202)和中央高校基本科研业务费专项资金资助项目.

2)司廷，特任教授，研究方向：实验流体力学、微纳尺度流动、流动稳定性.E-mail:tsi@ustc.edu.cn

李帅兵,杨睿,罗喜胜,司廷.气流作用下同轴带电射流的不稳定性研究.力学学报,2017,49(5):997-1007

Li Shuaibing,Yang Rui,Luo Xisheng,Si Ting.Instability study of an electri fi ed coaxial jet in a co fl owing gas stream.Chinese Journal of Theoretical and Applied Mechanics,2017,49(5):997-1007