开放空间高压氢气射流中点火爆炸的实验研究

马梦飞; 於星; 张爱凤; 张佳庆; 祝现礼; 王昌建

doi:10.11883/bzycj-2023-0037

开放空间高压氢气射流中点火爆炸的实验研究

doi: 10.11883/bzycj-2023-0037

马梦飞^{1, 2,},
於星^{1, 2},
张爱凤^{1, 2},
张佳庆^{3, 4},
祝现礼^{3, 4},
王昌建^{1, 5, ,}

1.
合肥工业大学土木与水利工程学院，安徽合肥 230009
2.
安徽国际氢能安全联合研究中心，安徽合肥 230009
3.
国家电网安徽省电力科学研究院，安徽合肥 230088
4.
电力火灾与安全防护安徽省重点实验室，安徽合肥 230088
5.
安全关键工业测控教育部工程研究中心，安徽合肥 230009

基金项目: 国家重点研发计划项目（2021YFB4001004）；国家电网公司总部科技项目（5419-202219075A-1-1-ZN）

详细信息

作者简介:
马梦飞（1998－　），男，硕士研究生，1625388658@qq.com

通讯作者:
王昌建（1975－　），男，博士，教授，博士生导师，chjwang@hfut.edu.cn

中图分类号: O383.1
计量
- 文章访问数: 190
- HTML全文浏览量: 78
- PDF下载量: 65
- 被引次数: 3
出版历程
- 收稿日期: 2023-02-10
- 修回日期: 2024-01-03
- 网络出版日期: 2024-06-18
- 刊出日期: 2024-06-18

An experimental study on ignition and explosion of high-pressure hydrogen jet in open space

MA Mengfei^{1, 2
,},
YU Xing^{1, 2},
ZHANG Aifeng^{1, 2},
ZHANG Jiaqing^{3, 4},
ZHU Xianli^{3, 4},
WANG Changjian^{1, 5
, ,}

1.
School of Civil Engineering, Hefei University of Technology, Hefei 230009, Anhui, China
2.
Anhui International Joint Research Center on Hydrogen Safety, Hefei 230009, Anhui, China
3.
State Grid Anhui Electric Power Research Institute, Hefei 230088, Anhui, China
4.
Anhui Province Key Laboratory for Electric Fire and Safety Protection, Hefei 230088, Anhui, China
5.
Engineering Research Center of Safety Critical Industrial Measurement and Control Technology, Ministry of Education, Hefei 230009, Anhui, China

摘要

摘要: 采用高速相机和压力传感器，对开放空间中稳态射流氢气点火爆炸初期的火焰行为和超压变化规律进行了实验研究。结果表明：在点火爆炸初期，火焰在点火电极处以球形向外扩散；爆炸后4~6 ms，火焰前锋达到最大位移，之后逐渐熄灭，最后形成射流火焰。火焰前锋位移主要受喷嘴直径影响，并随喷嘴直径的增大而增大。火焰宽度的变化规律与火焰前锋位移基本相似。整个爆炸过程仅出现1个超压峰值，正压维持时间约为1 ms。在同一点火距离处，峰值超压随氢气流量的增加而增大。在相同氢气流量下，峰值超压随点火距离的增大而减小。最大峰值超压与氢气流量成正比，与点火距离成反比。
- 氢气爆炸 /
- 火焰前锋位移 /
- 火焰宽度 /
- 最大峰值超压
Abstract: Experiments were carried out on the flame behaviors and overpressure evolutions at the initial stage of ignition and explosion of steady-state hydrogen jet in open space, while a high-speed camera and pressure transducers were employed to record the flame shape and overpressure. The results show that at the early stage of ignition and explosion, the flame propagated outward from the ignition electrode with a spherical shape. After 4-6 ms, the flame front reached its maximum displacement, gradually extinguished, and finally formed a jet flame. The displacement of the flame front was mainly affected by the nozzle diameter and increased with the nozzle diameter. The variations of flame width were basically similar to those of the flame front displacement. The entire explosion process only experienced one overpressure peak, with a positive pressure maintained for approximately 1 ms. At the same ignition distance, the peak overpressure increased with hydrogen flow rate. At the same hydrogen flow rate, the peak overpressure decreased with increasing ignition distance. The maximum peak overpressure was directly proportional to the hydrogen flow rate and inversely proportional to the ignition distance.
- hydrogen explosion /
- flame front displacement /
- flame width /
- maximum peak overpressure

HTML全文

弹性波散射理论一直以来都是弹性动力学中重要研究课题之一，巧妙地运用了一些数学物理方法解析地求解了一些复杂的波动问题，对地震工程、岩土工程及地下结构工程等相关技术的研究与应用有着重要价值。对弹性波在全空间中传播时遇缺陷发生散射的研究已日趋完善，相关的研究已扩展到半空间、四分之一空间等更复杂的情况，更多的界面模型被涉及^[1-8]。双相介质半空间中的缺陷对SH波的散射问题是在近几年才备受科研人员的重视，而且他们基本上是讨论入射SH波与缺陷在同一个介质中的情况^[2-6]，对于缺陷与入射波处于不同介质中的研究还非常少。本文中采用Green函数法、复变函数法、保角映射法、“镜像”法、极坐标移动技术以及“契合”的思想解析地求解双相介质弹性半空间内椭圆弹性夹杂对透射SH波的散射问题。并通过具体的算例得出在不同的入射角、SH波频率和介质性质的情况下椭圆夹杂周边环向动应力集中分布情况，以期获得一些具体的理论结果为相关科研及工程实际应用提供参考。

