陳長坤，張宇倫，趙小龍，雷 鵬

(1.中南大學(xué) 防災(zāi)科學(xué)與安全技術(shù)研究所，湖南 長沙 410075；2.中南大學(xué) 資源與安全工程學(xué)院，湖南 長沙 410083)

0 引言

隨著我國隧道數(shù)量日漸增多，易燃易爆危險(xiǎn)化學(xué)品運(yùn)輸量也在不斷增加，如甲醇、液化石油氣等[1]。這些危險(xiǎn)化學(xué)品在通過隧道時(shí)有發(fā)生爆炸的可能性，而一旦發(fā)生爆炸，往往會造成人員傷亡、隧道結(jié)構(gòu)破壞、交通運(yùn)輸系統(tǒng)癱瘓等嚴(yán)重后果[1-2]，如2017年5月河北浮圖峪5號隧道內(nèi)發(fā)生危險(xiǎn)品槽罐車爆炸事故，共造成15人死亡、21人受傷。在隧道內(nèi)可燃?xì)怏w爆炸過程中，障礙物的存在能夠?qū)Ρ▔毫龇植籍a(chǎn)生影響，因此，開展隧道內(nèi)障礙物對可燃?xì)怏w爆炸壓力場影響的研究具有必要性。

針對障礙物對于可燃?xì)怏w爆炸影響的研究是國內(nèi)外學(xué)者關(guān)注的熱點(diǎn)問題之一。Abdulmajid等[3-4]在圓柱形封閉管道中進(jìn)行甲烷-空氣爆炸實(shí)驗(yàn)，研究了阻塞比及障礙物間距對可燃?xì)怏w爆炸最大超壓值和火焰?zhèn)鞑ヌ匦缘挠绊?；Hall等[5]探究了障礙物數(shù)量對可燃?xì)怏w爆炸超壓的影響，發(fā)現(xiàn)爆炸超壓隨障礙物數(shù)量增加而增大，但存在上限，超過上限后，爆炸超壓隨障礙物數(shù)量增加而減?。籏indracki等[6]在密閉長管中進(jìn)行了甲烷-空氣爆炸實(shí)驗(yàn)，研究了點(diǎn)火源位置和阻塞比對可燃?xì)怏w爆炸發(fā)展過程的影響；Oh等[7]在封閉六面體容器中進(jìn)行液化石油氣-空氣爆炸實(shí)驗(yàn)，研究了容器內(nèi)置障礙物引起可燃?xì)怏w爆炸特性的變化情況； Masri等[8]進(jìn)行了障礙物形狀對火焰?zhèn)鞑ニ俣扔绊憣?shí)驗(yàn)，發(fā)現(xiàn)方形障礙物對火焰的加速作用比圓形和三角形障礙物更加明顯；尉存娟等[9-11]利用了水平管道式氣體爆炸實(shí)驗(yàn)系統(tǒng)，開展了置障下瓦斯爆炸實(shí)驗(yàn)，發(fā)現(xiàn)障礙物的存在對瓦斯爆炸超壓具有顯著的激勵(lì)作用；王海賓等[12-13]進(jìn)行了管道內(nèi)可燃?xì)怏w爆炸實(shí)驗(yàn)，研究了障礙物阻塞率和數(shù)量對可燃?xì)怏w爆炸超壓和火焰?zhèn)鞑ヌ匦缘挠绊?；郭丹彤等[14]運(yùn)用數(shù)值模擬手段，研究了阻塞率和障礙物結(jié)構(gòu)對可燃混合氣體爆炸壓力場的影響，發(fā)現(xiàn)同種障礙物結(jié)構(gòu)下，爆炸超壓隨阻塞率增加而呈現(xiàn)先增大后減小趨勢，相同阻塞率下，立體障礙物對爆炸壓力場的影響遠(yuǎn)大于平面障礙物。以上研究主要集中在密閉管道中障礙物的存在對可燃?xì)怏w爆炸壓力場分布的影響，為隧道內(nèi)存在障礙物情況下發(fā)生可燃?xì)怏w爆炸壓力場分布的研究提供了一定的理論依據(jù)和基礎(chǔ)。

