陸明飛，叢立新，周軍偉

(哈爾濱工業(yè)大學(xué) 海洋工程學(xué)院，山東 威海 264200)

瓦斯爆燃事故是礦井主要災(zāi)害之一，給煤炭行業(yè)造成了嚴(yán)重的人員財(cái)產(chǎn)損失。阻火器作為一種使火焰發(fā)生淬熄的阻燃裝置，能夠有效抑制礦井下火焰的傳播?；鹧嬖谄鋬?nèi)部微細(xì)通道的淬熄涉及到各種復(fù)雜的現(xiàn)象，如湍流流動(dòng)、傳熱效應(yīng)、流動(dòng)不穩(wěn)定性和邊界層等，一直是該領(lǐng)域研究的熱點(diǎn)和難點(diǎn)。

1932年，HOLM等采用燒嘴法首次試驗(yàn)測(cè)量了不同氣體火焰的淬熄直徑。1997年，周凱元提出了熄滅層理論來(lái)解釋淬熄，由于忽視了流體的黏性和內(nèi)部的化學(xué)反應(yīng)，在入口火焰速度較大時(shí)預(yù)測(cè)模型效果并不理想。2005年，宋占兵基于VC++建立了二維火焰在狹縫中的淬熄模型，研究專注火焰在狹縫中的流動(dòng)狀態(tài)和形狀變化，但模擬條件過(guò)于理想忽視了壁面溫度對(duì)淬熄的影響。2019年，ARIFF等試驗(yàn)探究了流動(dòng)狀態(tài)對(duì)火焰淬熄距離的影響。2020年，WAN等采用高速紋影攝影技術(shù)捕捉到了氫/空氣預(yù)混火焰在多孔板組阻火裝置內(nèi)的3種流動(dòng)現(xiàn)象：“通過(guò)”、“淬熄”和“接近極限”。但由于淬熄機(jī)理復(fù)雜，淬熄過(guò)程迅速，所得結(jié)果在細(xì)節(jié)方面仍有不足。

隨著計(jì)算機(jī)和仿真技術(shù)的發(fā)展，結(jié)果詳細(xì)、可控性強(qiáng)的商用CFD軟件已成為解決火焰淬熄問(wèn)題的主流方法。2014年，溫小萍等模擬了瓦斯爆燃火焰在狹縫中的淬熄得到了與周凱元試驗(yàn)相似的結(jié)果。2017年，孫少辰等基于波紋型阻火器內(nèi)火焰淬熄模擬提出了孔隙率和阻火單元厚度對(duì)火焰?zhèn)鞑サ挠绊憴C(jī)制，張建全等和梁勝龍等則著眼于阻火器的擴(kuò)張腔和擴(kuò)張角。然而模擬研究針對(duì)阻火芯整體，大都忽略了內(nèi)部阻火單元的具體形狀而將其抽象為孔隙率。BERLAD和POTTER很早就已發(fā)現(xiàn)阻火單元的截面形狀對(duì)阻火效果有著重要影響。喻健良、胡春明、李江濤等學(xué)者已對(duì)絲網(wǎng)、微型圓管和平板狹縫等結(jié)構(gòu)較為簡(jiǎn)單的阻火單元進(jìn)行了二維/三維數(shù)值模擬。而應(yīng)用廣泛的波紋型阻火器的阻火單元，研究相對(duì)缺失。且由于其結(jié)構(gòu)復(fù)雜，難以通過(guò)試驗(yàn)取得定量的結(jié)果。

因此，筆者針對(duì)以往研究的不足，選取波紋阻火器中單一阻火單元為研究對(duì)象，基于CFD流體動(dòng)力學(xué)軟件模擬不同工況下三維火焰在其內(nèi)的淬熄過(guò)程，進(jìn)而得出火焰的淬熄機(jī)理。研究主要集中在入口火焰速度和壁面溫度等對(duì)火焰淬熄發(fā)展的影響，包括溫度、化學(xué)反應(yīng)速率等。本研究有助于指導(dǎo)阻火器的安裝設(shè)計(jì)，為瓦斯的抑燃抑爆提供參考。