1. 问题模型

如图 1所示，由介质Ⅰ和介质Ⅱ组成的双相介质半空间内有一个椭圆夹杂，椭圆夹杂为介质Ⅲ，这3种介质均为连续、均匀、各向同性的弹性介质。椭圆夹杂长半轴和短半轴长度分别为a和b，夹杂中心到垂直边界B_V的距离为h，到水平边界B_H距离为d，建立如图 1所示直角坐标系xOy和x″O″y″。SH波从介质Ⅱ中入射，遇垂直边界B_V产生透射SH波进入介质Ⅰ中。主要研究在不同参数条件下介质Ⅰ中椭圆夹杂周边在透射SH波作用下产生的环向动应力集中情况，并对结果进行分析，突出反映透射SH波的危害性。

图 1 理论模型

Figure 1. Theoretical model

下载: 全尺寸图片幻灯片

2. 控制方程

二维平面SH波位移函数W(x, y)与时间的依赖关系为e^-iωt, 满足控制方程：

$\frac{{{\partial ^2}W}}{{\partial {x^2}}} + \frac{{{\partial ^2}W}}{{\partial {y^2}}} + {k^2}W = 0$

(1)

式中: $k = \frac{\omega }{{{c_{\rm{s}}}}}$ ， ${c_{\rm{s}}} = \sqrt {\frac{\mu }{\rho }}$ ，k为波数，ω为位移函数圆频率，c_s为介质的剪切波速，μ为介质的剪切模量，ρ为介质密度。引入复变量z=x+iy，z=x-iy，位移函数W(z, z)在复平面(z, z)上控制方程的表达形式为：

$\frac{{{\partial ^2}W}}{{\partial z\partial \bar z}} + \frac{1}{4}{k^2}W = 0$

(2)

引入保角映射函数：

$z = \omega \left( \eta \right) = R\left( {\eta + \frac{m}{\eta }} \right)$

(3)

式中: $\eta = \mathit{R}{\mathit{e}^{{\rm{i}}\theta }}, R = \frac{{a + b}}{2}, m = \frac{{a-b}}{{a + b}}$ ，a、b分别为椭圆长半轴和短半轴长。

控制方程在映射平面(η, η)上可以表示为：

$\frac{1}{{\omega '\left( \eta \right)\overline {\omega '\left( \eta \right)} }}\frac{{{\partial ^2}W}}{{\partial \eta \partial \bar \eta }} + \frac{1}{4}{k^2}W = 0$

(4)

与式(4)相应的应力表达式为：

${\tau _{rz}} = \frac{\mu }{{R\left| {\omega '\left( \eta \right)} \right|}}\left( {\eta \frac{{\partial W}}{{\partial \eta }} + \bar \eta \frac{{\partial W}}{{\partial \bar \eta }}} \right),\;\;\;\;\;\;\;\;{\tau _{\theta z}} = \frac{{{\rm{i}}\mu }}{{R\left| {\omega '\left( \eta \right)} \right|}}\left( {\eta \frac{{\partial W}}{{\partial \eta }} - \bar \eta \frac{{\partial W}}{{\partial \bar \eta }}} \right)$

(5)

3. Green函数

本文中在求解实际问题之前先构造问题的Green函数，其中Green函数Ⅰ具体为含椭圆弹性夹杂的四分之一空间在垂直边界上任一位置处的出平面点源荷载作用下产生的位移场，四分之一空间为介质Ⅰ，椭圆夹杂为介质Ⅲ，如图 2所示。采用“虚设点源”法构造Green函数入射波位移场表达式如下：

${G^{\left( {{\rm{in}}} \right)}}\left( {\eta ,\bar \eta } \right) = \frac{{\rm{i}}}{{2{\mu _1}}}H_0^{\left( 1 \right)}\left( {{k_1}\left| {\omega \left( \eta \right) - \omega \left( {{\eta _0}} \right)} \right|} \right) + \frac{{\rm{i}}}{{2{\mu _1}}}H_0^{\left( 1 \right)}\left( {{k_1}\left| {\omega \left( \eta \right) - \overline {\omega \left( {{\eta _0}} \right)} + 2{\rm{i}}d} \right|} \right)$

(6)

图 2 虚设点源模型

Figure 2. Model of dummy point source loads

下载: 全尺寸图片幻灯片

式中:i为虚数单位，H₀⁽¹⁾为零阶第一类Hankel函数，ω(η₀)为点源荷载的位置矢量。构造四分之一空间内椭圆弹性夹杂产生的散射波位移场表达式^[4]为：

${G^{\left( {\rm{s}} \right)}}\left( {\eta ,\bar \eta } \right) = \sum\limits_{n = - \infty }^\infty {{A_n}\sum\limits_{j = 1}^4 {S_n^{\left( j \right)}} }$

(7)

式中：

$S_n^{\left( 1 \right)} = H_n^{\left( 1 \right)}\left( {{k_1}\left| {\omega \left( \eta \right)} \right|} \right){\left[ {\frac{{\omega \left( \eta \right)}}{{\left| {\omega \left( \eta \right)} \right|}}} \right]^n}，$

$S_n^{\left( 2 \right)} = H_n^{\left( 1 \right)}\left( {{k_1}\left| {\omega \left( \eta \right) - 2h} \right|} \right){\left[ {\frac{{\omega \left( \eta \right) - 2h}}{{\left| {\omega \left( \eta \right) - 2h} \right|}}} \right]^{ - n}}，$

