高压氢气泄漏自燃形成喷射火的实验研究

闫伟阳; 潘旭海; 汪志雷; 华敏; 蒋益明; 王清源; 蒋军成

doi:10.11883/bzycj-2018-0394

高压氢气泄漏自燃形成喷射火的实验研究

doi: 10.11883/bzycj-2018-0394

闫伟阳^1,,
潘旭海^{1, 2, ,},
汪志雷¹,
华敏^{1, 2},
蒋益明¹,
王清源¹,
蒋军成^{1, 2}

1.
南京工业大学安全科学与工程学院，江苏南京 210009
2.
南京工业大学江苏省危险化学品本质安全与控制技术重点实验室，江苏南京 210009

基金项目: 国家重点研发计划（2017YFC0804700，2016YFC0800100）

详细信息

作者简介:
闫伟阳（1992－　），男，硕士研究生，455314754@qq.com

通讯作者:
潘旭海（1977－　），男，博士，教授，xuhaipan@njtech.edu.cn

中图分类号: O389
计量
- 文章访问数: 6465
- HTML全文浏览量: 1680
- PDF下载量: 84
- 被引次数: 1
出版历程
- 收稿日期: 2018-10-16
- 修回日期: 2019-02-12
- 刊出日期: 2019-11-01

Experimental investigation on spontaneous combustion of high-pressure hydrogen leakage to form jet fire

YAN Weiyang^1
,,
PAN Xuhai^{1, 2
, ,},
WANG Zhilei¹,
HUA Min^{1, 2},
JIANG Yiming¹,
WANG Qingyuan¹,
JIANG Juncheng^{1, 2}

1.
College of Safety Science and Engineering, Nanjing Tech University, Nanjing 210009, Jiangsu, China
2.
Jiangsu Key Laboratory of Intrinsic Safety and Control of Hazardous Chemicals, Nanjing Tech University, Nanjing 210009, Jiangsu, China

摘要

摘要: 为了探究高压氢气泄漏发生自燃时所需的临界初始释放压力随管道长度的变化规律，了解管内自燃火焰向管外喷射火焰转变的发展过程，本文利用压力、光电以及高速摄像等测试系统展开实验研究。实验结果表明：当管道长度相同，初始释放压力较低时，氢气泄漏不容易发生自燃；随着管道长度的增加，氢气发生自燃时的临界初始释放压力先缓慢减小后迅速增大；当管道长度一定时，初始释放压力越大，激波传播速度越快，氢气管内自燃的位置距离爆破片越近；气流通过激波马赫盘后，火焰燃烧加剧；随着时间的增加，火焰长度呈现先增大后逐渐减小的变化趋势，喷射火焰尖端的平均传播速度逐渐减小；火焰宽度呈现先增大后迅速减小至稳定值的变化规律。
- 氢气 /
- 自燃 /
- 管道 /
- 激波 /
- 喷射火焰
Abstract: In this paper, to investigate the variation of the critical initial release pressure with the pipeline length in high-pressure hydrogen leakage that leads to spontaneous combustion and the transition process from spontaneous combustion flame inside the tube to the jet flame outside the tube, we conducted experiments using a pressure gauge, a photoelectric and high-speed camera, etc. Our results showed that, at the same pipeline length and under a low initial release pressure, hydrogen is not apt to spontaneous combustion. The minimum initial release pressure of hydrogen spontaneous combustion decreases slowly and then increases rapidly with as the pipe length increases. At the same pipeline length, the greater the initial release pressure, the faster the shock wave propagation, and the closer the hydrogen self-ignition position inside the pipe to the rupture disc. It is found that the flame combustion is intensified after the airflow passed through the Mach disk. With the increase of time, the flame length increases first and then decreases gradually, the average propagation speed of the jet flame tip decreases gradually. The flame width increases first and then decreases rapidly to a stable value.
- hydrogen /
- spontaneous combustion /
- pipeline /
- shock wave /
- jet fire

HTML全文