本文針對甲醇具有易燃、易爆及在化工領(lǐng)域應(yīng)用廣泛等特點(diǎn)，選取其與空氣組成爆炸介質(zhì)，以常用的?；愤\(yùn)輸槽罐車建立障礙物模型。應(yīng)用CFD方法對隧道內(nèi)不同槽罐車分布下的甲醇蒸氣爆炸過程進(jìn)行數(shù)值模擬，在槽罐車尺寸和分布間距固定的情況下，研究隧道內(nèi)槽罐車數(shù)量及其相對點(diǎn)火源中心的位置對甲醇蒸氣爆炸壓力場分布規(guī)律的影響，以期為此類隧道內(nèi)爆炸事故的應(yīng)急救援提供一定的理論參考。

1 模型建立

1.1 物理模型及測點(diǎn)布置

以某隧道為原型建立二維物理模型，隧道長為110 m，高為6 m，如圖1所示。以常用的危化品運(yùn)輸槽罐車為原型建立二維障礙物模型，槽罐車長為10 m，高為3 m，布置間距固定為3 m，如圖2所示。運(yùn)用ICEM劃分結(jié)構(gòu)化網(wǎng)格，網(wǎng)格尺寸為0.1 m×0.1 m。隧道中設(shè)置29個(gè)測點(diǎn)，測點(diǎn)1～20布置在隧道頂部壁面上，間距為2 m，用以監(jiān)測甲醇蒸氣爆炸對隧道壁面的影響；測點(diǎn)A～I(xiàn)坐標(biāo)如圖1所示，A～F位于距地面2 m高度處，G～I(xiàn)位于距地面3.5 m高度處，用以監(jiān)測槽罐車周圍壓力，如圖1所示。

圖1 隧道模型及測點(diǎn)Fig.1 Tunnel model and monitors

圖2 初始條件及槽罐車布置Fig.2 Initial conditions and tank layout

1.2 工況及模擬條件設(shè)置

使用FLUENT軟件對隧道內(nèi)甲醇蒸氣爆炸過程進(jìn)行數(shù)值模擬，共設(shè)置6組工況。工況1不設(shè)置槽罐車；工況2，3，4分別設(shè)置1，2，3輛槽罐車；工況5，6只設(shè)置1輛槽罐車，其中心距點(diǎn)火源中心水平距離分別為23和36 m，如表1及圖2所示。可燃?xì)怏w濃度接近化學(xué)當(dāng)量比濃度時(shí)，爆炸產(chǎn)生的超壓最大[1]。假定隧道內(nèi)甲醇蒸氣在分布域內(nèi)的濃度為化學(xué)當(dāng)量比濃度，即甲醇蒸氣質(zhì)量分?jǐn)?shù)為13.44%，甲醇蒸氣分布在隧道中間10 m范圍內(nèi)。點(diǎn)火源中心設(shè)置在隧道的幾何中心(55 m，3 m)，點(diǎn)火源半徑為0.1 m，點(diǎn)火溫度為1 200 K。空氣域由氮?dú)夂脱鯕饨M成，體積分?jǐn)?shù)分別為79%和21%，初始溫度和壓力為300 K和101 325 kPa。

表1 工況設(shè)置Table1 Test conditions

采用壁面自適應(yīng)局部渦粘型(WALE)的大渦模擬(LES)模型描述湍流過程中混合氣體湍流變化，采用渦耗散概念模型(EDC)描述爆炸過程中混合氣體湍流流動與化學(xué)反應(yīng)的相互作用。參考文獻(xiàn)[14]，結(jié)合隧道內(nèi)甲醇蒸氣爆炸的預(yù)先模擬分析，確定模擬時(shí)間為0.12 s，時(shí)間步長則取為0.1 ms。同時(shí)為簡化分析，甲醇蒸氣爆炸的化學(xué)反應(yīng)考慮為一步不可逆反應(yīng)，且隧道壁面保持絕熱，并忽略熱輻射的影響。

1.3 網(wǎng)格獨(dú)立性驗(yàn)證

參考文獻(xiàn)[15]，對本文所使用的網(wǎng)格進(jìn)行獨(dú)立性驗(yàn)證，表2給出了3種不同尺寸網(wǎng)格下，t=50 ms時(shí)刻監(jiān)測點(diǎn)1的壓強(qiáng)及變化幅度。

表2 3種網(wǎng)格尺寸下的測點(diǎn)壓強(qiáng)Table2 Measured point pressure at 3 grid sizes