1 模型及參數(shù)設(shè)定

1.1 波紋型阻火器阻火單元物理模型

圖1為波紋型阻火器的結(jié)構(gòu)，阻火芯由2層平板夾1層波紋板纏繞而成，板間隙類似于一個(gè)個(gè)三棱柱狹縫，即波紋型阻火器的阻火單元，火焰的淬熄在這些狹縫內(nèi)發(fā)生。筆者將這些阻火單元統(tǒng)一簡(jiǎn)化為如圖2所示的三棱柱，其長(zhǎng)40 mm，截面為邊長(zhǎng)1 mm的等邊三角形，選擇波紋板的常用的鋁材作為阻火器材料。以當(dāng)量比為1的甲烷/空氣預(yù)混氣體完全燃燒后所得燃燒產(chǎn)物表示入口處進(jìn)入的火焰，溫度取絕熱燃燒溫度?；谝酝芯?，以1 700 K等溫面代表火焰面的位置。由于入口爆燃火焰速度極快，燃燒至熄滅的過(guò)程短暫，因此模型忽略了熱輻射。壁面設(shè)置為厚度為0的恒溫壁面，這也是研究火焰淬熄時(shí)采用較多的設(shè)定。

圖1 波紋型阻火器結(jié)構(gòu)Fig.1 Structure of crimped ribbon flame arrester

圖2 阻火單元幾何模型Fig.2 Geometrical model of crimped ribbon flame arrester elements

1.2 大渦模擬和燃燒模型

阻火單元內(nèi)火焰的傳播通常為高速可壓縮流。基于這一問(wèn)題，通過(guò)守恒方程經(jīng)過(guò)濾波函數(shù)過(guò)濾，結(jié)合本構(gòu)方程和狀態(tài)方程得到描述火焰淬熄過(guò)程的大渦模擬(LES)控制方程。在燃燒過(guò)程中密度變化較大，為使N-S方程較易封閉，對(duì)變量函數(shù)釆用密度加權(quán)(Favre)過(guò)濾。其中，表示過(guò)濾后的大尺度分量用上標(biāo)“-”表示，F(xiàn)avre濾波過(guò)濾處理用上標(biāo)“～”表示。

控制方程：

(1)

(2)

(3)

(4)

(5)

亞格子應(yīng)力項(xiàng)采用動(dòng)態(tài)Smagorinsky模型求解：

(6)

(7)

(8)

亞格子熱流和亞格子組分流采用梯度輸運(yùn)假設(shè)求解：

(9)

(10)

燃燒模型采用有限速率/渦耗散(Finite-Rate/Eddy-Dissipation)模型封閉組分守恒方程 (5)。該模型同時(shí)考慮到化學(xué)反應(yīng)動(dòng)力學(xué)與湍流混合速率對(duì)燃燒時(shí)化學(xué)反應(yīng)速率的影響，燃燒速率最終取2種機(jī)制影響較小的一方。

(11)

由湍流混合時(shí)間尺度控制的反應(yīng)速率如式(12)所示，取最小值：

(12)

(13)

式中，為指前因子；為溫度指數(shù)；為反應(yīng)活化能,J;′為反應(yīng)物的化學(xué)計(jì)量系數(shù)；為單位質(zhì)量湍動(dòng)能耗散率，m/s；為單元湍動(dòng)能，J/kg；″為生成物的化學(xué)計(jì)量系數(shù)；為物質(zhì)的分子量；為質(zhì)量分?jǐn)?shù)；下標(biāo)為自由標(biāo)，Rec為反應(yīng)物，Pro為生成物；和為經(jīng)驗(yàn)常數(shù)，取4.0，取0.5。

邊界條件和初始條件見(jiàn)表1，2。

表1 邊界條件參數(shù)

表2 初始條件參數(shù)