对爆炸加载壳体膨胀断裂的研究主要集中在膨胀断裂机理和破片尺寸及其速度分布2个方面。在破片尺寸及速度分布研究方面：Gurney^[1]提出了破片速度预估模型；Mott^[2]、Grady等^[3]、Hopson等^[4]、Zhou等^[5]和郑宇轩等^[6]基于一维理论提出了破片尺寸分布模型^[2-6]；对于二维和三维情况，理论预测的破片尺寸与实际差异较大^[7-8]。在壳体膨胀断裂机理方面，研究主要涉及柱壳的环向拉伸和剪切2种断裂模式，给出了断裂判据以及直观的裂纹扩展演化描述^[9-11]。其基本依据是壳体内的应力分布和材料断裂判据。从研究情况看，壳体膨胀断裂的萌生位置分别有内表面、外表面、内表面附近等不同认识，对应的扩展路径和方式也有不同认识^[9-14]。这些认识差异的起源，既有壳体内应力分布演化复杂，也有诊断信息不足造成分析难以收敛的原因。

对于壳体内的应力分布，早期研究一般采用二维轴对称模型的解析解来简化分析滑移爆轰加载柱壳，未考虑滑移爆轰带来的沿柱壳轴向的应力梯度^[9-11]；随着数值模拟技术的应用，研究者发现滑移爆轰加载柱壳内的空间应力梯度分布具有显著、不可忽略的影响^[15-17]。因此，将二维轴对称简化模型用于滑移爆轰加载的膨胀柱壳断裂分析，存在较大的认识偏差。对于数值模拟技术，虽然可以提供裂纹萌生前的应力分布描述，但是受材料模型参数获取难度的限制，其应用相对比较困难。因此，膨胀壳体断裂研究目前主要还是依赖实验诊断和简化应力分析开展。

近年来发展起来的爆炸丝起爆技术（简称线起爆技术）^[18-19]利用大电流通过金属丝并使其等离子体化，然后沿丝的长度方向同步起爆柱壳内部装填的炸药，从而能够在膨胀柱壳中实现一维轴对称应力的加载。在此状态下，就可以使用一维轴对称应力分析模型对线起爆膨胀柱壳的断裂问题进行有效分析。

在壳体膨胀断裂诊断方面，常用手段包括高速分幅照相^{[11, 16]}、高速狭缝扫描照相^[20]、X射线透射照相^[21]、DPS (Doppler detection system)测速^[16]。其中，高速分幅照相是常用诊断手段之一，可观测壳体表面发生的皱褶、冒烟（爆轰产物泄漏）以及外径尺寸等特征信息。以往研究一般采用冒烟特征作为壳体断裂诊断标准^[11]，也有一些研究者将皱褶特征作为壳体断裂的诊断标准^[22]。按照胡八一等^[11]的认识，从壳体表面皱褶发展到冒烟需要几微秒到十几微秒。高速狭缝扫描照相可观测空间固定位置上的壳体轮廓投影的不连续变化（如冒烟、宏观断裂等）和外径尺寸信息。X射线透射照相可获取壳体沿照相投影方向的密度分布，进而通过密度梯度变化识别壳体的裂纹宽度及分布、外径尺寸等特征信息。X射线透射照相的优点是不受撞击或者爆炸发光影响，缺点是图像的对比度和清晰度较差，定量分析精度低^[23]。DPS测速是近年来发展起来的常用诊断手段，用于测量壳体表面的垂直运动速度。DPS测速的优点是可以长时间精确测量。受柱壳弧度影响，上述照相诊断手段一般只能准确监测柱壳圆周部分角度范围内的断裂状态，因此诊断结果可能存在一定程度（或严重）的偏差。

对于均匀承载柱壳，断裂将使局部壳体承载失效，进而导致壳体的应力、应变和速度出现显著的非均匀分布。中低应变率下的传统断裂诊断方法是采用粘贴在裂纹附近的应变计监测应变扰动来判读断裂时刻^[24-25]。由于爆轰加载强度较高，不满足应变计的使用要求，因此可以采用一定数量、沿空间分布的DPS探头来监测壳体表面速度分布的演化，实现柱壳圆周范围内的断裂诊断。

1. 实　验

