李　蒙，杜　揚，李國慶，齊　圣，王世茂，韋世豪

(中國人民解放軍陸軍勤務(wù)學(xué)院 供油系，重慶 401331)

0　引言

可燃?xì)怏w爆炸是工業(yè)生產(chǎn)和生活領(lǐng)域爆炸災(zāi)害的主要形式之一。為了防止氣體爆炸帶來重大的經(jīng)濟(jì)損失和人員傷亡，人們通過實驗和模擬手段進(jìn)行了深入的研究。已有的研究成果顯示，可燃?xì)怏w所在的容器形狀[1-6]、初始濃度[1,3,7-8]、初始溫度與壓力、障礙物數(shù)量[2,11]都會對可燃?xì)怏w爆炸釋放的能量產(chǎn)生影響。

目前，對于全封閉管道內(nèi)油氣爆炸特性的研究較為廣泛，杜揚[9]對受限空間內(nèi)局部油氣爆燃超壓與火焰?zhèn)鞑ヌ匦赃M(jìn)行了研究，發(fā)現(xiàn)傳播過程可分為5個階段，且出現(xiàn)了tulip火焰，得到了最大爆燃超壓峰值出現(xiàn)在當(dāng)量比1.2處等結(jié)論；齊圣[10]對于受限空間油氣爆燃火焰形態(tài)進(jìn)行了探究，給出了區(qū)分3種火焰形態(tài)的臨界條件。對于半開口管道內(nèi)氣體爆炸特性的研究大多停留在對于內(nèi)場壓力變化的探究，對于外場壓力的變化以及火焰行為的研究偏少，并且采用的管道長徑比較小，一般不超過5[12-16]，研究成果普適性有限?；谏鲜鲈?，本文采用長徑比為10的有機(jī)玻璃全透明管道，通過改變初始油氣濃度，重點探究油氣濃度對油氣爆炸超壓峰值、升壓速率、火焰形態(tài)和火焰?zhèn)鞑ニ俣鹊鹊挠绊懸?guī)律。

1　實驗系統(tǒng)與方法

1.1　實驗系統(tǒng)裝置

實驗系統(tǒng)裝置圖如圖1所示。主實驗裝置為一段截面為正方形的長直管道(管道長L=1 000 mm，橫截面尺寸為100 mm×100 mm，容積V=10 L，管道長徑比L/D=10)。管道一端用鋼制的盲板密封，為了確保連接的密封良好，在盲板與管道法蘭之間使用密封墊圈密封；另一端開口使用聚乙烯薄膜密封，用以確保在管道充氣階段不會有可燃?xì)怏w外泄。

圖1　實驗系統(tǒng)示意Fig.1　The experimental system schematic diagram

實驗測試系統(tǒng)由壓力動態(tài)數(shù)據(jù)采集系統(tǒng)、高速攝影儀、配氣系統(tǒng)、碳?xì)錆舛葴y試系統(tǒng)、點火系統(tǒng)等組成。測試系統(tǒng)使用成都泰斯特公司的DAP 7.10，壓力傳感器使用寶雞市智星傳感器有限責(zé)任公司的ZXP660高頻瞬態(tài)壓力傳感器(量程0～200 kPa，精度誤差<0.3%)。高速攝影儀使用PHOTRON公司的FASTCAM-ultima 512，拍攝速度的設(shè)定為1 000 幀/秒。汽油蒸汽由配氣系統(tǒng)產(chǎn)生，并采用碳?xì)錅y試儀GXH-1050(使用環(huán)境0～40℃，重復(fù)性≤±1%，線性誤差≤2%F.S，量程為0.01%～100%)檢測油蒸汽的體積濃度，使之達(dá)到實驗所需的初始油氣濃度。實驗中產(chǎn)生油氣的方法如下：在真空泵的作用下，空氣在密閉管道和油瓶中循環(huán)，產(chǎn)生初始油氣混合物，此時閥門1，2，3，5打開，閥門4關(guān)閉(閥門編號見圖1)；當(dāng)油氣在密閉系統(tǒng)內(nèi)循環(huán)一段時間之后(根據(jù)實驗所需的初始油氣濃度決定時間長短)，再關(guān)閉閥門2和3，打開閥門4，讓系統(tǒng)內(nèi)混合氣體循環(huán)大約3 min，使氣體均勻混合。點火系統(tǒng)采用實驗室定制的抗干擾點火系統(tǒng)，點火范圍為2～20 J；點火頭放置于管道盲板中心位置，初始點火能量為6 J。