火焰的傳播和瞬態(tài)溫度密切相關(guān)，取阻火單元中心線上的瞬時(shí)溫度來(lái)驗(yàn)證網(wǎng)格的獨(dú)立性。如圖3所示，選用了10s作為時(shí)間步長(zhǎng)和301 536個(gè)節(jié)點(diǎn)的空間六面體網(wǎng)格，此時(shí)已滿足計(jì)算精度，繼續(xù)提升網(wǎng)格密度計(jì)算結(jié)果改善并不明顯。

圖3 0.5 ms時(shí)不同網(wǎng)格節(jié)點(diǎn)數(shù)和時(shí)間步長(zhǎng)下狹縫中心線上的溫度分布Fig.3 Temperature distribution along slit centerline at0.5 ms with different grid densities and time steps

2 模型驗(yàn)證

2.1 試驗(yàn)案例

數(shù)值計(jì)算結(jié)果與以往的平行板甲烷火焰淬熄試驗(yàn)結(jié)果進(jìn)行了比較。試驗(yàn)裝置布置如圖4所示，由長(zhǎng)4 m、寬0.11 m、高0.02 m的激波管和其他試驗(yàn)設(shè)備組成。在有機(jī)玻璃窗口G段設(shè)有長(zhǎng)192 mm、寬12 mm的平行板狹縫裝置，狹縫間距可調(diào)。本試驗(yàn)介質(zhì)為體積分?jǐn)?shù)為9.50%的甲烷/空氣預(yù)混氣。配氣方法是采用靜態(tài)分壓法，試驗(yàn)時(shí)向激波管內(nèi)充預(yù)混氣采用真空充氣的方法，以提高試驗(yàn)的精確度和可重復(fù)性?？扇?xì)庥秒娀鸹ㄔ诠芏它c(diǎn)燃，點(diǎn)火能量在0.1～25.0 J調(diào)節(jié)。從點(diǎn)火處到激波管窗口之間的距離也可以分段改變，通過(guò)改變這個(gè)距離和點(diǎn)火能量來(lái)獲得火焰進(jìn)入平行板狹縫時(shí)不同的傳播速度?；鹧孢M(jìn)入狹縫前的傳播速度用火焰?zhèn)鞲衅?、電荷放大器及瞬態(tài)波形存儲(chǔ)器測(cè)量獲得。不同狹縫間距下火焰淬熄距離經(jīng)相機(jī)B門拍攝觀測(cè)，考慮信號(hào)傳遞時(shí)間和觀測(cè)玻璃折射率修正后得出。

圖4 甲烷/空氣預(yù)混火焰淬熄試驗(yàn)裝置Fig.4 Experimental device of methane/air premixedflame quenching

2.2 模型驗(yàn)證

由于火焰在微細(xì)通道的傳播存在著共性，采用平行板試驗(yàn)數(shù)據(jù)同樣可以驗(yàn)證筆者數(shù)學(xué)模型的有效性。因此筆者根據(jù)試驗(yàn)設(shè)定和結(jié)果，建立長(zhǎng)40 mm、寬12 mm、狹縫間距0.4～1.2 mm不等的多組幾何模型，入口火焰速度和試驗(yàn)一樣都為90 m/s，內(nèi)部為體積分?jǐn)?shù)9.5%的甲烷/空氣混合氣，其他設(shè)定與本文所建立的數(shù)學(xué)模型相同。如圖5所示，以1 700 K等溫面觀測(cè)到的火焰淬熄距離數(shù)值模擬結(jié)果與試驗(yàn)觀測(cè)到的對(duì)比，結(jié)果基本一致，驗(yàn)證了模型的有效性?；鹧娲阆ň嚯x都隨著狹縫間距的增大而增大，且曲線的斜率也逐漸增大。模擬和試驗(yàn)結(jié)果間略有誤差，這可能是試驗(yàn)條件下火焰的不穩(wěn)定性和難以準(zhǔn)確控制火焰入口速度在90 m/s所導(dǎo)致的。