采用线起爆技术对装填粉末PETN (pentaerythritol tetranitrate)炸药的金属柱壳做一维柱面加载，装置结构见图1。线起爆金属丝安装在柱壳中轴位置。炸药与柱壳之间填充延展性较好的尼龙来约束粉末装药，同时抑制尼龙层碎裂而破坏加载均匀性。装填炸药的密度控制在(1.0～1.1)×10³ kg/m³范围内，直径为15 mm。柱壳内径为40 mm，外径为48 mm。装药、尼龙层和柱壳高度相同，均为160 mm。柱壳材料为304钢和45钢。每种柱壳材料各开展1次实验。

图 1 线起爆实验装置

Figure 1. The experimental device with linear initiation

下载: 全尺寸图片幻灯片

304钢柱壳实验仅采用DPS探头测量柱壳外壁的径向速度，以监测柱壳膨胀和断裂状态。在柱壳80 mm高度处，沿环向0°、90°、180°、270°方向各布置1个DPS探头，监测径向速度的对称性；基于柱壳沿高度方向的中心对称原则，在沿环向90°方向、80～140 mm高度范围内，按照20 mm等高度间隔布置了4个DPS探头，监测径向速度的一致性。

45钢柱壳实验除采用DPS探头测量柱壳外壁速度外，还增加了高速分幅照相检测内容，以直观观测柱壳膨胀和断裂的宏观形态。其中，照相光路占据沿环向0°方向，因此该方向上的DPS探头取消。另外，考虑柱壳端部边侧稀疏对径向速度的影响，取消了沿环向90°方向、140 mm高度处的DPS探头。其余DPS探头布局与304钢柱壳实验的相同。此外，在45钢柱壳外壁20～140 mm高度范围内，按照20 mm等高度间隔分别画了蓝色环形高度标线。

在时序控制方面，以金属丝通电起爆时刻为零时刻，并通过系统设置保障DPS测速和高速分幅照相与金属丝起爆的时基相同。两相邻高速分幅照片的时间间隔约为2 μs。

2. 结　果

304钢柱壳的径向速度监测结果见图2。受柱壳端部边侧稀疏卸载影响，140 mm高度处的速度明显低于柱壳中部其他高度处的速度。45钢柱壳的速度监测结果和高速分幅照相结果分别见图3和图4。图4(b)中5条裂纹带的最早出现时刻依次为39.2 μs （45°）、43.3 μs （0°）、47.4 μs （270°）、49.4 μs （315°）和51.5 μs （90°）。各相邻裂纹带的间距估算约为2 cm。受柱壳弧度影响，其中270°和90°裂纹带的最早出现时刻判读明显滞后。

图 2 304钢柱壳外壁的径向速度曲线

Figure 2. The radial velocity curves of the outer surface of the 304 steel cylinder

下载: 全尺寸图片幻灯片

图 3 45钢柱壳外壁的径向速度曲线

Figure 3. The radial velocity curves of the outer surface of the 45 steel cylinder

下载: 全尺寸图片幻灯片

图 4 45钢柱壳的高速分幅照相结果

Figure 4. The high-speed framing photography results of the 45 steel cylinder

下载: 全尺寸图片幻灯片

在图4(a)的有效段内，柱壳轮廓基本为直圆柱状。而在有效段外，圆柱轮廓有轻微的直径收缩。这一形态与图2中速度曲线积分获得的膨胀柱壳轮廓一致，即如图5中黑色曲线所示的中部平直凸出、两侧略滞后的鼓形轮廓。图6显示了图2中各速度曲线的对应位移状态。表1列出了柱壳有效段内各测点处的速度起跳时刻， t₁₁～t₁₃分别为沿环向90°方向，在80、100和120 mm高度处柱壳速度的起跳时刻； t₂₁～t₂₄分别为沿环向0°、90°、180°和270°方向，在80 mm高度处柱壳各测点的速度起跳时刻。

图 5 304钢柱壳的初始轮廓与运动35 μs后的轮廓

Figure 5. The contours of the outer-surface of the 304 steel cylinder at the initial state and the deformed state after 35 μs

下载: 全尺寸图片幻灯片

图 6 304钢柱壳沿环向90°方向、不同高度的外壁位移状态

Figure 6. The displacement of the outer surface of the 45 steel cylinder at various heights and the circle angle of 90°

下载: 全尺寸图片幻灯片

表 1 304钢和45钢柱壳外壁各测点的速度曲线起跳时刻

