EI、Scopus 收录
中文核心期刊

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

李帅兵, 杨睿, 罗喜胜, 司廷

李帅兵, 杨睿, 罗喜胜, 司廷. 气流作用下同轴带电射流的不稳定性研究[J]. 力学学报, 2017, 49(5): 997-1007. DOI: 10.6052/0459-1879-17-082
引用本文: 李帅兵, 杨睿, 罗喜胜, 司廷. 气流作用下同轴带电射流的不稳定性研究[J]. 力学学报, 2017, 49(5): 997-1007. DOI: 10.6052/0459-1879-17-082
Li Shuaibing, Yang Rui, Luo Xisheng, Si Ting. INSTABILITY STUDY OF AN ELECTRIFIED COAXIAL JET IN A COFLOWING GAS STREAM[J]. Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(5): 997-1007. DOI: 10.6052/0459-1879-17-082
Citation: Li Shuaibing, Yang Rui, Luo Xisheng, Si Ting. INSTABILITY STUDY OF AN ELECTRIFIED COAXIAL JET IN A COFLOWING GAS STREAM[J]. Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(5): 997-1007. DOI: 10.6052/0459-1879-17-082
李帅兵, 杨睿, 罗喜胜, 司廷. 气流作用下同轴带电射流的不稳定性研究[J]. 力学学报, 2017, 49(5): 997-1007. CSTR: 32045.14.0459-1879-17-082
引用本文: 李帅兵, 杨睿, 罗喜胜, 司廷. 气流作用下同轴带电射流的不稳定性研究[J]. 力学学报, 2017, 49(5): 997-1007. CSTR: 32045.14.0459-1879-17-082
Li Shuaibing, Yang Rui, Luo Xisheng, Si Ting. INSTABILITY STUDY OF AN ELECTRIFIED COAXIAL JET IN A COFLOWING GAS STREAM[J]. Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(5): 997-1007. CSTR: 32045.14.0459-1879-17-082
Citation: Li Shuaibing, Yang Rui, Luo Xisheng, Si Ting. INSTABILITY STUDY OF AN ELECTRIFIED COAXIAL JET IN A COFLOWING GAS STREAM[J]. Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(5): 997-1007. CSTR: 32045.14.0459-1879-17-082

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

基金项目: 

国家自然科学基金 11472270

国家自然科学基金 81327803

国家自然科学基金 11621202

中央高校基本科研业务费专项资金资助项目 

详细信息
    通讯作者:

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

  • 中图分类号: O361.5

