陳 鵬，李艷超，黃福軍，張玉濤

（1.中國礦業(yè)大學（北京）煤炭資源與安全開采國家重點實驗室，北京100083；2.中國礦業(yè)大學（北京）資源與安全工程學院，北京100083）

方孔障礙物對瓦斯火焰?zhèn)鞑ビ绊懙膶嶒炁c大渦模擬＊

陳 鵬1，2，李艷超2，黃福軍2，張玉濤2

（1.中國礦業(yè)大學（北京）煤炭資源與安全開采國家重點實驗室，北京100083；2.中國礦業(yè)大學（北京）資源與安全工程學院，北京100083）

為揭示置障管道內(nèi)甲烷／空氣預混火焰?zhèn)鞑ヌ匦?，運用高速攝影技術(shù)對甲烷／空氣預混火焰的形狀變化和火焰前鋒的速度特性進行實驗，并利用大渦模擬對管道內(nèi)的流場結(jié)構(gòu)進行數(shù)值分析。結(jié)果表明：置障管道內(nèi)依次出現(xiàn)了球形火焰、指尖形火焰及“蘑菇”狀火焰，且“蘑菇”狀火焰出現(xiàn)之后，火焰開始反向傳播；“蘑菇”狀火焰是雙渦旋結(jié)構(gòu)與火焰前鋒面相互作用的結(jié)果，而火焰的反向傳播是由流場中出現(xiàn)逆流結(jié)構(gòu)引起的；障礙物對火焰前鋒有明顯的加速作用；大渦模擬成功再現(xiàn)了實驗中觀察到的火焰形狀、火焰前鋒速度及流場結(jié)構(gòu)，說明大渦模擬適用于置障管道內(nèi)預混火焰?zhèn)鞑ヌ匦缘难芯俊?/p>

爆炸力學；火焰前鋒；大渦模擬；高速攝影技術(shù)；甲烷／空氣預混火焰；方孔障礙物

瓦斯爆炸是發(fā)生在煤礦井下最嚴重的災害之一。在實際礦井巷道中不可避免地存在各種障礙物，如通風設備、支護設施及各種管道等，當火焰與障礙物相互作用時，火焰的形態(tài)、火焰?zhèn)鞑ニ俣取⒘鲌鼋Y(jié)構(gòu)及爆炸壓力等都會發(fā)生顯著變化［1－4］。因此，研究火焰與障礙物的相互作用機制及由此引發(fā)的流場結(jié)構(gòu)變化對于有效預防瓦斯爆炸事故的發(fā)生及降低災害損失有著重大意義。

學者們對置障條件下的瓦斯爆炸機理、傳播特性以及火焰和流場相互作用進行了深入研究［5－6］。A.R.Masri等［7］進行了障礙物孔型對火焰?zhèn)鞑ニ俣扔绊憣嶒?，指出與三角形和圓形障礙物相比，方形障礙物對火焰的加速作用更明顯。R.Hall等［8］分析了障礙物數(shù)量對湍流預混火焰?zhèn)鞑ヌ匦缘挠绊懀l(fā)現(xiàn)爆炸超壓隨障礙物數(shù)量的增加而增大，但存在上限，超過上限后，爆炸超壓隨障礙物數(shù)量的增加而減小。P.Chen等［9］通過改變障礙物阻塞比證明了火焰?zhèn)鞑ニ俣群蛪毫υ鲩L速度隨阻塞比的增大而增加，并指出火焰穿越障礙物后速度的改變是火焰和湍流相互作用的結(jié)果。障礙物能顯著增強管道內(nèi)流場的湍流強度，提高火焰燃燒效率［10］。C.T.Johansen等［11］利用新型的紋影技術(shù)捕捉到了障礙物后的渦旋結(jié)構(gòu)，并指出火焰?zhèn)鞑ヌ匦宰兓母驹蛟谟诨鹧媾c渦旋的相互作用。V.D.Sarli等［12］利用粒子示蹤技術(shù)得到了障礙物后的渦旋結(jié)構(gòu)，實驗證明了火焰與渦旋的相互作用是火焰加速及火焰形態(tài)變化的主要原因。實驗手段很難得到氣體爆燃過程中火焰和流場的精細結(jié)構(gòu)，大渦模擬（large eddy simulation，LES）可以揭示火焰?zhèn)鞑ミ^程中的流動現(xiàn)象，有助于理解層流火焰向湍流火焰轉(zhuǎn)變、燃燒不穩(wěn)定性和火焰形態(tài)變化等［13］。本文中通過實驗和大渦模擬對含障礙物管道內(nèi)甲烷／空氣預混火焰?zhèn)鞑ヌ匦赃M行研究，研究置障管道內(nèi)火焰?zhèn)鞑サ膭恿W特性，通過對比實驗和模擬結(jié)果，驗證大渦模擬的有效性。