1.2　實驗內(nèi)容和方法

為充分探究油氣濃度對油氣泄壓的影響，本次實驗分別采用2.13%，1.87%，1.73%，1.5%，1.29% 和1.1%等6種不同的濃度進(jìn)行實驗。在對管道內(nèi)充入油氣之前，先使用聚乙烯薄膜將管道開口進(jìn)行端封堵以避免油氣外泄；實驗時，為記錄起爆過程管道內(nèi)的壓力—時間曲線，將5個壓力傳感器Pa，Pb，Pc，Pd，Pe均勻的布置在管道上表面，同時使用高速攝影儀記錄了油氣從起爆到熄滅的整個過程中火焰形態(tài)和火焰鋒面位置的變化過程。為了保證各個測試系統(tǒng)工作的同步性，實驗中采用同步觸發(fā)控制裝置來保證點火系統(tǒng)、高速攝影系統(tǒng)、壓力采集系統(tǒng)的同步觸發(fā)。為了保證管道內(nèi)可燃?xì)怏w與空氣均勻混合，每次點火前管道內(nèi)的氣體均會靜置約30 s。實驗預(yù)混氣體初始溫度為實驗溫度，預(yù)混氣體的初始壓力為當(dāng)?shù)卮髿鈮?。為保證數(shù)據(jù)的準(zhǔn)確性，減小偶然誤差對實驗造成的影響，每組實驗都進(jìn)行了至少3次重復(fù)試驗。

2　實驗結(jié)果與討論

2.1　油氣爆炸超壓隨時間的變化規(guī)律

從圖2可以明顯觀察到2個峰值：Pv與Pmax，其中，Pv是因為管道開口處薄膜破裂瞬間產(chǎn)生的，大多數(shù)的文獻(xiàn)[13]將其定義為“泄壓峰值”；Pmax則是由火焰鋒面穿過開口處聚集的高壓油氣團(tuán)發(fā)生爆炸產(chǎn)生的沖擊波導(dǎo)致[20]。由此，可以將超壓—時間變化曲線劃分為4個階段：

第1階段，定容爆炸。在開口處薄膜尚未破裂時，管道內(nèi)為密閉空間，當(dāng)空間內(nèi)的油氣混合物被點燃后，管道內(nèi)的溫度迅速上升，氣體受熱膨脹，管道內(nèi)的壓力隨之迅速增加。

圖2　油氣爆炸超壓特性曲線Fig.2　The overpressure characteristic curve of gasoline-air vapor explosion

第2階段，開口泄流。當(dāng)管道內(nèi)壓力達(dá)到泄壓峰值時，管口處的薄膜破裂，管道內(nèi)高度壓縮的氣體泄出，管道內(nèi)的壓力也隨之下降。

第3階段，外部爆炸。當(dāng)管道內(nèi)的壓縮氣體泄出后，在管口附近形成較高濃度的聚集油氣團(tuán)，火焰鋒面穿過油氣團(tuán)時，隨即引發(fā)劇烈的爆炸，氣體受熱膨脹，壓力急劇上升，產(chǎn)生外部超壓。以濃度為1.5%時為例，對壓力曲線(圖2、圖3)以及火焰位置圖線(圖7)進(jìn)行對比，可以知道，在火焰離開管道之后，管道內(nèi)超壓達(dá)到最大值，對于外部超壓曲線來說，其壓力在火焰離開管道之前沒有較大的變化，在火焰離開管道之后，出現(xiàn)了壓力波動并達(dá)到超壓峰值，由此可以判斷，外部氣體燃燒爆炸導(dǎo)致了外部的超壓。對于這一現(xiàn)象，相關(guān)文獻(xiàn)[2，22-23]也給出了類似解釋。