Table 1. Jump-up times in velocity curves of the outside surfaces of the 304 steel and 45 steel cylinders

柱壳材料	t₁₁/μs	t₁₂/μs	t₁₃/μs	t₂₁/μs	t₂₂/μs	t₂₃/μs	t₂₄/μs
304钢	10.27	10.17	10.09	10.27	10.27	10.14	10.25
45钢	12.13	11.97	11.87		12.13	12.11	12.06

下载: 导出CSV

| 显示表格

3. 加载状态

加载状态分析用于确认加载满足一维柱面状态，从而为承载结构的应力分析提供简单的状态基础。从图2(a)和3(a)来看，80～120 mm高度范围的各速度曲线在柱壳断裂前基本保持重合，说明加载响应处于良好的柱面状态。而从图2(b)和图3(b)来看，沿环向各角度的速度曲线也基本保持重合，说明加载响应的轴对称性良好。因此，可以定性判断2种钢柱壳均处于良好的一维柱面状态。

在定量判断方面，考虑加载状态偏离一维柱面将造成波形前沿出现一定程度的空间离散。定义空间离散度为单位长度波形前沿的最大距离差。为了便于实验监测，假定加载波在柱壳内稳定传播，将波形前沿的最大距离差转化为最大时刻差，得到空间离散度：

$\psi = \frac{{{{\Delta }}s}}{{{{\Delta }}L}} = \frac{{c{{ \Delta }}t}}{{{{\Delta }}L}}$

(1)

式中：ψ为空间离散度，Δs为波形前沿最大距离差，ΔL为Δs的统计长度，Δt为波形前沿最大时刻差（或速度起跳时刻晃动），c为应力波在柱壳中沿径向传播的速度。

按照一般小量考虑，可接受的空间离散度应不超过5%。因此，由式(1)可得速度起跳时刻晃动的要求如下：

${{\Delta }}t {\text{≤}} \frac{{{{\Delta }}L}}{c} \times 5\text{%}$

(2)

考虑柱壳在自由膨胀阶段的静水压相对较低，加载波的传播速度可简单近似为材料声速（计为5 km/s）。对于外壁半径为24 mm、有效段长80 mm的柱壳段，分别按照沿外壁环绕一周、沿轴向80 mm的长度统计柱壳有效段的波形前沿最大距离差，则可以通过式(2)计算出满足一维柱面状态的径向速度起跳时刻晃动上限：沿环向的时刻晃动Δt≤0.24 μs，沿轴向的时刻晃动Δt≤0.80 μs。从表1的数据来看，2种钢柱壳沿环向不同角度的速度起跳时刻晃动（0.13、0.07 μs）和沿轴向不同高度处的速度起跳时刻晃动（0.18、0.26 μs）均小于上述2个限值，表明柱壳近似处于一维柱面加载状态。

4. 加载过程及断裂诊断分析

图2(b)、3(b)显示，柱壳外壁的速度曲线呈现多峰值加载特征。根据应力波分析，在各峰值前的加速阶段，尼龙层对柱壳做碰撞加载；在峰值后的减速阶段，柱壳与尼龙层脱离、自由膨胀并在环向拉伸应力作用下近似线性减速。减速阶段的柱壳自由膨胀特征表现为速度峰值前后的起跳时刻间隔大于应力波在尼龙层内/外壁来回反射一次的时间（约11 μs）以及各减速阶段相近的加速度值（−(4.0～5.0) ×10⁶ m/s²）。

在图3(b)的第2个峰值后的减速阶段，45钢柱壳在初始约5 μs时段内保持均匀速度分布，并且环向承载水平与第1个减速阶段的近似相等（a′≈a₁）。从37 μs时刻开始，沿环向90°方向的速度曲线呈现出较长时间的迅速下降趋势。假设此时段的柱壳仍为完整承载体，则该方向上的环向承载水平明显高于其余2个方向的（a″≈2a′），而其余2个方向的环向承载水平未发生变化。考虑到柱壳处于自由膨胀状态，其环向承载水平仅由柱壳承载能力控制。因此，由a″≈2a′关系可进一步推断该方向的柱壳承载力出现大幅度的阶跃强化。该推断结果与柱壳材料流动应力连续演化和柱壳质量沿环向均匀分布或连续演化（对应颈缩状态）的认识相悖，因此说明假设错误。柱壳在沿环向90°方向附近应发生了初始断裂并丧失环向承载能力，导致该方向的柱壳径向速度不受环向拉伸应力控制，而其余2个方向的柱壳则保持完整、环向承载水平未发生变化。