表2可以看出，隨著網(wǎng)格尺寸的增大，監(jiān)測點(diǎn)1壓強(qiáng)的變化幅度也逐漸增大，當(dāng)網(wǎng)格尺寸達(dá)到0.10 m×0.10 m時(shí)，模擬結(jié)果與0.08 m×0.08 m相比，其變化幅度為3.195%，計(jì)算誤差在可接受的范圍之內(nèi)，綜合考慮計(jì)算量等因素，選取主體網(wǎng)格尺寸為0.10 m×0.10 m。同時(shí)考慮到近壁面處湍流特性更為復(fù)雜，對近壁面位置的網(wǎng)格也做了一定程度的細(xì)化。

2 結(jié)果及分析

2.1 槽罐車數(shù)量對甲醇蒸氣爆炸壓力場的影響分析

圖3給出了不同槽罐車數(shù)量工況的4個(gè)時(shí)刻甲醇蒸氣爆炸壓強(qiáng)云圖。當(dāng)隧道內(nèi)沒有槽罐車時(shí)(如圖3(a)所示)，甲醇蒸氣爆炸產(chǎn)生的壓強(qiáng)相比于其他3組工況較小，壓力波整體表現(xiàn)為沿隧道縱向傳播。當(dāng)隧道內(nèi)存在1輛槽罐車時(shí)(如圖3(b)所示)，甲醇蒸氣爆炸產(chǎn)生的壓強(qiáng)明顯增大，槽罐車的存在使混合氣體流動更加紊亂，湍流程度得到加強(qiáng)。36 ms時(shí)，在壓力波沿隧道縱向傳播方向，槽罐車前端形成壓力渦團(tuán)，槽罐車前端測點(diǎn)A(60 m，2 m)壓力增大；54 ms時(shí)，混合氣體流動的紊亂程度加劇，湍流程度繼續(xù)加強(qiáng)，壓力波逐漸越過槽罐車，槽罐車上端測點(diǎn)G(65 m，3.5 m)及隧道上壁面測點(diǎn)壓強(qiáng)增大；100 ms時(shí)，壓力波繼續(xù)向隧道兩端傳播，壓強(qiáng)逐漸減小。當(dāng)隧道內(nèi)存在2輛槽罐車時(shí)(如圖3(c)所示)，甲醇蒸氣爆炸產(chǎn)生的壓強(qiáng)繼續(xù)增大，因槽罐車的阻礙使混合氣體流動的紊亂程度加強(qiáng)而形成更多數(shù)量的壓力渦團(tuán)，槽罐車周圍壓力激增，壓力波越過槽罐車后繼續(xù)向隧道兩端傳播，強(qiáng)度不斷衰減。當(dāng)隧道內(nèi)存在3輛槽罐車時(shí)(如圖3(d)所示)，甲醇蒸氣爆炸產(chǎn)生的壓強(qiáng)達(dá)到4組工況中的最大值，壓力波傳播更加紊亂，槽罐車局部形成壓力渦團(tuán)數(shù)量更多，槽罐車周圍壓強(qiáng)明顯增大。

圖3 不同數(shù)量槽罐車工況甲醇蒸氣爆炸壓強(qiáng)云圖Fig.3 Pressure contours of methanol steam explosion in different quantity tank truck conditions

圖4給出了不同槽罐車數(shù)量工況，隧道頂端測點(diǎn)1(55 m，6 m)壓強(qiáng)隨時(shí)間變化曲線。圖4(a)為隧道內(nèi)沒有槽罐車的工況，測點(diǎn)1的超壓峰值為81.4 kPa，而隧道內(nèi)分別存在1，2，3輛槽罐車的工況(如圖4(b)，(c)，(d)所示)，測點(diǎn)1的超壓峰值分別達(dá)到284.9，369.2，389.5 kPa，爆炸超壓峰值明顯增大。對比圖4時(shí)間坐標(biāo)可以看出，當(dāng)隧道內(nèi)沒有槽罐車時(shí)，甲醇蒸氣爆炸達(dá)到最大超壓需要52.9 ms；隧道內(nèi)存在槽罐車時(shí)，工況2，3，4中甲醇蒸氣爆炸達(dá)到最大超壓分別需要44.5，44.4，44.6 ms，爆炸達(dá)到最大超壓所用時(shí)間更短。槽罐車的存在對隧道內(nèi)甲醇蒸氣爆炸反應(yīng)有一定的激勵(lì)作用。