圖3　管道外部超壓曲線Fig.3　External overpressure characteristic curve

第4階段，衰減震蕩。爆炸發(fā)生后，爆炸生成物與周圍的空氣進(jìn)行劇烈的熱交換，油氣團(tuán)內(nèi)部溫度下降，壓力急劇衰減。管道內(nèi)部因為壓力下降，產(chǎn)生負(fù)壓形成局部的真空區(qū)，將部分未燃燒油氣混合物倒吸入管道，與管道內(nèi)的火焰鋒面碰撞時再次發(fā)生爆炸，壓力再次上升；爆炸后管道內(nèi)的空氣再次被擠出管道，管道內(nèi)壓力下降，隨之反復(fù)進(jìn)行這2個過程，壓力振蕩衰減。

2.2　油氣濃度對油氣爆炸最大爆炸超壓峰值變化規(guī)律的影響

2.2.1管道內(nèi)部超壓

圖4為不同油氣濃度下油氣管道內(nèi)部超壓峰值曲線，實驗圖像可以看出，在不同油氣濃度下，長直管道內(nèi)油氣爆炸超壓隨濃度變化呈現(xiàn)先增加后減小的趨勢。

圖4　管道內(nèi)部超壓特性曲線Fig.4　Internal overpressure characteristic curve

在濃度為1.73%時，最大爆炸超壓峰值取到最大值，大小約為10.416 kPa，對于這一現(xiàn)象，在已發(fā)表的文獻(xiàn)中有過論述[10,17]，其作用機(jī)理為：在油氣混合物處于較低濃度時，氧氣濃度較為充足，油氣混合物含量較少，可以實現(xiàn)充分的反應(yīng)，此時油氣爆炸超壓的大小取決于油氣混合物初始的濃度；隨著油氣混合物濃度的增加，反應(yīng)越來越劇烈，當(dāng)反應(yīng)物剛好可以消耗完氧氣時，爆炸超壓達(dá)到最大值，此時油氣混合物的濃度被稱為化學(xué)當(dāng)量濃度；當(dāng)濃度繼續(xù)增大時，氧氣含量相對較少，化學(xué)反應(yīng)進(jìn)行不完全，反應(yīng)產(chǎn)物中生成大量的CO，CO被O2氧化為CO2，這個過程中大量能量沒有釋放出來，同時CO2含量的增加抑制了爆炸的強(qiáng)度，所以隨著油氣混合物濃度繼續(xù)的增大，油氣爆炸超壓值持續(xù)下降。

2.2.2管道外部超壓

由圖5可以看出，在不同油氣濃度下，長直管道外部壓力隨著油氣濃度的變化呈現(xiàn)一直增大的趨勢。這是由于隨著油氣濃度的增大，薄膜破裂后泄放到外部的油氣混合物逐漸增加，在外部發(fā)生的爆炸超壓液隨之增大。由于空氣的稀釋作用，在此次的實驗中，外部的油氣混合物濃度并沒有達(dá)到當(dāng)量比濃度，所以呈現(xiàn)一直上升的趨勢。在濃度為1.73%時，外部超壓變化速率陡升，外部爆炸超壓發(fā)生躍升，說明此時外部油氣濃度達(dá)到了較為合適的濃度。

圖5　管道外部超壓特性曲線Fig.5　 External overpressure characteristic curve

2.2.3油氣濃度對平均升壓速率的影響

定義平均升壓速度=最大爆炸超壓峰值/達(dá)到最大爆炸超壓峰值所用的時間。升壓速度的大小直觀的反映了爆炸發(fā)生時能量釋放速度的快慢，升壓速度越大，能量釋放速度越快。由圖6曲線可以看出，初始濃度為1.29%時，升壓速度較小，隨著濃度的增加，升壓速度加快，并在濃度為1.73%時達(dá)到最大，約為333.45 kPa/s，當(dāng)濃度繼續(xù)增加時，升壓速率隨之下降。同時可以從曲線看出，升壓速率增加的變化率隨著濃度的增加而減小，其減小的變化率也隨著濃度的增加而減小。

