Study on mixing induced by Richtmyer-Meshkov instability by using buoyancy-drag model
-
摘要: 采用浮阻力模型对激波管低压缩和激光加载高压缩情况下的Richtmyer-Meshkov不稳定性诱导混合现象进行了研究。通过与实验和理论分析结果进行比较发现:为了达到好的吻合, Richtmyer-Meshkov不稳定性情况下阻力系数的取值范围(2.0~5.36)比Rayleigh-Taylor不稳定性情况下的值(3.3~4.0)宽得多; 而在Richtmyer-Meshkov不稳定性情况下, 高压缩时阻力系数的不确定度(约为3.36)明显高于低压缩时的值(约为1.46), 模型的进一步完善还有待于更精确实验的验证。研究显示:指数律经验公式中指数随工况的不同而显著变化, 目前工程设计中采用指数律经验公式是粗糙的。
-
关键词:
- 流体力学 /
- Richtmyer-Meshkov不稳定性 /
- 浮阻力模型 /
- 扰动界面 /
- 压缩比
Abstract: The mixing induced by Richtmyer-Meshkov instability under lower and high compression ratios is studied by using buoyancy-drag model.It is found by the comparison between the experimental and theoretical analytic results that in order to achieve good agreement, the range of the drag coefficient value chosen in Richtmyer-Meshkov instability is much wider than that in Rayleigh-Taylor instability; while in Richtmyer-Meshkov instability the uncertainty of the drag coefficient under high compression is larger than that under lower compression.It is pointed that the further improvement of the model needs to be validated by more accurate experiment.Moreover the study reveals that the exponent in empirical expression varies with calculation conditions remarkably and the use of the empirical expression in current engineering design is rough. -
根据中国城市建设统计年鉴[1]的记载,截至2019年,我国可燃气体管道在城市地下空间的长度达76.7万公里。然而可燃气体管道时常因阀门失灵、管道破损、操作失误等原因而泄漏,一旦可燃气体在城市地下浅埋管沟内积聚,在遇到点火点时,极易发生严重的可燃气体爆炸事故,同时造成恶劣的社会影响。如2013年11月,山东青岛一条输油管道在地下暗渠中发生了严重的可燃气体爆炸事故,造成了62人死亡、136人受伤[2]。2014年7月,台湾高雄一条输送丙烯液体的地下管道发生了严重的蒸气云爆炸事故,造成了32人遇难、321人受伤[3]。因此,对城市地下浅埋管沟内可燃气体爆炸的荷载进行研究,以期减轻爆炸灾害效应,具有重要的经济价值和社会意义。
学者们对长直空间内可燃气体爆炸已开展了不同类型的试验研究。王东武等[4]、司荣军[5]、Ma等[6]分别在大尺寸的巷道、圆形管道和天然气管道内进行了可燃气体爆炸试验,研究爆炸冲击波在装置内的传播规律,得到了冲击波在受限空间内的超压峰值分布特征和衰减过程。在小尺寸可燃气体爆炸试验方面,Ciccarelli等[7]在边长7.6 cm、长3.6 m的管道中进行了不同障碍物拥塞率的多组天然气爆炸试验;Na’inna等[8]在直径为16.2 cm、长4.5 m的管道中进行了多组体积分数为10%的甲烷爆炸试验,发现了超压和火焰速度与障碍物尺寸符合平方关系;孙庆文[9]对一段尺寸为0.11 m×0.11 m×6 m、顶部设置了若干个泄爆口的小比例综合管廊进行了试验研究,结果表明:设置泄爆口往往能降低管沟内50%以上的超压和冲量值。另外,Hou等[10]结合了管沟的实际情况,自主设计了一种Y形结构的分支管沟,还原了山东青岛“11•22”爆炸事故的破坏场景,分析了管沟内外的超压时程曲线,划分了管沟内超压变化的不同阶段,研究表明管沟盖板会同时受到前驱波和燃烧波的影响。目前,试验研究大多是在长径比较大的小截面管道中进行的。受环境条件、试验操作危险性和外界其他复杂因素的约束,关于大尺寸原型试验研究的报道较少。
数值模拟分析具有经济性和高效性的优点,也是研究可燃气体爆炸的一种重要手段。目前,全球通用的燃气爆炸效应分析数值模拟软件主要有FLACS、AutoReaGas、Fluent、CFX等,其中FLACS是一款被众多国际组织和标准认证的气体爆炸分析专业软件。宫广东等[11]应用AutoReaGas软件定量研究了管道中瓦斯爆炸的特性,分析了不同初始压力、初始温度、障碍物形状、尺寸及空间位置等条件对管道中瓦斯爆炸超压及火焰速度的影响。龚燚[12]利用Fluent软件对规范中防火分区内一段长为200 m的综合管廊燃气仓进行了可燃气体爆炸数值模拟分析,结果表明:当泄爆口设置在投料仓和通风口之间时,泄爆效果较好。董浩宇[13]、刘洋等[14]、陈晓坤等[15]、王涛[16]、Hisken等[17]和Zhang等[18]利用FLACS软件分别对综合管廊、密闭方形管道、独头巷道、圆形玻璃管道、带有障碍物的方形管道、隧道进行了可燃气体爆炸数值模拟,研究了爆炸冲击波在不同的长直空间结构内的传播特性。目前的研究多是通过建立长直空间的数值模型分析空间内部的爆炸荷载,而对长直空间内的荷载通过泄爆口向地面传播的规律及特征的研究较少。
本文中,将基于大尺寸管沟可燃气体爆炸试验,利用FLACS软件对管沟内可燃气体爆炸产生的冲击波通过泄爆口向地面传播的机理、超压时程曲线特征及超压峰值分布规律进行数值模拟研究。
1. 模型验证
1.1 试验
试验系统的主体是管沟装置,总长度为30 m,由10节长3 m、截面尺寸为0.8 m×0.8 m的长方体组装而成,如图1(a)所示。试验系统还包含浓度仪、点火系统、燃气装置、气体混合风机、传感器及数据采集装置等[19],如图1(b)所示。其中管沟端部中间位置预留小孔用以安装点火装置,管沟尾部设置有尺寸为0.8 m×0.8 m的泄爆口,管沟顶部也设置有尺寸为0.6 m×0.6 m的泄爆口,泄爆口的数量和尺寸可以根据实际工况灵活调整。
试验以常见的可燃气体甲烷作为研究对象,分别设置了不同的泄爆口方式和甲烷体积分数,具体工况见表1。点火点位于管沟端部中央位置,仅在管沟端部这一节充入甲烷,即可燃气云尺寸为0.8 m×0.8 m×3 m。
表 1 管沟可燃气体爆炸工况记录Table 1. Working condition record of combustible gas explosion in pipe trench工况 泄爆口位置 顶部泄爆口数目 甲烷体积分数/% 1-A 尾部+顶部 3 7.5 1-B 尾部+顶部 3 8.5 1-C 尾部+顶部 3 9.5 1-D 尾部+顶部 3 10.5 1-E 尾部+顶部 3 11.5 2-C 密闭,无泄爆口 0 9.5 3-C 尾部 0 9.5 1.2 模型验证及网格敏感性分析
基于上述管沟试验装置,利用FLACS软件建立了相应的数值模型。通过分析以往城市地下管沟燃气爆炸事故可知,爆炸荷载如果通过管沟顶部或尾部的孔口直接传播到地面,通常会造成比较大的危害,因此选取表1中试验工况1-C来验证,点火点设置在管沟端部中间位置,管沟内中线上设置10个压力测点(P1~P10),第1个测点距端部水平距离为1.1 m,其余相邻测点之间的水平间距均为3 m,如图2所示。
为了验证数值模型的网格敏感性,根据管沟截面尺寸以及泄爆口大小,选取了0.05、0.10和0.20 m共3种网格尺寸进行数值模拟,部分测点(P1、P5、P10)的数值模拟结果与试验数据的对比如图3所示。可以看到,当网格尺寸为0.20 m时,数值模拟的曲线变化趋势与试验结果基本吻合,可以较好地呈现测点的超压随时间的变化规律,而当网格尺寸为0.10和0.05 m时,模拟结果均大于试验数据。这主要是由于预混合燃烧的亚网格模型在网格单元尺寸较小时适用性较差,对燃烧速度的计算往往存在较大偏差,模拟结果会高估爆炸超压。综合考虑计算效率及准确性,选择0.20 m的网格尺寸来进行数值模拟。
2. 地面可燃气体爆炸冲击波传播规律分析
2.1 物理模型
城市管沟通常浅埋在土壤里,由于铺设的长度较长,为了便于检修和后期维护,每隔一定距离都会设置一个检查井。根据国家标准规定,当管沟高度为0.8~1.0 m时,检查井的最大间距为90 m[20]。因此在实际试验中,建立一段尺寸为180 m×1 m×1 m的城市浅埋管沟模型,如图4所示,其中检查井之间的距离为90 m,检查井上方的盖板为无约束泄爆口,为了研究管沟内的爆炸冲击波通过泄爆口向地面传播的分布规律,沿检查井X方向两侧各预留长度为45 m的管沟,沿管沟Y方向两侧各预留长度为10 m的地面,管沟顶部与地面平齐、地面上方空气域的高度为20 m(见图5),由于本文的边界条件设为平面波,因此可以避免反射回来的波对结果的影响。
考虑极端的爆炸情况,假设AB两端检查井之间(90 m×1 m×1 m)填充满体积分数为9.5%的甲烷,点火点位于管沟内2个检查井的中心。环境温度与大气压强分别设置为为20 ℃、100 kPa,边界条件设置为 平面波,湍流模型采用k-ε模型。为排除其他外界因素对管沟可燃气体爆炸的影响,假设:(1)爆炸产生的能量不会与管沟进行传热而损失;(2)管沟内没有其他障碍物;(3)壁面绝对光滑;(4)在可燃气体爆炸作用下管沟不发生形变。
在土壤及空气域所在的区域设置核心网格,网格尺寸为0.2 m,核心区域外伸展系数为1.2,如图5所示。在泄爆口上方的空气域中,沿X、Y、Z轴方向上分别设置14个测点,用于记录超压时程曲线的变化,由于在泄爆口附近的荷载较大、超压变化幅度较大,因此在尾部开孔附近的空气域中测点布置得较密集,在远离开孔的空气域中测点布置得较稀疏,具体设置如表2和图6所示。
表 2 测点的坐标Table 2. Coordinates of measuring pointsX 轴方向 Y 轴方向 Z 轴方向 测点 坐标/m 测点 坐标/m 测点 坐标/m X1 (25,0,1) Y1 (44.5,−10,1) Z1 (44.5,0,1) X2 (30,0,1) Y2 (44.5,−7,1) Z2 (44.5,0,2) X3 (35,0,1) Y3 (44.5,−5,1) Z3 (44.5,0,3) X4 (40,0,1) Y4 (44.5,−4,1) Z4 (44.5,0,4) X5 (41,0,1) Y5 (44.5,−3,1) Z5 (44.5,0,5) X6 (42,0,1) Y6 (44.5,−2,1) Z6 (44.5,0,6) X7 (43,0,1) Y7 (44.5,−1,1) Z7 (44.5,0,7) X8 (46,0,1) Y8 (44.5,1,1) Z8 (44.5,0,8) X9 (47,0,1) Y9 (44.5,2,1) Z9 (44.5,0,9) X10 (48,0,1) Y10 (44.5,3,1) Z10 (44.5,0,10) X11 (49,0,1) Y11 (44.5,4,1) Z11 (44.5,0,11) X12 (54,0,1) Y12 (44.5,5,1) Z12 (44.5,0,13) X13 (59,0,1) Y12 (44.5,7,1) Z12 (44.5,0,16) X14 (64,0,1) Y12 (44.5,10,1) Z12 (44.5,0,20) 2.2 管沟外冲击波传播过程及机理分析
通过对管沟外部X、Y、Z轴等3个方向不同测点处超压时程曲线特征的分析,将其划分为3个阶段:稳定段、超压峰值Δp1段、超压峰值Δp2段,如图7~8所示。随后研究了超压、燃料、燃烧产物以及速度矢量的动态变化过程,并由此阐述了不同阶段爆炸荷载产生的机理。
阶段Ⅰ—稳定段(
t <200 ms):管沟内的超压值不断提高,地面上测点的超压值几乎可以忽略且基本不变。图9为燃料和燃烧产物的体积分数二维分布云图。管沟内的可燃气体遇点火后发生燃烧反应,而此时产生的燃烧产物仅在管沟内部传播,未传播到泄爆口和地面(见图9(b)),同时管沟内的燃料还未通过泄爆口被排放到空气域中(见图9(a)),因此测点未出现超压峰值。阶段Ⅱ—超压峰值Δp1段(200 ms<
t <500 ms):空气域经历第1个超压峰值Δp1期间。随着燃烧反应的不断加剧,管沟内的火焰不断朝泄爆口方向传播,推动管沟内残余燃料通过泄爆口喷向空气域中(见图9(a))。当t =300 ms时,空气域中的残余燃料在流场的作用下(见图10)具有一定的运动速度,波的最大速度为数十米每秒,因此燃料会压缩分界面处的空气,形成一个有限幅度的压力扰动,导致空气域中产生压缩波(见图8)。当空气压缩波的超压作用结束后,测点就受稀疏波的影响,超压值迅速减小乃至低于大气压力。当t=350 ms时,空气域中的残余燃料不断减少,压缩波继续向周围空气域传播,其强度也逐渐变小直至为零。阶段Ⅲ—超压峰值Δp2段(
t >500 ms):管沟外空气域经历第2个超压峰值Δp2及之后的小幅振荡直至消失。在管沟中由于有泄爆口的存在,当波传播到泄爆口附近时,管沟内部气体与外部大气之间的存在明显的压力差,管沟内部压力能够得到释放,从而会出现负波峰。当t =825 ms时,火焰沿管沟的方向到达泄爆口附近,同时受管沟尾部壁面约束的作用而向空气域内传播,火焰波阵面产生的燃烧产物具有高温特性(如图11所示),并与空气域发生热量的交换,使得燃烧产物迅速膨胀,压缩周围空气介质,并推动周围空气向外运动,形成了密度突跃的阵面(如图12所示),即该强扰动以冲击波形式在空气中传播,因此出现第2个超压峰值Δp2。此时波的传播速度如图13所示,在该阶段中波的最大速度能够达到数百米每秒。由于火焰波的强度要远大于阶段Ⅱ中压缩波的能量,同时波的传播速度也远大于阶段Ⅱ中的速度,因此第Ⅲ阶段的最大超压峰值Δp2要远大于第Ⅱ阶段的超压峰值Δp1。随着燃烧产物的膨胀,其压力、密度和运动速度不断下降,随之能量密度就逐渐减小,当燃烧产物内的压力下降到大气压力时,火焰波就脱离燃烧产物向前运动并在尾部形成稀疏区。当t =850 ms时,火焰波不断衰减并大致呈不规则的半球形在空气域中向外传播,在泄爆口两侧沿Y方向大致呈对称分布,但是在泄爆口两侧沿X轴方向并不呈均匀对称分布(见图14),其可能的原因是受空气流场的影响,当火焰传到泄爆口后,偏向于沿火焰传播的方向分布。同时泄爆口附近还有一部分强度较小的残余火焰会向空气域中传播能量较低的火焰波,导致该阶段会出现多峰值波动状,但是随着振荡的进行,超压峰值在逐渐减小。2.3 超压时程曲线特征分析
在三维坐标系中,选取在X、Y、Z轴方向上具有代表性的全过程超压时程曲线,如图15所示。由图15(a)可以看出,沿X轴方向,各个测点上的第1个超压峰值Δp1相差较小,第2个超压峰值Δp2相差较大。在沿X轴方向上,最大超压峰值Δp2位于测点X8,约为13 kPa,且在泄爆口两侧,测点离泄爆口的距离越远,超压峰值Δp1越小。同时,也可以观察到测点X1~X7一侧的超压峰值小于测点X8~X14一侧的超压峰值。
由图15(b)可以看出,在沿Y轴方向上,最大峰值超压Δp2位于距离泄爆口最近的测点Y7和Y8,其值约为7 kPa,并且在泄爆口两侧超压峰值Δp2随着距离的增加呈对称衰减的趋势。
由图15(c)可以看出,在沿Z轴方向上,第2个超压峰值Δp2的最大值位于测点Z1,其值约为14.5 kPa,且随着测点高度的上升,火焰波在空气域中逐渐衰减导致超压峰值Δp2也逐渐减小。
2.4 超压峰值分布规律
通过对超压时程曲线的分析,可以看出,城市地下浅埋管沟内可燃气体爆炸冲击波通过泄爆口到达地面后有2个超压峰值Δp1和Δp2,图16为这2个超压峰值在X、Y、Z方向的分布。以泄爆口位置为参照点,管沟内冲击波传播的方向为X的正方向,管沟内冲击波传播的反方向为X的负方向,Z和Y方向则与坐标轴设置保持一致。由图16可知,与超压峰值Δp2相比,超压峰值Δp1很小,且各个测点之间波动幅度不大。对于超压峰值Δp2,其值随泄爆口距离的增加而衰减。同时,在Y方向上的超压峰值会明显小于X正方向和Z方向上的超压峰值,这是由于爆炸冲击波的传播具有方向性而引起的。
为进一步分析超压峰值Δp2的衰减规律,参照空气域中的超压峰值分布的数据,可以对其进行非线性曲线拟合,得到测点的超压峰值Δp2和泄爆口之间的距离d大致满足指数函数型的关系:
Δp2=a−bcd (1) 式中:a、b、c为常数。
根据已有测点的数据,通过拟合分别可以得到沿X、Y、Z方向上的Δp2-d关系式:
Δp2={0.574+21.301×0.624d沿X正方向0.131+6.001×0.589d沿X负方向0.518+12.775×0.492d沿Y方向0.865+19.364×0.685d沿Z方向 (2) 式中:Δp2为超压峰值,kPa;d为与泄爆口之间的距离,m。该拟合公式适用于甲烷体积分数为9.5%、盖板为无约束的情况。
将式(2)中超压峰值Δp2与距离d的指数函数型方程关系作拟合曲线,如图17所示。可以得到:曲线的拟合质量较理想,拟合度均在98.8%以上;Δp2拟合曲线在沿Z和X正方向上的超压峰值最大,但同时其衰减速度(斜率)也是最大的,同时Δp2拟合曲线在沿X负方向上的超压峰值最小,其衰减速度(斜率)也是最小的,当测点与泄爆口的距离足够远时,在Z、X、Y方向上的超压峰值大致相等。
3. 结 论
利用FLACS软件建立了城市地下浅埋管沟内燃气爆炸的数值模型,基于管沟可燃气体爆炸的试验数据验证了模型的准确性,随后分析了爆炸冲击波通过泄爆口向地面传播的机理和超压峰值分布规律,具体结论如下。
(1) 爆炸荷载通过泄爆口向地面传播过程可分为3个阶段:阶段Ⅰ,即稳定段,该阶段测点的超压值很小,基本未受影响;阶段Ⅱ,即超压峰值Δp1段,经历第1个超压峰值Δp1的期间;阶段Ⅲ,即超压峰值Δp2段,经历第2个超压峰值Δp2及之后的小幅振荡直至消失。
(2) 阶段Ⅱ的Δp1较小,主要是由残余燃料挤压空气产生压缩波引起的;而阶段Ⅲ的Δp2较Δp1更危险,其主要是由燃料燃烧产生火焰波引起的。
(3) 超压峰值Δp2和与泄爆口之间距离d的关系大致满足
Δp2=a−bcd 的关系式,且拟合效果较好。 -
表 1 实验中采用的流体和脉冲加速度性质参数
Table 1. Fluid combinations and characteristics for impusive accerleration experiments
No. 流体1 流体2 ρ1/(g·cm-3) ρ2/(g·cm-3) R A We Re 1 H2O CCl2F2 1.000 1.57 1.57 0.22 4 000 2 600 2 SF6 C4H10 0.067 0.81 12.10 0.85 1 100 8 000 3 SF6 CCl2F2 0.067 1.57 23.40 0.92 11 000 23 000 4 SF6 CCl2F2 0.032 1.57 49.10 0.96 6 000 25 000 -
[1] Richtmyer R D. Taylor instability in shock acceleration of compressible fluids[J]. Communicational Pure Applied Mathematics, 1960, 13(1): 297-319. [2] Meshkov E E. Instability of the interface of two gases accelerated by a shock wave[J]. Soviet Fluid Dynamics, 1969, 4(1): 101-104. doi: 10.1007/BF01015969 [3] Dimonte G, Remington B. Richtmyer-Meshkov experiments on the Nova laser at high compression[J]. Physical Review Letters, 1993, 70(12): 1806-1809. http://www.ncbi.nlm.nih.gov/pubmed/10053391?dopt=Abstract [4] Dimonte G, Frerking C E, Schneider M. Richtmyer-Meshkov instability in the turbulent regime[J]. Physical Review Letters, 1995, 74(24): 4855-4858. http://europepmc.org/abstract/MED/10058616 [5] Vetter M, Sturtevant B. Experiments on the Richtmyer-Meshkov instability of an air/SF6interface[J]. Shock Waves, 1995, 4(5): 247-252. doi: 10.1007/BF01416035 [6] Jourdan G, Houas L, Haas J F, et al. Thickness and volume measurements of a Richtmyer-Meshkov instability induced mixing zone in a square shock tube[J]. Journal of Fluid Mechanics, 1997, 349: 67-94. http://www.ingentaconnect.com/content/cupr/00221120/1997/00000349/00000001/art00003 [7] Alon U, Hecht J, Ofer D, et al. Power laws and similarity of Rayleigh-Taylor and Richtmyer-Meshkov mixing fronts at all density ratios[J]. Physical Review Letters, 1995, 74(4): 534-537. http://www.ncbi.nlm.nih.gov/pubmed/10058782 [8] Aglitskiy Y, Velikovich A L, Karasik M, et al. Basic hydrodynamics of Richtmyer-Meshkov-type growth and os-cillations in the inertial confinement fusion-relevant conditions[J]. Philosophical Transactions of Royal Society A, 2010, 368(1916): 1739-1768. http://www.ncbi.nlm.nih.gov/pubmed/20211882 [9] Llor A. Statistical hydrodynamic models for developed mixing instability flows[S]. Springer, 2005. [10] 杨玟, 王丽丽, 张树道. Rayleigh-Taylor不稳定性诱导湍流混合的数值模拟[J].工程力学, 2011, 28(6): 236-241. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gclx201106035Yang Min, Wang Li-li, Zhang Shu-dao. Numerical simulation of turbulent mixing induced by Rayleigh-Taylor instability[J]. Engineering Mechanics, 2011, 28(6): 236-241. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gclx201106035 [11] 杨玟, 王丽丽, 张树道, 等.用湍流模型研究Richtmyer-Meshkov不稳定性诱导的湍流混合[J].空气动力学学报, 2010, 28(1): 119-123. http://www.cqvip.com/Main/Detail.aspx?id=32960053Yang Min, Wang Li-li, Zhang Shu-dao, et al. The study of turbulent mixing induced by Richtmyer-Meshkov instability using turbulence model[J]. Acta Aerodynamica Sinica, 2010, 28(1): 119-123. http://www.cqvip.com/Main/Detail.aspx?id=32960053 [12] Cheng B. Modeling chaotic mixing[J]. Nuclear Weapons Journal, 2010, 1: 8-17. [13] Layzer D. On the gravitational instability of two superposed fluids in a gravitational field[J]. Astrophysics Journal, 1955, 122(1): 1-12. http://www.ams.org/mathscinet-getitem?mr=71198 [14] 杨玟, 王丽丽, 周海兵, 张树道.用浮阻力模型研究Rayleigh-Taylor不稳定性诱发混合现象[J].工程力学, 2013, 30(4): 385-391. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gclx201304057Yang Min, Wang Li-li, Zhang Shu-dao. Study on mixing induced by Rayleigh-Taylor instability using buoyancydrag model[J]. Engineering Mechanics, 2013, 30(4): 385-391. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gclx201304057 [15] Dimonte G, Schneider M. Density ratio dependence of Rayleigh-Taylor mixing for sustained and impulsive acceleration histories[J]. Physics of Fluids, 2000, 12(2): 304-321. doi: 10.1063/1.870309 [16] Dimonte G, Schneider M. Turbulent Richtmyer-Meshkov instability experiments with strong radiatively driven shocks[J]. Physics of Plasmas, 1997, 4(12): 4347-4357. [17] Dimonte G. Nonlinear evolution of the Rayleigh-Taylor and Richtmyer-Meshkov instabilities[J]. Physics of Plasmas, 1999, 6(5): 2009-2015. doi: 10.1063/1.873491 期刊类型引用(2)
1. 杨石刚,蔡炯炜,杨亚,孙文盛,门敬敏. 城市地下浅埋管沟可燃气体爆炸的灾害效应(Ⅱ):影响因素分析及后果评估. 爆炸与冲击. 2023(01): 157-168 . 本站查看
2. 杨峰,翟红波,苏健军,李尚青,肖洋,刘伟. 浅埋爆炸空气冲击波峰值压力计算方法研究. 火工品. 2023(01): 42-45 . 百度学术
其他类型引用(2)
-