C2H4/N2O预混气体的爆轰性能与火焰淬熄特性

李玉艳 蒋榕培 李智鹏 徐森 潘峰 解立峰

引用本文:
Citation:

C2H4/N2O预混气体的爆轰性能与火焰淬熄特性

    作者简介: 李玉艳(1985-),女,博士研究生,主要从事推进剂安全性能研究. E-mail:lyy.piao@163.com;
    通讯作者: 徐森, pfhycc@163.com
  • 中图分类号: O643.2; TJ55

Detonation and Quenching Characteristics of Premixed C2H4/N2O

    Corresponding author: XU Sen, pfhycc@163.com ;
  • CLC number: O643.2; TJ55

  • 摘要: 采用自制的燃爆实验装置对C2H4/N2O预混气体的爆轰性能与火焰淬熄特性进行了实验研究。结果表明:在大直径有机玻璃管中预混气体均经历了爆燃转爆轰过程,点火初期火焰速度及加速度在内径为5、10和15 mm的管道中依次减小;预混气体中加入CO2(2.4%,质量分数)后,火焰加速进程明显延缓,点火初期处于稳定燃烧阶段;预混气体的稳定爆速为2 207 m/s,爆压为3.92 MPa,与理论值一致;常压下预混火焰在小直径不锈钢管中的临界淬熄管径为0.5~0.7 mm,预混气体火焰传播速度越大,管径越大,淬熄越困难。依据淬熄管径、湍流火焰速度和淬熄管道长度的关系,可计算防回火管道的有效长度,从而为防回火装置设计提供参考。
  • 图 1  有机玻璃管C4示意图(单位:mm)

    Figure 1.  Schematic of the PMMA channel C4 (Unit:mm)

    图 2  组合管道示意图(单位:mm)

    Figure 2.  Schematic of the combination channel (Unit:mm)

    图 3  不同管径管道中火焰阵面传播速度随时间变化曲线

    Figure 3.  Flame speed as a function of time in channels with various diameters

    图 4  不同管径管道中火焰加速度随时间的变化曲线

    Figure 4.  Flame acceleration as a function of time in channels with various diameters

    图 5  C1管中C2H4/N2O和C2H4/N2O/CO2的火焰传播速度随时间变化曲线

    Figure 5.  Flame speed of C2H4/N2O and C2H4/N2O/CO2 as a function of time in the channel C1

    图 6  管道C4中火焰前端速度随时间变化曲线

    Figure 6.  Flame speed as a function of time in the channel C4

    图 7  管道C4中8个压力传感器记录的压力曲线

    Figure 7.  Pressure profiles versus time obtained by eight pressure gauges in the channel C4

    图 8  压力峰值变化曲线

    Figure 8.  Plots of maximum overpressure obtained in the channel C4

    图 9  火焰传播速度、冲击波速度和C-J速度曲线

    Figure 9.  Flame speed, shock wave velocity and C-J velocity curves

    图 10  组合管道中火焰传播实测图像

    Figure 10.  Image of the flame propagation in combination channels

    图 11  组合管道中预混火焰淬熄情况

    Figure 11.  Diagram of the flame behaviors in combination channels

    图 12  加速管长度不同时火焰传播速度随时间变化曲线

    Figure 12.  Flame acceleration process as a function of time in accelerating channels with different lengths

    表 1  有机玻璃管道尺寸

    Table 1.  Geometrical characteristics of PMMA channels

    Channell/mmd/mml/dV/mL
    C1, smooth1 400 5280 27.5
    C2, smooth1 40010140109.9
    C3, smooth1 40015 93247.3
    C4, rough2 00015133112.5
     Note: l is the length of the PMMA channel; d is the inner diameter of the PMMA channel; V is the volume of the channel.
    下载: 导出CSV
  • [1] 朱成财, 韩伟, 于忻立, 等. 氧化亚氮基单元复合推进剂技术研究述评 [J]. 火箭推进, 2016, 42(2): 79–85. doi: 10.3969/j.issn.1672-9374.2016.02.015
    ZHU C C, HAN W, YU X L, et al. Review of nitrous-oxide-based composite monopropellants technology [J]. Journal of Rocket Propulsion, 2016, 42(2): 79–85. doi: 10.3969/j.issn.1672-9374.2016.02.015
    [2] 宋长青, 徐万武, 张家奇, 等. 氧化亚氮推进技术研究进展 [J]. 火箭推进, 2014, 40(2): 7–15. doi: 10.3969/j.issn.1672-9374.2014.02.002
    SONG C Q, XU W W, ZHANG J Q, et al. Research progress of nitrous oxide propulsion technology [J]. Journal of Rocket Propulsion, 2014, 40(2): 7–15. doi: 10.3969/j.issn.1672-9374.2014.02.002
    [3] MUNGAS G, VOZOFF M, RISHIKOF B. NOFBXTM: a new non-toxic, “green” propulsion technology with high performance and low cost [C]//63rd International Astronautical Congress. Naples, Italy, 2012.
    [4] GOHARDANI A S, STANOJEV J, DEMAIRE A, et al. Green space propulsion: opportunities and prospects [J]. Progress in Aerospace Sciences, 2014, 71: 128–149. doi: 10.1016/j.paerosci.2014.08.001
    [5] ROY G D, FROLOV S M, BORISOV A A. Pulse detonation propulsion: challenges, current status, and future perspective [J]. Progress in Energy and Combustion Science, 2004, 30(6): 545–672. doi: 10.1016/j.pecs.2004.05.001
    [6] WERLINGY L, LAUCK F, FREUDENMANN D, et al. Experimental investigation of the ignition, flame propagation and flashback behavior of a premixed green propellant consisting of N2O and C2H4 [J]. Journal of Energy and Power Engineering, 2017, 11: 735–752.
    [7] MUSCAT V I, DENTON M B, SUDDENDORF R F. A Flashback-resistant burner for use with the nitrous oxide-acetylene flame [J]. Spectroscopy Letters, 1973, 6(9): 563–567. doi: 10.1080/00387017308060845
    [8] GIBBON D, BAKER A, NICOLINI D, et al. The design, development and in-flight performance of a low power resistojet thruster [C]//39th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit. Huntsville, Alabama, 2003.
    [9] VENKATESH P B, D’ENTREMONT J, MEYER S E, et al. High-pressure combustion and deflagration-to-detonation transition in ethylene/nitrous oxide mixtures [C]//8th U.S. National Combustion Meeting. Park City, Utah, 2013: 158–165.
    [10] VENKATESH P B, GRAZIANO T J, BANE S P M, et al. Deflagration-to-detonation transition in nitrous oxide-ethylene mixtures and its application to pulsed propulsion systems [C]//55th AIAA Aerospace Sciences Meeting. Grapevine, TX, 2017: 0372.
    [11] WERLING L K, HOCHHEIMER B, BARAL A L, et al. Experimental and numerical analysis of the heat flux occurring in a nitrous oxide/ethene green propellant combustion demonstrator [J]. Journal of the American College of Surgeons, 2013, 186(5): 562–569.
    [12] HOCHHEIMER B, PERAKIS N, WERLING L, et al. Test facilities to assess properties of a nitrous oxide/ethene premixed bipropellant for satellite propulsion system [C]//5th CEAS Air & Space Conference. Delft, The Netherlands, 2014.
    [13] ZHANG B, LIU H, WANG C. Detonation velocity behavior and scaling analysis for ethylene-nitrous oxide mixture [J]. Applied Thermal Engineering, 2017, 127: 671–678. doi: 10.1016/j.applthermaleng.2017.08.016
    [14] MOVILEANU C, RAZUS D, MITU M, et al. Explosion of C2H4-N2O-N2 in elongated closed vessels [C]//7th European Combustion Meeting. Budapest, Hungary, 2015.
    [15] POWELL O A, PAPAS P, DREYER C. Laminar burning velocities for hydrogen-, methane-, acetylene-, and propane-nitrous oxide flames [J]. Combustion Science & Technology, 2009, 181(7): 917–936.
    [16] 李智鹏, 孙海云, 蒋榕培, 等. 乙烯-氧化亚氮层流预混燃烧过程研究 [J]. 火箭推进, 2018, 44(5): 37–42. doi: 10.3969/j.issn.1672-9374.2018.05.006
    LI Z P, SUN H Y, JIANG R P, et al. Study on premixed laminar combustion process of ethylene/nitrous oxide mixture [J]. Journal of Rocket Propulsion, 2018, 44(5): 37–42. doi: 10.3969/j.issn.1672-9374.2018.05.006
    [17] NEWMAN-LEHMAN T, GRANA R, SESHADRI K, et al. The structure and extinction of nonpremixed methane/nitrous oxide and ethane/nitrous oxide flames [J]. Proceedings of the Combustion Institute, 2013, 34(2): 2147–2153. doi: 10.1016/j.proci.2012.05.102
    [18] 程关兵, 李俊仙, 李书明, 等. 氢气/丙烷/空气预混气体爆轰性能的实验研究 [J]. 爆炸与冲击, 2015, 35(2): 249–252. doi: 10.11883/1001-1455(2015)02-0249-06
    CHENG G B, LI J X, LI S M, et al. An experimental study on detonation characteristics of binary fuels hydrogen/propane-air mixtures [J]. Explosion and Shock Waves, 2015, 35(2): 249–252. doi: 10.11883/1001-1455(2015)02-0249-06
    [19] 张博, 白春华. H2-O2/Air直接起爆形成爆轰临界能量的预测模型 [J]. 高压物理学报, 2013, 27(5): 719–724. doi: 10.11858/gywlxb.2013.05.010
    ZHANG B, BAI C H. Theoretical prediction model of critical energy for direct detonation initiation in H2-O2/air mixtures [J]. Chinese Journal of High Pressure Physics, 2013, 27(5): 719–724. doi: 10.11858/gywlxb.2013.05.010
    [20] 王鲁庆, 马宏昊, 王波, 等. 氢气/甲烷-空气爆轰波在含环形障碍物圆管内传播的试验研究 [J]. 高压物理学报, 2018, 32(3): 035203. doi: 10.11858/gywlxb.20170687
    WANG L Q, MA H H, WANG B, et al. Detonation propagation in hydrogen/methane-air mixtures in a round tube filled with orifice plates [J]. Chinese Journal of High Pressure Physics, 2018, 32(3): 035203. doi: 10.11858/gywlxb.20170687
    [21] WU M H, BURKE M P, SON S F, et al. Flame acceleration and the transition to detonation of stoichiometric ethylene/oxygen in microscale tubes [J]. Proceedings of the Combustion Institute, 2007, 31(2): 2429–2436. doi: 10.1016/j.proci.2006.08.098
    [22] WANG C, HUANG F L, ADDAI E K, et al. Effect of concentration and obstacles on flame velocity and overpressure of methane-air mixture [J]. Journal of Loss Prevention in the Process Industries, 2016, 43: 302–310. doi: 10.1016/j.jlp.2016.05.021
    [23] ALIOU S, ASHWIN C, ABDELLAH H. Mean structure of one-dimensional unstable detonations with friction [J]. Journal of Fluid Mechanics, 2014, 743(3): 503–533.
    [24] 路长, 李毅, 潘荣锟. 管道截面对氢气/空气预混爆炸影响的实验研究 [J]. 火灾科学, 2015, 24(2): 68–74. doi: 10.3969/j.issn.1004-5309.2015.02.02
    LU C, LI Y, PAN R K. Experimental study of the duct cross section effects on the hydrogen/air premixed explosion [J]. Fire Safety Science, 2015, 24(2): 68–74. doi: 10.3969/j.issn.1004-5309.2015.02.02
    [25] 王成, 韩文虎, 宁建国. 边界层和障碍物对湍流火焰加速机理的研究 [C]//第十五届全国激波与激波管学术会议. 杭州, 2012.
    [26] LIU F, GUO H, SMALLWOOD G J. The chemical effect of CO2 replacement of N2 in air on the burning velocity of CH4 and H2 premixed flames [J]. Combustion and Flame, 2003, 133(4): 495–497. doi: 10.1016/S0010-2180(03)00019-1
    [27] PARK J, LEE K, LEE E. Effects of CO2 addition on flame structure in counter flow diffusion flame of H2/CO2/N2 fuel [J]. International Journal of Hydrogen Energy, 2001, 25(6): 469–485.
    [28] PARK J, HWANG D, CHOI J, et al. Chemical effects of CO2 addition to oxidizer and fuel streams on flame structure in H2-O2 counter flow diffusion flames [J]. International Journal of Energy Research, 2003, 27(13): 1205–1220. doi: 10.1002/er.946
    [29] POWELL O, PAPAS P. Flame structure measurements of nitric oxide in hydrocarbon-nitrous-oxide flames [J]. Journal of Propulsion & Power, 2015, 28(5): 1052–1059.
    [30] 李猛, 王宏, 陈雪莉. 复杂化学平衡应用计算程序 [J]. 兵器装备工程学报, 2010, 31(9): 132–134. doi: 10.3969/j.issn.1006-0707.2010.09.043
    [31] BABKIN V S. Filtrational combustion of gases. Present state of affairs and prospects [J]. Pure & Applied Chemistry, 1993, 65(2): 335–344.
    [32] KELLENBERGER M, CICCARELLI G. Advancements on the propagation mechanism of a detonation wave in an obstructed channel [J]. Combustion and Flame, 2018, 191: 195–209. doi: 10.1016/j.combustflame.2017.12.023
    [33] LIBERMAN M A, IVANOV M F, KIVERIN A D, et al. Deflagration-to-detonation transition in highly reactive combustible mixtures [J]. Acta Astronautica, 2010, 67(7/8): 688–701.
    [34] CHAKRAVARTHII D M K, DEVARAJAN M, SUBRAMANI S. Experimental and numerical investigation of pressure drop and heat transfer coeffcient in converging-diverging microchannel heat sink [J]. Heat and Mass Transfer, 2017, 53(7): 2265–2277. doi: 10.1007/s00231-017-1978-7
    [35] FAN A, WAN J, LIU Y, et al. Effect of bluff body shape on the blow-off limit of hydrogen/air flame in a planar micro-combustor [J]. Applied Thermal Engineering, 2014, 62(1): 13–19. doi: 10.1016/j.applthermaleng.2013.09.010
    [36] YANG S, SHY S. Global quenching of premixed CH4/air flames: effects of turbulent straining, equivalence ratio, and radiative heat loss [J]. Proceedings of the Combustion Institute, 2002, 29(2): 1841–1847. doi: 10.1016/S1540-7489(02)80223-1
    [37] YANG W, DENG C, ZHOU J, et al. Experimental and numerical investigations of hydrogen-air premixed combustion in a converging-diverging micro tube [J]. International Journal of Hydrogen Energy, 2014, 39(7): 3469–3476. doi: 10.1016/j.ijhydene.2013.12.102
  • [1] 毛浩清黄炜超李斌解立峰 . 强点火条件下RP-3航空煤油燃爆特性实验研究. 高压物理学报, 2018, 32(2): 025201-1-025201-8. doi: 10.11858/gywlxb.20170583
    [2] 桑圣军郭浩哲李斌王永旭王志平解立峰 . 爆速对纳米铝粉爆炸烧结性能的影响. 高压物理学报, 2018, 32(2): 023401-1-023401-7. doi: 10.11858/gywlxb.20170581
    [3] 韩志伟解立峰邓吉平倪欧琪解一超陈际洋 . 乳化炸药爆轰合成纳米氧化铈及表征. 高压物理学报, 2014, 28(3): 319-323. doi: 10.11858/gywlxb.2014.03.009
    [4] 赵国民张若棋陈刚王诚洪 . 铅壳柔爆索爆炸特性实验研究. 高压物理学报, 2001, 15(2): 91-96 . doi: 10.11858/gywlxb.2001.02.003
    [5] 李成兵吴国栋经福谦 . CH4-O2-N2预混气体爆炸二维数值模拟. 高压物理学报, 2009, 23(5): 367-376 . doi: 10.11858/gywlxb.2009.05.008
    [6] 晏成立于川李良忠仝延锦 . 爆炸磁通量压缩发生器金属管爆炸试验研究. 高压物理学报, 1999, 13(1): 76-80 . doi: 10.11858/gywlxb.1999.01.014
    [7] 张博白春华 . C2H2-O2-Ar 混合气体爆轰特征参数研究. 高压物理学报, 2013, 27(2): 287-291. doi: 10.11858/gywlxb.2013.02.017
    [8] 宁建国李健王成赵慧 . 基于基元反应模型的H2-O2-N2爆轰数值模拟. 高压物理学报, 2011, 25(5): 395-400 . doi: 10.11858/gywlxb.2011.05.003
    [9] 张博白春华 . H2-O2/Air直接起爆形成爆轰临界能量的预测模型. 高压物理学报, 2013, 27(5): 719-724. doi: 10.11858/gywlxb.2013.05.010
    [10] 王小红李晓杰张越举曲艳东孙贵磊谢兴华 . 爆轰法制备纳米MnFe2O4粉体的实验研究. 高压物理学报, 2007, 21(2): 173-177 . doi: 10.11858/gywlxb.2007.02.009
    [11] 张博Lee John H S白春华 . 高浓度氩气稀释对C2H2-2.5O2气体直接起爆临界能量影响的实验研究. 高压物理学报, 2012, 26(1): 55-62. doi: 10.11858/gywlxb.2012.01.008
    [12] 曲艳东李晓杰赵铮欧阳欣 . 爆轰和燃烧法合成SrAl2O4:Eu2+,Dy3+纳米发光粉. 高压物理学报, 2008, 22(2): 175-180 . doi: 10.11858/gywlxb.2008.02.011
    [13] 李大红 . 凝聚炸药中超压爆轰的实验研究. 高压物理学报, 1987, 1(1): 81-87 . doi: 10.11858/gywlxb.1987.01.011
    [14] 曲艳东闫鸿浩李晓杰刘华新 . 裸爆和特制半球形结构内爆超压对比实验研究. 高压物理学报, 2012, 26(1): 95-101. doi: 10.11858/gywlxb.2012.01.014
    [15] 李科斌李晓杰王小红闫鸿浩曹景祥 . 基于连续压导探针的水箱法测量爆压. 高压物理学报, 2019, 33(2): 023201-1-023201-7. doi: 10.11858/gywlxb.20180652
    [16] 韩杰邹广田崔田刘冰冰潘耀武高春晓杨海滨 . Sr2Bi4Ti5O18的压致结构相变. 高压物理学报, 1998, 12(3): 233-235 . doi: 10.11858/gywlxb.1998.03.012
    [17] 曾代朋陈军谭多望 . 超压爆轰产物冲击绝热线的实验研究. 高压物理学报, 2010, 24(1): 76-80 . doi: 10.11858/gywlxb.2010.01.014
    [18] 李银成 . 爆轰产物的冷比内能与冷压的函数表达式. 高压物理学报, 2005, 19(1): 71-79 . doi: 10.11858/gywlxb.2005.01.013
    [19] 陈军曾代朋谭多望王荣波 . 超压爆轰产物声速以及热力学CJ点的实验研究. 高压物理学报, 2010, 24(3): 201-205 . doi: 10.11858/gywlxb.2010.03.007
    [20] 马秋菊邵俊程王众山刘家平 . 氢气比例和点火能量对CH4-H2混合气体爆炸强度影响的实验研究. 高压物理学报, 2020, 34(1): 015201-1-015201-6. doi: 10.11858/gywlxb.20190803
  • 加载中