INSTABILITY STUDY OF AN ELECTRIFIED COAXIAL JET IN A COFLOWING GAS STREAM

  • 摘要: 通过对气体驱动同轴电流动聚焦的实验模型进行简化,开展了电场力和惯性力共同作用下同轴带电射流的不稳定性理论研究.在流动为无黏、不可压缩、无旋的假设下,建立了三层流体带电射流物理模型并得到了扰动在时间域内发展演化的解析形式色散关系,利用正则模方法求解色散方程发现了流动的不稳定模态,进而分析了主要控制参数对不稳定模态的影响.结果表明,在参考状态下轴对称模态的最不稳定增长率最大,因此轴对称扰动控制整个流场.外层气流速度越高,气体惯性力越大,射流的界面越容易失稳.内外层液-液同轴射流之间的速度差越大,射流越不稳定.表面张力对射流不稳定性起到促进作用.轴向电场对射流不稳定性具有双重影响:当加载电场强度较小时,射流不稳定性被抑制;当施加电压大于某一临界值时,轴向电场会促进射流失稳.临界电压的大小与界面上自由电荷密度和射流表面扰动发展关系密切.这些结果与已有的实验现象吻合,能够对实验的过程控制提供理论指导.
    Abstract: Instability study of electrified coaxial jet coupling the electric and inertial forces is performed based on the simplified experimental model of gas-driven coaxial electro-flow focusing. Under the assumption that the fluids are inviscid, incompressible and irrotational, a triple-layer electrified fluid 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 modes of the flow are calculated and the effects of mainly controllable parameters on the unstable modes are analyzed. The results indicate that the axisymmetric mode dominates the complete flow as the maximum growth rate of the axisymmetric mode is the largest among all unstable modes. As the velocity of outer gas stream increases, the inertial force can definitely promote the jet instability. The jet will become more unstable as the velocity difference between the inner and outer liquid jets increases. The surface tension also promotes the jet instability. The axial electric field has two-fold influence 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 field can suppress the jet instability when its intensity is smaller than the critical value; otherwise, the applied axial electric field 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.
  • 电流体力学(electrohydrodynamics, EHD)是流体力学和电动力学之间的交叉学科.随着纳米技术、微机电系统和生物技术的迅猛发展,EHD研究在纳米材料、微流动、生物医学等领域显示了巨大的应用价值.电雾化(electrospray,ES)是EHD的一个重要分支,可描述为利用单轴或复合结构的毛细管,在电场力作用下促使流体界面拉伸成微细射流[1-3].在ES中,带电液体射流周围的气体是静止的.不同于ES,近年来另一种有效的微细液体射流产生方法被提出,称之为流动聚焦(flow focusing,FF),也就是从毛细管流出的液体由高速运动的气体驱动,经小孔聚焦后形成稳定的锥射流结构[3-5].可以看出,在FF中液体射流周围的气体是高速运动的. FF和ES形成的液体射流将破碎成微纳米量级的液滴、颗粒或胶囊,具有相似的射流不稳定性机理. Gañán-Calvo等[6]和司廷等[7]将FF与ES相结合,提出了电流动聚焦(electro-flow focusing,EFF)技术,能够使雾化更稳定,制备的产物尺寸更小,因而具有重要的应用前景.

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

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

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

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

    图  1  气流作用下液-液同轴带电射流的简化模型,其中1, 2, 3分别表示内层、中间层和外层流体
    Figure  1.  Simplified model of electrified coaxial liquid-liquid jet in a coflowing gas stream. 1, 2 and 3 stand for the inner, middle and outer fluids

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

    下面从量纲分析的角度解释忽略重力、体电荷以及流体黏性的合理性.在流体力学中无量纲参数 $Fr$ 数(惯性力/重力)以及 $Bo$ 数(重力/表面张力)与重力相关.取一组同轴射流实验[26]的参数作为参考状态:速度 $U_{0}=2 $ m/s;长度 $L_{0}=30 \mu $ m;密度 $\rho_{0}=1 000$ kg/m $^{3}$ ;黏度 $\mu_{0}=10^{-3}$ kg/(s $\cdot$ m);气体压力差 $\Delta p_{\rm g}=5$ kPa,液体流量速度 $Q_{0}=50$ mL/h,表面张力系数 $\gamma _{0}=7.17×10^{-2}$ N/m, 重力加速度 $g=9.8$ m/s $^{2}$ .可得

    $$Fr = \dfrac{U_0^2 }{gL_0 } = 1.36×10^4 , \ \ Bo = \dfrac{\rho _0 L_0^2 g}{\gamma _0 } = 1.23×10^{ - 4} , \\ Re = \dfrac{2}{\mu _0 }\left( {\dfrac{\rho _0^3 \Delta p_g }{8\pi ^2}} \right)^{0.25} Q_0^{0.5} \approx 132$$

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

    水的电导率为 $\sigma_{\rm w}=1.4×10^{-4}$ S/m,介电常数 $\varepsilon _{\rm w}\varepsilon_{0} \approx 7.1 × 10^{ - 10}$ C $^{2}$ /(N $\cdot$ m $^{2}$ ).可得电松弛时间以及对流运动特征时间为

    $$\tau _{\rm e} = \dfrac{\varepsilon _{\rm w} \varepsilon _0 }{\sigma _{\rm w} } = 5.07×10^{ - 6} {\rm s} , \ \ \tau _{\rm F} = \dfrac{L_0 }{U_0 } = 1.5×10^{ - 5} {\rm s}$$

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

    由于假设流体为不可压流体,射流满足连续性方程为: $\nabla \cdot {\pmb u}_i = 0$ .其中 $i=1$ , 2, 3分别表示内层、中间层以及最外层流体.因为流动无旋,故可引入速度势函数 $\varphi_{i}$

    $$\Delta \varphi _i = 0$$ (1)

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

    $$\rho _i \left( {\dfrac{\partial {\pmb u}_i }{\partial t} + {\pmb u}_i \cdot \nabla {\pmb u}_i } \right) = - \nabla p_i$$ (2)

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

    $$\nabla \cdot {\pmb D}_i = 0$$ (3)
    $$\nabla ×{\pmb E}_i ={\bf 0}$$ (4)

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

    (1) 射流对称轴上的扰动量有限以及无穷远处的扰动为零

    $$\left. {\varphi'_1 } \right|_{r \to 0} < \infty , \ \ \left. {\varphi'_3 } \right|_{r \to \infty } = 0$$ (5)
    $$\left. {V'_1 } \right|_{r \to 0} < \infty , \ \ \left. {V'_3 } \right|_{r \to \infty } = 0$$ (6)

    其中 $\varphi'_{i}$ 为速度势扰动量, $V'_{i}$ 为电势扰动量.

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

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

    $$\dfrac{{\rm{d}}\left( {r_i - \eta _i \left( {\theta, z, t} \right)} \right)}{ {\rm{d}} t} = 0$$

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

    $$u'_{ir} = \dfrac{\partial \eta_i }{\partial t} + U_{iz} \dfrac{\partial \eta _i }{\partial z}$$ (7)

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

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

    $$\gamma \left( {\nabla \cdot {\pmb n}} \right){\pmb n} = \left( {{\pmb T}_2 - {\pmb T}_1 } \right) \cdot {\pmb n}$$ (8)

    同样在界面2上有

    $$\gamma \left( {\nabla \cdot {\pmb n}} \right){\pmb n} = \left( {{\pmb T}_3 - {\pmb T}_2 } \right) \cdot {\pmb n}$$ (9)

    其中 ${\pmb n}$ 为界面上法向单位矢量, ${\pmb T}$ 为水动力学张量 ${\pmb T}^{h}$ 与电应力张量 ${\pmb T}^{e}$ 之和 ${\pmb T}={\pmb T}^{h}+{\pmb T}^{e}$

    $${\pmb n} = {\pmb e}_r - \dfrac{1}{r}\dfrac{\partial \eta }{\partial \theta }{\pmb e}_\theta - \dfrac{\partial \eta }{\partial z}{\pmb e}_z \\ {\pmb T}^h = - p{\pmb \delta } + \mu \left[{\nabla {\pmb u} + \left( {\nabla {\pmb u}} \right)^{\rm T}} \right] \\ T_{ik}^e = \varepsilon E_k E_i - \dfrac{\varepsilon \delta _{ki} }{2}E_m E_m$$

    其中 $\varepsilon $ 为介电常数, $\delta $ 为Kronecker函数.

    将 ${\pmb T}^{h}$ 以及 ${\pmb T}^{e}$ 代入平衡方程(8) 和式(9) 得

    $$\gamma \left( {\nabla \cdot {\pmb n}_1 } \right) = p_1 - p_2 + {\pmb n}_1 \cdot \left( {{\pmb T}_2^e - {\pmb T}_1^e } \right) \cdot {\pmb n}_1$$ (10)
    $$\gamma \left( {\nabla \cdot {\pmb n}_2 } \right) = p_2 - p_3 + {\pmb n}_2 \cdot \left( {{\pmb T}_3^e - {\pmb T}_2^e } \right) \cdot {\pmb n}_2$$ (11)

    (4) 电场边界条件

    ${\pmb E}_{i}$ ( $i=1$ , 2, 3) 分别表示内层、中间层以及外层流体中的电场强度, $\varepsilon_{i}$ 表示介质的绝对介电常数, ${\pmb D}_{i}=\varepsilon_{i}{\pmb E}_{i}$ 表示电位移矢量, $q_{\rm s}$ 为界面1上的自由电荷密度, $\sigma_{i}$ 为物质的电导率, ${\pmb J}_{i} =\sigma_{i}{\pmb E}_{i}$ 表示电流密度.忽略磁场作用,由法拉第感应定律得

    $${\pmb n}_1 ×\left( {{\pmb e}_2 - {\pmb e}_1 } \right) = {\bf 0} , \quad {\pmb n}_2 ×\left( {{\pmb e}_3 - {\pmb e}_2 } \right) = {\bf 0}$$ (12)

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

    由高斯定律得

    $${\pmb n}_1 \cdot \left( {{\pmb D}_2 - {\pmb D}_1 } \right) = q_{\rm s} , \quad {\pmb n}_2 \cdot \left( {{\pmb D}_3 - {\pmb D}_2 } \right) = 0$$ (13)

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

    $${\pmb n}_1 \cdot \left( {{\pmb J}_2 - {\pmb J}_1 } \right) = 0 , \quad {\pmb n}_2 \cdot \left( {{\pmb J}_3 - {\pmb J}_2 } \right) = 0$$ (14)

    假设流动无黏,采用Kelvin-Helmholtz速度型, ${\pmb U}_{i}(r, \theta, z)=(0, 0 U_{z})$ .

    基本电场为轴向均匀电场: ${\pmb E}_{i}= E_{0}{\pmb e}_{z}$ , $i=1$ , 2, 3.

    当界面受到任意小扰动时,它们将离开平衡位置发生微小的位移,新位置可以表示为: $r_i = R_i + \eta _i $ ( $i=1$ , 2表示界面1和界面2).

    内层、中间层以及外层流体中,流体速度、压力、电势等物理量可以分解成基本量加扰动量的形式

    $${\pmb u}_i = \left( {u_{ir}, u_{i\theta }, u_{iz} } \right) = \left( {u'_{ir}, u'_{i\theta }, U_{iz} + u'_{iz} } \right) \\ \left( {p_i, V_i, \varphi _i, q_s } \right) = \left( {P_{i0} + p'_i, \varPsi _{i0} + V'_i, \varPhi _i + \varphi'_i, Q_{is} + q'_s } \right)$$

    式中 $u_{ir}$ , $u_{i\theta }$ , $u_{iz}$ , $p_{i}$ , $V_{i}$ , $\varphi_{i}$ 分别表示流体 $i$ ( $i=1$ , 2, 3分别表示内层、中间层和外层流体)中的径向速度分量、周向速度分量、轴向速度分量、压力、电势以及速度势. $u'_{ir}$ , $u'_{i\theta }$ , $u'_{iz}$ , $p'_{i}$ , $V'_{i}$ , $\varphi'_{i}$ , $ q'_{is}$ 为相应物理量的扰动量; $U_{iz}$ , $P_{i0}$ , $V_{i0}$ , $\varPhi_{i}$ , $Q_{is}$ 分别表示无扰动情况下的轴向速度分量、压力、电势、速度势、自由电荷密度.

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

    $$\left( {\eta _i, u'_{ir}, u'_{i\theta }, u'_{iz}, p'_i, V'_i, \varphi'_i } \right) = \\ \qquad \left( {\overline {\eta _i }, \overline {u_{ir} }, \overline {u_{i\theta } }, \overline {u_{iz} }, \overline {p_i }, \overline {V_i }, \overline {\varphi _i } } \right) {\rm e}^\Re$$

    其中 $\Re = \omega t + {\rm i}\left( {kz + n\theta } \right)$ .扰动量是 $r$ 的函数.

    由 $\varphi_{i} = \varPhi_{i} + \varphi'_{i} $ 以及方程(1) 和式(5) 得

    $${\pmb u}_1 = U_{1z} {\pmb e}_z + \nabla \left[{C_1 I_n \left( {kr} \right) {\rm e}^\Re } \right] \\ {\pmb u}_2 = U_{2z} {\pmb e}_z + \nabla \left[{\left( {C_2 I_n \left( {kr} \right) + C_3 K_n \left( {kr} \right)} \right) {\rm e}^\Re } \right] \\ {\pmb u}_3 = U_{3z} {\pmb e}_z + \nabla \left[{\left( {C_4 I_n \left( {kr} \right) + C_5 K_n \left( {kr} \right)} \right) {\rm e}^\Re } \right]$$

    其中

    $$C_1 = \dfrac{{\rm i}kU_{1z} + \omega }{kI'_n \left( {kR_1 } \right)}\overline {\eta _1 } , \ C_4 = 0 , \ C_5 = \dfrac{{\rm i}kU_{3z} + \omega }{kK'_n \left( {kR_2 } \right)}\overline {\eta _2 } \\ C_2 = \dfrac{{\rm i}kU_{2z} + \omega }{k}\dfrac{K'_n \left( {kR_2 } \right)\overline {\eta _1 } - K'_n \left( {kR_1 } \right)\overline {\eta _2 } }{I'_n \left( {kR_1 } \right)K'_n \left( {kR_2 } \right) - I'_n \left( {kR_2 } \right)K'_n \left( {kR_1 } \right)} \\ C_3 = \dfrac{{\rm i}kU_{2z} + \omega }{k}\dfrac{I'_n \left( {kR_1 } \right)\overline {\eta _2 } - I'_n \left( {kR_2 } \right)\overline {\eta _1 } }{I'_n \left( {kR_1 } \right)K'_n \left( {kR_2 } \right) - I'_n \left( {kR_2 } \right)K'_n \left( {kR_1 } \right)}$$

    将 ${\pmb u}_{i}$ 代入动量方程(2) 得

    $$p_1 = - \rho _1 \left[{\dfrac{1}{2}U_{1z}^2 + C_1 \left( {{\rm i}kU_{1z} + \omega } \right)I_n \left( {kr} \right) {\rm e}^\Re } \right] \\ p_2 = - \rho _2 \left[{\dfrac{U_{2z}^2 }{2} + \left( {{\rm i}kU_{2z} + \omega } \right)\left( {C_2 I_n \left( {kr} \right) + C_3 K_n \left( {kr} \right)} \right) {\rm e}^\Re } \right] \\ p_3 = - \rho _3 \left[{\dfrac{1}{2}U_{3z}^2 + C_5 \left( {{\rm i}kU_{3z} + \omega } \right)K_n \left( {kr} \right) {\rm e}^\Re } \right]$$

    由式(4) 可引入势函数 $V_{i}$ : ${\pmb E}_i = - \nabla V_i \ \left( {i = 1, 2, 3} \right)$ 又由式(3) 可得

    $$\Delta V_i = 0 \ \ \left( {i = 1, 2, 3} \right)$$ (15)

    求解式 $V_i = \varPsi _{i0} + \overline {V_i } {\rm e}^\Re $ 及式(6) 式(15) 可得

    $$\overline {V_1 } \left( r \right) = G_1 I_n \left( {kr} \right)$$ (16)
    $$\overline {V_2 } \left( r \right) = G_2 I_n \left( {kr} \right) + G_3 K_n \left( {kr} \right)$$ (17)
    $$\overline {V_3 } \left( r \right) = G_5 K_n \left( {kr} \right)$$ (18)

    将 $\overline {V_i } \left( r \right)$ 代入势方程 ${\pmb E}_i = - \nabla V_i $ ,可得扰动后的电场

    $${\pmb E}_1 = E_0 {\pmb e}_z - \nabla \left[{G_1 I_n \left( {kr} \right){\rm e}^\Re } \right]$$ (19)
    $${\pmb e}_2 = E_0 {\pmb e}_z - \nabla \left[{\left( {G_2 I_n \left( {kr} \right) + G_3 K_n \left( {kr} \right)} \right) {\rm e}^\Re } \right]$$ (20)
    $${\pmb e}_3 = E_0 {\pmb e}_z - \nabla \left[{G_5 K_n \left( {kr} \right) {\rm e}^\Re } \right]$$ (21)

    由式(12) ~式(14) 得

    $$G_1 = - \dfrac{{\rm i}\overline {\eta _1 } E_0 }{I'_n \left( {kR_1 } \right)}, G_6 = \dfrac{K'_n \left( {kR_2 } \right)}{K_n \left( {kR_2 } \right)}$$
    $$G_5 = \dfrac{\varepsilon _2 }{\varepsilon _3 }\left( {\dfrac{G_2 I'_n \left( {kR_2 } \right)}{K'_n \left( {kR_2 } \right)} + \dfrac{G_3 K'_n \left( {kR_2 } \right)}{K'_n \left( {kR_2 } \right)}} \right) - \dfrac{{\rm i}k\overline {\eta _2 } \left( {\varepsilon _3 - \varepsilon _2 } \right)E_0 }{k\varepsilon _3 K'_n \left( {kR_2 } \right)}$$
    $$G_2 = \Bigg [-\dfrac{{\rm i}\overline {\eta _1 } E_0 I_n \left( {kR_1 } \right)}{I'_n \left( {kR_1 } \right)}\left( {\dfrac{\varepsilon _2 }{\varepsilon _3 }-1} \right)-$$
    $$\qquad {\rm i}\left( {1 - \dfrac{\varepsilon _2 }{\varepsilon _3 }} \right)\overline {\eta _2 } E_0 \dfrac{K_n \left( {kR_1 } \right)}{ K'_n \left( {kR_2 } \right)}\Bigg] \Bigg / $$
    $$\qquad \Bigg [I_n \left( {kR_1 } \right)\left( {\dfrac{\varepsilon _2 }{\varepsilon _3 }-1} \right)-\dfrac{K_n \left( {kR_1 } \right)}{K'_n \left( {kR_2 } \right)}\cdot $$
    $$\qquad \left( {\dfrac{\varepsilon _2 }{\varepsilon _3 }I'_n \left( {kR_2 } \right)-I_n \left( {kR_2 } \right)G_6 } \right) \Bigg] $$
    $$ G_3 = \Bigg [{\rm i}\left( {1-\dfrac{\varepsilon _2 }{\varepsilon _3 }} \right)\overline {\eta _2 } E_0 + $$
    $$\qquad \left( {\dfrac{\varepsilon _2 }{\varepsilon _3 }I'_n \left( {kR_2 } \right)-I_n \left( {kR_2 } \right)G_6 } \right)\dfrac{{\rm i}\overline {\eta _1 } E_0 }{I'_n \left( {kR_1 } \right)}\Bigg] \Bigg /$$
    $$\qquad \Bigg [\left( {\dfrac{\varepsilon _2 }{\varepsilon _3 }-1} \right)K'_n \left( {kR_2 } \right)-\dfrac{K_n \left( {kR_1 } \right)}{I_n \left( {kR_1 } \right)}\cdot $$
    $$\qquad \left( {\dfrac{\varepsilon _2 }{\varepsilon _3 }I'_n \left( {kR_2 } \right)-I_n \left( {kR_2 } \right)G_6 } \right) \Bigg] $$

    选择中间层射流半径 $R_{2}$ ,密度 $\rho_{2}$ ,无扰动速度 $U_{2}$ ,表面张力系数 $\gamma _{2}$ 以及空气介电常数 $\varepsilon_{3}$ 作为特征尺度.设定无量纲参数: $\alpha =R_{2}k$ 为无量纲波数; $\beta =\omega R_{2}/U_{2}$ 为无量纲复频率; $\rho _{12}=\rho _{1}/\rho_{2}$ 为内层流体与中间层流体密度比; $\rho_{32}=\rho _{3}/\rho _{2}$ 为外层流体与中间层流体密度比; $U_{12}=U_{1z}/U_{2z}$ 为内层流体与中间层流体的速度比; $U_{32}= U_{3z}/U_{2z}$ 为外层流体与中间层流体的速度比; $a=R_{1}/R_{2}$ 为无扰动时界面1与界面2的位置比; $\varGamma=\gamma_{1}/\gamma_{2}$ 为界面1与界面2的张力系数比; $\varepsilon_{13}=\varepsilon_{1}/\varepsilon_{3}$ 为内层流体与外层流体的介电常数比; $\varepsilon _{23}=\varepsilon_{2}/\varepsilon_{3}$ 为中间层流体与外层流体的介电常数比; $We=\rho_{2}U_{2}^{2}R_{2}/\gamma_{2}$ 为Weber数; $Eu=\varepsilon _{3}E_{0}^{2}/\rho_{2}U_{2}^{2}$ 为电Euler数.

    将 $p_{1}$ , $p_{2}$ , $p_{3}$ , ${\pmb E}_{1}$ , ${\pmb E}_{2}$ , ${\pmb E}_{3}$ 代入动力学边界条件(10) 和(11),经过运算,最终得到无量纲化的色散方程

    $$- \left( {{\rm i}\alpha + \beta } \right)^4{H}_1 - \left( {{\rm i}\alpha U_{12} + \beta } \right)^2\left( {{\rm i}\alpha U_{32} + \beta } \right)^2{H}_2 + \\ \qquad \left( {{\rm i}\alpha + \beta } \right)^2\left[{\left( {{\rm i}\alpha U_{12} + \beta } \right)^2{H}_3 + \left( {{\rm i}\alpha U_{32} + \beta } \right)^2{H}_4 } \right] - \\ \qquad ( {\rm i}\alpha U_{12} + \beta )^2{H}_5 + \left( {{\rm i}\alpha U_{32} + \beta } \right)^2{H}_6 + \\ \qquad \left( {{\rm i}\alpha + \beta } \right)^2{H}_7 + {H}_8 = 0$$ (22)

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

    $$\dfrac{\overline {\eta _1 } }{\overline {\eta _2 } } = \dfrac{\aleph _4 + \aleph _2 }{\aleph _5 + \aleph _3 } = \left| {\dfrac{\overline {\eta _1 } }{\overline {\eta _2 } }} \right|{\rm e}^{{\rm i}\left( {\phi _1 - \phi _2 } \right)}$$ (23)

    其中 $\phi _1 - \phi _2 $ 表示界面运动的相位差, $\phi _1 - \phi _2 = 0$ 表示界面同向运动, $\phi _1 - \phi _2 = \pi $ 表示界面反向运动.

    $${H}_1 = \dfrac{\varDelta _3 \varDelta _4 }{\varDelta _5 } + \dfrac{\varDelta _1 \varDelta _6 }{\varDelta _5 } , \ \ {H}_2 = \dfrac{\rho _{12} \rho _{32} \varDelta _5 }{\varDelta _7 \varDelta _8 } $$
    $${H}_3 = \dfrac{\rho _{12} \varDelta _3 \varDelta _5 }{\varDelta _5 \varDelta _7 } , \ \ {H}_4 = \dfrac{\rho _{32} \varDelta _4 }{\varDelta _8 } , \ \ {H}_5 = \dfrac{\rho _{12} \varDelta _5 \aleph _2 }{\varDelta _7 } $$
    $${H}_6 = \dfrac{\rho _{32} \aleph _1 }{\varDelta _8 } , \ \ {H}_7 = \varDelta _4 \aleph _2 - \dfrac{\varDelta _3 }{\varDelta _5 }\aleph _1 - \varDelta _1 \aleph _3 $$
    $${H}_8 = \aleph _1 \aleph _2 , \ \ \aleph _3 = \dfrac{\alpha ^2E_u \left( {\varepsilon _{23} - 1} \right)\varDelta _6 \varepsilon _{23} }{\left( {\varepsilon _{23} \varDelta _4 - \varDelta _8 \varDelta _2 } \right)\varDelta _7 }$$
    $$\aleph _1 = \dfrac{\varGamma \alpha }{Wea^2}\left[{1-n^2-\left( {\alpha a} \right)^2} \right]\varDelta _5 + \dfrac{\alpha ^2E_u }{\varDelta _7 }\left( {\varepsilon _{23} - \varepsilon _{13} } \right)\varDelta _5 \\ \aleph _2 = \dfrac{\alpha ^2\varDelta _2 E_u \left( {\varepsilon _{23} - 1} \right)^2}{\varepsilon _{23} \varDelta _4 - \varDelta _8 \varDelta _2 } + \alpha \dfrac {1 - n^2 - \alpha ^2 }{We} \\ \aleph _4 = \dfrac{\rho _{32} \left( { {\rm i}\alpha U_{32} + \beta } \right)^2}{\varDelta _8 } - \dfrac{\varDelta _3 }{\varDelta _5 }\left( {{\rm i}\alpha + \beta } \right)^2 \\ \aleph _5 = \left( {{\rm i}\alpha + \beta } \right)^2\dfrac{\varDelta _6 }{\varDelta _5 }$$
    $$\varDelta _1 = I_n \left( {a\alpha } \right)K'_n \left( {a\alpha } \right) - I'_n \left( {a\alpha } \right)K_n \left( {a\alpha } \right)\\ \varDelta _2 = I_n \left( {a\alpha } \right)K_n \left( \alpha \right) - I_n \left( \alpha \right)K_n \left( {a\alpha } \right)\\ \varDelta _3 = I'_n \left( {a\alpha } \right)K_n \left( \alpha \right) - I_n \left( \alpha \right)K'_n \left( {a\alpha } \right)\\ \varDelta _4 = I_n \left( {a\alpha } \right)K'_n \left( \alpha \right) - I'_n \left( \alpha \right)K_n \left( {a\alpha } \right)\\ \varDelta _5 = I'_n \left( {a\alpha } \right)K'_n \left( \alpha \right) - I'_n \left( \alpha \right)K'_n \left( {a\alpha } \right)\\ \varDelta _6 = I_n \left( \alpha \right)K'_n \left( \alpha \right) - I'_n \left( \alpha \right)K_n \left( \alpha \right)\\ \varDelta _7 = \dfrac{I'_n \left( {a\alpha } \right)}{I_n \left( {a\alpha } \right)} , \ \ \varDelta _8 = \dfrac{K'_n \left( \alpha \right)}{K_n \left( \alpha \right)}$$

    当 $E_{0}= U_{3z}=0$ 时色散方程转化为

    $$\dfrac{\overline {\eta _1 } }{\overline {\eta _2 } } = \Bigg [\dfrac{\rho _{32} \left( {\beta + {\rm i}\alpha U_{32} } \right)^2}{\varDelta _8 }-\dfrac{\left( {\beta + {\rm i}\alpha } \right)^2\varDelta _3 }{\varDelta _5 } +\\ \qquad \dfrac{\alpha }{We_2 }\left( {1-n^2-\alpha ^2} \right)\Bigg] \Bigg / \Bigg [ \dfrac{\left( {\beta + {\rm i}\alpha } \right)^2\varDelta _6 }{\Delta _5 } \Bigg] \\ \dfrac{\overline {\eta _2 } }{\overline {\eta _1 } } = \Bigg [\dfrac{\rho _{12} \left( {\beta + {\rm i}\alpha U_{12} } \right)^2}{\varDelta _7 }-\dfrac{\alpha \left( {1-n^2-a^2\alpha ^2} \right)}{We_1 a^2} - \\ \qquad \left( {\beta + {\rm i}\alpha } \right)^2\dfrac{\varDelta _4 }{\varDelta _5 } \Bigg] \Bigg / \Bigg [-\left( {\beta + i\alpha } \right)^2 \dfrac{\varDelta _1 }{\varDelta _5 } \Bigg]$$

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

    无量纲色散关系式(22) 是关于无量纲复频率β的四次代数方程,给定一组( $n, \alpha $ )可得两对四个复特征值β.由 $\Re = \omega t +{\rm i}\left( {kz + n\theta } \right)$ 及 $\beta =\omega R_{2}/U_{2}$ 可知,β的实部 $\beta_{r}$ 为正值对应于射流的不稳定状态.

    López-Herrera等[27]以水作为导体以及葵花籽油作为电介质开展了同轴ES实验.相应试剂的物理参数见表 1.葵花籽油与空气之间界面表面张力系数 $\gamma_{2}$ 为32.8 mN/m,葵花籽油与水之间界面张力系数 $\gamma_{1}$ 为16.6 mN/m.实验装置包括内外直径为200 $\mu$ m和400 $\mu$ m的内部针管以及内外直径为420 μm和800 μm的外部针管.这样内层流体为半径约100 μm的圆柱射流,中间层流体为半径约10 μm的环形射流.作为中间层流体的葵花籽油流速 $U_{2} \approx 1$ m/s,又因 $R_{2} \approx 110\mu $ m.加在同轴射流上的电压约为1 kV,电极距离约为1 cm[19],可得 $E_{0} \approx 10^{5}$ V/m.空气的绝对介电常数为 $\varepsilon_{\rm air}=8.85×10^{-12}$ C $^{2}$ /(N $\cdot$ m $^{2})$ .根据实验参数,本文拟设定无量纲参数的参考值近似为 $\rho _{12}=1.19$ , $\rho_{32}=0.001$ , $\varepsilon _{13}=80$ , $\varepsilon _{23}=3.4$ , $\varGamma=0.51$ , $a=0.9$ , $U_{12}=1.25$ , $U_{32}=2$ , $We=3$ , $Eu=10^{ - 4}$ .

    表  1  实验试剂的物理特性
    Table  1.  Physical property of the tested material
    下载: 导出CSV 
    | 显示表格

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

    图  2  周向波数 $n$ 对扰动增长率 $\beta_{r}$ 的影响
    Figure  2.  Effect of the azimuthal wave number $n$ on the temporal growth rate $\beta_{r}$ of perturbations

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

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

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

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

    图  4  电欧拉数 $Eu$ 对扰动增长率 $\beta_{r}$ 的影响
    Figure  4.  Effect of the electric Euler number $Eu$ on the temporal growth rate $\beta_{r}$

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

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

    图  5  外层流体与中间层流体流速比 $U_{32}$ 对增长率 $\beta_{r}$ 的影响
    Figure  5.  Effect of the velocity ratio $U_{32}$ between the outer layer fluid and middle layer fluid on the growth rate $\beta_{r}$

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

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

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

    图  6  同轴射流间速度比 $U_{12}$ 对增长率 $\beta_{r}$ 的影响
    Figure  6.  Effect of the velocity ratio $U_{12}$ between the coaxial jet on the growth rate $\beta_{r}$

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

    当内层流体流速相对于中间层流体流速较大时( $U_{12}>1$ ),在长波区出现类对称模,短波区出现过渡模.如图 6(c)所示,随着流速 $U_{12}$ 的增大,过渡模的影响愈来愈显著,最大扰动增长率 $\beta_{r\max}$ 变大,射流不稳定性由长波不稳定性向短波不稳定性转化.

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

    韦伯数 $We = \rho_{2}U_{2}^{2}R_{2}/\gamma_{2}$ 与增长率 $\beta_{r}$ 之间的关系反映了中间层流体与外层流体之间界面(即界面2) 上表面张力对射流稳定性影响.

    设界面1上的韦伯数 $We_{1}=\rho_{1}U_{1}^{2}R_{1}/\gamma_{1}$ ,又因界面2上韦伯数为

    $$We_1 = We\dfrac{aU_{12}^2 \rho _{12} }{\varGamma }$$ (25)

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

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

    图  7  韦伯数 $We$ 对增长率 $\beta_{r}$ 的影响
    Figure  7.  Effect of the Weber number $We$ on the growth rate $\beta_{r}$

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

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

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

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

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

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

  • 图  1   气流作用下液-液同轴带电射流的简化模型,其中1, 2, 3分别表示内层、中间层和外层流体

    Figure  1.   Simplified model of electrified coaxial liquid-liquid jet in a coflowing gas stream. 1, 2 and 3 stand for the inner, middle and outer fluids

    图  2   周向波数 $n$ 对扰动增长率 $\beta_{r}$ 的影响

    Figure  2.   Effect of the azimuthal wave number $n$ on the temporal growth rate $\beta_{r}$ of perturbations

    图  3   (a)类反对称模, (b)类对称模, (c)过渡模

    Figure  3.   (a) parasinuous mode, (b) paravaricose mode, (c) transitional mode

    图  4   电欧拉数 $Eu$ 对扰动增长率 $\beta_{r}$ 的影响

    Figure  4.   Effect of the electric Euler number $Eu$ on the temporal growth rate $\beta_{r}$

    图  5   外层流体与中间层流体流速比 $U_{32}$ 对增长率 $\beta_{r}$ 的影响

    Figure  5.   Effect of the velocity ratio $U_{32}$ between the outer layer fluid and middle layer fluid on the growth rate $\beta_{r}$

    图  6   同轴射流间速度比 $U_{12}$ 对增长率 $\beta_{r}$ 的影响

    Figure  6.   Effect of the velocity ratio $U_{12}$ between the coaxial jet on the growth rate $\beta_{r}$

    图  7   韦伯数 $We$ 对增长率 $\beta_{r}$ 的影响

    Figure  7.   Effect of the Weber number $We$ on the growth rate $\beta_{r}$

    表  1   实验试剂的物理特性

    Table  1   Physical property of the tested material

    下载: 导出CSV
  • [1]

    Fernández de la Mora J. The fluid dynamics of Taylor cones. Annu Rev Fluid Mech, 2007, 39:217-243 doi: 10.1146/annurev.fluid.39.050905.110159

    [2]

    Melcher JR, Taylor GI. Electrohydrodynamics:a review of the role of interfacial shear stresses. J Fluid Mech, 1969, 1:111-146 doi: 10.1146/annurev.fl.01.010169.000551

    [3]

    Barrero A, Loscertales IG. Micro-and nanoparticles via capillary flows. Annu Rev Fluid Mech, 2007, 39:89-106 doi: 10.1146/annurev.fluid.39.050905.110245

    [4]

    Gañán-Calvo AM. Generation of steady liquid microthreads and micron-sized monodisperse sprays in gas streams. Phys Rev Lett, 1998, 80:285-288 doi: 10.1103/PhysRevLett.80.285

    [5] 司廷, 李广滨, 尹协振.流动聚焦及射流不稳定性.力学进展, 2017, 47:201706 http://youxian.cnki.com.cn/yxdetail.aspx?filename=LXJZ20170113001&dbname=CAPJ2015

    Si Ting, Li Guangbin, Yin Xiezhen. Flow focusing and jet instability. Advances in Mechanics, 2017, 47:201706 (in Chinese) http://youxian.cnki.com.cn/yxdetail.aspx?filename=LXJZ20170113001&dbname=CAPJ2015

    [6]

    Gañán-Calvo AM, López-Herrera JM, Riesco-Chueca P. The combination of electrospray and flow focusing. J Fluid Mech, 2006, 566:421-445 doi: 10.1017/S0022112006002102

    [7] 司廷, 李广滨, 田瑞军等.电场作用下流动聚焦的实验研究.力学学报, 2011, 43(6):1030-1036 doi: 10.6052/0459-1879-2011-6-lxxb2011-156

    Si Ting, Li Guangbin, Tian Ruijun, et al. Experimental study of the flow focusing under an electric field. Chinese J Theo Appl Mech, 2011, 43(6):1030-1036 (in Chinese) doi: 10.6052/0459-1879-2011-6-lxxb2011-156

    [8] 司廷. 流动聚焦的实验和理论研究. [博士论文]. 合肥: 中国科学技术大学, 2009

    Si Ting. Experimental and theoretical investigation on flow focusing.[PhD Thesis]. Hefei:University of Science and Technology of China, 2009 (in Chinese)

    [9]

    Gañán-Calvo AM, Montanero JM, Martín-Banderas L, et al. Building functional materials for health care and pharmacy from microfluidic principles and flow focusing. Adv Drug Deliv Rev, 2013, 65:1447-1469 doi: 10.1016/j.addr.2013.08.003

    [10]

    Si T, Li F, Yin XY, et al. Modes in flow focusing and instability of coaxial liquid-gas jets. J Fluid Mech, 2009, 629:1-23 doi: 10.1017/S0022112009006211

    [11]

    Si T, Li F, Yin XY, et al. Spatial instability of coflowing liquid-gas jets in capillary flow focusing. Phys Fluids, 2010, 22:112105 doi: 10.1063/1.3490066

    [12]

    Li GB, Luo XS, Si T, et al. Temporal instability of coflowing liquidgas jets under an electric field. Phys Fluids, 2014, 26(5):054101 doi: 10.1063/1.4875109

    [13]

    Si T, Feng HX, Luo XS, et al. Formation of steady compound conejet modes and multilayered droplets in a tri-axial capillary flow focusing device. Microfluid. Nanofluid, 2015, 18:967-977 doi: 10.1007/s10404-014-1486-8

    [14]

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

    [15]

    Si T, Yin CS, Gao P, et al. Steady cone-jet mode in compound fluidic electro-flow focusing for fabricating multicompartment microcapsules. Appl Phys Lett, 2016, 108:021601 doi: 10.1063/1.4939632

    [16] 李芳. 同轴带电射流的稳定性研究. [博士论文]. 合肥: 中国科学技术大学, 2007

    Li Fang. Stability investigation on coaxial electrified jets.[PhD Thesis]. Hefei:University of Science and Technology of China, 2007 (in Chinese)

    [17]

    Gañán-Calvo AM, González-Prieto R, Riesco-Chueca P, et al. Focusing capillary jets close to the continuum limit. Nat Phys, 2007, 3(10):737-742 doi: 10.1038/nphys710

    [18]

    Li F, Yin XY, Yin XZ. Instability analysis of a coaxial jet under a radial electric field in the nonequipotential case. Phys Fluids, 2006, 18:037101 doi: 10.1063/1.2181604

    [19]

    Li F, Yin XY, Yin XZ. Linear instability of a coflowing jet under an axial electric field. Phys Rev E, 2006, 74:036304 doi: 10.1103/PhysRevE.74.036304

    [20]

    Loscertales IG, Barrero A, Guerrero I, et al. Micro/nano encapsulation via electrified coaxial liquid jets. Science, 2002, 295:1695-1698 doi: 10.1126/science.1067595

    [21]

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

    [22]

    Martín-Banderas L, Flores-Mosquera M, Riesco-Chueca P, et al. Flow focusing:a versatile technology to produce size-controlled and specific morphology microparticles. Small, 2005, 7:688-692 http://www.nsfc.gov.cn/Portals/0/fj/fj20160106_01.xls

    [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 doi: 10.1017/S0022112010000443

    [24]

    Li F, Yin XY, Yin XZ. Instability of a viscous coflowing jet in a radial electric field. J Fluid Mech, 2008, 596:285-311

    [25] 李广滨, 司廷, 尹协振.电场作用下无黏聚焦射流的时间不稳定性研究.力学学报, 2012, 44(5):876-883 doi: 10.6052/0459-1879-12-032

    Li Guangbin, Si Ting, Yin Xiezhen. Temporal instability study of an inviscid focused jet under an electric field. Chinese J Theo Appl Mech, 2012, 44(5):876-883 (in Chinese) doi: 10.6052/0459-1879-12-032

    [26] 李广滨. 复合流动聚焦的实验和理论研究. [博士论文]. 合肥: 中国科学技术大学, 2016

    Li Guangbin. Experimental and theoretical investigation on co-flow focusing.[PhD Thesis]. Hefei:University of Science and Technology of China, 2016 (in Chinese)

    [27]

    López-Herrera JM, Barrero A, López A, et al. Coaxial jets generated from electrified Taylor cones. Scaling laws. J Aerosol Sci, 2003, 34:535-552 doi: 10.1016/S0021-8502(03)00021-1

图(7)  /  表(1)
计量
  • 文章访问数:  1504
  • HTML全文浏览量:  337
  • PDF下载量:  767
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-03-11
  • 网络出版日期:  2017-06-04
  • 刊出日期:  2017-09-17

目录

/

返回文章
返回