圖6　不同濃度下升壓速率曲線Fig.6　The curve of rate of pressure rise under different concentrations

2.3　油氣濃度對火焰行為的影響

2.3.1油氣濃度對火焰鋒面位置的影響

圖7顯示的是在5種不同油氣濃度(YCH)工況下火焰鋒面位置隨時間的變化關(guān)系。火焰鋒面的具體位置是通過測量管道底部與火焰鋒面之間在管道軸線方向上的最大距離來獲得[13]。從圖中可以看出，濃度的變化對于火焰鋒面位置的變化有著較為顯著地影響。當(dāng)濃度較小時(以YCH等于1.29%為例)，火焰鋒面位置變化速度較慢，達(dá)到同一位置所需要花費的時間最多，當(dāng)濃度上升時(以YCH等于1.5%為例)，火焰鋒面位置變化速度有所提升。較為巧合的是，當(dāng)濃度達(dá)到1.73%和1.87%時，二者的火焰鋒面移動速度十分相似，這說明在這個濃度范圍內(nèi)，參與反應(yīng)的反應(yīng)物當(dāng)量比較為合適，化學(xué)反應(yīng)速度達(dá)到了最大。當(dāng)濃度繼續(xù)增大時，火焰鋒面位置變化速度下降。同時可以從圖中觀察，火焰鋒面達(dá)到出口位置時，各工況在濃度由小到大的順序下分別花費了54，42，34，34和50 ms，達(dá)到最大火焰鋒面位置所花時間分別是57，50，42，44，53 ms左右。可見濃度的變化對火焰鋒面位置的傳播速度有著較大的影響，在當(dāng)量濃度比下，火焰鋒面位置的傳播速度最大，并且在當(dāng)量濃度比下火焰鋒面?zhèn)鞑サ淖畲缶嚯x也更遠(yuǎn)。

圖7　5種不同濃度工況下火焰鋒面位置曲線Fig.7　The flame front position under five different concentrations

2.3.2火焰形態(tài)隨時間變化規(guī)律

圖8為高速攝影儀記錄下的在2種不同濃度下管道內(nèi)油氣爆炸火焰?zhèn)鞑バ螒B(tài)隨時間變化的圖像。根據(jù)火焰鋒面位置曲線，選取2種位置曲線較為一致的工況進(jìn)行火焰形態(tài)的分析。從圖8中可見，管道內(nèi)火焰的瞬變行為可以分為以下4個階段：

圖8　油氣濃度為1.73%(結(jié)構(gòu)2)和1.87%(結(jié)構(gòu)1)時火焰?zhèn)鞑バ螒B(tài)Fig.8　The flame propagation form with gasoline-air concentration of 1.73%(configuration 1) and 1.87%(configuration 2)

第1階段，0時刻至30 ms，層流拉伸階段。當(dāng)管道內(nèi)氣體被點爆時，點火端發(fā)出明亮的火焰，以點火端為中心，火焰陣面呈現(xiàn)規(guī)則的球面形狀向四周擴(kuò)散，當(dāng)火焰陣面分別到達(dá)右壁面、上下壁面后，火焰被擠壓發(fā)生變形，以橢圓形形狀向開口處運動，火焰陣面球面曲率逐漸減小，火焰陣面逐漸轉(zhuǎn)變?yōu)椤爸讣庑汀?。同時可以明顯的看出，燃燒的火焰區(qū)域有著明顯的分區(qū)現(xiàn)象，火焰前沿有一個近球形的黃色火焰燃燒區(qū)域，火焰外圍為藍(lán)色火焰鋒面區(qū)域，這種分區(qū)現(xiàn)象是因為在藍(lán)色火焰鋒面區(qū)域，油氣混合物與空氣接觸發(fā)生了完全的燃燒；而在黃色火焰燃燒區(qū)域，由于高溫的誘導(dǎo)，同時因為部分燃燒不完全,產(chǎn)生了一定量的碳粒子，故而發(fā)出了明亮的黃色。