$S_n^{\left( 3 \right)} = {\left( { - 1} \right)^n}H_n^{\left( 1 \right)}\left( {{k_1}\left| {\omega \left( \eta \right) - 2{\rm{i}}d} \right|} \right){\left[ {\frac{{\omega \left( \eta \right) - 2{\rm{i}}d}}{{\left| {\omega \left( \eta \right) - 2{\rm{i}}d} \right|}}} \right]^{ - n}}，$

$S_n^{\left( 4 \right)} = {\left( { - 1} \right)^n}H_n^{\left( 1 \right)}\left( {{k_1}\left| {\omega \left( \eta \right) - 2{\rm{i}}d - 2h} \right|} \right){\left[ {\frac{{\omega \left( \eta \right) - 2{\rm{i}}d - 2h}}{{\left| {\omega \left( \eta \right) - 2{\rm{i}}d - 2h} \right|}}} \right]^n}。$

椭圆弹性夹杂内部驻波表达式如下：

${G^{\left( {\rm{t}} \right)}}\left( {\eta ,\bar \eta } \right) = \sum\limits_{{\rm{n}} = - \infty }^\infty {{{\rm{B}}_n}{J_n}\left( {{k_3}\left| {\omega \left( \eta \right)} \right|} \right){{\left[ {\frac{{\omega \left( \eta \right)}}{{\left| {\omega \left( \eta \right)} \right|}}} \right]}^n}}$

(8)

然后依据椭圆夹杂边界上应力和位移连续条件可建立如下方程组：