圖5為統(tǒng)計(jì)隧道內(nèi)29個(gè)測點(diǎn)在以上4組工況甲醇蒸氣爆炸過程中分別出現(xiàn)的最大超壓峰值繪制的爆炸最大超壓峰值隨槽罐車數(shù)量變化曲線。當(dāng)隧道內(nèi)沒有槽罐車時(shí)(工況1)，甲醇蒸氣爆炸最大超壓峰值89.4 kPa遠(yuǎn)小于其他3組工況；當(dāng)隧道內(nèi)存在1輛槽罐車時(shí)(工況2)，爆炸最大超壓峰值明顯增大，達(dá)到337.5 kPa。當(dāng)隧道內(nèi)存在2輛槽罐車時(shí)(工況3)，爆炸最大超壓峰值為369.2 kPa。當(dāng)隧道內(nèi)存在3輛槽罐車時(shí)(工況4)，爆炸最大超壓峰值繼續(xù)增大，達(dá)到559.4 kPa，為不存在槽罐車情況下的6.3倍?？梢钥闯?，甲醇蒸氣爆炸的最大超壓峰值隨隧道內(nèi)槽罐車數(shù)量的增加而明顯增大。槽罐車的存在使爆炸過程中混合氣體流動更加紊亂，加強(qiáng)了混合氣體的湍流程度，使爆炸反應(yīng)更完全，爆炸反應(yīng)產(chǎn)物迅速膨脹，產(chǎn)生更大的爆炸壓力。同時(shí)，溫度升高和壓力增大會促進(jìn)混合氣體的湍流流動，繼續(xù)加強(qiáng)混合氣體的湍流程度。爆炸反應(yīng)與湍流流動之間相互作用使隧道內(nèi)甲醇蒸氣爆炸最大超壓峰值隨槽罐車數(shù)量增加而顯著增大。

圖5 甲醇蒸氣爆炸最大超壓峰值變化曲線Fig.5 Maximum overpressure peak variation curve in methanol vapor explosion

圖6 隧道頂端測點(diǎn)超壓峰值變化曲線Fig.6 Curve of overpressure peak at the top of tunnel

圖6給出了工況1～4隧道頂端測點(diǎn)1～20超壓峰值變化曲線?？梢钥闯?，隧道內(nèi)存在槽罐車工況的隧道頂端測點(diǎn)超壓峰值明顯大于沒有槽罐車工況的超壓峰值。隧道頂端測點(diǎn)超壓峰值隨距點(diǎn)火源中心距離增加而減小，靠近隧道出口位置測點(diǎn)超壓峰值最小，說明隧道開口端具有一定的泄壓作用。同時(shí)，存在槽罐車工況的隧道壁面壓強(qiáng)雖然整體呈現(xiàn)衰減趨勢，但在槽罐車上方隧道壁面測點(diǎn)壓強(qiáng)增大，這是由于槽罐車的存在使壓力波經(jīng)過槽罐車時(shí)表現(xiàn)更為紊亂，在槽罐車周圍形成壓力渦團(tuán)，使局部壓力增大。

2.2 槽罐車位置對甲醇蒸氣爆炸壓力場的影響分析

圖7給出了工況2，5，6(只設(shè)1輛槽罐車，其中心距點(diǎn)火源水平距離分別為10，23，36 m)的爆炸壓強(qiáng)云圖。當(dāng)槽罐車設(shè)置在距點(diǎn)火源水平距離5 m處時(shí)，甲醇蒸氣爆炸壓力波整體沿隧道縱向傳播，混合氣體在運(yùn)動過程中受槽罐車的阻礙而更加紊亂，混合氣體湍流程度得到加強(qiáng)，爆炸反應(yīng)更加劇烈，槽罐車前端形成大量壓力渦團(tuán)，使槽罐車前端測點(diǎn)A(60 m，2 m)壓力增大，54 ms時(shí)，混合氣體流動紊亂加劇，湍流程度繼續(xù)加強(qiáng)，壓力波逐漸越過槽罐車，槽罐車上端測點(diǎn)G(65 m，3.5 m)及隧道上壁面測點(diǎn)壓強(qiáng)增大。100 ms時(shí)，壓力波繼續(xù)向隧道兩端傳播，強(qiáng)度衰減。當(dāng)槽罐車分別設(shè)在距點(diǎn)火源水平距離為23和36 m處時(shí)(如圖7(b)，(c)所示)，甲醇蒸氣爆炸壓力場與工況2表現(xiàn)無明顯的差別，壓力波整體表現(xiàn)為沿隧道縱向傳播，混合氣體流動遇到槽罐車阻礙后，混合氣體湍流程度得到加強(qiáng)，爆炸反應(yīng)更加劇烈，在槽罐車周圍形成壓力渦團(tuán)，使槽罐車周圍壓力增大，壓力波越過槽罐車后繼續(xù)向隧道兩端傳播，強(qiáng)度衰減。