第2階段，>30～40 ms，破壞變形階段。隨著管道內(nèi)壓力的震蕩增加，管道端口處薄膜破裂，在壓力波的驅(qū)動下管道破口處產(chǎn)生強(qiáng)烈的泄流現(xiàn)象，破口處的擾動增強(qiáng)了管道內(nèi)的湍流度。根據(jù)Rayleigh-Taylor不穩(wěn)定可以知道，當(dāng)?shù)兔芏鹊娜紵蜌猱a(chǎn)物向高密度未燃油氣反應(yīng)物加速時，其過程是不穩(wěn)定的，其表現(xiàn)為火焰陣面的扭曲與變形。同時，火焰陣面前方的未燃?xì)怏w受熱膨脹，使得火焰陣面處的密度場發(fā)生較為復(fù)雜的變化。未燃燒區(qū)域產(chǎn)生較小的渦流，火焰陣面進(jìn)一步發(fā)生拉伸變形，這被稱為亥姆霍茲不穩(wěn)定性?？梢詮膱D中看到破口處火焰陣面與管道內(nèi)層流拉伸時形態(tài)的不同，此時火焰陣面呈現(xiàn)毛刷狀湍流火焰的形態(tài)。

第3階段，>40～45 ms，管道外爆炸階段。在壓力波與火焰的驅(qū)動下，大量高速運動的未燃油氣從破口處泄出，發(fā)生迅速的膨脹與形變，管道外部流場發(fā)生深刻復(fù)雜的變化?；鹧姘l(fā)生反轉(zhuǎn)、拉伸與分層的現(xiàn)象，形成類似蘑菇云的形狀。在圖8(a)中時間為45 ms時，可以看到火焰上層為藍(lán)色的火焰、下層為橙黃色火焰。這是因為，在上層火焰中，反應(yīng)物為以甲烷為主的低碳原子數(shù)分子，其反應(yīng)呈現(xiàn)淡藍(lán)色火焰；在下層火焰中，反應(yīng)物為高碳原子數(shù)分子，其反應(yīng)呈現(xiàn)的是亮黃色的顏色。

第4階段，45 ms之后為回流燃燒階段。當(dāng)管道外部油氣爆炸反應(yīng)結(jié)束后，管道內(nèi)部氧氣被消耗殆盡，剩余大量低密度油氣產(chǎn)物，在氣壓作用下，外部空氣及未燃燒完的油氣一同倒吸入管道內(nèi)，此時的燃燒為不完全燃燒，火焰呈現(xiàn)橙紅色的外觀。

2.3.3油氣濃度對油氣爆炸火焰?zhèn)鞑ニ俣鹊挠绊?/p>

火焰?zhèn)鞑ニ俣仁敲枋鲇蜌獗ǔ潭鹊闹匾獏?shù)，圖9反應(yīng)的是在不同油氣濃度下，火焰?zhèn)鞑ニ俣入S時間的變化曲線。

圖9　不同濃度下火焰的傳播速度Fig.9　The flame propagation speed under different concentrations

定義火焰?zhèn)鞑ニ俣葹椋?/p>

式中：xn與xn-1分別代表第n時刻與第n-1時刻火焰鋒面的位置；△t代表這2個時刻的時間間隔；因為時間間隔非常短，所以v近似看做中間時刻的瞬時速度。

從圖9可以看出，在所有濃度下，火焰?zhèn)鞑ニ俣榷际浅尸F(xiàn)先增大，后減小的趨勢。不同濃度下最大火焰?zhèn)鞑ニ俣炔煌?，其中，濃度?.73%時最大火焰?zhèn)鞑ニ俣茸畲?，達(dá)到了102.93 m/s；隨著濃度的增加，火焰?zhèn)鞑ニ俣鹊募铀俣认仍龃?，后減小，濃度為1.73%和1.87%時，火焰?zhèn)鞑ニ俣壬仙羁臁Ｍ瑫r，還可以從圖中看出，在火焰?zhèn)鞑ニ俣冗_(dá)到峰值前，上升較為平穩(wěn)，接近于線性。這是因為火焰在管道內(nèi)傳播時，為層流拉伸狀態(tài)，火焰保持勻速傳播；而在管道外空間爆炸后，能量迅速釋放，火焰的傳播速度隨即迅速下降。