图(12)表(1)
计量
  • 文章访问数:  287
  • 阅读全文浏览量:  229
  • PDF下载量:  2
出版历程
  • 收稿日期:  2019-12-17
  • 录用日期:  2020-01-05
  • 刊出日期:  2020-08-05

C2H4/N2O预混气体的爆轰性能与火焰淬熄特性

    作者简介:李玉艳(1985-),女,博士研究生,主要从事推进剂安全性能研究. E-mail:lyy.piao@163.com
    通讯作者: 徐森, pfhycc@163.com
  • 1. 南京理工大学化工学院,江苏 南京 210094
  • 2. 国家民用爆破器材质量监督检验中心,江苏 南京 210094
  • 3. 北京航天实验技术研究所,北京 100074

摘要: 采用自制的燃爆实验装置对C2H4/N2O预混气体的爆轰性能与火焰淬熄特性进行了实验研究。结果表明:在大直径有机玻璃管中预混气体均经历了爆燃转爆轰过程,点火初期火焰速度及加速度在内径为5、10和15 mm的管道中依次减小;预混气体中加入CO2(2.4%,质量分数)后,火焰加速进程明显延缓,点火初期处于稳定燃烧阶段;预混气体的稳定爆速为2 207 m/s,爆压为3.92 MPa,与理论值一致;常压下预混火焰在小直径不锈钢管中的临界淬熄管径为0.5~0.7 mm,预混气体火焰传播速度越大,管径越大,淬熄越困难。依据淬熄管径、湍流火焰速度和淬熄管道长度的关系,可计算防回火管道的有效长度,从而为防回火装置设计提供参考。

