高级检索

  • ISSN 1006-3080
  • CN 31-1691/TQ
引用本文:
Citation:

液体射流内部气泡对黏性流体气流式雾化的影响

    作者简介: 吴兆伟(1991-),男,安徽安庆人,博士生,研究方向为多相流气液雾化。E-mail:Y20150076@mail.ecust.edu.cn;
    通讯作者: 赵辉, zhaohui@ecust.edu.cn
  • 中图分类号: O359+.1

Effects of Bubbles in the Liquid Jet on theAir-Blast Atomization of Viscous Liquidd

    Corresponding author: Hui ZHAO, zhaohui@ecust.edu.cn ;
  • CLC number: O359+.1

  • 摘要: 使用马尔文激光粒度仪研究了同轴双通道喷嘴液体射流内部气泡对黏性流体气流式雾化的影响。以甘油水溶液(昂色格数Oh为2.14×10-3~0.407)作为液相,空气作为气相,研究了不同环隙气流速度(ug为94.8~142 m/s)和液体射流速度(ul,0为0.29~1.21 m/s)下雾化液滴直径的变化规律。实验发现:当Oh较小时,低液体射流速度条件下雾化液滴直径随气液体积比的增大而减小,高液体射流速度条件下雾化液滴直径随气液体积比的增大而增大;当Oh较大时,液滴直径随气液体积比的增大而增大。建立了一个物理模型用来描述黏性流体的雾化液滴直径随气液体积比的变化规律,模型计算值与实验结果之间的偏差小于±5%。
  • 图 1  实验装置示意图 1–Syringe pump; 2, 8, 9, 12–Valve; 3–Experimental atomizer; 4, 14–Gas cylinder; 5–Liquid tube rotameter; 6–Glass tube rotameter flowmeter; 7–Gauge; 10, 11–Malvern laser particle size analyzer; 13–Air blower; 15–Pressurized liquid vessel; 16–Pump; 17–Liquid storage; 18–Computer

    Figure 1.  Schematic of the experimental setup

    图 2  喷嘴示意图

    Figure 2.  Aerated air-blast atomizer geometry

    图 3  甘油水溶液的雾化液滴直径随环隙气流速度变化(ul,0=0.29 m/s)

    Figure 3.  Dependence of the droplet size D32 on the co-flow air stream velocity(ul,0=0.29 m/s)

    图 4  D32计算值和实验测量值之间的对比

    Figure 4.  Comparison of the experimental D32 to the calculated D32

    图 5  气泡直径随气液体积比的变化

    Figure 5.  Bubble sizes at varying V

    图 6  含气泡黏性流体雾化液滴直径随时间的变化图 (图中竖实线对应停止通气的时刻,ug=118 m/s,ul,0=0.29 m/s)

    Figure 6.  Time evolution of D32 of viscous liquid at different operation conditions. The vertical solid line is corresponding to the time stopping aerating the gas (ug=118 m/s, ul,0=0.29 m/s)

    图 7  D32随气液体积比的变化 (ug=118 m/s,ul,0=0.29 m/s)

    Figure 7.  Effects of changing V on D32 (ug=118 m/s, ul,0=0.29 m/s)

    图 8  D32随气液体积比V的变化 (ug=118 m/s,ul,0=0.29 m/s)

    Figure 8.  Dimensionless droplet size D32 at varying V (ug=118 m/s, ul,0=0.29 m/s)

    图 9  量纲为—液滴直径D随气液体积比的变化 (ul,0=0.29 m/s)

    Figure 9.  Dependence of dimensionless droplet size D on V (ul,0=0.29 m/s)

    图 10  量纲为一液滴直径D随气液体积比的变化 (ug=142 m/s)

    Figure 10.  Dimensionless droplet size D at varying V (ug=142 m/s)

    图 11  D模型计算值与实验测量值的对比

    Figure 11.  Comparison of the measured D to the calculated D

    表 1  不同液体物性参数表

    Table 1.  Properties of the different liquids

    Liquidρl / (kg·m−3)σ / (mN·m−1)μl / (mPa·s)Oh
    Water+Glycerol (w=60%)1 16066.916.73.46×10-2
    Water+Glycerol (w =75%)1 20066.556.30.115
    Water+Glycerol (w =80%)1 21066.01040.212
    Water+Glycerol (w =85%)1 23065.62000.407
    下载: 导出CSV

    表 2  甘油水溶液的雾化液滴直径随液体射流速度的变化(ug=142 m/s)

    Table 2.  Dependence of the droplet size D32 on liquid apparent velocity (ug=142 m/s)

    ul / (m·s−1)Measure value/μmPredicted value/μm
    Oh =3.46×10−2Oh =0.407Oh =3.46×10−2Oh =0.407
    0.2925.5±0.8830.6±0.5625.531.5
    0.3927.9±0.5033.2±0.5226.532.7
    0.5029.8±0.3634.0±0.3327.433.9
    下载: 导出CSV
  • [1] MARMOTTANT P, VILLERMAUX E. On spray formation[J]. Journal of Fluid Mechanics, 2004, 498: 73-111. doi: 10.1017/S0022112003006529
    [2] ZHAO H, WU ZW, LI WF, et al. Nonmonotonic effects of aerodynamic force on droplet size of prefilming air-blast atomization[J]. Industrial & Engineering Chemistry Research, 2018, 57(5): 1726-1732.
    [3] GAD H M, IBRAHIM I A, ABDEL-BAKY M E, et al. Experimental study of diesel fuel atomization performance of air blast atomizer[J]. Experimental Thermal and Fluid Science, 2018, 99: 211-218. doi: 10.1016/j.expthermflusci.2018.07.006
    [4] 石仲璟, 王学生, 陈琴珠. 乙氧基化反应喷嘴雾化和反应性能[J]. 华东理工大学学报(自然科学版), 2017, 43(2): 273-279.
    [5] CANTON J, AUTERI F, CARINI M. Linear global stability of two incompressible coaxial jets[J]. Journal of Fluid Mechanics, 2017, 824: 886-911. doi: 10.1017/jfm.2017.290
    [6] MATAS JP, DELON A, CARTELLIER A. Shear instability of an axisymmetric air–water coaxial jet[J]. Journal of Fluid Mechanics, 2018, 843: 575-600. doi: 10.1017/jfm.2018.167
    [7] URBÁN A, MALÝ M, JÓZSA V, et al. Effect of liquid preheating on high-velocity airblast atomization: From water to crude rapeseed oil[J]. Experimental Thermal and Fluid Science, 2019, 102: 137-151. doi: 10.1016/j.expthermflusci.2018.11.006
    [8] KHAN M A, GADGIL H, KUMAR S. Influence of liquid properties on atomization characteristics of flow-blurring injector at ultra-low flow rates[J]. Energy, 2019, 171: 1-13. doi: 10.1016/j.energy.2019.01.006
    [9] JIANG DJ, LIU HF, LI WF, et al. Modeling atomization of a round water jet by a high-speed annular air jet based on the self-similarity of droplet breakup[J]. Chemical Engineering Research and Design, 2012, 90(2): 185-192. doi: 10.1016/j.cherd.2011.07.006
    [10] BOIKO V M, NESTEROV A Y, POPLAVSKI S V. Liquid atomization in a high-speed coaxial gas jet[J]. Thermophysics and Aeromechanics, 2019, 26(3): 385-398. doi: 10.1134/S0869864319030077
    [11] CHEN L, LIU Z, SUN P, et al. Formulation of a fuel spray smd model at atmospheric pressure using design of experiments (doe)[J]. Fuel, 2015, 153: 355-360. doi: 10.1016/j.fuel.2015.03.013
    [12] SHEN J, LI X. Instability of an annular viscous liquid jet[J]. Acta Mechanica, 1996, 114: 163-183.
    [13] ALISEDA A, HOPFINGER E J, LASHERAS J C, et al. Atomization of viscous and non-newtonian liquids by a coaxial, high-speed gas jet. Experiments and droplet size modeling[J]. International Journal of Multiphase Flow, 2008, 34(2): 161-175. doi: 10.1016/j.ijmultiphaseflow.2007.09.003
    [14] JOSEPH D D, BEAVERS G S, FUNADA T. Rayleigh–taylor instability of viscoelastic drops at high weber numbers[J]. Journal of Fluid Mechanics, 2002, 453: 109-132. doi: 10.1017/S0022112001006802
    [15] VARGA C M, LASHERAS J C, HOPFINGER E J. Initial breakup of a small-diameter liquid jet by a high-speed gas stream[J]. Journal of Fluid Mechanics, 2003, 497: 405-434. doi: 10.1017/S0022112003006724
    [16] ZHAO H, LIU HF, XU JL, et al. Breakup and atomization of a round coal water slurry jet by an annular air jet[J]. Chemical Engineering Science, 2012, 78: 63-74. doi: 10.1016/j.ces.2012.05.007
    [17] URBáN A, ZAREMBA M, MALý M, et al. Droplet dynamics and size characterization of high-velocity airblast atomization[J]. International Journal of Multiphase Flow, 2017, 95: 1-11. doi: 10.1016/j.ijmultiphaseflow.2017.02.001
    [18] GAñáN-CALVO A M. Enhanced liquid atomization: From flow-focusing to flow-blurring[J]. Applied Physics Letters, 2005, 86(21): 214101. doi: 10.1063/1.1931057
    [19] SIMMONS B M, AGRAWAL A K. Spray characteristics of a flow-blurring atomizer[J]. Atomization and Sprays, 2010, 20(9): 821-835. doi: 10.1615/AtomizSpr.v20.i9.60
    [20] CUSTER J R, RIZKF N K. Influence of design concept and liquid properties on fuel injector performance[J]. Journal of Propulsion and Power, 1988, 4(4): 378-384. doi: 10.2514/3.23077
    [21] KONSTANTINOV D, MARSH R, BOWEN P, et al. Effervescent atomization for industrial energy-technology review[J]. Atomization and Sprays, 2010, 20(6): 525-552. doi: 10.1615/AtomizSpr.v20.i6.40
    [22] KOURMATZIS A, LOWE A, MASRI A R. Combined effervescent and airblast atomization of a liquid jet[J]. Experimental Thermal and Fluid Science, 2016, 75: 66-76. doi: 10.1016/j.expthermflusci.2016.02.002
    [23] WU ZW, SHI ZH, ZHAO H, et al. Effects of bubbles in the liquid jet on the air-blast atomization[J]. Fuel, 2020, 266: 117117. doi: 10.1016/j.fuel.2020.117117
    [24] LEFEBVRE A H. Airblast atomization[J]. Progress in Energy and Combustion Science, 1980, 6: 233-261. doi: 10.1016/0360-1285(80)90017-9
    [25] RIZK N K ,LEFEBVRE. A H Spray characteristics of plain-jet airblast atomizers[J]. Journal of Engineering for Gas Turbines and Power, 1983, 106(3): 634-640.
    [26] JASUJA A K. Atomization of crude and residual fuel oils[J]. J Eng Power, 1979, 101: 250-258.
    [27] JASUJA A K. Airblast atomization of alternative liquid petroleum fuels under high pressure conditions[J]. Journal of Engineering for Gas Turbines and Power, 1982, 103: 514-518.
    [28] LASHERAS J C, VILLERMAUX E, HOPFINGER E J. Break-up and atomization of a round water jet by a high-speed annular air jet[J]. Journal of Fluid Mechanics, 1998, 357: 351-379. doi: 10.1017/S0022112097008070
    [29] MURAI Y, SHIRATORI T, KUMAGAI I, et al. Effective viscosity measurement of interfacial bubble and particle layers at high volume fraction[J]. Flow Measurement and Instrumentation, 2015, 41: 121-128. doi: 10.1016/j.flowmeasinst.2014.10.006
    [30] LLEWELLIN E W, MANGA M. Bubble suspension rheology and implications for conduit flow[J]. Journal of Volcanology and Geothermal Research, 2005, 143: 205-217.
    [31] TASAKA Y, KIMURA T, MURAI Y. Estimating the effective viscosity of bubble suspensions in oscillatory shear flows by means of ultrasonic spinning rheometry[J]. Experiments in Fluids, 2014, 56(1): 1867.
    [32] MURAI Y, OIWA H. Increase of effective viscosity in bubbly liquids from transient bubble deformation[J]. Fluid Dynamics Research, 2008, 40(7-8): 565-575. doi: 10.1016/j.fluiddyn.2007.12.009
  • [1] 牛永健杨澜朱诗杰陈雪莉 . 多喷嘴对置式水煤浆气化细灰的润湿特性. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20200328001
    [2] 许叶龙刘迎圆惠虎於晔鸿 . 固液混合过程的数值模拟及实验研究. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180709001
    [3] 杨家鹏熊文莉安琦 . 滚子直径误差对双列调心轴承疲劳寿命的影响. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180411001
    [4] 熊祥军邵方明张祖渊 . 直径限定可靠性计算的冗余边的检测算法. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.
    [5] 胡志梁于新海赖焕新 . 拟塑性流体在六直叶搅拌器中的模拟与比较. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180517006
    [6] 肖桐王千红沈盈莺惠卿陈晗程振民 . 颗粒床内固液过滤的三维CFD-DEM模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190223001
    [7] 肖卫强周国俊蒋健胡安福詹望成郭杨龙 . 加热卷烟制品传热与烟气流动过程数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20191029004
    [8] 刘状查如俊邵斐凌昊 . 石油树脂钻井液用页岩抑制剂的研制. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20191209003
    [9] 周志豪沈中杰梁钦锋许建良刘海峰 . 神府煤掺石油焦气流床气化细渣的CO2气化反应特性. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20200331002
    [10] 张圆飞闫宇强宗原曹发海 . 以鲍尔环为内构件的鼓泡塔反应器的流体力学性能. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180427001
  • 加载中