1 實驗系統(tǒng)

圖1所示為甲烷／空氣預混火焰?zhèn)鞑嶒炏到y(tǒng)示意圖。該系統(tǒng)由定容燃燒管道、高速攝影系統(tǒng)、混氣系統(tǒng)、高壓點火系統(tǒng)和同步控制系統(tǒng)等組成。實驗管道為一方形直管，長度為500mm，橫斷面尺寸為110mm×80mm。障礙物（厚度4mm）與點火點相距200mm，方孔截面積為50mm×50mm，障礙物的材質(zhì)為阻燃有機玻璃板，實驗前后，方孔玻璃板無結(jié)構(gòu)損壞和形變，障礙物的上下端與固定在管道內(nèi)的凹槽相連接。為提供高速攝影所需光路，管道前后兩側(cè)由石英玻璃板制成。實驗中的點火器和高速攝像系統(tǒng)由同步控制器進行同步控制。高速攝像機的拍攝速度設置為960s－1。實驗氣體的甲烷的體積分數(shù)為7%，管道的初始溫度和壓力分別為298K和0.1MPa?？扇細獗桓吣茳c火器點燃之前，靜置60s，點火電極安裝在管道左端的法蘭上，高能點火器點火能可調(diào)范圍是3～20J。

圖1 實驗系統(tǒng)示意圖Fig.1 Sketch of experimental system

2 數(shù)學模型

2.1 控制方程

隨著計算機性能的提高，計算流體力學在燃燒科學領域應用日益廣泛［14－15］。本文中大渦模擬所采用的控制方程為濾波后的三維瞬態(tài)質(zhì)量守恒方程、動量守恒方程、能量守恒方程和反應進程變量方程［3，16］。大渦模擬放棄對全尺度范圍上渦結(jié)構(gòu)瞬時運動的模擬，只將比網(wǎng)格尺度大的湍流運動通過瞬時N－S方程直接計算，而小尺度渦對大渦運動的影響則通過亞格子模型（sub－gridscale model）進行建模處理。其中，反應進程變量方程控制著甲烷／空氣反應物的化學反應狀態(tài)，通過控制反應進程變量捕捉火焰燃燒面，大渦模型的控制方程為［17］：

式中：LES濾波的參量標注為橫線上標，質(zhì)量權(quán)重濾波的參量標注為波浪線上標；ρ為密度，p為壓力，ui、uj為速度分量，t為時間，應力張量σij由分子黏度μ決定，可表示為，τij為亞網(wǎng)格尺度應力，將其定義為為顯焓，λ為熱導率，亞網(wǎng)格熱焓通量可通過梯度假設近似為：，其中：μSGS為亞網(wǎng)格黏度，PrSGS為亞網(wǎng)格的普朗特數(shù)，cp為定壓比熱，T為溫度。

2.2 初始條件和邊界條件

點火前，管道內(nèi)甲烷／空氣混合氣保持靜止，初溫和初壓分別為298K和0.1MPa。設反應過程變量c為量綱一化的產(chǎn)物質(zhì)量分數(shù)，相應地，c＝0對應為未燃反應物，c＝1對應為已燃產(chǎn)物。將初始化學反應過程變量在整個流體區(qū)域設置為c＝0，并通過將c設置為0.1來實現(xiàn)點火，點火半徑為10mm。管道邊界條件為無滑移、絕熱邊界條件。網(wǎng)格為結(jié)構(gòu)網(wǎng)格，單元特征尺寸為2mm×2mm×2mm。應用ANSYS Fluent軟件作為計算平臺，壓力速度耦合采用SIMPLE算法，對流項采用二階逆風格式，擴散項采用中心差分格式。