圖5 試驗(yàn)和模擬所得不同狹縫間距下的淬熄距離Fig.5 Quenching distances of experiment andsimulation with different slit sizes

3 甲烷/空氣預(yù)混氣體火焰淬熄過(guò)程

本部分重點(diǎn)考查了爆燃火焰在波紋型阻火器中的傳播特征，以入口火焰速度60 m/s的甲烷/空氣預(yù)混火焰為例，探究其在邊長(zhǎng)1 mm正三角形截面的阻火單元中的傳播規(guī)律，壁面溫度為300 K。

圖6為波紋型阻火器阻火元件中線面上的溫度分布隨時(shí)間的變化。從圖6可以看出，初始時(shí)刻，溫度場(chǎng)主要受入口火焰的影響，預(yù)混氣體被點(diǎn)燃后還未及時(shí)與火焰面前未燃?xì)怏w充分換熱，因此在阻火單元徑向上出現(xiàn)分層，而在縱向上除靠近壁面處溫度基本相同。隨著入口處火焰的進(jìn)入，熱量通過(guò)熱傳導(dǎo)和活性中心的擴(kuò)散等傳遞到火焰面前方預(yù)混氣體，靠近火焰表面的氣體達(dá)到著火點(diǎn)后首先被點(diǎn)燃，燃燒釋放的熱量繼續(xù)傳遞，預(yù)熱更前方的氣體達(dá)到燃點(diǎn)。就這樣，未反應(yīng)氣體不斷被加熱進(jìn)入火焰面，燃燒產(chǎn)物不斷離開火焰面，循環(huán)往復(fù)，使得火焰向前移動(dòng)。而由于壁面的存在，火焰熱損失隨著火焰的傳播不斷增加，在靠近壁面處形成了低溫的熄滅層。熄滅層厚度隨著熱損失的增加逐漸增加，當(dāng)火焰發(fā)展到一定程度時(shí)交匯，此時(shí)火焰發(fā)生淬熄。

圖6 中線面溫度隨時(shí)間的變化云圖Fig.6 Change of temperature cloud chart onthe midline plane with time

圖7為在阻火元件中心線1 mm處，0.3～0.8 ms時(shí)的縱向溫度變化?？梢钥吹诫S著火焰的傳播，溫度場(chǎng)略有萎縮，這是因?yàn)榛鹧鎰傞_始進(jìn)入阻火單元時(shí)，單元內(nèi)部存在著未燃的預(yù)混氣體，它們?cè)诨鹧孢M(jìn)入后被點(diǎn)燃并釋放熱量。而火焰進(jìn)入三角單元初始時(shí)刻速度較快，時(shí)間較短，因此還未來(lái)得及與壁面進(jìn)行換熱，火焰并不穩(wěn)定。隨著單元內(nèi)預(yù)混氣體燃燒殆盡，壁面的散熱量增加，火焰逐漸穩(wěn)定下來(lái)，火焰從高溫向低溫的過(guò)渡最終趨于穩(wěn)定。其中左側(cè)溫度降低較為明顯；而右側(cè)溫度在一開始略有升高，之后才隨時(shí)間的變化而下降。這是由波紋板阻火器阻火單元的獨(dú)特形狀導(dǎo)致的，溫度場(chǎng)的左端貼近截面三角形的底面；右端則貼近三角形的頂端，同時(shí)受到兩側(cè)壁面夾擊，其通過(guò)壁面散失的熱量大于右側(cè)。且頂部由于碰撞銷毀的自由基數(shù)量更多，如圖8所示，隨著火焰的傳播，三角形底面更加貼近活化中心因而獲得更多燃燒釋放的熱量；同時(shí)頂部空間較小，可供反應(yīng)的可燃?xì)怏w也相對(duì)底部較少。因此右端淬熄效果優(yōu)于左端，溫度更快趨于穩(wěn)定。