从图2(b)的第3个峰值时刻开始，沿环向0°方向的304钢速度曲线呈现出较长时间的水平演化趋势（加速度值为零），而其他3个方向的速度曲线则呈现自由膨胀的减速趋势（a′≈a₁≈a₂）。同样假设304钢柱壳为完整承载体，则沿环向0°方向的环向承载水平衰减至零，而其余3个方向的环向承载水平与之前各减速阶段的近似相等。由此推断，柱壳在沿环向0°方向附近发生断裂、丧失环向承载能力，而其余各方向上的柱壳则保持完整；如假设不成立，同样可以得出上述结论，因此说明该结论是唯一、真实的。

为了直观展示分布式表面速度诊断方法的效果，将高速分幅照相诊断的45钢柱壳断裂情况与速度诊断结果进行了对比。在空间关系上，速度诊断的初始断裂位置与图4(b)中沿环向90°方向的裂纹带位置重叠，说明速度诊断的断裂位置是正确的；在时间关系上，速度诊断的初始断裂时刻（37.0 μs）明显早于同方向上照相诊断的裂纹最早出现时刻（约51.5 μs），并稍早于照相诊断的初始断裂时刻（39.2 μs），说明速度诊断方法可消除柱壳弧度影响，速度诊断的初始断裂时刻为柱壳的初始断裂时刻。

表2是基于图2(b)和3(b)中的速度曲线诊断得到的304钢和45钢柱壳的初始断裂参数。其中，断裂时刻根据速度曲线簇的分叉（或演化趋势变化）时刻判读；断裂应变通过断裂时刻前的速度曲线积分获得的位移量计算；平均应变率为柱壳外表速度起跳至断裂时刻的应变率平均值，用来近似表征柱壳开始承载变形至断裂时刻的应变率水平。从表2的断裂参数看，45钢和304钢柱壳的一维柱面动态拉伸断裂性能存在明显差异。45钢柱壳的断裂应变相对较小，其延展性弱于304钢柱壳。

表 2 钢柱壳的初始断裂参数

Table 2. Initial fracture parameters of steel cylinders

材料	断裂时刻/μs	断裂应变/%	平均应变率/s⁻¹
304钢	51	37	0.9×10⁴
45钢	37	24	1.0×10⁴

下载: 导出CSV

| 显示表格

5. 结　论

采用电爆炸丝起爆PETN粉末炸药并驱动尼龙对金属柱壳加载的方式，结合分布式表面速度监测和高速分幅照相监测，开展了2种钢柱壳的膨胀断裂实验。通过实验研究，实现了金属柱壳的一维柱面均匀加载，建立了可监测柱壳圆周范围内初始断裂信息的分布式表面速度诊断方法，并获得了45钢和304钢柱壳的初始断裂参数（含断裂应变、平均应变率），具体结论如下。

(1)基于线起爆粉末炸药并驱动尼龙加载的方式，可在金属柱壳中部、半柱高范围内实现一维柱面均匀加载。柱壳外表速度起跳时刻的晃动满足加载状态判据要求。

(2)基于均匀承载壳体断裂引起的局部承载失效将导致均匀分布的速度曲线簇出现分叉（或演化趋势变化）的原理，分布式表面速度诊断方法可准确获取柱壳圆周范围内的初始断裂信息。标记沿柱壳环向分布的径向速度曲线簇出现分叉（或演化趋势变化）的时刻为断裂时刻，标记与速度曲线簇分叉的曲线对应区域为断裂位置。

(3)在相同的一维柱面动态膨胀加载条件下，45钢柱壳的断裂应变（或延展性）低于304钢柱壳的。

感谢金山、但加坤、郭昭亮、陈浩玉、罗振雄、张振涛、莫俊杰、李军、赵延安、谢明强、王远在实验设计、实施及结果分析方面给予的帮助。