3 結(jié)果與討論

3.1 火焰前鋒發(fā)展的高速攝影圖片

圖2所示為置障條件下甲烷／空氣預混火焰?zhèn)鞑サ母咚贁z影圖像。t＝10ms時，由于未受到障礙物和管道側(cè)壁影響，火焰陣面以球形向前傳播；之后，管道側(cè)壁開始阻礙火焰的自由發(fā)展，火焰形狀由球形轉(zhuǎn)變成指尖形，且火焰裙邊在t＝32ms時接觸管道側(cè)壁，火焰燃燒表面積開始減??；t＝42ms時，障礙物開始作用于火焰前鋒，平滑的弧面火焰前沿開始變尖；t＝44ms時，火焰恰好穿過障礙物的方孔；t＝48ms時，穿過方孔的火焰再次發(fā)生變形，并在t＝50ms時形成“蘑菇”狀的火焰前端；50ms之后“蘑菇”狀火焰開始破碎，值得注意的是，“蘑菇”狀靠近上下壁面的部分開始反向傳播，并在t＝53ms時反向穿過方孔結(jié)構(gòu)。從火焰形態(tài)特征而言，置障條件下管道內(nèi)的預混火焰經(jīng)歷了層流向湍流燃燒的轉(zhuǎn)變。

圖2 置障管道內(nèi)甲烷／空氣預混火焰?zhèn)鞑サ母咚贁z影圖像Fig.2 Sequences of high－speed images of premixed methane／air flame propagating in an obstructed duct

3.2 火焰前鋒發(fā)展的大渦模擬

圖3所示為置障條件下甲烷／空氣預混火焰?zhèn)鞑サ拇鬁u模擬，其中燃燒面的反應進程變量c＝0.1。

圖3 甲烷／空氣預混火焰?zhèn)鞑サ拇鬁u模擬Fig.3 Large eddy simulation of premixed methane／air flame propagating in obstructed duct

從圖3中可以發(fā)現(xiàn)，基于大渦模型的數(shù)值模擬成功再現(xiàn)了實驗觀察到的火焰形態(tài)及傳播現(xiàn)象，如球形火焰、指尖形火焰、“蘑菇”狀火焰及其反向流動等。這也證明了大渦模型對于模擬置障條件下甲烷／空氣預混火焰?zhèn)鞑サ目尚行浴?/p>

3.3 流場結(jié)構(gòu)

圖4所示為置障條件下甲烷／空氣預混火焰?zhèn)鞑サ牧鲌鼋Y(jié)構(gòu)。數(shù)值模擬中紅色代表已燃氣體，藍色代表未燃氣體，白色箭頭線代表速度矢量。

圖4 置障管道內(nèi)甲烷／空氣預混燃燒的流場結(jié)構(gòu)Fig.4 Flow field of premixed methane／air flame propagating in an obstructed duct

在t＝42ms之前，管道內(nèi)的氣體流動是層流狀態(tài)，流線分布均勻無交叉；t＝49ms時，未燃區(qū)域出現(xiàn)雙渦旋結(jié)構(gòu)，且渦旋范圍不斷擴大，渦旋中心逐漸向障礙物移動；t＝50ms時，已燃區(qū)的流動方向由原來的正向流動反轉(zhuǎn)成逆向；t＝51ms時，未燃區(qū)域方孔障礙物附近也出現(xiàn)逆流。可以推斷，“蘑菇”狀火焰的形成與雙渦旋結(jié)構(gòu)密切相關(guān)，而燃燒后期的火焰反向傳播是由流場的逆流結(jié)構(gòu)造成的。

3.4 火焰前鋒的動力學變化

圖5所示為火焰前鋒位置隨時間的變化關(guān)系的數(shù)值模擬和實驗結(jié)果?？梢?，t＝42ms時刻是火焰前鋒位置隨著時間變化的轉(zhuǎn)折點，在42ms之前，變化速率較為緩慢，在42ms之后，變化速率突然增大。圖6所示為火焰前鋒速度隨時間的變化特性。從圖6中可以看出，在32ms之前，火焰前鋒速度隨著燃燒表面積的增大而緩慢上升；當t＝32ms時，火焰裙邊觸壁，燃燒表面積突然減小，進而引起火焰?zhèn)鞑ニ俣鹊男》葴p??；當t＝44ms時，火焰前鋒速度達到最大值17m／s，此時火焰前鋒正在穿越障礙物的方孔。