圖7 中線面1 mm處各時(shí)刻的溫度分布Fig.7 Temperature distribution of differenttime at 1 mm of the midline plane

圖8 中線面反應(yīng)速率隨時(shí)間的變化云圖Fig.8 Change of reaction rate cloud chart onthe midline plane with time

圖8為在狹縫中線面上的反應(yīng)速率云圖。燃燒化學(xué)反應(yīng)采用單步不可逆甲烷與氧氣燃燒完全反應(yīng)模型進(jìn)行描述。從圖8可以看到，隨著火焰進(jìn)入阻火單元距離的增長(zhǎng)，壁面對(duì)火焰的淬熄效果明顯，靠近壁面位置的燃燒受到抑制，而入口進(jìn)入的高溫已燃?xì)怏w影響逐漸下降，火焰鋒面附近溫度的下降直接導(dǎo)致了燃燒的減弱，反應(yīng)區(qū)域逐漸由一開始與入口等寬的平面逐漸向阻火單元中心收縮，變?yōu)橐粓A曲面，反應(yīng)速率也隨之降低。

火焰的傳播速度由火焰的燃燒速度和流體質(zhì)點(diǎn)的速度組成，燃燒的增強(qiáng)與減弱直接影響火焰鋒面的速度變化。由于黏性作用等，入口火焰在傳播一段距離后流體質(zhì)點(diǎn)速度降低，此時(shí)火焰的傳播速度主要依靠單元內(nèi)本身所含有可燃?xì)怏w燃燒所釋放的熱量來(lái)維持。阻火單元內(nèi)氣體是有限的，燃燒釋放的熱量遠(yuǎn)低于壁面處的熱量損失。隨著火焰的傳播，更少的可燃?xì)怏w被點(diǎn)燃，釋放的熱量更少，燃燒化學(xué)反應(yīng)速率也更低，反應(yīng)區(qū)域逐漸減小，直至火焰熄滅，火焰區(qū)域逐漸萎縮至圓錐形。減小燃燒化學(xué)反應(yīng)區(qū)域有利于火焰的淬熄，而當(dāng)阻火單元截面足夠大即燃燒化學(xué)反應(yīng)區(qū)域足夠大時(shí)，甚至?xí)霈F(xiàn)火焰分離的現(xiàn)象，即混合氣燃燒產(chǎn)生的火焰脫離入口處火焰，此時(shí)火焰很有可能通過(guò)阻火器導(dǎo)致阻火器失效。

4 不同工況下火焰淬熄過(guò)程

4.1 入口火焰速度影響下火焰淬熄過(guò)程

火焰在管道內(nèi)的傳播是一個(gè)加速、加壓的過(guò)程，在管內(nèi)火焰?zhèn)鞑ブ磷杌鹌鲿r(shí)往往伴隨有極快的速度。圖9為不同速度火焰在壁溫300 K的1 mm邊長(zhǎng)正三角形截面阻火單元內(nèi)傳播，0.5 ms時(shí)刻的溫度分布云圖以及溫度沿狹縫中心線的變化?？梢钥吹饺肟谒俣仍酱螅邷貐^(qū)域向前推進(jìn)的距離越大，火焰?zhèn)鞑サ母h(yuǎn)，此時(shí)火焰?zhèn)鞑ゾ嚯x近似與速度大小成正比。且入口速度較大的火焰，高溫向低溫區(qū)的過(guò)渡區(qū)域要稍微厚于入口速度較小的火焰。

圖9 0.5 ms時(shí)不同入口火焰速度下中線面上的溫度云圖及狹縫中心線上的溫度分布Fig.9 Temperature cloud chart on the midline plane anddistribution along slit centerline at0.5 ms with different inlet flame velocities

