尚拓強，楊斯博

（北京航天石化技術(shù)裝備工程有限公司，北京100176）

可燃氣體廣泛存在于石油、化工、煉化等眾多行業(yè)中， 因其泄漏或在非預期條件下燃燒導致的災害事故也是人類生產(chǎn)生活中常見的事故類型之一。 可燃氣體燃燒導致的火災事故，通常又會發(fā)展成爆燃、爆炸等嚴重災害，造成巨大的人身傷害及財產(chǎn)損失。如2007 年7 月11 日， 發(fā)生在山東某化工集團的火災爆炸事故， 因閥門故障導致大量氮氫氣體外泄并著火， 造成9 人死亡，1 人受傷， 直接經(jīng)濟損失450萬元的慘烈后果。又如2003 年9 月發(fā)生在某石化公司常減壓裝置的閃爆事故， 因操作人員誤操作導致瓦斯氣聚集后被點燃而引發(fā)， 事故造成3 人死亡，1人重傷。

工業(yè)生產(chǎn)中的火災爆炸事故經(jīng)常因可燃氣體泄漏、聚集后被引燃導致。 因此，了解可燃氣體火焰在空間內(nèi)形成、發(fā)展、擴散的過程，對抑制火災爆炸的發(fā)展和降低火災爆炸事故帶來的損失具有重要意義。 國內(nèi)外學者也已經(jīng)對由可燃氣體引發(fā)的此類災害進行了廣泛而深入的研究， 旨在盡可能地預防其發(fā)生或減少其危害。

何學超，孫金華等[1]對密閉管道內(nèi)甲烷氣預混火焰的發(fā)展過程進行了實驗研究，采用高速攝影系統(tǒng)，得到了火焰?zhèn)鞑リ嚸娴牡湫吞卣鳎?并分析得出了其原因機理。 WESTBROOK C K[2]等對火焰中烴類燃料層流火焰模型的氧化反應機理進行了研究， 確定了各種反應速率對層流火焰性能的影響。 DENIS V,LUC V[3]等討論了預混湍流燃燒的平均燃燒速率數(shù)值模型， 說明了將數(shù)值模擬作為研究工具用于預混湍流燃燒傳播的方法。盧捷等[4]對煤氣預混氣在封閉管道內(nèi)的火焰加速現(xiàn)象進行了研究， 認為火焰加速是由于火焰前端未燃氣體被前驅(qū)壓縮波加熱和產(chǎn)生的湍流的正反饋造成的。

1 數(shù)值模擬

1.1 物理模型

為簡化計算，建立了長為500 mm，寬為80 mm的封閉二維矩形管道物理模型，如圖1 所示。管道內(nèi)充滿混合均勻的甲烷與空氣， 其中甲烷濃度為10%（vol%），模擬計算初始（t=0 時）在圖示位置處點燃混合氣體。

圖1 物理模型圖

1.2 數(shù)學模型

甲烷氣體燃燒及傳播是一個復雜的過程， 不僅涉及到相關(guān)化學反應，傳播過程更與層流流動、湍流流動及兩種流態(tài)的轉(zhuǎn)變過程密不可分。 為了簡化計算過程， 突出燃燒火焰形態(tài)及火焰發(fā)展過程的研究重點，因此需對該過程的數(shù)學模型進行部分簡化。 本文中采用了可壓縮形式的N-S 方程，Simple 格式求解有限差分方程。 同時采用標準k-ε 方程來計算湍流流動， 其中湍動能k 和耗散速率ε 由標準k-ε 傳遞方程得到：

式中：k—湍動能，J

ε—湍動能耗散率，%

t—時間，s

x—坐標位置，m

u—速度，m/s

μ—黏度，Pa·s

μt—湍流速度，m/s

Gk—由層流速度梯度產(chǎn)生的湍動能，J

Gb—由浮力產(chǎn)生的湍動能，J

YM—過度擴散產(chǎn)生的波動

C—常量，試驗得出

σk—k 方程（式（1））的湍流普朗特數(shù)

σε—ε 方程（式（2））的湍流普朗特數(shù)

Sk，Sε—自定義項

其中，湍流速度 μt由式（3）確定：

式中，Cμ為常量。

湍流燃燒中的燃燒反應速率， 同時受到化學反應速率和流動參數(shù)的影響。 由于Fluent 自帶的渦耗散概念（EDC）模型可以包括詳細的化學反應機理，因此本文采用EDC 模型用以反應甲烷在燃燒過程中的詳細化學反應機理， 反應時間由小渦生存時間以及化學反應本身需要的時間共同控制。 小渦的空間分辨率ε*可由下式得出：

式中：Cξ—容積比率常數(shù)，Cξ=2.137 7

v—運動黏度，m2/s

τ*為細微結(jié)構(gòu)的特征時間尺度，物質(zhì)經(jīng)過τ*開始反應，τ*由下式給出：

式中：Cτ—時間尺度常數(shù)，Cτ=0.408 2

組分守恒方程中的化學反應源項由下式給出：

Yi—化學反應之前的組分i 的質(zhì)量分數(shù)，%

ρ—反應物密度，kg/m3

1.3 數(shù)值計算初始條件

設定管道物理模型四個邊界均為絕熱壁面邊界，且壁面無滑移。

設定所有氣體為理想氣體，氣體比熱隨溫度變化。

設定t0時刻為點火時刻， 此時管道內(nèi)未引燃區(qū)域 T（t0）=300 K，P（t0）=101 kPa，u（t0）=0，v（t0）=0，甲烷的體積分數(shù)為10%。 t0時刻管道內(nèi)已點燃區(qū)域T（t0）=2 000 K，P（t0）=101 kPa，u（t0）=0，v（t0）=0。