圖5 火焰前鋒位置隨時間的變化特性Fig.5 Histories of flame front surface position

圖6 火焰前鋒速度隨時間的變化特性Fig.6 Histories of flame front surface velocity

4 結(jié) 論

本文中通過實驗和數(shù)值模擬研究了方孔障礙物對管道內(nèi)甲烷／空氣預混火焰?zhèn)鞑ヌ匦缘挠绊?。采用高速攝影技術(shù)揭示了火焰過程中的形狀變化和火焰?zhèn)鞑ニ俣忍匦裕么鬁u模擬對火焰?zhèn)鞑バ袨樽隽诉M一步分析。得到結(jié)論如下：

（1）由于未受到方孔障礙物的影響，前期火焰經(jīng)歷了球形火焰和指尖形火焰兩個階段，指形火焰觸壁之前，火焰燃燒面積快速增大；觸壁之后，火焰燃燒面積迅速減小。臨近方孔時，平滑的弧面火焰前沿開始變尖，并且火焰前沿在穿過方孔障礙物后形成“蘑菇”狀火焰；火焰前沿接近管道尾端時，“蘑菇”狀火焰的靠近上下壁面的部分開始反向傳播，并在t＝53ms時反向穿過方孔結(jié)構(gòu)。

（2）火焰前沿形狀的變化與管道內(nèi)的流場結(jié)構(gòu)密切相關(guān)。通過分析流場結(jié)構(gòu)可以發(fā)現(xiàn)，“蘑菇”狀火焰的形成與雙渦旋結(jié)構(gòu)有直接聯(lián)系，而燃燒后期的火焰反向傳播是由流場中出現(xiàn)的逆流造成的。

（3）方孔障礙物對預混火焰的加速作用顯著，火焰未穿過障礙物之前，火焰裙邊觸壁時刻，火焰?zhèn)鞑ニ俣冗_到5.72m／s；當火焰穿越方孔時，火焰?zhèn)鞑ニ俣仍俅芜_到最大值17m／s。

（4）大渦模擬成功再現(xiàn)了實驗中觀察到的火焰形狀、火焰前鋒速度及流場結(jié)構(gòu)。說明大渦模擬對于模擬置障管道內(nèi)甲烷／空氣預混火焰?zhèn)鞑サ目尚行浴?/p>

［1］Dorofeev S B.Flame acceleration and explosion safety applications［J］.Proceedings of the Combustion Institute，2011，33（2）：2161－2175.

［2］Alharbi A，Masri A R，Ibrahim S S.Turbulent premixed flames of CNG，LPG，and H2propagating past repeated obstacles［J］.Experimental Thermal and Fluid Science，2014，56（7）：2－8.

［3］Johansen C T，Ciccarelli G.Modeling the initial flame acceleration in an obstructed channel using large eddy simulation［J］.Journal of Loss Prevention in the Process Industries，2013，26（4）：571－585.

［4］陳志華，葉經(jīng)方，范寶春，等.方形管內(nèi)楔形障礙物對火焰結(jié)構(gòu)與傳播的影響［J］.爆炸與沖擊，2006，26（3）：208－213.Chen Zhihua，Ye Jingfang，F(xiàn)an Baochun，et al.Effects of a wedge obstacle on flame propagation and its structure［J］.Explosion and Shock Waves，2006，26（3）：208－213.

［5］Ciccarelli G，Dorofeev S.Flame acceleration and transition to detonation in ducts［J］.Progress in Energy ＆Combustion Science，2008，34（4）：499－550.

［6］Kundu S，Zanganeh J，Moghtaderi B.A review on understanding explosions from methane－air mixture［J］.Journal of Loss Prevention in the Process Industries，2016（40）：507－523.