圖10為0.5 ms時(shí)不同入口火焰速度下，反應(yīng)速率沿狹縫中線的變化?？梢钥吹椒磻?yīng)的活化中心基本與溫度的過(guò)渡層位置一致，都隨著入口速度的增加傳播得更遠(yuǎn)。同時(shí)反應(yīng)速率和反應(yīng)區(qū)域的厚度隨著入口速度的增大而增大。這是由于入口速度增大，不僅加快了高溫區(qū)向前傳播，還使得火焰拉伸越加明顯，同時(shí)火焰面更容易發(fā)生褶皺，導(dǎo)致火焰鋒面的面積更大，與氧氣接觸發(fā)生反應(yīng)的面積也就更大。這加劇了燃燒反應(yīng)的劇烈程度，使得反應(yīng)速率增大；反應(yīng)速率越大，反應(yīng)產(chǎn)生的熱量也就更多，熱量通過(guò)熱傳導(dǎo)等傳遞給前方未反應(yīng)氣體，使其更快達(dá)到著火點(diǎn)，加快了火焰面的傳播速度。可以看出，火焰速度與反應(yīng)速率之間存在一種相互促進(jìn)的關(guān)系，一方面較大的火焰速度促使反應(yīng)速率增大，而反過(guò)來(lái)更大的反應(yīng)速率也加速了火焰的發(fā)展。一般來(lái)說(shuō)，狹窄通道中火焰速度越快，單位時(shí)間傳播距離越遠(yuǎn)，對(duì)前方氣體所做推進(jìn)功也就越大，火焰的動(dòng)能損失也就越大。在通道截面一定時(shí)，其動(dòng)能損失與火焰?zhèn)鞑ゾ嚯x成正比。根據(jù)動(dòng)能定理，隨著入口處火焰速度的增加，火焰?zhèn)鞑ゾ嚯x的增長(zhǎng)速度應(yīng)逐漸放緩。但由于這種動(dòng)態(tài)的相互促進(jìn)關(guān)系的存在，使得火焰燃燒速率提高，加快了火焰的傳播速度，最終呈現(xiàn)出上述傳播距離與入口火焰速度成線性關(guān)系的現(xiàn)象。

圖10 0.5 ms時(shí)不同入口火焰速度下狹縫中心線上的反應(yīng)速率Fig.10 Reaction rate along slit centerline at0.5 ms with different inlet flame velocities

圖11為不同壁面溫度下淬熄距離隨入口火焰速度的變化。入口火焰速度越大，火焰的淬熄距離越遠(yuǎn)。雖然由于火焰速度與燃燒速率間的相互促進(jìn)作用，補(bǔ)足了火焰部分動(dòng)能損失，但隨著反應(yīng)速率的逐漸減小，入口速度較大的火焰動(dòng)能損失相較速度較低的火焰仍然更多。最終使得圖中曲線的斜率逐漸降低，由入口火焰速度的增加導(dǎo)致的淬熄距離增加逐漸放緩。

圖11 不同壁面溫度下淬熄距離隨入口火焰速度的變化Fig.11 Quenching distance of flame with different inlet flamevelocities at different wall temperatures

4.2 壁面溫度影響下火焰淬熄過(guò)程

阻火器的壁面溫度會(huì)隨著其使用環(huán)境和使用次數(shù)的變化而變化，因此有必要研究壁溫對(duì)爆燃火焰在阻火單元內(nèi)傳播的影響。圖12為入口速度60 m/s火焰在不同壁面溫度 (300，450，600，750 K) 的1 mm邊長(zhǎng)正三角形截面阻火單元內(nèi)傳播，0.5 ms時(shí)刻的溫度分布云圖以及溫度沿狹縫中心線的變化?？梢钥吹奖诿鏈囟仍酱蟊诿嫣幍臏囟忍荻仍讲幻黠@，同時(shí)高溫區(qū)域的寬度和火焰向前推進(jìn)的距離越大，說(shuō)明壁面溫度越大越有利于火焰在狹縫中的傳播。與圖9(b)中不同入口火焰速度下觀察到的高溫段不同，圖12(b)中火焰的高溫段幾乎完全重合，僅在溫度梯度變化處出現(xiàn)分離。說(shuō)明不同壁面溫度下，火焰?zhèn)鞑ゾ嚯x的增加主要是高溫向低溫的過(guò)渡區(qū)域厚度的增加。

