Correlation between the critical tube diameter and annular interval for detonation wave in high-concentration argon diluted mixtures
-
摘要: 建立圆管及环形管道系统研究临近极限下爆轰波在管道内传播失效机理。选用C2H2+2.5O2+70%Ar气体,采用光纤探针测量爆轰波在管道内传播速度,用烟迹法记录管道内爆轰波胞格结构。结果表明:初始压力远大于爆轰极限压力时,爆轰波在管道内以稳定速度传播;随着初始压力的减小,爆轰波速度逐渐降低;当初始压力一定时,爆轰波速度随着管道尺寸的减小而逐渐减小;当初始压力达到临界压力时,爆轰波在进入到管道内后其速度会逐渐衰减直至爆轰波完全失效。对于不同几何尺寸的圆管与环管,通过引入无量纲参数d/λ及w /λ(d为圆管管径,w为环管间距,λ为爆轰胞格尺寸)得出,爆轰波在管道内传播的临界圆管直径为环形间距的2倍,与理论模型结果相吻合,验证了稳态气体基于爆轰波波面曲率的失效机理。Abstract: Detonation tube including driver section and test section was built to investigate the failure mechanism of detonation wave near the limits. The mixture of C2H2+2.5O2+70%Ar was investigated experimentally. Fiber optics was used to measure detonation velocity. Smocked foils were used to record the detonation cellular structure. The results show that, with the initial pressure far lager than the critical pressure, detonation wave propagates at a constant value in the tubes. Detonation velocity decreases with the decreasing initial pressure. With a given initial pressure, the detonation velocity decreased as the tube diameter (or channel interval) decreased. Under the critical pressure, the detonation velocity propagated a short distance in the tubes and then decreased gradually until complete failure. For different geometries tubes and channels, by introducing dimensionless parameter d /λ and w/λ (d the dameter of the round tube, w the interval of the annular channel and λ the size of detonation cellular), the results show that the critical thickness is half of the critical diameter. Good agreement is found between the experimental measurements in both geometries which supports this conclusion and theoretical mode. The failure mechanisms based on the detonation front curvature for stable detonation in mixtures that are highly argon diluted are well defined.
-
爆轰波在管道内传播时需要满足一定的初始条件(初始压力、当量比等)及边界条件(管道尺寸、几何形状), 若上述条件无法满足, 则爆轰波不能够在管道内自持稳定传播, 并最终失效, 此时达到爆轰极限[1]。爆轰极限的研究对深入理解爆轰波在管道内的传播机理具有重要意义。爆轰波波面是瞬态三维非稳定结构, 且在传播过程中会与管道壁面产生相互扰动[2-4], 并导致爆轰波热量、动量及质量的损失[5]。当条件接近爆轰极限时, 爆轰波的传播变得非常复杂且不稳定, 在介于自持稳定传播与失效之间出现一系列的非稳定传播现象, 如:单头螺旋爆轰、驰振爆轰、低速爆轰等[6-7]。近期研究表明, 对于高浓度氩气稀释气体, 其化学反应区具有层流结构且胞格结构非常规则, 可用ZND模型描述其爆轰波结构[8-10]。对于此类气体, 爆轰波波面非稳定性对于其传播影响很小, 因此, 在研究此类气体的爆轰极限时, 可以忽略爆轰波波面非稳定性因素而仅从边界条件对其传播影响着手研究爆轰极限[11]。根据Fay-Dabora理论模型[12], 管道壁面的边界层效应导致爆轰波化学反应区中的流场扩散, 使得爆轰波波阵面产生弯曲, 导致热量及动量的损失, 从而使得爆轰速度衰减。当爆轰波波面曲率达到某个临界值时, 爆轰波波速特征值不再存在, 此时达到爆轰极限, 爆轰波失效。
本文中研究高浓度氩气稀释气体在2种不同几何形状(圆管、环管)管道内的传播现象及爆轰极限, 对比分析临界管径与临界间距之间的关系, 进一步验证和阐述高浓度氩气稀释气体爆轰波的失效机理。
1. 实验装置
实验装置见图 1。圆管爆轰管道由驱动段和测试段组成(图 1(a)), 驱动段为长1.3 m, 内径为70 mm, 壁厚4.5 mm的钢管。测试段为透明的聚碳酸酯塑料管。环形爆轰管道由驱动段、过渡段及测试段组成(图 1(b)), 其中间距(w)为3.2及5.9 mm的环形通道通过在内径为50.8 mm的聚碳酸酯塑料管内分别插入外径为44.4与39.0 mm光滑铜管构成。圆管及环管几何尺寸见表 1~2, 其中:L为管道长度, t为管壁厚度, d为圆管内径。由于环管间距相对于环管直径很小, 故爆轰波在环管内传播时, 其径向弯曲曲率对其传播的影响可以忽略, 因此可以将其视为二维通道。采用环管而不采用矩形管道的优点在于, 环管能够避免矩形管道的边缘效应[11, 13]。采用高压电火花起爆, 并在驱动段内靠近起爆端添加一长200 mm的Shchelkin螺旋用以增加对爆炸波的扰动, 从而更易在驱动段内形成稳态的爆轰。采用分压法配置C2H2+2.5O2+70%Ar气体。研究表明C2H2+2.5O2+70%Ar混合气体爆轰时具有规则胞格结构, 且化学反应区具有层流结构[1, 4, 9], 因此被称为稳定气体。
表 1 圆管几何尺寸Table 1. Round tube sized/mm L/mm t/mm 1.5 2 438 0.8 3.2 2 438 1.6 12.7 4 118 3.2 31.7 4 118 3.2 50.8 3 048 3.2 表 2 环管几何尺寸Table 2. Annular channel sizew/mm L/mm t/mm 3.2 1 828 3.2 5.9 1 828 3.2 实验前采用真空泵对管道抽真空, 使管道内的压力低于10 Pa, 然后通过阀门控制加入混合气体, 并通过电子压力表控制进气量。在驱动段上安装PCB压力传感器(型号:113B24, 响应频率:500 kHz, 量程:0~6.9 MPa, 输出电压:0~5 V), 用以测定爆轰波的压力和到达各压力传感器的时间, 从而计算在驱动段内是否形成并达到稳态的爆轰波。在测试段管道上等间距布置光纤探针, 通过光纤探针测量测试段内爆轰波到达各测点的时间, 并计算出爆轰波在相邻两测点间的实际传播速度。通过CEA程序[14]计算混合气体在该初始状态下的理论爆轰CJ速度。采用烟迹法记录爆轰波在管道内胞格变化规律并测量对应条件下胞格尺寸。
2. 实验结果及讨论
2.1 爆轰极限临界压力的确定
图 2为不同压力下, C2H2+2.5O2+70%Ar气体在内径为50.8 mm圆形管道内的爆轰波传播速度, x为距管道入口的距离。初始压力远离爆轰极限临界压力时(p0=9 kPa), 爆轰波进入到管道内后, 其平均速度接近于理论爆轰速度vCJ, 并以该速度在管道内稳定传播; 随着初始压力的降低, 爆轰波在管道内传播速度逐渐降低; 初始压力临近极限压力时(p0=1.1 kPa), 爆轰波在进入到管道后仍能够以一稳定速度传播, 此时爆轰波的平均速度
≈0.8vCJ; 若进一步降低预混气体的初始压力(p0=1 kPa), 爆轰波在进入到测试段管道后速度将逐渐衰减, 在距离管道入口大约1 400 mm处, 速度降至0.3vCJ, 此后光学探针已不能检测到火焰信号, 表明爆轰波已经完全失效。因此, 内径为50.8 mm圆管对应的爆轰极限临界压力约为1.1 kPa。对于环形管道, 爆轰波在其内部传播与在圆管管道内类似, 如图 3所示。因此可以同样采用逐步降低预混气体初始压力方式获得对应管道内爆轰极限临界压力。
2.2 初始压力及管道几何尺寸对爆轰速度的影响
图 4为不同管道内, 初始压力与爆轰速度关系。随着初始压力的降低, 爆轰波速度逐渐减小, 这是由于C2H2+2.5O2+70%Ar稳定气体化学反应区具有典型的ZND结构。当混合气体初始压力降低时, 爆轰波诱导区长度增加, 边界层扩散导致的能量损失逐渐增大, 从而导致爆轰波的化学反应变缓, 并使得爆轰速度降低[12]。在临界压力处, 当爆轰波失效时, 最小爆轰速度与管道直径、环形间距无关, 其值约为0.8vCJ。同时由图 4得到, 当初始压力一定时, 爆轰速度随管道尺寸减小而降低。这是由于在小内径管道中, 边界层所占比例增加, 边界层效应增大, 从而导致爆轰速度的降低。
对于圆形管, 当管径(d)为1.5、3.2、12.7、31.7、50.8 mm时, 对应的爆轰极限极限压力(pc)分别为20、9.8、3.2、1.5、1.1 kPa; 对于环形管, 当间距(w)为3.2和5.9 mm时, 对应的爆轰极限极限压力分别为7.0和3.5 kPa。通过比较不同几何尺寸管道所对应的临界压力可以看出, 爆轰极限临界压力随着管径或间距的减小而显著增大, 说明随着管径或间距的减小, 管道内壁边界层所占的比例不断增加, 在小尺寸管道内管道边界层效应对爆轰波在管道内的传播起到主导作用。
2.3 临界圆管内径与环管间距的关系
当混合气初始压力一定时, 圆管临界管径尺寸为环管临界间距的2倍[12], 这是因为爆轰波在环管内传播时, 由于环管可视为二维通道, 爆轰波波面弯曲只发生在一个平面上, 而对于圆管, 爆轰波波面弯曲则发生在两个相互垂直的平面上。本实验中共采用了5种不同内径的圆管及2种间距的环管, 并得到其对应的临界压力。由于临界压力值不同, 为了消除初始压力不同对结果的影响, 引入无量纲参数dc/λ及wc/λ, 通过无量纲化的方法消除初始压力对结果的影响, 并得到一个普适规律。无量纲参数dc/λ及wc/λ中, dc、wc和λ分别为临界管道直径、临界环形间距及胞格尺寸, 其中dc和wc可反映爆轰波在管道内传播的边界条件, 而胞格尺寸作为爆轰动态参数[1], 可反映爆轰气体的敏感度。对于给定的气体, 胞格尺寸大小与初始压力相关, 因此胞格尺寸可间接反映爆轰气体初始压力。
图 5为典型的爆轰波胞格结构随初始压力变化规律(圆管, d=50.8 mm)。由图可知, C2H2+2.5O2+70%Ar具有规则的胞格结构, 且随着初始压力的逐渐降低, 爆轰波胞格尺寸逐渐增加。通过实验测定了不同管道及初始压力下爆轰波胞格尺寸(如图 6所示), 并采用了最小二乘法拟合出胞格尺寸与初始压力之间的关系式:
λ=320.82p−1.33580 (1) 式中:λ单位为mm, p0单位为kPa。通过(1)式, 可计算得到各管道在临界压力处对应的胞格尺寸, 并计算出相应的dc/λ、wc/λ, 结果见表 3~4。由表 3~4可以得出, dc/λ和wc/λ的平均值分别为0.185和0.095 4。对于不同内径的圆管及环管, 在极限压力处, dc/λ与wc/λ的平均比值为1.95, 这与理论结果(dc/wc=2)相吻合, 证明对于稳态气体, 可采用理论分析预测其爆轰极限。
表 3 圆管实验结果Table 3. Experimental results of round tubesd/mm pc/kPa λ dc/λ 50.8 1.1 306.47 0.166 31.7 1.5 202.51 0.157 12.7 3.2 73.60 0.173 3.2 9.8 16.50 0.194 1.5 20.0 6.36 0.236 表 4 环管实验结果Table 4. Experimental results of annular channelsw/mm pc/kPa λ wc/λ 5.9 3.5 65.30 0.0903 3.2 6.0 31.79 0.1006 对于高浓度氩气稀释的稳态气体, 其胞格结构非常规则, 因此可采用经典的ZND模型进行分析, 故除了采用胞格尺寸所建立无量纲参数外, 亦可采用ZND诱导区长度(ΔZND)建立无量纲参数, 以此分析临界管径与临界间距间关系。ZND诱导区长度可通过Chemkin程序[15]和San Diego化学反应动力学原理[16]计算得到, 结果见表 5和6。B.Bazathazajan等[17]证明了San Diego化学反应机理对于计算C2H2+2.5O2+70%Ar气体爆轰波ZND诱导区长度的有效性。由表 5和6可知, dc/ΔZND与wc/ΔZND的平均值分别为24.07和12.25。dc/wc的平均值为1.96。当采用ZND诱导区长度建立无量纲参数时, 结果同样与理论结果吻合。进一步证明了采用爆轰波波面曲率失效机理鱼刺爆轰极限的正确性。
表 5 圆管ZND化学反应区长度Table 5. ZND reaction zone length of round tubesd/mm pc/kPa ΔZND/mm dc/ΔZND 50.8 1.1 1.818 27.94 31.7 1.5 1.273 24.94 12.7 3.2 0.553 22.96 3.2 9.8 0.148 21.54 1.5 20.0 0.065 23.08 表 6 环管ZND化学反应区长度Table 6. ZND reaction zone length of annular channelsw/mm pc/kPa ΔZND/mm wc/ΔZND 5.9 3.5 0.482 12.2 3.2 6.0 0.259 12.3 3. 结论
(1) 当初始压力远大于爆轰极限临界压力时, 爆轰波在管道内以一稳定速度传播。随着预混气初始压力的减小, 爆轰波在管道内传播的速度逐渐降低。当初始压力达到临界压力时, 爆轰波在进入到管道内并传播一段距离后, 其速度将逐渐衰减直至爆轰波衰变成爆燃波并完全失效。通过比较不同几何尺寸管道内爆轰波传播速度可知, 爆轰波速度随着管道几何尺寸的减小而减小, 表明爆轰波在管道内传播时受到边界条件及初始条件的影响。
(2) 对于不同几何尺寸管道, 通过无量纲参数分析研究得出, 爆轰波在管道内传播时临界圆管直径为临界环管间距的2倍, 这与理论预测相符合。表明高浓度氩气稀释气体爆轰波在管道内传播时, 管道壁面的边界层效应所导致的爆轰波化学反应区内的动量、热量损失是导致此种高浓度氩气稀释气体爆轰波在管道内传播失效的主要原因。通过对爆轰波在圆形及环形管道内的失效分析, 进一步验证了此种高浓度氩气稀释气体基于爆轰波波面曲率的失效机理。
-
表 1 圆管几何尺寸
Table 1. Round tube size
d/mm L/mm t/mm 1.5 2 438 0.8 3.2 2 438 1.6 12.7 4 118 3.2 31.7 4 118 3.2 50.8 3 048 3.2 表 2 环管几何尺寸
Table 2. Annular channel size
w/mm L/mm t/mm 3.2 1 828 3.2 5.9 1 828 3.2 表 3 圆管实验结果
Table 3. Experimental results of round tubes
d/mm pc/kPa λ dc/λ 50.8 1.1 306.47 0.166 31.7 1.5 202.51 0.157 12.7 3.2 73.60 0.173 3.2 9.8 16.50 0.194 1.5 20.0 6.36 0.236 表 4 环管实验结果
Table 4. Experimental results of annular channels
w/mm pc/kPa λ wc/λ 5.9 3.5 65.30 0.0903 3.2 6.0 31.79 0.1006 表 5 圆管ZND化学反应区长度
Table 5. ZND reaction zone length of round tubes
d/mm pc/kPa ΔZND/mm dc/ΔZND 50.8 1.1 1.818 27.94 31.7 1.5 1.273 24.94 12.7 3.2 0.553 22.96 3.2 9.8 0.148 21.54 1.5 20.0 0.065 23.08 表 6 环管ZND化学反应区长度
Table 6. ZND reaction zone length of annular channels
w/mm pc/kPa ΔZND/mm wc/ΔZND 5.9 3.5 0.482 12.2 3.2 6.0 0.259 12.3 -
[1] Lee J H S. The detonation phenomenon[M]. Cambrige, UK: Cambridge University Press, 2008. [2] Dupre G, Peraldi O, Lee J H S, et al. Propagation of detonation waves in an acoustic absorbing-walled tube[J]. Progress in Astronautics and Aeronautics, 1988, 114: 248-263. doi: 10.2514/5.9781600865886.0248.0263 [3] Teodoczyk A, Lee J H S. Detonation attenuation by foams and wire meshes lining the walls[J]. Shock Waves, 1995, 4(4): 225-236. doi: 10.1007/BF01414988 [4] Raduledscu M I, Lee J H S. The failure mechanism of gaseous detonations: Experiments in porous wall tubes[J]. Combustion and Flame, 2002, 131: 29-46. http://www.sciencedirect.com/science/article/pii/S0010218002003905 [5] Zeldovich Y B. On the theory of the propagation of detonation in gaseous systems[R]. Soviet Union: Soviet Physics-JETP, 1940. [6] Lee J J, Dupre G, Knystautas R, et al. Doppler interferometer study of unstable detonations[J]. Shock Waves, 1995, 5(3): 175-181. [7] Camargo A, Ng H D, Chao J, et al. Propagation of near-limit gaseous detonations in small diameter tubes[J]. Shock Waves, 2010, 20(6): 499-508. doi: 10.1007/s00193-010-0253-3 [8] Radulescu M I. The propagation and failure mechanism of gaseous detonations: Experiments in porous-walled tubes[D]. Montreal, Canada: McGill University, 2003. [9] Radulescu M I, Ng H D, Lee J H S, et al. The effect of argon dilution on the stability of acetylene-oxygen detonations[J]. Proceedings of the Combustion Institute, 2002, 29(2): 2825-2831. https://www.sciencedirect.com/science/article/pii/S1540748902803455 [10] Ng H D, Radulescu M I, Higgins A J, et al. Numerical investigation of the instability for one-dimensional Chapman-Jouguet detonations with chain-branching kinetics[J]. Combustion Theory and Modeling, 2005, 9: 385-401. doi: 10.1080/13647830500307758 [11] Chao J, Ng H D, Lee J H S. Detonation limits in thin annular channels[J]. Proceedings of the Combustion Institute, 2009, 32(2): 2349-2354. http://www.sciencedirect.com/science/article/pii/S1540748908003660 [12] Fay J A. Two-dimensional gaseous detonations: Velocity deficit[J]. Physics of Fluids, 1959, 2(3): 283-289. doi: 10.1063/1.1705924 [13] Meredith J, Ng H D, Lee J H S. Detonation diffraction from an annular channel[J]. Shock Waves, 2010, 20(6): 449-455. doi: 10.1007/s00193-010-0256-0 [14] Mcbride B J, Gordon S. Computer program for calculation of complex chemical equilibrium compositions and applications[R]. NASA, 1996. [15] Kee R J, Rupley F M, Millerja. Chemkin-Ⅱ: A fortran chemical kinetics package for the analysis of gas-phase chemical kinetics[R]. Sandia National Laboratories, 1989. [16] The San Diego Mechanism. Chemical-kinetic mechanisms for combustion applications[EB/OL]. http://web.eng.ucsd.edu/mae/groups/combustion/mechanism.html. [17] Varatharajan B, Williams F A. Chemical-kinetic descriptions of high-temperature ignition and detonation of acetylene-oxygen diluents systems[J]. Combustion and Flame, 2001, 124(4): 624-645. https://www.sciencedirect.com/science/article/pii/S0010218000002352 期刊类型引用(4)
1. 韩文虎,张博,王成. 气相爆轰波起爆与传播机理研究进展. 爆炸与冲击. 2021(12): 65-96 . 本站查看
2. 赵焕娟,伯玉兰,张英华,严屹然. Ar稀释C_2H_2+2.5O_2预混气高频爆轰的端面结构. 爆炸与冲击. 2018(05): 1121-1129 . 本站查看
3. 赵焕娟,John H.S.Lee,张英华,严屹然. 爆轰波三波点擦除烟迹表面积碳机制. 工程科学学报. 2017(03): 335-341 . 百度学术
4. 喻健良,管清韦,闫兴清,姚福桐,张东,张师哲. 初始条件对管道内爆轰波传播特性影响研究. 科学技术与工程. 2017(04): 126-131 . 百度学术
其他类型引用(3)
-