圖7 不同槽罐車位置工況甲醇蒸氣爆炸壓強(qiáng)云圖Fig.7 Pressure contours of methanol steam explosion in different tanker locations

圖8 隧道頂端測點(diǎn)超壓峰值變化曲線Fig.8 Peak overpressure variation curves at the top of the tunnel

圖8給出了工況1，2，5，6的隧道頂端測點(diǎn)1～20超壓峰值變化曲線。當(dāng)僅設(shè)置一個(gè)槽罐車(工況2，5，6)時(shí)，隧道頂端測點(diǎn)超壓峰值明顯大于沒有槽罐車工況的超壓峰值。槽罐車的存在加強(qiáng)了爆炸過程中混合氣體的湍流程度，使爆炸反應(yīng)更加劇烈，導(dǎo)致產(chǎn)生更大壓力，釋放更多能量，而溫度和壓力的提升通過爆炸反應(yīng)與湍流流動之間的相互作用，繼續(xù)加強(qiáng)了混合氣體的湍流程度，進(jìn)一步激勵(lì)爆炸反應(yīng)產(chǎn)生更大壓強(qiáng)，從而使設(shè)置槽罐車工況的隧道頂端壁面測點(diǎn)超壓峰值明顯大于沒有槽罐車工況下的超壓峰值。同時(shí)，工況1(未設(shè)置槽罐車)的隧道頂端測點(diǎn)壓強(qiáng)呈現(xiàn)出隨距點(diǎn)火源中心水平距離增加而衰減的趨勢；而工況2，5，6(只設(shè)置1輛槽罐車，其中心距點(diǎn)火源中心水平距離分別為10，23，36 m)的隧道頂端測點(diǎn)壓強(qiáng)只是在整體上呈現(xiàn)衰減趨勢，在爆炸過程中槽罐車周圍會積聚大量壓力渦團(tuán)，使槽罐車周圍的隧道頂端壁面測點(diǎn)壓強(qiáng)局部增大。隧道內(nèi)槽罐車相對點(diǎn)火源中心的位置對甲醇蒸氣爆炸壓力場的分布有一定的影響。

3 結(jié)論

1)在特定條件下隧道內(nèi)槽罐車的存在對甲醇蒸氣爆炸有一定的激勵(lì)作用，可以使爆炸反應(yīng)更加劇烈，壓力波傳播更加紊亂，爆炸達(dá)到最大超壓的時(shí)間減少。

2)在甲醇蒸氣爆炸過程中，槽罐車的存在可提高混合氣體的湍流程度，激勵(lì)了爆炸反應(yīng)，進(jìn)而又促進(jìn)了混合氣體的湍流流動，二者之間的相互作用使甲醇蒸氣爆炸超壓隨槽罐車數(shù)量的增加而顯著增大，模擬過程中槽罐車數(shù)量從0增加到4時(shí)，該隧道內(nèi)甲醇蒸氣爆炸超壓峰值由89.4 kPa增大到559.4 kPa。

3)當(dāng)僅存在1輛槽罐車時(shí)，槽罐車相對點(diǎn)火源中心的位置可以影響爆炸過程中槽罐車周圍空間壓力場的分布，使槽罐車周圍壓強(qiáng)增大。

[1] 陳長坤,雷鵬,史聰靈,等. 公路隧道內(nèi)氫氣和丙烷爆炸數(shù)值模擬對比分析[J]. 中國安全生產(chǎn)科學(xué)技術(shù), 2017, 13(1): 151-155.

CHEN Changkun, LEI Peng, SHI Congling, et al. Comparative numerical simulation analysis on explosion of hydrogen and propane in highway tunnel [J]. Journal of Safety Science and Technology, 2017, 13(1): 151-155.

[2] 王海賓. 甲烷爆炸沖擊波對動物損傷研究[D]. 太原：中北大學(xué), 2015.