图(11)表(2)
计量
  • 文章访问数:  94
  • HTML全文浏览量:  64
  • PDF下载量:  1
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-04-05
  • 网络出版日期:  2020-07-16

液体射流内部气泡对黏性流体气流式雾化的影响

    作者简介:吴兆伟(1991-),男,安徽安庆人,博士生,研究方向为多相流气液雾化。E-mail:Y20150076@mail.ecust.edu.cn
    通讯作者: 赵辉, zhaohui@ecust.edu.cn
  • 1. 华东理工大学洁净煤技术研究所,上海煤气化工程技术研究中心,上海 200237
  • 2. 上海理工大学能源与动力工程学院,上海 200093

摘要: 使用马尔文激光粒度仪研究了同轴双通道喷嘴液体射流内部气泡对黏性流体气流式雾化的影响。以甘油水溶液(昂色格数Oh为2.14×10-3~0.407)作为液相,空气作为气相,研究了不同环隙气流速度(ug为94.8~142 m/s)和液体射流速度(ul,0为0.29~1.21 m/s)下雾化液滴直径的变化规律。实验发现:当Oh较小时,低液体射流速度条件下雾化液滴直径随气液体积比的增大而减小,高液体射流速度条件下雾化液滴直径随气液体积比的增大而增大;当Oh较大时,液滴直径随气液体积比的增大而增大。建立了一个物理模型用来描述黏性流体的雾化液滴直径随气液体积比的变化规律,模型计算值与实验结果之间的偏差小于±5%。