2 計算結(jié)果及分析

根據(jù)所建立模型的模擬計算結(jié)果， 得到t0時刻后二維內(nèi)火焰發(fā)展過程如圖2 所示?？梢钥闯觯鸪跻讶紖^(qū)域為高溫氣體區(qū)域， 隨著熱量向未燃燒區(qū)域傳播， 未燃區(qū)域內(nèi)氣體溫度逐漸達到著火溫度，甲烷-空氣混合氣體隨即被引燃。 可以看到20 ms 內(nèi)，當火焰以球面形式向四周傳播時， 已燃區(qū)軸線方向火焰?zhèn)鞑ニ俣瓤?，其他位置傳播速度慢?此外，由于初期火焰陣面表面積較小， 與未燃氣體接觸面積也較小，因此火焰?zhèn)鞑ニ俣日w較慢，火焰為層流燃燒狀態(tài)。 隨著火焰在管道中發(fā)展，燃燒區(qū)域逐漸擴大，火焰陣面表面積也逐漸增大， 火焰?zhèn)鞑ニ俣纫搽S之越來越快。

圖2 火焰發(fā)展過程圖

因密閉管道中火焰發(fā)展空間受限， 管道上下壁面的約束不僅使得火焰形態(tài)發(fā)生變化， 成為不規(guī)則的橢圓形（10 ms 時）。 22 ms 時可以看出，靠近壁面的火焰發(fā)展方向出現(xiàn)偏轉(zhuǎn)， 火焰與已燃燒區(qū)域內(nèi)產(chǎn)物發(fā)生作用，改變此處的流動狀態(tài)，火焰開始轉(zhuǎn)變?yōu)橥牧鳡顟B(tài)，并形成兩個氣體回流區(qū)域。 受此影響，靠近壁面的火焰?zhèn)鞑ニ俣燃涌欤?而火焰軸線處的傳播速度也因受到回流影響而減緩， 火焰?zhèn)鞑ッ娴那蕼p小，逐漸在傳播方向上形成一個近似平面（22 ms時），在28 ms 時，軸線處火焰?zhèn)鞑ニ俣纫呀?jīng)明顯滯后于軸線兩側(cè)火焰。 此時，火焰陣面分層逐漸消失，燃燒整體轉(zhuǎn)變湍流狀態(tài)。在35 ms 時，火焰陣面的凹陷程度達到最大，隨后開始逐漸減小，直至60 ms 時成為一個平面繼續(xù)傳播至管道末端。

管道內(nèi)速度矢量隨時間變化如圖3 所示。由圖可以看出，20 ms 以前， 火焰陣面附近的速度的方向基本一致，火焰以橢圓形向四周傳播。在22 ms 時，靠近上下壁面兩側(cè)的速度方向與靠近軸線方向的速度方向不再相同，形成渦流區(qū)域，導致靠近軸線處傳播慢，兩側(cè)傳播快，火焰陣面逐漸變成平面。 隨著渦流區(qū)域面積增大，兩側(cè)渦流逐漸靠近火焰陣面，進一步發(fā)展后靠近軸線處火焰陣面內(nèi)陷，至35 ms 左右凹陷程度達到最大。 隨后兩個渦流區(qū)開始逐漸遠離火焰陣面，渦流區(qū)面積也逐漸減小，至60 ms 左右火焰陣面變?yōu)橐粋€平面后保持該形態(tài)發(fā)展至管道末端。

圖3 管道內(nèi)的速度矢量圖

可以看出， 火焰受壁面壓迫后產(chǎn)生的渦流是火焰陣面形態(tài)發(fā)生改變的重要原因之一。同時，火焰陣面作為甲烷燃燒的化學反應區(qū)， 與其相關(guān)的可能對火焰發(fā)展產(chǎn)生影響的因素不在本文探討范圍。

3 結(jié)論

本文以可燃氣體甲烷在封閉管道中的發(fā)展過程作為研究對象， 選取適當?shù)奈锢砗蛿?shù)學模型對該過程進行簡化處理，并進行數(shù)值模擬計算，主要得出以下結(jié)論：

（1）管道內(nèi)氣體被引燃后，火焰首先以球面形式向四周傳播， 由于兩側(cè)壁面限制了火焰徑向傳播速度，火焰陣面逐漸被壓縮成橢圓形，至火焰?zhèn)鞑サ缴舷卤诿婧螅诿媪黧w發(fā)展成為湍流狀態(tài)，傳播速度加快，而軸線處火焰受到反向氣流的作用，傳播速度減慢，使火焰陣面前端逐漸由橢圓形變?yōu)槠矫?，進而發(fā)展為凹陷狀。在繼續(xù)傳播一段時間后，火焰陣面恢復為一個平面并持續(xù)傳播至管道末端。

（2）火焰?zhèn)鞑ミ^程中，火焰陣面后方上下兩側(cè)出現(xiàn)兩個渦流區(qū)域，渦流區(qū)逐漸增大并靠近火焰陣面，進一步加大近壁面和軸線處的傳播速度差， 火焰陣面凹陷程度也隨之增大，當凹陷達到最大后，渦流區(qū)面積減小并開始遠離火焰陣面， 凹陷也隨之減小并最終成為一個平面。 說明火焰?zhèn)鞑ブ挟a(chǎn)生的渦流區(qū)域是火焰陣面形態(tài)變化的重要原因之一。

（3）預混的甲烷氣體火焰在封閉管道內(nèi)的傳播過程中會由起初的層流流動逐步轉(zhuǎn)化為湍流流動。