图 1 高压氢气泄漏自燃的扩散点火理论

Figure 1. Diffusion ignition theory of high-pressure hydrogen self-ignition

下载: 全尺寸图片幻灯片

图 2 高压氢气泄漏模拟实验装置示意图

1. Pressure reducing valve；2. High-pressure pneumatic valve；3. Controller；4. Pressure gauge；5. Tanks；6. Vacuum pump；7. Rupture disc；8. Data processor；9. High-speed camera；10. Computer

Figure 2. Illustration of the experimental platform

下载: 全尺寸图片幻灯片

图 3 部分零件实物图

Figure 3. Some actual of parts

下载: 全尺寸图片幻灯片

图 4 氢气自燃时管道长度与初始释放压力的关系

Figure 4. Pipe lengths vs. initial release pressure of hydrogen self-ignition

下载: 全尺寸图片幻灯片

图 5 压力传感器P₁监测到的压力波形变化

Figure 5. Pressure waveform variation monitored by pressure sensor P₁

下载: 全尺寸图片幻灯片

图 6 管道长度1 700 mm下的2组实验压力图

Figure 6. Pressure of pipe length 1 700 mm

下载: 全尺寸图片幻灯片

图 7 管道长度1 700 mm下的2组实验光电图

Figure 7. Photo-electricity of pipe length 1 700 mm

下载: 全尺寸图片幻灯片

图 8 管长700 mm不同初始释放压力下的激波传播速度图

Figure 8. Shock wave propagation speed of pipe length 700 mm under different initial release pressures

下载: 全尺寸图片幻灯片

图 9 管内火焰向管外喷射火焰转变的过程图

Figure 9. Transition of flame inside pipe to the jet flame outside pipe

下载: 全尺寸图片幻灯片

图 10 管长700 mm初始释放压力7.08 MPa下管外火焰的发展图

Figure 10. Flame development process of pipe length 700 mm and initial release pressure 7.08 MPa

下载: 全尺寸图片幻灯片

图 11 管长300 mm和700 mm初始释放压力相近的管外火焰形态变化

Figure 11. Flame development process with similar initial release pressure of 300 mm and 700 mm

下载: 全尺寸图片幻灯片

图 12 管长1 200、1 700和2 200 mm初始释放压力相近的管外火焰形态变化

Figure 12. Flame development process with similar initial release pressure of 1 200, 1 700 and 2 200 mm

下载: 全尺寸图片幻灯片

参考文献(22)