$\left\{ \begin{array}{l} {G^{\left( {{\rm{in}}} \right)}}\left( {\eta ,\bar \eta } \right) + {G^{\left( {\rm{s}} \right)}}\left( {\eta ,\bar \eta } \right) = {G^{\left( {\rm{t}} \right)}}\left( {\eta ,\bar \eta } \right)\\ \tau _{rz}^{\left( {{\rm{in}}} \right)}\left( {\eta ,\bar \eta } \right) + \tau _{rz}^{\left( {\rm{s}} \right)}\left( {\eta ,\bar \eta } \right) = \tau _{rz}^{\left( {\rm{t}} \right)}\left( {\eta ,\bar \eta } \right) \end{array} \right.$

(9)

在方程两边同时乘以e^-imθ，然后在(-π, π)区间内积分并进行有限项截断求出系数A_n和B_n，本文中取n=8。Green函数Ⅰ为:

${G_1}\left( {\eta ,\bar \eta } \right) = G_{\rm{I}}^{\left( {{\rm{in}}} \right)}\left( {\eta ,\bar \eta } \right) + G_{\rm{I}}^{\left( {\rm{s}} \right)}\left( {\eta ,\bar \eta } \right)$

(10)

Green函数Ⅱ为四分之一空间中无椭圆夹杂时的位移场，即:

${G_2}\left( {\eta ,\bar \eta } \right) = G_{{\rm{II}}}^{\left( {{\rm{in}}} \right)}\left( {\eta ,\bar \eta } \right)$

(11)

4. 入射波、反射波、透射波和散射波

采用“镜像”法将半空间双相介质问题转化为全空间双相介质问题，入射波、反射波和透射波可分别表示为：

$\begin{array}{l} {W^{\left( {{\rm{in}}} \right)}} = {W_0}\exp \left\{ {\frac{{{\rm{i}}{k_2}}}{2}\left[ {\left( {\omega \left( \eta \right) + h} \right){{\rm{e}}^{ - {\rm{i}}\left( {{\rm{ \mathsf{ π} }} - {\alpha _0}} \right)}} - \left( {\overline {\omega \left( \eta \right)} + h} \right){{\rm{e}}^{ - {\rm{i}}{\alpha _0}}}} \right]} \right\} + \\ {W_0}\exp \left\{ {\frac{{{\rm{i}}{k_2}}}{2}\left[ {\left( {\omega \left( \eta \right) + h - 2{\rm{i}}d} \right){{\rm{e}}^{ - {\rm{i}}\left( {{\rm{ \mathsf{ π} }} - {\alpha _0}} \right)}} - \left( {\overline {\omega \left( \eta \right)} + h + 2{\rm{i}}d} \right){{\rm{e}}^{ - {\rm{i}}{\alpha _0}}}} \right]} \right\} \end{array}$

(12)

$\begin{array}{l} {W^{\left( {{\rm{re}}} \right)}} = {W_2}\exp \left\{ {\frac{{{\rm{i}}{k_2}}}{2}\left[ {\left( {\omega \left( \eta \right) + h} \right){{\rm{e}}^{{\rm{i}}\left( {{\rm{ \mathsf{ π} }} - {\alpha _0}} \right)}} - \left( {\overline {\omega \left( \eta \right)} + h} \right){{\rm{e}}^{{\rm{i}}{\alpha _0}}}} \right]} \right\} + \\ {W_2}\exp \left\{ {\frac{{{\rm{i}}{k_2}}}{2}\left[ {\left( {\omega \left( \eta \right) + h - 2{\rm{i}}d} \right){{\rm{e}}^{ - {\rm{i}}\left( {{\rm{ \mathsf{ π} }} - {\alpha _0}} \right)}} - \left( {\overline {\omega \left( \eta \right)} + h + 2{\rm{i}}d} \right){{\rm{e}}^{ - {\rm{i}}{\alpha _0}}}} \right]} \right\} \end{array}$

(13)

$\begin{array}{l} {W^{\left( {{\rm{tr}}} \right)}} = {W_4}\exp \left\{ {\frac{{{\rm{i}}{k_1}}}{2}\left[ {\left( {\omega \left( \eta \right) + h} \right){{\rm{e}}^{ - {\rm{i}}\left( {{\rm{ \mathsf{ π} }} - {\alpha _4}} \right)}} - \left( {\overline {\omega \left( \eta \right)} + h} \right){{\rm{e}}^{ - {\rm{i}}{\alpha _4}}}} \right]} \right\} + \\ {W_4}\exp \left\{ {\frac{{{\rm{i}}{k_1}}}{2}\left[ {\left( {\omega \left( \eta \right) + h - 2{\rm{i}}d} \right){{\rm{e}}^{{\rm{i}}\left( {{\rm{ \mathsf{ π} }} - {\alpha _4}} \right)}} - \left( {\overline {\omega \left( \eta \right)} + h + 2{\rm{i}}d} \right){{\rm{e}}^{{\rm{i}}{\alpha _4}}}} \right]} \right\} \end{array}$

(14)

运用Snell定律可以得到入射波、反射波和透射波位移幅值之间的关系式如下：

${W_2} = {W_0}\frac{{\cos {\alpha _0} - \frac{{{c_2}}}{{{c_1}}}\frac{{{\mu _1}}}{{{\mu _2}}}\cos {\alpha _4}}}{{\cos {\alpha _0} + \frac{{{c_2}}}{{{c_1}}}\frac{{{\mu _1}}}{{{\mu _2}}}\cos {\alpha _4}}},\;\;\;\;\;\;\;\;\;\;{W_4} = {W_0}\frac{{2\cos {\alpha _0}}}{{\cos {\alpha _0} + \frac{{{c_2}}}{{{c_1}}}\frac{{{\mu _1}}}{{{\mu _2}}}\cos {\alpha _4}}}$

(15)

由于入射波的作用，同样在介质Ⅰ中会产生散射波，在椭圆夹杂内部会产生驻波，其过程及相关系数的求解参考构造Green函数的过程，在此不再叙述。

5. 定解积分方程组

如图 3所示，本文中将双相介质沿垂直界面剖开，在剖分面上作用连续分布的入平面和出平面点源荷载，分别表示为f₁(r″₀, θ″₀)和f₂(r″₀, θ″₀)，在垂直边界上满足应力连续条件，即:

$\tau _{\theta ''z''}^{\left( {\rm{I}} \right)}\sin {{\theta ''}_0} + {f_1}\left( {{{r''}_0} + {{\theta ''}_0}} \right) = \tau _{\theta ''z''}^{\left( {{\rm{II}}} \right)}\sin {{\theta ''}_0} + {f_2}\left( {{{r''}_0} + {{\theta ''}_0}} \right)$

(16)

图 3 剖分面模型

Figure 3. Model of cutaway interface

下载: 全尺寸图片幻灯片

介质Ⅰ中存在透射波和散射波，即:

${W^{\left( {\rm{I}} \right)}} = {W^{\left( {{\rm{tr}}} \right)}} + {W^{\left( {\rm{s}} \right)}},\;\;\;\;\;\;\;\tau _{\theta z}^{\left( {\rm{I}} \right)} = \tau _{\theta z}^{\left( {{\rm{tr}}} \right)}$

(17)

介质Ⅱ中存在入射波和反射波，即:

${W^{\left( {{\rm{II}}} \right)}} = {W^{\left( {{\rm{in}}} \right)}} + {W^{\left( {{\rm{re}}} \right)}},\;\;\;\;\;\;\;\tau _{\theta z}^{\left( {{\rm{II}}} \right)} = \tau _{\theta z}^{\left( {{\rm{in}}} \right)} + \tau _{\theta z}^{\left( {{\rm{re}}} \right)}$

(18)

然后利用“契合”的思想将介质Ⅰ和介质Ⅱ“契合”在一起，在垂直界面处τ_θz⁽ⁱⁿ⁾+τ_θz^(re)=τ_θz^(tr)，于是可得f₁(r″₀, θ″₀)=f₂(r″₀, θ″₀)。根据前面构造的Green函数，通过积分的方法可以得到外力系f₁(r″₀, θ″₀)与f₂(r″₀, θ″₀)相应的位移表达式，即:

${W^{\left( {{f_1}} \right)}}\left( {r'',\theta ''} \right) = \int_0^h {{f_1}\left( {{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right){G_1}\left( {r'',\theta '',{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right){\rm{d}}{{r''}_0}} + \\ \int_0^\infty {{f_1}\left( {{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right){G_1}\left( {r'',\theta '',{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right){\rm{d}}{{r''}_0}}$

(19)

${W^{\left( {{f_2}} \right)}}\left( {r'',\theta ''} \right) = - \int_0^h {{f_2}\left( {{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right){G_2}\left( {r'',\theta '',{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right){\rm{d}}{{r''}_0}} - \\ \int_0^\infty {{f_2}\left( {{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right){G_2}\left( {r'',\theta '',{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right){\rm{d}}{{r''}_0}}$

(20)

利用垂直边界处的关系：W^(Ⅰ)+W^(f₁)=W^(Ⅱ)+W^(f₂)，W⁽ⁱⁿ⁾+W^(re)=W^(tr)可得到如下定解积分方程组：

$\begin{array}{*{20}{c}} {\int_0^h {{f_1}\left( {{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right)\left[ {{G_1}\left( {r'',\frac{{\rm{ \mathsf{ π} }}}{2},{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right) + {G_2}\left( {r'',\frac{{\rm{ \mathsf{ π} }}}{2},{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right)} \right]{\rm{d}}{{r''}_0}} + }\\ {\int_0^\infty {{f_1}\left( {{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right)\left[ {{G_1}\left( {r'',\frac{{\rm{ \mathsf{ π} }}}{2},{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right) + {G_2}\left( {r'',\frac{{\rm{ \mathsf{ π} }}}{2},{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right)} \right]{\rm{d}}{{r''}_0}} = \\ - {W^{\left( {\rm{s}} \right)}}\left| {_{{{\theta ''}_0} = \frac{{\rm{ \mathsf{ π} }}}{2}}} \right.} \end{array}$

(21)

$\begin{array}{*{20}{c}} {\int_0^h {{f_1}\left( {{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right)\left[ {{G_1}\left( {r'',\frac{{\rm{ \mathsf{ π} }}}{2},{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right) + {G_2}\left( {r'',\frac{{\rm{ \mathsf{ π} }}}{2},{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right)} \right]{\rm{d}}{{r''}_0}} + }\\ {\int_0^\infty {{f_1}\left( {{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right)\left[ {{G_1}\left( {r'',\frac{{\rm{ \mathsf{ π} }}}{2},{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right) + {G_2}\left( {r'',\frac{{\rm{ \mathsf{ π} }}}{2},{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right)} \right]{\rm{d}}{{r''}_0}} = \\ - {W^{\left( {\rm{s}} \right)}}\left| {_{{{\theta ''}_0} = - \frac{{\rm{ \mathsf{ π} }}}{2}}} \right.} \end{array}$

(22)

式中:G₁与G₂分别为之前构造的在介质Ⅰ和介质Ⅱ中的Green函数，利用散射波的衰减性并进行有限项截断求得未知附加力系。

6. 动应力集中因子

本文中主要是运用弹性波的散射理论来研究弹性波的绕射与动应力集中问题。环向动应力集中因子可写成如下形式：

$\tau _{\theta z}^ * = \left| {\frac{\tau }{{{\tau _0}}}} \right|$

(23)

式中:τ₀=μ₂k₂W₄为半空间透射SH波应力的最大幅值，τ为椭圆夹杂周边环向总应力。

在椭圆夹杂周边环向应力表达式：

${\tau _{\theta z}} = \tau _{\theta z}^{\left( {{\rm{tr}}} \right)} + \tau _{\theta z}^{\left( {\rm{s}} \right)} + \int_0^h {{f_1}\left( {{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right)\frac{{{\mu _1}}}{r}\frac{{\partial {G_1}\left( {r,\theta ,{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right)}}{{\partial \theta }}{\rm{d}}{{r''}_0}} + \\ \int_0^\infty {{f_1}\left( {{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right)\frac{{{\mu _1}}}{r}\frac{{\partial {G_1}\left( {r,\theta ,{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right)}}{{\partial \theta }}{\rm{d}}{{r''}_0}}$

(24)

7. 算例及结果分析

通过引入量纲一参数并对其赋值得到了一些具体算例的结果，给出了椭圆夹杂周边环向动应力集中因子在不同参数情况下的分布情况。取量纲一参数 ${\mu _{21}} = \frac{{{\mu _2}}}{{{\mu _1}}},{\mu _{31}} = \frac{{\mu 3}}{{{\mu _1}}},{k_{21}} = \frac{{{k_2}}}{{{k_1}}},{k_{31}} = \frac{{k3}}{{{k_1}}}$ ；透射波量纲一波数为k₁b；椭圆夹杂位置坐标量纲一参数为 $\frac{d}{b}$ 和 $\frac{h}{b}$ 。

图 4所示为椭圆夹杂周边环向动应力集中因子τ_θz^*随不同入射角α₀的分布情况。此处取k₂₁=1.0，k₃₁=1.0，μ₂₁=2.0，μ₃₁=0.5， $\frac{a}{b}$ =1.25, $\frac{h}{b} = \frac{d}{b} = 12.0$ 。可以看出，当SH波水平向入射时，椭圆夹杂周边环向动应力集中因子明显大于其他入射角度时的相应值，且当入射角α₀=45°时的相应值最小，α₀=30°和α₀=60°时的相应值较为接近且大小处于居中位置。当透射波波数k₁b=0.1时，即在“准静态”情况下，SH波水平入射时椭圆夹杂周边环向动应力集中因子最大值是入射角α₀=45°时相应最大值的1.7倍，比其他透射波波数情况下相应值大。由此可以得出，当SH波水平入射产生透射波时，椭圆夹杂周边环向动应力集中情况最严重，且不同的透射波波数对动应力集中的分布有一定的影响。

图 4 SH波以不同角度入射时动应力集中因子的分布情况

Figure 4. Distribution of dynamic stress concentration factor with different incident angles disturbed by SH wave

下载: 全尺寸图片幻灯片

图 5所示为SH波水平入射产生透射波时，椭圆夹杂周边环向动应力集中因子τ_θz^*随不同波数比k₂₁的分布情况。其中量纲一参数k₃₁=1.2，μ₃₁=1.0，μ₂₁=0.2， $\frac{a}{b}$ =1.25， $\frac{h}{b}$ =2.0， $\frac{d}{b}$ =12.0。可以看出，当透射波波数k₁b=0.1时，即在“准静态”情况下，波数比k₂₁对夹杂周边环向动应力集中因子的分布几乎没有影响，在其他透射波波数情况下，夹杂周边环向动应力集中因子最大值随k₂₁增大而略微减小，整体上变化并不大。当透射波波数发生变化时，夹杂周边环向动应力集中因子最大值的位置发生了明显变化。由此可以得出，垂直界面右侧介质的性质对椭圆夹杂周边动应力集中因子的分布影响较小，但是透射波频率对椭圆夹杂周边动应力集中因子最大值的位置影响较大。

图 5 SH波水平向入射时动应力集中因子随k₂₁的分布情况

Figure 5. Distribution of dynamic stress concentration factor around the elliptic inclusion edge with k₂₁ disturbed by SH wave horizontally

下载: 全尺寸图片幻灯片

图 6所示为SH波水平向入射产生透射波时，椭圆弹性夹杂周边环向动应力集中因子τ_θz^*随波数比k₃₁变化的分布情况。其中量纲一参数k₂₁=1.0，μ₂₁=2.0，μ₃₁=1.0， $\frac{a}{b}$ =1.25， $\frac{h}{b} = \frac{d}{b} = 12.0$ 。可以看出，当透射波波数k₁b=0.1时，即在“准静态”情况下，波数比k₃₁对夹杂周边环向动应力集中因子的分布影响较小。在中频和高频透射波情况下，夹杂周边环向动应力集中因子变化非常明显。当k₁b=0.5，k₃₁=4.0时，环向动应力集中因子τ_θz^*达到了极值|τ_θz^*|=7.63，比k₁b=0.5，k₃₁=0.5时的相应极值增大了3.06倍。由此可得知，在中频和高频透射波情况下，波数比k₃₁对椭圆夹杂周边环向动应力集中因子影响非常大，且极值的位置会随透射波频率的变化而发生改变。

图 6 SH波水平向入射时动应力集中因子随k₃₁的分布情况

Figure 6. Distribution of dynamic stress concentration factor around the elliptic inclusion edge with k₃₁ disturbed by SH wave horizontally

下载: 全尺寸图片幻灯片

图 7所示为SH波水平入射产生透射波时，椭圆弹性夹杂周边在θ=90°处环向动应力集中因子τ_θz^*随 $\frac{h}{b}$ 变化的分布情况。其中量纲一参k₂₁=1.0，k₃₁=1.0，μ₂₁=2.0，μ₃₁=2.0， $\frac{a}{b}$ =1.25， $\frac{d}{b}$ =12.0。可以看出，椭圆夹杂周边在θ=90°处环向动应力集中因子τ_θz^*随 $\frac{h}{b}$ 的增大而呈周期性变化，透射波频率越大，周期越小，且振幅呈衰减趋势，衰减到一定程度后逐渐趋于稳定。当透射波波数k₁b=0.1时，即在“准静态”情况下接近为一条直线。可以得知，椭圆弹性夹杂距离垂直边界一定距离后，该距离对夹杂周边环向动应力集中因子的分布影响可以忽略，且这个距离值随透射波频率增大而增大。

图 7 SH波水平向入射时动应力集中因子随h/b的分布情况

Figure 7. Distribution of dynamic stress concentration factor around the elliptic inclusion edge with h/b disturbed by SH wave horizontally

下载: 全尺寸图片幻灯片

8. 结论

SH波入射的角度、介质的性质、透射波的频率以及椭圆夹杂与垂直边界的距离对椭圆夹杂周边环向动应力集中因子的分布均有不同程度的影响。当SH波水平入射产生透射波时，椭圆夹杂周边环向动应力集中程度较大。在一定范围内，垂直界面右侧的介质性质对动应力集中因子的极值影响较小，但在不同透射波频率下该极值的位置会发生明显变化。在中频和高频透射波情况下，椭圆夹杂的性质对其周边动应力集中的分布影响较大，且在一定条件下动应力集中程度会非常严重。当椭圆夹杂距离垂直边界一定距离后，该距离的影响可以忽略不计。总之，在实际工程中应该重视双相介质中透射波对结构可能造成的不利影响。

图 1 实验装置示意图

Figure 1. Layout of the experimental apparatus

下载: 全尺寸图片幻灯片

图 2 氢气火焰图像处理过程（工况18，p₀=0.5 MPa, d=4 mm）

Figure 2. Image processing of hydrogen flame (case 18, p₀=0.5 MPa, d=4 mm)

下载: 全尺寸图片幻灯片

图 3 初始阶段火焰演变过程（工况18，p₀=0.5 MPa, d=4 mm）

Figure 3. Evolutions of the flame at the initial stage (case 18, p₀=0.5 MPa, d=4 mm)

下载: 全尺寸图片幻灯片

图 4 平均体积分数随距离的变化（工况18，p₀=0.5 MPa, d=4 mm）

Figure 4. Variation of mean concentration with distance (case 18, p₀=0.5 MPa, d=4 mm)

下载: 全尺寸图片幻灯片

图 5 3 ms时氢气火焰前锋的位移和宽度

Figure 5. Displacement and width of hydrogen flame front at 3 ms

下载: 全尺寸图片幻灯片

图 6 火焰前锋位移和宽度随点火距离的变化

Figure 6. Variations of flame front displacement and width with ignition distance

下载: 全尺寸图片幻灯片

图 7 超压随时间变化趋势（工况17, p₀=0.5 MPa, d=4 mm）

Figure 7. Evolution of overpressure with time (case 17, p₀=0.5 MPa, d=4 mm)

下载: 全尺寸图片幻灯片

图 8 不同流量下峰值超压随点火距离的变化

Figure 8. Evolution of peak overpressure with ignition distance at different flow rates

下载: 全尺寸图片幻灯片

图 9 不同点火距离下最大峰值超压随流量的变化

Figure 9. Relation between the flow rate and maximum peak overpressure at different ignition distances

下载: 全尺寸图片幻灯片

图 10 最大峰值超压与点火距离的关系

Figure 10. Relationship between the maximum peak overpressure and ignition distance

下载: 全尺寸图片幻灯片

表 1 氢气点火实验参数

Table 1. Experimental parameters of hydrogen ignition

工况	d/mm	点火距离/cm	p₀/MPa	m/（g·s⁻¹）	最大火焰前锋位移/mm	$\bar \varphi$
1	2	16	0.3	0.78	229.6	0.314
2	2	20	0.3	0.78	199.5	0.254
3	2	30	0.3	0.78	73.5	0.171
4	2	37	0.3	0.78
5	2	16	0.5	1.17	339.5	0.325
6	2	20	0.5	1.17	267.4	0.269
7	2	30	0.5	1.17	147.0	0.188
8	2	37	0.5	1.17	145.6	0.156
9	2	16	0.7	1.56	387.8	0.322
10	2	20	0.7	1.56	338.8	0.282
11	2	30	0.7	1.56	229.6	0.196
12	2	37	0.7	1.56	88.2	0.164
13	3	16	0.5	2.64	350.0	0.449
14	3	20	0.5	2.64	338.1	0.377
15	3	30	0.5	2.64	348.6	0.269
16	3	37	0.5	2.64	264.6	0.224
17	4	16	0.5	4.69	576.5	0.574
18	4	20	0.5	4.69	592.2	0.471
19	4	30	0.5	4.69	321.3	0.342
20	4	37	0.5	4.69	342.3	0.287

下载: 导出CSV

参考文献(25)

[1]	LI H W, TANG Z S, LI J L, et al. Investigation of vented hydrogen-air deflagrations in a congested vessel [J]. Process Safety and Environmental Protection, 2019, 129: 196–201. DOI: 10.1016/j.psep.2019.07.009.
[2]	郝腾腾, 王昌建, 颜王吉, 等. 氢气泄爆作用下结构动力响应特性研究 [J]. 爆炸与冲击, 2020, 40(6): 065401. DOI: 10.11883/bzycj-2019-0412. HAO T T, WANG C J, YAN W J, et al. Structural dynamical characteristics induced by vented hydrogen explosion [J]. Explosion and Shock Waves, 2020, 40(6): 065401. DOI: 10.11883/bzycj-2019-0412.
[3]	RUI S C, WANG C J, GUO S S, et al. Hydrogen-air explosion with concentration gradients in a cubic enclosure [J]. Process Safety and Environmental Protection, 2021, 151: 141–150. DOI: 10.1016/j.psep.2021.05.003.
[4]	WANG L Q, LI T, MA H H. Explosion behaviors of hydrogen-nitrous oxide mixtures at reduced initial pressures [J]. Process Safety and Environmental Protection, 2021, 153: 11–18. DOI: 10.1016/j.psep.2021.07.010.
[5]	曹勇, 郭进, 胡坤伦, 等. 点火位置对氢气-空气预混气体泄爆过程的影响 [J]. 爆炸与冲击, 2016, 36(6): 847–852. DOI: 10.11883/1001-1455(2016)06-0847-06. CAO Y, GUO J, HU K L, et al. Effect of ignition locations on vented explosion of premixed hydrogen-air mixtures [J]. Explosion and Shock Waves, 2016, 36(6): 847–852. DOI: 10.11883/1001-1455(2016)06-0847-06.
[6]	余明高, 阳旭峰, 郑凯, 等. 障碍物对甲烷/氢气爆炸特性的影响 [J]. 爆炸与冲击, 2018, 38(1): 19–27. DOI: 10.11883/bzycj-2017-0172. YU M G, YANG X F, ZHENG K, et al. Effect of obstacles on explosion characteristics of methane/hydrogen [J]. Explosion and Shock Waves, 2018, 38(1): 19–27. DOI: 10.11883/bzycj-2017-0172.
[7]	陈昊, 郭进, 王金贵, 等. 破膜压力对氢气-甲烷-空气泄爆的影响 [J]. 爆炸与冲击, 2022, 42(11): 115401. DOI: 10.11883/bzycj-2021-0418. CHEN H, GUO J, WANG J G, et al. Effects of vent burst pressure on hydrogen-methane-air deflagration in a vented duct [J]. Explosion and Shock Waves, 2022, 42(11): 115401. DOI: 10.11883/bzycj-2021-0418.
[8]	SWAIN M R, FILOSO P A, SWAIN M N. An experimental investigation into the ignition of leaking hydrogen [J]. International Journal of Hydrogen Energy, 2007, 32(2): 287–295. DOI: 10.1016/j.ijhydene.2006.06.041.
[9]	PANDA P P, HECHT E S. Ignition and flame characteristics of cryogenic hydrogen releases [J]. International Journal of Hydrogen Energy, 2017, 42(1): 775–785. DOI: 10.1016/j.ijhydene.2016.08.051.
[10]	FRIEDRICH A, BREITUNG W, STERN G, et al. Ignition and heat radiation of cryogenic hydrogen jets [J]. International Journal of Hydrogen Energy, 2012, 37(22): 17589–17598. DOI: 10.1016/j.ijhydene.2012.07.070.
[11]	KOBAYASHI H, MUTO D, DAIMON Y, et al. Experimental study on cryo-compressed hydrogen ignition and flame [J]. International Journal of Hydrogen Energy, 2020, 45(7): 5098–5109. DOI: 10.1016/j.ijhydene.2019.12.091.
[12]	MOGI T, HORIGUCHI S. Experimental study on the hazards of high-pressure hydrogen jet diffusion flames [J]. Journal of Loss Prevention in the Process Industries, 2009, 22(1): 45–51. DOI: 10.1016/j.jlp.2008.08.006.
[13]	GUO J, WANG C J, LI M H, et al. Experimental study on spark ignition of non-premixed hydrogen jets [J]. Journal of Thermal Science, 2014, 23(4): 354–358. DOI: 10.1007/s11630-014-0717-3.
[14]	闫伟阳, 潘旭海, 汪志雷, 等. 高压氢气泄漏自燃形成喷射火的实验研究 [J]. 爆炸与冲击, 2019, 39(11): 115402. DOI: 10.11883/bzycj-2018-0394. YAN W Y, PAN X H, WANG Z L, et al. Experimental investigation on spontaneous combustion of high-pressure hydrogen leakage to form jet fire [J]. Explosion and Shock Waves, 2019, 39(11): 115401. DOI: 10.11883/bzycj-2018-0394.
[15]	AHMED S F, MASTORAKOS E. Spark ignition of lifted turbulent jet flames [J]. Combustion and Flame, 2006, 146(1/2): 215–231. DOI: 10.1016/j.combustflame.2006.03.007.
[16]	AHMED S F, BALACHANDRAN R, MASTORAKOS E. Measurements of ignition probability in turbulent non-premixed counterflow flames [J]. Proceedings of the Combustion Institute, 2007, 31(1): 1507–1513. DOI: 10.1016/j.proci.2006.07.089.
[17]	GRUNE J, SEMPERT K, KUZNETSOV M, et al. Experimental study of ignited unsteady hydrogen jets into air [J]. International Journal of Hydrogen Energy, 2011, 36(3): 2497–2504. DOI: 10.1016/j.ijhydene.2010.04.152.
[18]	ROYLE M, WILLOUGHBY D B. Consequences of catastrophic releases of ignited and unignited hydrogen jet releases [J]. International Journal of Hydrogen Energy, 2011, 36(3): 2688–2692. DOI: 10.1016/j.ijhydene.2010.03.141.
[19]	GRUNE J, SEMPERT K, KUZNETSOV M, et al. Experimental study of ignited unsteady hydrogen releases from a high pressure reservoir [J]. International Journal of Hydrogen Energy, 2014, 39(11): 6176–6183. DOI: 10.1016/j.ijhydene.2013.08.076.
[20]	BIRCH A D, BROWN D R, DODSON M G, et al. The structure and concentration decay of high pressure jets of natural gas [J]. Combustion Science and Technology, 1984, 36(5/6): 249–261. DOI: 10.1080/00102208408923739.
[21]	BIRCH A D, HUGHES D J, SWAFFIELD F. Velocity decay of high pressure jets [J]. Combustion Science and Technology, 1987, 52(1): 161–171. DOI: 10.1080/00102208708952575.
[22]	SCHEFER R W, HOUF W G, WILLIAMS T C, et al. Characterization of high-pressure, underexpanded hydrogen-jet flames [J]. International Journal of Hydrogen Energy, 2007, 32(12): 2081–2093. DOI: 10.1016/j.ijhydene.2006.08.037.
[23]	DUAN Q L, XIAO H H, GAO W, et al. Experimental study on spontaneous ignition and flame propagation of high-pressure hydrogen release via a tube into air [J]. Fuel, 2016, 181: 811–819. DOI: 10.1016/j.fuel.2016.05.066.
[24]	TAKENO K, OKABAYASHI K, KOUCHI A, et al. Dispersion and explosion field tests for 40 MPa pressurized hydrogen [J]. International Journal of Hydrogen Energy, 2007, 32(13): 2144–2153. DOI: 10.1016/j.ijhydene.2007.04.018.
[25]	KIKUKAWA S, YAMAGA F, MITSUHASHI H. Risk assessment of Hydrogen fueling stations for 70 MPa FCVs [J]. International Journal of Hydrogen Energy, 2008, 33(23): 7129–7136. DOI: 10.1016/j.ijhydene.2008.08.063.

施引文献

期刊类型引用(2)

1.	屈恩相，张景颢，王伟业，胡兴森，王浩宁，徐静怡，李聪. 镜像法求解SH波散射问题的应用综述. 安徽建筑. 2024(07): 115-117 . 百度学术
2.	屈恩相，齐辉，郭晶，杨杰，郑易. “分区”与“契合”思想求解半空间含凸起地形对SH波散射问题的研究进展. 黑龙江工业学院学报(综合版). 2022(04): 74-82 . 百度学术

其他类型引用(1)

资源附件(0)

访问统计

图(10) / 表(1)

计量

文章访问数: 190
HTML全文浏览量: 78
PDF下载量: 65
被引次数: 3

1. 问题模型
2. 控制方程
3. Green函数
4. 入射波、反射波、透射波和散射波
5. 定解积分方程组
6. 动应力集中因子
7. 算例及结果分析
8. 结论

开放空间高压氢气射流中点火爆炸的实验研究

doi: 10.11883/bzycj-2023-0037

作者简介: 马梦飞（1998－ ），男，硕士研究生，1625388658@qq.com

通讯作者: 王昌建（1975－ ），男，博士，教授，博士生导师，chjwang@hfut.edu.cn

计量

出版历程