English Abstract

  • C2H4/N2O作为一种新型绿色单组元推进剂,具有高比冲、低冰点、易贮存、推力系统简单等优良性能,近年来受到广泛关注[1-3]。虽然其理化性质稳定,可于常温下带压贮存,但在某些外界条件激发下仍会被点燃,瞬间释放大量热量和气体,发生爆炸[4-5];然而,当它被点燃后,在火焰传播过程中会沿管道发生回火,引发事故[6-8]

    目前关于C2H4/N2O的研究已有文献报道。Venkatesh等[9-10]研究了高压条件下C2H4/N2O在不锈钢合金管中的燃烧特性,并测量了C2H4/N2O在高压及长径比为68的燃烧管中的爆轰速度、爆轰压力和爆燃转爆轰成长距离。德国航空航天中心[11-12]对C2H4/N2O单元复合推进剂的点火、燃烧及传热特性进行了实验研究,设计了适用于该推进系统的发动机,同时开展了发动机防回火研究。Zhang等[13]建立了不同初始压力和管道长径比条件下C2H4/N2O预混气体的爆轰行为,得到了爆速与初始压力、管径及C-J爆轰速度等参数的关系。Movileanu等[14]在初始压力为50~150 kPa,60% N2稀释条件下,针对C2H4/N2O在圆柱形管道中的爆炸压力和压力上升速率进行了实验研究,得到了最大爆炸压力与初始压力的关系。Powell等[15]对不同压力及当量比时C2H4/N2O的层流燃烧速度进行了试验及模拟研究。李智鹏等[16]采用预混模块得到了不同压力和氧/燃比时C2H4/N2O的层流火焰速度、火焰温度和燃烧质量流率变化规律。Newman-Lehman等[17]对CH4/N2O和C2H6/N2O的火焰结构和临界淬熄条件进行了实验研究,结果表明在变径管道中N2O比O2更容易使火焰熄灭。

    本工作拟在常温常压下对C2H4/N2O预混气体的爆燃规律及火焰淬熄特性进行实验研究,利用高速摄影、压力测试系统,采用不同尺寸的有机玻璃管和不锈钢管,分析管径、惰性气体对火焰传播特性的影响;通过实验获得预混气体的爆速、爆压、淬熄管径等参数,探讨不同湍流火焰传播速度和管径条件下预混火焰在小尺寸管道中的传播规律,为提高推进剂安全性能及发动机防回火设计提供一定的理论基础。

    • 实验装置主要由燃烧管道、高速摄像、配气系统、点火系统、压力测试系统及同步控制系统组成。燃烧管道分为两种:单一的有机玻璃燃烧管道和组合管道,其中有机玻璃燃烧管用于研究预混气体的爆燃特性,组合管道用于研究预混火焰的淬熄特性。

      有机玻璃管的尺寸如表1所示。燃烧管C1、C2和C3的长度均为1 400 mm,内径分别为5、10和15 mm,壁厚均为5 mm。燃烧管C4的长度为2 000 mm,内径15 mm,点火端内置长300 mm的螺旋环,螺旋环的螺距为5 mm,外径为15 mm;距离点火点100、400、600、800、1 000、1 200、1 400和1 600 mm处装有压力传感器,分别标记为P1、P2、P3、P4、P5、P6、P7和P8,传感器截面刚好与玻璃管内壁齐平;压力传感器型号为PCB-113B24,测量的最大误差为2.3%;根据传感器间隔距离和压力波经过传感器的时刻,可计算冲击波速度[18-19]。燃烧管C4示意图如图1所示。组合管道由不锈钢细管和有机玻璃管组成,其中:不锈钢细管长度为1 000 mm,位于中间位置,两端通过转换接头与有机玻璃管连接;前端玻璃管用于预混火焰加速,称为加速管;后端玻璃管用于观察火焰是否通过中间不锈钢细管,称为验证管。组合管道示意图如图2所示。

      Channell/mmd/mml/dV/mL
      C1, smooth1 400 5280 27.5
      C2, smooth1 40010140109.9
      C3, smooth1 40015 93247.3
      C4, rough2 00015133112.5
       Note: l is the length of the PMMA channel; d is the inner diameter of the PMMA channel; V is the volume of the channel.

      表 1  有机玻璃管道尺寸

      Table 1.  Geometrical characteristics of PMMA channels

      图  1  有机玻璃管C4示意图(单位:mm)

      Figure 1.  Schematic of the PMMA channel C4 (Unit:mm)

      图  2  组合管道示意图(单位:mm)

      Figure 2.  Schematic of the combination channel (Unit:mm)

      配气时,通过高精度质量流量计控制气体质量。采用堀场汇博隆精密仪器有限公司的S48 300HMT流量计,最大量程有500 mL/min和200 mL/min两种规格,误差为 ± 0.5%。配气前,先将气袋(10 L)抽至真空,再依次充入C2H4和N2O,充气完毕后关闭气袋阀门,静置24 h,以便气体混合均匀[20]。采用电阻丝点火,点火位置位于燃烧管一端,电阻丝直径0.6 mm,长100 mm,点火时通电电压为20 V,功率为20 W。采用日本PHOTRON公司的“FASTCAM”系列高速摄像机拍摄,拍摄速度为104帧每秒。

      N2O、C2H4、CO2的纯度均为99.99%,N2O/C2H4的质量比为9∶1,接近化学当量比,氧化剂略过量。

    • 充气前将真空泵、气袋和燃烧管道通过三通管道连接起来;在室温下采用真空泵先将燃烧管道抽至真空,接着用预先配好的气体(气袋中的气体)清洗管路3次,然后将混合气体充入燃烧管道中至压力为0.1 MPa,关闭进气阀门;点火前通过同步控制系统将管道末端的电磁阀打开,电磁阀打开1 s后接通点火电源,开始点火,同时启动高速摄像和压力测试系统,记录数据。每种实验条件下测试3次。

    • 图3为火焰阵面传播速度随时间的变化曲线,处理高速摄影图像时以出现可见火焰作为零时刻。由图3可知:C1(5 mm)、C2(10 mm)和C3(15 mm)3种管道中预混火焰从点火端传播至末端的时间分别为1.2、1.4和1.6 ms。点火初期(0~0.8 ms阶段),5 mm管中火焰传播速度明显大于10和15 mm管。5 mm管中预混气体被引燃之后,火焰传播速度开始迅速增大,0.2 ms时速度为250 m/s,0.4 ms时已增至1 044 m/s,0.8 ms时速度进一步增至1 661 m/s;10 mm管中,0.2~0.4 ms时间段内火焰传播速度缓慢增大,约为250 m/s,0.4~0.6 ms时间段内,火焰加速变快,0.6 ms时速度增至约600 m/s,0.8 ms时火焰速度为1 541 m/s;15 mm管中,0.2~0.8 ms之间火焰传播速度增大缓慢,速度由238 m/s增至457 m/s,此阶段预混气体处于相对平稳的燃烧阶段。因此当管径大于临界爆轰直径时,管径尺寸越大,点火初期同一时刻对应的火焰速度越小,与文献报道的结果[21]相一致。同时,从图3可以看出,5、10和15 mm 3种管径的管道中预混气体火焰最大速度分别为1 871、2 021和2 341 m/s,对应的时刻分别为1.0、1.0和1.2 ms。随着管径的增大,预混火焰传播最大速度呈现增大趋势,这主要是因为管径越小,热损失越严重,爆轰速度越小。

      图  3  不同管径管道中火焰阵面传播速度随时间变化曲线

      Figure 3.  Flame speed as a function of time in channels with various diameters

      图4为火焰加速度随时间的变化曲线。由图4可知:点火后0~0.4 ms时间段内,火焰加速度随着管径的增大依次减小。预混气体被点燃后,5 mm管中火焰加速度随时间几乎呈直线下降,0.2 ms时加速度值最大,说明5 mm管中的预混气体被点燃后很快进入爆燃转爆轰阶段;10和15 mm管中火焰加速度呈现先增大后减小的趋势,点火初期(0~0.4 ms),10和15 mm管中气体处于相对平稳燃烧阶段,火焰加速度较小,但后期随着燃烧加剧,进入爆燃转爆轰阶段后,火焰加速度明显增大,超过5 mm管中的加速度值,当火焰速度接近该管径的爆轰极限值时,加速速率开始下降。

      图  4  不同管径管道中火焰加速度随时间的变化曲线

      Figure 4.  Flame acceleration as a function of time in channels with various diameters

      综上所述,点火初期火焰速度及加速度均随管径的增大而减小。主要是由于:火焰与管壁之间的黏性效应促使火焰在管道中拉伸,形成经典的“Tulip”火焰;初始阶段火焰速度呈指数增长,然后压缩波对未反应气体进行预热,并与边界层相互作用,使火焰进一步加速;当通道尺度变小时,黏性效应加剧,边界层与内部火焰更容易产生速度梯度,火焰表面更容易产生褶皱,从而进一步促进“Tulip”火焰的形成[22-23]。研究表明:起爆首先发生在阻力较大的边界层失稳区,湍流效应及漩涡运动诱导形成“Tulip”火焰的同时加快了燃烧反应速率,尺寸越小,湍流加速形成的“Tulip”火焰越早,管内温度也越高,火焰加速速率越快[24-25]。因此,点火初期同一时刻对应的火焰速度和加速度均随着管径的增大而减小。

    • 图5为C1(5 mm)管中CO2稀释前后火焰传播速度随时间的变化曲线。结果表明:初始阶段两者差异较大,C2H4/N2O体系点火后很快进入爆燃转爆轰阶段,0.4 ms时火焰已加速至1 044 m/s;对于C2H4/N2O/CO2体系,0~1.2 ms阶段火焰处于缓慢加速阶段,1.2 ms时火焰传播速度为433 m/s。因此CO2明显降低了初始阶段火焰的加速速率,同时也降低了火焰在管道中的最大传播速度,两者的最大传播速度分别为1 871 m/s和1 781 m/s。研究表明:CO2作为稀释剂时,热力学效应起主导作用,即CO2通过改变混合物的热容,降低了火焰温度和化学反应速率,使得火焰传播速度降低[26];同时CO2也参与了动力学反应[27-28],反应步骤可表示为

      图  5  C1管中C2H4/N2O和C2H4/N2O/CO2的火焰传播速度随时间变化曲线

      Figure 5.  Flame speed of C2H4/N2O and C2H4/N2O/CO2 as a function of time in the channel C1

      $\rm {N_2}O + \left( M \right)\rightleftharpoons{\;\;\;\;\;\;\;\;} {N_2} + O + \left( M \right)$

      $\rm {N_2}O + H \rightleftharpoons{\;\;\;\;\;\;\;\;} {N_2} + OH$

      $\rm {N_2}O + H \rightleftharpoons{\;\;\;\;\;\;\;\;} NH + NO$

      $\rm C{O_2} + H \rightleftharpoons{\;\;\;\;\;\;\;\;} CO + OH$

      C2H4/N2O为自由基反应,反应步骤(1)为链引发反应,反应步骤(2)和反应步骤(3)为链传递过程,以上几步基元反应是消耗N2O的主要反应,该反应步骤对火焰传播速度的影响很大,CO2通过基元反应(4)转化为CO,减少了H自由基与N2O的反应几率,从而使主体化学反应速率降低,预混气体燃烧速率随之降低,火焰传播速度随之减小[29]

    • 采用NASA的CEA化学平衡软件[30]估算了N2O与C2H4反应的C-J爆轰参数。若系统内各组分的摩尔数没有自发的变化趋势,则称为达到化学平衡。在质量守恒的约束条件下,采用吉布斯自由能最小化为化学平衡判据,建立体系的化学平衡方程

      $G = {{g}} + \sum\limits_{i = 1}^i {{\lambda _i}} {b_i} - b_i^0$

      式中:g为每千克混合物的吉布斯自由能,$b_i^0$为每千克反应物中元素i的量,bi为每千克混合物中元素i的量,λi为拉格朗日乘子。根据变分原理,对方程两边取变分,得到

      ${\text{δ}} G = \sum\limits_{j = 1}^{{\rm{NS}}} {\left( {{\mu _j} + \sum\limits_{i = 1}^t {{\lambda _i}} {a_{ij}}} \right)} {\text{δ}} {n_j} + \sum\limits_{i = 1}^t {\left( {{b_i} - b_i^0} \right)}{\text{δ}}{\lambda _i} = 0$

      式中:µj为每千克物质j的化学势,aij为物质j中元素i的量。由于任何不等于零的小量${\text{δ}} $nj${\text{δ}} $λi是相互独立的,因此可得

      ${\mu _j} + \sum\limits_{i = 1}^t {{\lambda _i}} {a_{ij}} = 0\;\;\;\;\;(j = 1, \ldots,{\rm{NS}})$

      ${b_i} - b_i^0 = 0\;\;\;\;\;(i = 1, \ldots,t)$

      计算时初始温度(T0)和压力(p0)分别设定为298 K和0.1 MPa;反应物的量按照N2O与C2H4反应的理想方程式设定,即

      $\rm 6{N_2}O + {C_2}{H_4} =\!\!=\!\!= 6{N_2} + 2C{O_2} + 2{H_2}O$

      则反应物的C-J爆轰速度(vCJ)为2 200 m/s,爆压(pCJ)为3.90 MPa。

      图6为C4管中火焰阵面传播速度随时间的变化曲线。点火后预混气体经历了点火初期、爆燃转爆轰、稳定爆轰和爆轰衰减4个阶段,爆轰阶段的平均爆速为2 207 m/s,理论C-J爆速为2 200 m/s,实验与理论计算值基本一致,误差为0.3%。点火初期,由于加速环的扰动作用使得燃烧迅速转变为湍流燃烧,随着湍流强度的增大,湍流氛围内较小的涡团能量增强,促进了火焰的正向拉伸,使得火焰传播速度加快,同时在火焰前方形成一个高压流场区域[31],因此预混气体在点火后很快进入爆燃转爆轰阶段。螺旋环的加速过程遵循“热点”机理,预混气体被点燃后,压缩波与障碍物相互作用产生局部热点,随着反应的进行,“热点”逐渐演变为微爆轰,产生的局部爆轰越频繁,火焰传播速度越大[5, 32]图7显示了C4管中压力演变过程。由图7可知,100 mm处压力峰值已达到最大值6.46 MPa,400 mm处压力峰值下降至3.46 MPa,此时已形成典型的爆轰波形,预混气体已完成爆燃转爆轰过程,此后超压峰值逐渐趋于稳定。图8为不同位置压力峰值变化曲线。由图8可知,在1 000~1 600 mm区间内,压力峰值处于相对稳定阶段,波动较小,平均压力峰值为3.92 MPa,理论计算值为3.90 MPa,误差为0.5%。

      图  6  管道C4中火焰前端速度随时间变化曲线

      Figure 6.  Flame speed as a function of time in the channel C4

      图  7  管道C4中8个压力传感器记录的压力曲线

      Figure 7.  Pressure profiles versus time obtained by eight pressure gauges in the channel C4

      图  8  压力峰值变化曲线

      Figure 8.  Plots of maximum overpressure obtained in the channel C4

      冲击波和火焰传播速度随位移的变化关系如图9所示。由图9可知:在100~1 000 mm之间,冲击波速度振荡较大,总体上冲击波速度大于火焰传播速度;在100~400 mm之间,冲击波平均速度已上升至2 325 m/s,在此区间内火焰最大速度约为1 500 m/s;在800~1 000 mm之间,冲击波速度达到最大值3 146 m/s。这是由于气体爆燃属于非理想爆炸,其爆炸能量不是瞬时释放,气体受热膨胀产生前驱压力波,其速度比火焰速度高得多,行进在火焰阵面前方,火焰在压缩过的未燃气中传播。在爆燃转爆轰过程中前驱冲击波阵面仍然行进在火焰面的前方,但二者的间距在减小[33]。在1 000~1 600 mm之间,冲击波振荡幅度变小,冲击波平均速度与火焰速度基本保持一致,燃烧管内形成稳定爆轰,爆轰以2 149~2 461 m/s向未燃区传播,此区间内冲击波和火焰平均速度分别为2 277和2 235 m/s,与C-J理论爆轰速度的相对偏差分别为3.5%和1.6%,实验结果与理论计算基本一致。

      图  9  火焰传播速度、冲击波速度和C-J速度曲线

      Figure 9.  Flame speed, shock wave velocity and C-J velocity curves

    • 图10为组合管道中火焰传播实测图像。右端为长1 000 mm、内径15 mm的有机玻璃加速管,中间为长1 000 mm、内径1.2 mm的不锈钢燃烧管,左端为长500 mm、内径15 mm的有机玻璃验证管。从图10可知,该条件下火焰从右端燃烧管经过中间的不锈钢管道传播至左端的验证管中。

      图  10  组合管道中火焰传播实测图像

      Figure 10.  Image of the flame propagation in combination channels

      图11为预混火焰在组合管道中的传播情况,“×”表示火焰未通过中间不锈钢管,“■”表示通过,l表示点火端玻璃管长度,d表示中间不锈钢管道内径。由图11可知,在本实验条件下,C2H4/N2O的淬火管径(dq)为0.5 mm < dq < 0.7 mm。当d < 0.5 mm时,预混火焰均无法通过中间管道到达验证管;当d为1.2和2.0 mm时,火焰均可通过。由2.1节可知,在大尺寸管道中,边界效应导致预混气体火焰加速度随管径的减小而增大,但当燃烧管道窄至微米级时,管径越小,火焰越容易淬熄。火焰发生淬熄主要由热损失和火焰过度拉伸导致[34]。管径越小,管道表面积与容积的比值越大,火焰通过管壁的热损失越大,当热损失率大于燃烧反应的放热速率时,火焰熄灭[35];过度拉伸使得火焰表面发生扭曲变形,大量未燃气体进入反应区,热扩散作用导致反应区热损失严重,火焰发生熄灭[36]。尤其在变径管道中,火焰更容易因热损失和火焰过度拉伸而淬熄[37]

      图  11  组合管道中预混火焰淬熄情况

      Figure 11.  Diagram of the flame behaviors in combination channels

      d = 0.7 mm,前端加速管l = 100 mm时,验证管中未观察到火焰,但当l ≥ 500 mm时,火焰均可通过不锈钢管到达验证管,因此,增加前端加速管有利于火焰通过不锈钢窄管,即淬熄管径减小。图12为预混火焰在加速管中的传播速度随时间变化曲线。由图12可知,当l为100、500、1 000和2 000 mm时,火焰阵面到达不锈钢管入口处的速度分别为303、671、2 204和2 246 m/s。因此,对于湍流火焰来讲,传播速度越大,淬熄越困难。根据设计阻火器的经验公式[6]

      图  12  加速管长度不同时火焰传播速度随时间变化曲线

      Figure 12.  Flame acceleration process as a function of time in accelerating channels with different lengths

      $ L > 2{S\rm_t} \cdot d\rm_q^2 $

      式中:L为淬熄管道长度,cm;St为湍流火焰传播速度,m·s−1dq为淬熄管径,cm。当dq = 0.5 mm时,加速管中湍流火焰的最大速度为2 533 m/s,根据式(10)可得中间不锈钢管长度应大于12.67 cm,才能使火焰完全熄灭;当dq = 0.3 mm时,L > 4.56 cm。

    • (1)在大直径半封闭有机玻璃管中,预混火焰经历了缓慢加速阶段、爆燃转爆轰阶段和爆轰阶段;点火初期,火焰速度和加速度随着管径的减小而增大,内径为5 mm的管道中火焰加速最快,10和15 mm管道中火焰加速依次减小。CO2会明显降低预混火焰的加速进程和最大爆轰速度,在不改变推进性能的前提下,加入CO2有利于提高预混气体的安全性能。C2H4/N2O的稳定爆速为2 207 m/s,稳定爆轰压力为3.92 MPa,与理论计算值一致。

      (2)实验得到常压下预混气体C2H4/N2O在小直径不锈钢管中淬熄管径为0.5 mm < dq < 0.7 mm;湍流火焰传播速度越大,管径越大,火焰淬熄越困难。同时依据淬熄管径、湍流火焰速度和淬熄长度的关系计算出防回火管道的长度,为阻火器的设计提供参考。

参考文献 (37)

目录

    /

    返回文章
    返回