[1]	NICOLETTI G, ARCURI N, NICOLETTI G, et al. A technical and environmental comparison between hydrogen and some fossil fuels [J]. Energy Conversion and Management, 2015, 89(89): 205–213.
[2]	ONO R, ODA T. Spark ignition of hydrogen-air mixture [C] // 2008: 012003. DOI: https://doi.org/10.1088/1742-6596/142/1/012003.
[3]	XU B P, WEN J X, DEMBELE S, et al. The effect of pressure boundary rupture rate on spontaneous ignition of pressurized hydrogen release [J]. Journal of Loss Prevention in the Process Industries, 2009, 22(3): 279–287. DOI: 10.1016/j.jlp.2008.07.007.
[4]	WOLINSKI M, WOIANSKI P. Investigation into the mechanism of the diffusion ignition of a combustible gas flowing into an oxidizing atmosphere [C] // Proceedings of the 14th symposium on combustion, 1973: 1217−1223.
[5]	FREDERICK L D, MARCOS C, ZHENWEI Z, et al. Spontaneous ignition of pressurized releases of hydrogen and natural gas into air [J]. Combustion Science and Technology, 2007, 179(4): 663–694.
[6]	MOGI T, WADA Y, OGATA Y, et al. Self-ignition and flame propagation of high-pressure hydrogen jet during sudden discharge from a pipe [J]. International Journal of Hydrogen Energy, 2009, 34(14): 5810–5816. DOI: 10.1016/j.ijhydene.2009.04.079.
[7]	MOGI T, KIM D, SHIINA H, et al. Self-ignition and explosion during discharge of high-pressure hydrogen [J]. Journal of Loss Prevention in the Process Industries, 2008, 21(2): 199–204. DOI: 10.1016/j.jlp.2007.06.008.
[8]	LEE H J, KIM Y R, KIM S H, et al. Experimental investigation on the self-ignition of pressurized hydrogen released by the failure of a rupture disk through tubes [J]. Proceedings of the Combustion Institute, 2011, 33(2): 2351–2358. DOI: 10.1016/j.proci.2010.06.040.
[9]	KITABAYASHI N, WADA Y, MOGI T, et al. Experimental study on high pressure hydrogen jets coming out of tubes of 0.1-4.2 m in length [J]. International Journal of Hydrogen Energy, 2013, 38(19): 8100–8107. DOI: 10.1016/j.ijhydene.2012.10.040.
[10]	段强领. 高压氢气泄漏自燃机理及其火焰传播特性实验研究[D]. 合肥: 中国科学技术大学, 2016. DUAN Qiangling. Experimental study of spontaneous ignition and subsequent flame propagation of high-pressure hydrogen release [D]. Heifei: University of Science and Technology of China, 2016.
[11]	GOLUB V V, BAKLANOV D I, BAZHENOVA T V, et al. Shock-induced ignition of hydrogen gas during accidental or technical opening of high-pressure tanks [J]. Journal of Loss Prevention in the Process Industries, 2007, 20(4): 439–446.
[12]	GOLUB V V, BAKLANOV D I, GOLOVASTOV S V, et al. Mechanisms of high-pressure hydrogen gas self-ignition in tubes [J]. Journal of Loss Prevention in the Process Industries, 2008, 21(2): 185–198. DOI: 10.1016/j.jlp.2007.06.012.
[13]	KIM Y R, LEE H J, KIM S, et al. A flow visualization study on self-ignition of high pressure hydrogen gas released into a tube [J]. Proceedings of the Combustion Institute, 2013, 34(2): 2057–2064. DOI: 10.1016/j.proci.2012.07.020.
[14]	GRUNE J, SEMPERT K, KUZNETSOV M, et al. Experimental investigation of flame and pressure dynamics after spontaneous ignition in tube geometry [J]. International Journal of Hydrogen Energy, 2014, 39(35): 20396–20403. DOI: 10.1016/j.ijhydene.2014.05.046.
[15]	KANEKO W, SHII K. Effects of diaphragm rupturing conditions on self-ignition of high-pressure hydrogen [J]. International Journal of Hydrogen Energy, 2016, 41(25): 10969–10975. DOI: 10.1016/j.ijhydene.2016.04.211.
[16]	KANEKO W, ISHII K. An experimental study on the mechanism of self-ignition of high-pressure hydrogen [J]. International Journal of Hydrogen Energy, 2017, 42(11): 7374–7379. DOI: 10.1016/j.ijhydene.2016.06.046.
[17]	XU B P, WEN J X. Numerical study of spontaneous ignition in pressurized hydrogen release through a length of tube with local contraction [J]. International Journal of Hydrogen Energy, 2012, 37(22): 17571–17579. DOI: 10.1016/j.ijhydene.2012.04.150.
[18]	陈强. 激波管流动的理论和实验技术[D]. 合肥: 中国科学技术大学, 1979. CHEN Qiangling. Theory and experimental techniques of shock tube flows [D]. Heifei: University of Science and Technology of China, 1979.
[19]	SPENCE D A, WOODS B A. A review of theoretical treatments of shock-tube attenuation [J]. Journal of Fluid Mechanics, 2006, 19(2): 161–174.
[20]	LEE B J, JEUNG I S. Numerical study of spontaneous ignition of pressurized hydrogen released by the failure of a rupture disk into a tube [J]. International Journal of Hydrogen Energy, 2009, 34(20): 8763–8769. DOI: 10.1016/j.ijhydene.2009.08.034.
[21]	SHEN X, SUN J. Numerical simulation on the spontaneous ignition of leaking high pressure hydrogen from terminal unit [J]. Physics Procedia, 2012, 33(6): 1833–1841.
[22]	STUDER E, JAMOIS D, JALLAIS S, et al. Properties of large-scale methane/hydrogen jet fires [J]. International Journal of Hydrogen Energy, 2009, 34(23): 9611–9619. DOI: 10.1016/j.ijhydene.2009.09.024.