3　結(jié)論

1)半開口管道內(nèi)油氣爆炸可以分為以下4個階段：定容爆炸階段、開口泄流階段、外部爆炸階段和衰減震蕩階段。管道油氣爆炸超壓峰值呈現(xiàn)隨濃度增加先增大后減小的趨勢。升壓速率隨著濃度的增加也呈現(xiàn)先增大后減小的趨勢。在當(dāng)量濃度比下的油氣爆炸升壓速率最快，同時超壓峰值也最大。

2)初始油氣濃度對火焰鋒面?zhèn)鞑ニ俣扔兄@著影響。在當(dāng)量濃度比下，火焰鋒面位置的傳播速度最大，并且在當(dāng)量濃度比下火焰鋒面?zhèn)鞑サ淖畲缶嚯x也更遠(yuǎn)。

3)管道內(nèi)的火焰瞬變行為可以分為4個階段：層流拉伸階段、破壞變形階段、管道外爆炸階段以及回流燃燒階段?；鹧孀畛跻郧蛐坞S后以“指尖型”形態(tài)在管道內(nèi)傳播，經(jīng)過管道端面破口處的擾動后，火焰呈現(xiàn)毛刷狀湍流火焰的形態(tài)?；鹧媸芡牧鲾_動以及Rayleigh-Taylor不穩(wěn)定等的影響，在管道外翻轉(zhuǎn)拉伸形成蘑菇云。

4)初始油氣濃度對火焰?zhèn)鞑バ螒B(tài)以及傳播速度有明顯的影響。對火焰?zhèn)鞑バ螒B(tài)的影響主要體現(xiàn)在破壞變形以及管道外爆炸階段，隨著濃度增加，爆炸半徑先增大后減小，火焰?zhèn)鞑ニ俣瘸尸F(xiàn)相同的變化規(guī)律。

[1]Tomlin G, Johnson D M, Cronin P, et al. The effect of vent size and congestion in large-scale vented natural gas/air explosions[J]. Journal of Loss Prevention in the Process Industries, 2015, 35: 169-181.

[2]Fakandu B M, Andrews G E, Phylaktou H N. Vent burst pressure effects on vented gas explosion reduced pressure[J]. Journal of Loss Prevention in the Process Industries, 2015, 36: 429-438.

[3]Bao Q, Fang Q, Zhang Y, et al. Effects of gas concentration and venting pressure on overpressure transients during vented explosion of methane-air mixtures[J]. Fuel, 2016, 175:40-48.

[4]Ambrosi G, Bosco L, Camosso M E. Experimental and numerical study of methane-air deflagrations in a vented enclosure[J]. Fire Safety Science, 2008, 9(2):1043-1054.

[5]Keenan J J, Makarov D V, Molkov V V. Rayleigh-Taylor instability: modelling and effect on coherent deflagrations[J]. International Journal of Hydrogen Energy, 2014, 39(35):20467-20473.

[6]QUILLATRE P, VERMOREL O, POINSOT T, et al. Large eddy simulation of vented deflagration[J].Industrial & Engineering Chemistry Research, 2013,52(33):11414-11423.

[7]Guo J, Sun X, Rui S, et al. Effect of ignition position on vented hydrogen-air explosions[J]. Explosion & Shock Waves, 2016, 40(45):15780-15788.

[8]Guo J, Wang C, Liu X. Experimental study on duct-vented explosion of hydrogen-air mixtures in a wide range of equivalence ratio[J]. Industrial & Engineering Chemistry Research, 2016, 55(35): 9518-9523.

[9]杜揚, 李陽超, 齊圣,等. 受限空間內(nèi)局部油氣爆燃超壓與火焰?zhèn)鞑ヌ匦訹J]. 后勤工程學(xué)院學(xué)報, 2016, 32(5):33-38.

DU Yang, LI Yangchao, QI Sheng, et al. Characteristics of partial gasoline-air mixture deflagration’s overpressure and flame propagation in confined space[J]. Journal of Logistical Engineering University, 2016, 32(5):33-38.