English Abstract

  • 液体雾化广泛应用于能源、化工、航空航天以及交通运输等领域[1-4]。气流式雾化通过气液之间的相互作用,将高速气体的能量传递给低速液体从而实现液体射流的破裂和雾化。随着科技和社会的发展,节能减排已经成为社会共识,同轴双通道雾化喷嘴由于其结构简单,在较低的射流压力下能实现较好的雾化和气液混合效果,减小副产物的生成和污染物的排放,一直受到工业和科学界的青睐[3, 5]

    气流式雾化的过程和机理比较复杂,受很多因素的影响。喷嘴结构[6]、气液性质[7-8]、气液质量流量比[9]、气液速度[10]等都会对雾化的最终结果产生影响。在众多气流式雾化应用场景中,液体黏度的影响不可忽略。Chen等[11]利用响应面方法研究了气液射流压力和液体黏度对雾化液滴直径的影响,发现黏度对雾化液滴直径的影响仅次于液体射流压力。Shen和Li[12]通过线性稳定性分析了气流式雾化过程,发现液体粘度对雾化存在双重影响。Aliseda等[13]在Joseph等[14]的通用色散关系式的基础上,将Varga等[15]的瑞利泰勒模型拓展到黏性和非牛顿流体,得到了雾化液滴直径和韦伯数 (We) 以及昂色格数 (Oh) 之间的关系。

    雾化效果会显著影响后续化学反应或者燃烧的进行,因此对气流式雾化效果进行表征非常重要。常见的雾化表征参数很多,有射流断裂长度[16]、雾化角[3]、液滴直径[17]等。其中,液滴索特尔平均直径D32经常被用来描述雾化液滴的大小。雾化液滴D32的减小能够提高比表面积,从而加快反应或燃烧的进行。Calvo[18]通过调整气体通道和液体通道的相对位置,发明了一种新的雾化方式——分散流雾化。与传统的气流式雾化相比,分散流雾化能产生更小的雾化液滴[19]。自1988年[20]气泡雾化被提出以来,由于其能显著提高雾化效果,一直是雾化研究的热点[21]。Kourmatzis[22]等发现将气泡雾化与气流式雾化相结合,在气液相对速度较小时可以削弱雾化波动性同时可以提高雾化的弥散特性。Wu等[23]发现向水射流中注入气泡后,气液混合物在环隙气流作用下雾化产生的液滴随注入气体量的增大呈现先减小后增大的趋势。

    虽然前人的研究表明,向液体射流中注入气体能够提高气流式雾化的效果,但对于含气泡粘性流体的气流式雾化尚缺少系统性的研究。本文在前人工作的基础上,实验研究了液体射流内部气泡对粘性流体气流式雾化的影响,旨在为含气泡流体气流式雾化的工业应用提供参考。

    • 实验装置如图1所示,来自储罐的液体在加压罐中被加压,通过液体流量计(精度±1.5%)进入喷嘴,并与来自气体注射泵(精度±0.35%)的空气在喷嘴中心通道混合形成气液混合物。气液混合物在环隙高速气流的作用下发生雾化,形成大量小液滴。实验使用马尔文激光粒度仪(Malvern Instruments Ltd)测量雾化液滴直径从而评判雾化的效果。马尔文激光粒度仪的激光直径设置为10 mm,激光波长为632.8 nm,镜头焦距300 mm,采样周期1 s,液滴直径测量范围为0.10~900 μm。使用气体玻璃转子流量计(精度±1.5%)来测量环隙气流的流量。实验使用的雾化喷嘴如图2所示。d0=0.50 mm是充气通道的直径,来自注射泵的气体通过此通道进入液体中心通道(d1=2.96 mm)。d2=5.10 mm、d3=13.2 mm分别是喷嘴外侧环隙通道的内径和外径。d4=25.00 mm是环隙气流入口内径,环隙气流由此进入喷嘴。通过配置不同浓度的甘油水溶液得到不同粘度的雾化液体,实验中采取的甘油水溶液的物性参数见表1Oh数用来定义液体黏度对表面张力的比值,即

      图  1  实验装置示意图 1–Syringe pump; 2, 8, 9, 12–Valve; 3–Experimental atomizer; 4, 14–Gas cylinder; 5–Liquid tube rotameter; 6–Glass tube rotameter flowmeter; 7–Gauge; 10, 11–Malvern laser particle size analyzer; 13–Air blower; 15–Pressurized liquid vessel; 16–Pump; 17–Liquid storage; 18–Computer

      Figure 1.  Schematic of the experimental setup

      Liquidρl / (kg·m−3)σ / (mN·m−1)μl / (mPa·s)Oh
      Water+Glycerol (w=60%)1 16066.916.73.46×10-2
      Water+Glycerol (w =75%)1 20066.556.30.115
      Water+Glycerol (w =80%)1 21066.01040.212
      Water+Glycerol (w =85%)1 23065.62000.407

      表 1  不同液体物性参数表

      Table 1.  Properties of the different liquids

      图  2  喷嘴示意图

      Figure 2.  Aerated air-blast atomizer geometry

      其中,${\mu _{\rm{l}}}$${\rho _{\rm{l}}}$$\sigma $分别是液体黏度,密度和表面张力。

    • 首先研究了不同黏度甘油水溶液的气流式雾化过程。甘油水溶液通过喷嘴后,在环隙气流的作用下雾化,形成大量小液滴。实验中,通过保持加压罐中的压力不变获得稳定的液体流量Ql。在液体射流内部不存在气泡时,液体射流平均出口速度为:

      通过调节鼓风机阀门获得稳定的环隙气体流量Qcg。因此,环隙气流的速度可以表示为:

      液体气流式雾化的效果不仅受液体黏度的影响,而且受气液射流速度的影响。因此,研究了不同黏度的流体在不同环隙气流和液体射流速度下的雾化液滴直径D32的变化情况,如图3表2所示。从图3可以看出,当液体射流速度保持不变时,提高环隙气流速度会产生更细小的液滴。这是因为环隙气流速度的提高使得气液之间的相对速度变大,气液间剪切作用变强。当环隙气流速度保持不变时,提高液体黏度,雾化液滴尺寸逐渐变大。这种增大的趋势随着环隙气流速度的提高而逐渐变得不明显。这是由于液体黏度会抑制射流表面的不稳定性,从而更难产生细小的液滴。但随着环隙气流速度的提高,气液间相互作用逐渐增强。相对于气液间相互作用,黏度对射流不稳定波的抑制作用逐渐变得不显著。不同黏性流体之间的雾化效果区别不明显。表2是环隙气流速度ug=142 m/s时,不同液体射流速度下,液滴直径D32的变化趋势。从表中可以看出,随着液体射流速度的增大,液滴直径D32逐渐增大。这是由于液体射流速度的提高减小了气液间速度差,降低了气液间的剪切作用。从表中还可以看出,随着液体黏度的增大,液体黏度对表面不稳定波的抑制作用逐渐增强,液滴直径D32逐渐增大。

      ul / (m·s−1)Measure value/μmPredicted value/μm
      Oh =3.46×10−2Oh =0.407Oh =3.46×10−2Oh =0.407
      0.2925.5±0.8830.6±0.5625.531.5
      0.3927.9±0.5033.2±0.5226.532.7
      0.5029.8±0.3634.0±0.3327.433.9

      表 2  甘油水溶液的雾化液滴直径随液体射流速度的变化(ug=142 m/s)

      Table 2.  Dependence of the droplet size D32 on liquid apparent velocity (ug=142 m/s)

      图  3  甘油水溶液的雾化液滴直径随环隙气流速度变化(ul,0=0.29 m/s)

      Figure 3.  Dependence of the droplet size D32 on the co-flow air stream velocity(ul,0=0.29 m/s)

      Aliseda等[13]假设黏度对D32的影响是具有加法性质的,通过理论推导得到了和实验结果吻合较好的D32计算模型。Lefebvre[24]、Rizk和Lefebvre[25]、Jasuja等[26-27]通过实验结果拟合得到的D32关系式也呈现相似的形式。作者在之前的研究中已经得到了相似工况下水的气流式雾化的D32计算模型[23],因此黏性流体的D32拟合关系式可以写作:

      其中ABC是和实验条件相关的拟合系数。通过拟合可以得到A=0.0397,B=0.383,C=0.591。拟合结果如图3表2所示。图4是甘油水溶液雾化液滴直径模型计算值与实验测量值之间的比较。从表2图3图4可以看出,模型计算值与实验测量值之间的误差在±5%以内,拟合结果与实验结果吻合较好。

      图  4  D32计算值和实验测量值之间的对比

      Figure 4.  Comparison of the experimental D32 to the calculated D32

    • 通过注射泵进入喷嘴的气体在表面张力的作用下形成气泡,并与液体混合形成气液混合物。气液混合物随后在环隙气流的作用下发生雾化。为了定量描述液体射流内部气泡对气流式雾化的影响,定义了气液体积比V

      其中Qag是注入喷嘴的气体体积流量。

      实验测得了不同液体射流速度下下气泡直径db随气液体积比V的变化,如图5所示。气泡直径测量是在没有环隙气流存在下进行的,每个工况下至少进行了7次测量并取其平均值作为气泡平均直径。从图中可以看出,在液体射流速度保持不变时,随着气液体积比的增大,气泡直径逐渐增大并逐渐逼近喷嘴中心通道直径。因此,气泡直径可以表示为气液体积比的函数。

      图  5  气泡直径随气液体积比的变化

      Figure 5.  Bubble sizes at varying V

      实验使用马尔文激光粒度仪测量了不同工况下雾化液滴直径。图6ug=118 m/s,ul,0=0.29 m/s时,Oh=0.115和0.407的黏性流体含气泡雾化液滴直径随气液体积比V的变化趋势图。图中的竖实线对应的是停止通气的时刻,为t=30 s。从图6中可以看出,通气前后雾化液滴直径存在明显的不同。当Oh=0.115时,如图6(a)所示,内部气体的加入减小了雾化液滴直径。随着气液体积比的提高,雾化液滴直径逐渐减小。当Oh=0.407时,雾化液滴直径的变化趋势与Oh=0.115时并不相同。雾化液滴直径随气液体积比的增大逐渐增大。在实验过程中,计算了停止通气前30 s的雾化液滴直径的平均值作为对应气液体积比下的含气泡黏性流体的雾化液滴平均直径。图7ug=118 m/s,ul,0=0.29 m/s时,不同Oh下含气泡黏性流体雾化液滴直径随气液体积比的变化趋势图。从图中可以看出,在Oh数较小时,随着气液体积比的逐渐提高,雾化液滴直径逐渐减小。随着Oh数的增加,这种减小的趋势逐渐减弱。当Oh=0.212时,这种减小的趋势几乎变得不可见。进一步提高Oh数至0.407,气液体积比的增加反而会使得液滴直径逐渐增加。

      图  6  含气泡黏性流体雾化液滴直径随时间的变化图 (图中竖实线对应停止通气的时刻,ug=118 m/s,ul,0=0.29 m/s)

      Figure 6.  Time evolution of D32 of viscous liquid at different operation conditions. The vertical solid line is corresponding to the time stopping aerating the gas (ug=118 m/s, ul,0=0.29 m/s)

      图  7  D32随气液体积比的变化 (ug=118 m/s,ul,0=0.29 m/s)

      Figure 7.  Effects of changing V on D32 (ug=118 m/s, ul,0=0.29 m/s)

      为了表征液体射流内部气泡对气流式雾化的影响,定义雾化量纲为—液滴直径D

      其中(D32)0是不存在内部气泡时粘性流体的雾化液滴直径,(D32)b是对应工况下液体射流内部存在气泡时的雾化液滴直径。

      图8ug=118 m/s,ul=0.29 m/s时,存在内部气泡的黏性流体雾化液滴直径D随气液体积比的变化趋势图。如图8所示,液体射流内部存在气体时,量纲为—液滴直径D的变化趋势和液体Oh有关。Oh较小时,提高气液体积比会导致较小的量纲为—液滴直径D。但Oh较大时,提高气液体积比会增大液滴的量纲为—直径D

      图  8  D32随气液体积比V的变化 (ug=118 m/s,ul,0=0.29 m/s)

      Figure 8.  Dimensionless droplet size D32 at varying V (ug=118 m/s, ul,0=0.29 m/s)

      图  9  量纲为—液滴直径D随气液体积比的变化 (ul,0=0.29 m/s)

      Figure 9.  Dependence of dimensionless droplet size D on V (ul,0=0.29 m/s)

      图  10  量纲为一液滴直径D随气液体积比的变化 (ug=142 m/s)

      Figure 10.  Dimensionless droplet size D at varying V (ug=142 m/s)

      为了探究环隙气流速度对含气泡黏性流体雾化的影响,实验研究了ug=94.8~142 m/s,ul,0=0.29 m/s时含气泡黏性流体的雾化过程。图9是不同黏度下,量纲为—液滴直径D随气液体积比V的变化趋势图。从图中可以看出,当气液体积比保持不变时,提高环隙气流速度会导致量纲为—液滴直径D的减小。当Oh=3.46×10-2时,量纲为—液滴直径D随着气液体积比的增加而逐渐减小。当Oh=0.115时,量纲为—液滴直径D随着气液体积比的增加而逐渐减小,但减小的趋势随着气液体积比的增加和环隙流速度的减小而逐渐减小。当Oh=0.212时,随着气液体积比的增加,量纲为—液滴直径D呈现先减小而后增加的趋势。环隙气流速度越大,这种增加的趋势就愈加的不明显。当Oh=0.407时,提高气液体积比会导致量纲为—液滴直径D的逐渐增大。

      为了探究不同液体射流速度下,液体射流内部气泡对粘性流体气流式雾化的影响,研究了ug=142 m/s,ul=0.29~1.21 m/s时不同黏性流体的雾化过程。图10是不同液体射流速度下,雾化液滴直径D随着气液体积比的变化。从图中可以看出,当气液体积比保持不变时,随着液体射流速度的增大,量纲为—液滴直径逐渐增大。当Oh=3.46×10-2时,低液体射流速度下的量纲为—液滴直径D随着气液体积比的增大而逐渐减小;高液体射流速度下量纲为—液滴直径D随着气液体积比的增大而逐渐增大。实验也研究了ul=1.21 m/s时水的含气泡雾化(图10(a))。从图中可以看出,水的雾化液滴直径变化趋势和黏性流体并不相同,这可能是流体物性不同导致的。当Oh=0.407时,随着气液体积比的增大,量纲为—液滴直径逐渐增大。

      液体黏度会阻碍射流表面不稳定波的发展,从而产生更大的液滴。因此在其他条件不变时,提高液体黏度会产生更大的液滴。对本文中特定的工况,液体的体积流量保持不变,因此当气体被注入液体射流内部时,液体表观速度会随着注入气体量的增加而增加。当环隙气流速度保持不变时,液体射流内部气体的增加会导致环隙气流和液体射流之间速度差的减小,从而削弱气液间剪切作用,形成更大的雾化液滴。降低环隙气流速度或者提高液体射流速度都会减小环隙气流和液体射流之间的速度差,从而产生更大的雾化液滴。前人研究发现,减小射流液体和环隙气流之间的密度比$\frac{{{\rho _{\rm{l}}}}}{{{\rho _{\rm{g}}}}}$会缩短液体射流表面瑞利-泰勒不稳定波的波长,从而产生更小的雾化液滴直径[1, 15, 28]。而气体的注入会导致液体射流密度的减小,从而减小密度比$\frac{{{\rho _{\rm{l}}}}}{{{\rho _{\rm{g}}}}}$,产生更小的液滴。含气泡流体最终的雾化液滴直径依赖于这两种作用的相对大小。当Oh较小时,气液密度比的影响较为显著,雾化液滴直径随着气液体积比的增大而逐渐减小。当Oh较大时,气液相对速度的影响较显著,提高气液体积比会产生更大的液滴。

    • 本研究通过保持加压罐中的压力不变,获得了稳定的液体体积流量。被注射泵注入喷嘴的气体会与液体射流混合形成气液混合物。为了简化分析,假设此时的气液混合物是均匀的。因此,气液混合物的速度可以表示为:

      其中k值是考虑到实际气液混合物对均匀气液混合物的偏差而设置的修正系数,通过作者之前的研究,可以得到k=24.5[23]。此时韦伯数Weg,b可以表示为:

      气液混合物的密度${\rho _{\rm{m}}}$可以表示为:

      前人研究发现,溶体中的颗粒或者气泡会影响溶液的黏度[29-32]。Llewellin和Mang[30]发现在液体中加入气泡可能会减小溶液的黏度,也可能会提高液体的黏度,最终结果取决于量纲为一参数CaCd

      其中db$\overset{\centerdot }{\mathop{\gamma }}\,$$\overset{\centerdot \centerdot }{\mathop{\gamma }}\,$ 分别是气泡直径,剪切应变率,剪切应变率的变化率。通过计算可以得到在本实验条件下Ca<1,Cd<1。因此,气液混合物的粘度可以表示为:

      其中,${\mu _{\rm{l,b}}}$是气液混合物的黏度。k1为修正系数,且${k_1} \geqslant 1$

      因此,含气泡黏性流体的Ohb可以表示为:

      因此,结合方程式(4)~式(6), 式(10)、(11)和(15)可以得到量纲无一直径D的表达式为:

      通过拟合实验数据可以得到,k1=1.42。拟合曲线见图8~图10中实线。实验测量值与模型计算值之间的比较见图11。从图8-11可以看出,拟合曲线和实验值吻合较好,量纲为一液滴直径模型计算值与实验测量值之间的偏差在±5%以内。

      图  11  D模型计算值与实验测量值的对比

      Figure 11.  Comparison of the measured D to the calculated D

    • (1) 液体射流内部气泡对黏性流体的气流式雾化存在多方面的影响。首先,液体射流内部气泡的存在会提高液体射流的表观速度,从而降低环隙气流和液体射流之间的相对速度,因此有利于产生更大的雾化液滴;其次,液体射流内部气泡会降低液体密度,缩短射流表面的瑞利-泰勒不稳定波,从而有利于产生更小的液滴。第三,气泡的加入会提高气液混合物的黏度,从而有利于产生更大的液滴。含气泡黏性流体最终的雾化液滴直径取决于这三者的竞争作用。

      (2) 当Oh较小时,在液体射流速度较小的情况下D随气液体积比的增大而减小;在液体射流速度较大的情况下D随气液体积比的增大而增大。当Oh较大时,D随气液体积比的增大而增大。

      (3) 本文建立了一个理论模型用来描述粘性流体雾化液滴直径随气液体积比的变化规律,雾化液滴直径模型计算值与实验测量值之间的偏差在±5%以内。

(11)  表(2) 参考文献 (32) 相关文章 (10)

目录

    /

    返回文章