［7］Masri A R，Ibrahim S S，Nehzat N，et al.Experimental study of premixed flame propagation over various solid obstructions［J］.Experimental Thermal and Fluid Science，2000，21（1／2／3）：109－116.

［8］Hall R，Masri A R，Yaroshchyk P，et al.Effects of position and frequency of obstacles on turbulent premixed propagating flames［J］.Combustion and Flame，2009，156（2）：439－446.

［9］Chen P，Li Y C，Huang F J，et al.Experimental and LES investigation of premixed methane／air flame propagating in a chamber for three obstacle BR configurations［J］.Journal of Loss Prevention in the Process Industries，2016，41（5）：48－54.

［10］Wen X P，Yu M G，Liu Z C，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（25）：730－738.

［11］Johansen C T，Ciccarelli G.Visualization of the unburned gas flow field ahead of an accelerating flame in an obstructed square channel［J］.Combustion and Flame，2009，156（2）：405－416.

［12］Sarli V D，Benedetto A D，Russo G，et al.Large eddy simulation and PIV measurements of unsteady flames accel－erated by obstacles［J］.Flow Turbulence and Combustion，2009，83（2）：227－250.

［13］Sarli V D，Benedetto A D，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／2／3）：71－78.

［14］Ibrahim S S，Gubba S R，Malalasekera W，et al.Calculations of explosion deflagration flames using a dynamic flame surface density model［J］.Combustion，Explosion，and Shock Waves，2012，48（4）：393－405.

［15］馬秋菊，張奇，龐磊.巷道壁面與瓦斯爆炸相互作用的數(shù)值模擬［J］.爆炸與沖擊，2014，34（1）：23－27.Ma Qiuju，Zhang Qi，Pang Lei.Numerical simulation on interaction between laneway surface and methane explosion［J］.Explosion and Shock Waves，2014，34（1）：23－27.

［16］Gubba S R，Ibrahim S S，Malalasekera W，et al.Measurements and LES calculations of turbulent premixed flame propagation past repeated obstacles［J］.Combustion and Flame，2011，158（12）：2465－2481.

［17］Zimont V L，Battaglia V.Joint RANS／LES approach to premixed flame modelling in the context of the TFC combustion model［J］.Flow Turbulence and Combustion，2006，77（1）：305－331.

LES approach to premixed methane／air flame propagating in the closed duct with a square－h(huán)ole obstacle

Chen Peng1，2，Li Yanchao2，Huang Fujun2，Zhang Yutao2
（1.State Key Laboratory of Coal Resources and Safety Mining，China University of Mining ＆Technology，Beijing100083，China；2.School of Resource and Safety Engineering，China University of Mining ＆Technology，Beijing100083，China）

Aiming at revealing the characteristics of premixed methane／air flame propagating in an obstructed duct.A 4mm thick obstacle with a square hole of 50mm×50mm was equipped at 210mm from the ignition source.In the experiment，the high－speed video photography was used to study the flame shape changes and flame front speed.In the numerical simulation，the large eddy simulation（LES）was applied to investigating the flow structure.The results demonstrate that the flame－tip successively takes on a spherical，finger and mushroom－like shape，and the flame begins to propagate in the backward direction after the mushroom－like flame appears.The mushroom－like flame can be explained by the interaction of the flame with two vortexes，and the reverse flow emerged in the flow field leads to the backward motion of the premixed flame.The flame speed is accelerated significantly due to the obstacle and the flame tip speed reaches the maximum value of 17m／s when the flame passes through the square hole of the obstacle.The flame shape changes and the flow structure observed in the experiments can be well reproduced in the numerical simulations using the LES model.It is indicated that the LES model can be used to predict the premixed flame propagating in an obstructed duct.

mechanics of explosion；flame front surface；large eddy simulation；high－speed video photography；premixed methane／air flame；square－h(huán)ole obstacle

O381；TD712國標學科代碼：13035

A

10.11883／1001－1455（2017）01－0021－06

（責任編輯 王易難）

2015－05－20；

2015－08－25

國家自然科學基金項目（51274205）；煤炭資源與安全開采國家重點實驗室開放課題項目（SKLCRSM10KFB13）

陳 鵬（1971— ），男，博士，副教授，chenpeng＠cumtb.edu.cn。