[10]齊圣，杜揚，梁建軍，等.受限空間油氣爆燃火焰形態(tài)[J]. 爆炸與沖擊, 2016,36(6): 832-838.

QI Sheng, DU Yang, LIANG Jianjun, et al. Flame patterns of gasoline-air mixture deflagration in a confined space[J]. Explosion and Shock Waves, 2016, 36(6): 832-838.

[11]杜揚, 李國慶, 王世茂, 等. 障礙物數(shù)量對油氣泄壓爆炸特性的影響[J]. 化工學(xué)報, 2017, 68(7): 2946-2955.

DU Yang, LI Guoqing, WANG Shimao, et al. Effects of obstacle number on characteristics of vented gasoline-air mixture explosions[J]. CIESC Journal, 2017, 68(7): 2946-2955.

[12]Wen X, Yu M, Liu Z, et al. Large eddy simulation of methane-air deflagration in an obstructed chamber using different combustion models[J]. Journal of Loss Prevention in the Process Industries, 2012, 25(4):730-738.

[13]Wen X, Yu M, Ji W, et al. Methane-air explosion characteristics with different obstacle configurations[J]. International Journal of Mining Science and Technology, 2015, 25(2):213-218.

[14]PARK D J,LEE Y S,GREEN A R.Experiments on the effects of multiple obstacles in vented explosion chambers[J]. Journal of Hazardous Materials,2008,153(1-2):340-350.

[15]Park D J, Lee Y S, Green A R. Prediction for vented explosions in chambers with multiple obstacles[J]. Journal of Hazardous Materials, 2008, 155(1-2):183-192.

[16]Di S V, Di B A, Russo G. Sub-grid scale combustion models for large eddy simulation of unsteady premixed flame propagation around obstacles[J]. Journal of Hazardous Materials, 2010, 180(1-3):71-78.

[17]王世茂, 杜揚, 李國慶, 等. 局部開口受限空間油氣爆燃的超壓瞬變與火焰行為[J]. 化工學(xué)報, 2017, 68(8): 3310-3318.

WANG Shimao, DU Yang, LI Gguoqing, et al. Overpressure transients and flame behaviors of gasoline-air mixture deflagration in confined space with local opening[J]. CIESC Journal, 2017, 68(8): 3310-3318.

[18]杜揚，王世茂，齊圣,等. 油氣在頂部含弱約束結(jié)構(gòu)受限空間內(nèi)的爆炸特性[J]. 爆炸與沖擊, 2017, 37(1): 53-60.

DU Yang, WANG Shimao, QI Sheng, et al. Explosion of gasoline/air mixture in confined space with weakly constrained structure at the top[J]. Explosion and Shock Waves, 2017, 37(1): 53-60.

[19]Pedersen H H, Tomlin G, Middha P, et al. Modelling large-scale vented gas explosions in a twin-compartment enclosure[J]. Journal of Loss Prevention in the Process Industries, 2013, 26(6):1604-1615.

[20]溫小萍, 武建軍, 解茂昭.瓦斯爆炸火焰結(jié)構(gòu)與壓力波的耦合規(guī)律[J]. 化工學(xué)報, 2013,64(10):3871-3877.

WEN Xiaoping, WU Jianjun, XIE Maozhao. Coupled relationship between flame structure and pressure wave of gas explosion[J].CIESC Journal, 2013, 64(10):3871-3877.

[21]Hisken H, Enstad G A, Middha P, et al. Investigation of concentration effects on the flame acceleration in vented channels[J]. Journal of Loss Prevention in the Process Industries, 2015, 36(1):447-459.

[22]Fakandu B M, Mbam C J, Andrews G E, et al. Gas explosion venting: External explosion turbulent flame speeds that control the overpressure[J].Chemical Engineering Transactions, 2016(53):1-6.

[23]Chao J, Bauwens C R, Dorofeev S B. An analysis of peak overpressures in vented gaseous explosions[J]. Proceedings of the Combustion Institute, 2011, 33(2):2367-2374.