圖12 0.5 ms時(shí)不同壁面溫度下中線面上的溫度云圖和狹縫中心線上的溫度分布Fig.12 Temperature cloud chart on the midline plane anddistribution along slit centerline at 0.5 ms with differentwall temperatures

圖13為不同壁面溫度下的反應(yīng)速率云圖。隨著壁溫的增加，反應(yīng)活化中心更加靠前。反應(yīng)活化中心的最大反應(yīng)速率基本都在0.8 kmol/(m·s)左右，且反應(yīng)速率在0.35 kmol/(m·s)以上的區(qū)域面積基本一致，靠近壁面兩側(cè)的反應(yīng)速率增加較為明顯。這導(dǎo)致火焰?zhèn)鞑コ跗诒诿鏈囟葘?duì)火焰的燃燒速度的影響并不大。

圖13 0.5 ms時(shí)不同壁面溫度下中線面上的反應(yīng)速率云圖Fig.13 Reaction rate cloud chart on the midlineplane at 0.5 ms with different wall temperatures

圖14為狹縫中線1 mm處0.5 ms時(shí)不同壁面溫度下的縱向溫度變化?？梢钥吹?，雖然狹縫兩側(cè)受壁面影響，溫度各不相同，但在縱向上0.2～0.6 mm的部分溫度基本保持一致。這也就解釋了為何壁面溫度在火焰?zhèn)鞑ミ^(guò)程中對(duì)反應(yīng)活化中心的影響并不大。雖然壁面溫度的增加有利于未燃反應(yīng)氣體的預(yù)熱，但其值仍然遠(yuǎn)小于傳播過(guò)程中甲烷/空氣爆燃火焰的溫度，預(yù)熱效果有限。而高溫的已燃?xì)怏w的熱傳導(dǎo)和燃燒釋放的熱量巨大，這是造成未燃?xì)怏w預(yù)熱和火焰?zhèn)鞑サ闹饕颉?/p>

圖14 0.5 ms時(shí)不同壁面溫度下中線面1 mm處溫度分布Fig.14 Temperature distribution at 1 mm of the midlineplane with different wall temperatures at 0.5 ms

圖15 不同入口火焰速度下淬熄距離隨壁面溫度的變化Fig.15 Quenching distance of flame with different walltemperatures at different inlet flame velocities

圖15為不同入口火焰速度下淬熄距離隨壁面溫度的變化。淬熄距離隨著壁面溫度的增加而增加，且入口火焰速度越大增加越明顯。較高的壁面溫度降低了壁面的散熱能量，不利于通過(guò)壁面將燃燒釋放的熱量傳遞出去，這會(huì)導(dǎo)致熄滅層厚度減小，火焰的淬熄距離增加，最終導(dǎo)致火焰穿過(guò)阻火器，造成阻火器失效。

5 結(jié) 論

(1)火焰在波紋型阻火器阻火單元內(nèi)的傳播受到入口處火焰和內(nèi)部燃燒化學(xué)反應(yīng)的共同作用。隨著火焰的傳播，燃燒化學(xué)反應(yīng)逐漸占據(jù)主導(dǎo)。

(2)火焰速度和燃燒化學(xué)反應(yīng)速率間存在相互促進(jìn)關(guān)系，入口火焰速度150 m/s時(shí)的最大燃燒反應(yīng)速率是30 m/s時(shí)的2倍。入口火焰速度越小，淬熄距離越短，但影響逐漸降低。

(3)壁面溫度改變對(duì)化學(xué)反應(yīng)速率影響較小,實(shí)驗(yàn)條件下不同溫度下最大燃燒反應(yīng)速率都在最大反應(yīng)速率基本都為0.8 kmol/(m·s)。提高壁面溫度將阻礙壁面散熱，使得熄滅層厚度減小，淬熄過(guò)程延遲。