[3] NA’INNA A M, PHYLAKTOU H N, ANDREWS G E. The acceleration of flames in tube explosions with two obstacles as a function of the obstacle separation distance [J]. Journal of Loss Prevention in the Process Industries, 2013, 26(6): 1597-1603.

[4] NA’INNA A M, PHYLAKTOU H N, ANDREWS G E. Effects of obstacle separation distance on gas explosions: the influence of obstacle blockage ratio [J]. Procedia Engineering, 2014(84): 306-319.

[5] HALL R, MASRI A R, YAROSHEHYK P, et al. Effects of position and frequency of obstacles on turbulent premixed propagating flames [J]. Combustion & Flame, 2009, 156(2): 439-446.

[6] KINDRACKI J, KOBIERA A, RARATA G, et al. Influence of ignition position and obstacles on explosion development in methane-air mixture in closed vessels [J]. Journal of Loss Prevention in the Process Industries, 2007, 20(4-6): 551-561.

[7] OH K H, KIM H, KIM J B, et al. A study on the obstacle-induced variation of the gas explosion characteristics [J]. Journal of Loss Prevention in the Process Industries, 2001, 14(6): 597-602.

[8] MASRI A R, IBRAHIM S S, NEHZAT N, et al. Experimental study of premixed flame propagation over various solid obstructions [J]. Experimental Thermal & Fluid Science, 2000, 21(1): 109-116.

[9] 尉存娟,譚迎新,郝永飛,等. 密閉管道內(nèi)障礙物對瓦斯爆炸壓力的影響[J]. 煤炭科學(xué)技術(shù), 2013, 41(7): 86-88, 128.

YU Cunjuan, TAN Yingxin, HAO Yongfei, et al. Influence on gas explosion pressure by obstacles in the closed tube [J]. Coal Science and Technology, 2013, 41(7): 86-88,128.

[10] 陳成,胡雙啟,譚迎新. 螺旋形障礙物對瓦斯爆炸影響的研究[J]. 中國安全生產(chǎn)科學(xué)技術(shù), 2012, 8(6): 44-47.

CHEN Cheng, HU Shuangqi, TAN Yingxin. Study on influence of spiral obstacle to gas explosion [J]. Journal of Safety Science and Technology, 2012, 8(6): 44-47.

[11] 秦澗,譚迎新,王志青,等. 管道內(nèi)障礙物形狀對瓦斯爆炸影響的試驗(yàn)研究[J]. 煤炭科學(xué)技術(shù), 2012, 40(2): 60-62.

QIN Jian, TAN Yingxin, WANG Zhiqing, et al. Experiment study on obstacle shape in pipeline affected to gas explosion [J]. Coal Science and Technology, 2012, 40(2): 60-62.

[12] 王海賓,尉存娟,譚迎新. 水平管道內(nèi)障礙物對氣體爆炸壓力影響的研究[J]. 廣州化工, 2011, 39(24): 43-46.

WANG Haibin, YU Cunjuan, TAN Yingxin. Effects of pressure of gas explosion in the pipeline of obstacles [J]. Guangzhou Chemical Industry, 2011, 39(24): 43-46.

[13] 許航,譚迎新,孟璐,等. 長管中立體障礙物對瓦斯爆炸特性影響的研究[J]. 中國安全生產(chǎn)科學(xué)技術(shù), 2012, 8(9): 17-21.

XU Hang, TAN Yingxin, MENG Lu, et al. Study on effects of stereo obstacles on methane explosion in long pipe [J]. Journal of Safety Science and Technology, 2012, 8(9): 17-21.

[14] 郭丹彤,呂淑然,楊凱. 障礙物布置對氣體爆炸壓力場的影響效果研究[J]. 中國安全生產(chǎn)科學(xué)技術(shù), 2015, 11(9): 88-93.

GUO Dantong, LYU Shuran, YANG Kai. Research on impact effect of obstacle arrangement to pressure field of gas explosion [J]. Journal of Safety Science and Technology, 2015, 11(9): 88-93.

[15] 祝釗. 管道瓦斯爆炸流場及其影響因素?cái)?shù)值模擬研究[J]. 中國礦業(yè)大學(xué)學(xué)報(bào), 2017, 46(2): 300-305.

ZHU Zhao. Numerical simulation of gas explosion flow field characteristics and influencing factors in pipelines [J]. Journal of China University of Mining & Technology, 2017, 46(2): 300-305.