李書明，戴文元，程關(guān)兵

（中國民航大學航空工程學院，天津 300300）

管道中丙烷/空氣預(yù)混火焰?zhèn)鞑ヌ匦缘臄?shù)值模擬

李書明，戴文元，程關(guān)兵

（中國民航大學航空工程學院，天津 300300）

選用標準k-ε湍流模型,采用4步反應(yīng)機制來?；榈难趸?利用Fluent軟件研究了常溫常壓下丙烷/空氣預(yù)混氣體在光滑的方管道中火焰的傳播特性,獲得火焰?zhèn)鞑ニ俣取⒒鹧娼Y(jié)構(gòu)、表面積、火焰到達位置和時刻沿管道變化情況。結(jié)果表明:火焰?zhèn)鞑ニ俣认仍黾雍笾饾u減小。在離點火端50 mm處到達最大值,之后逐漸降低。在t=21 ms時,火焰開始逐漸呈現(xiàn)典型的郁金香結(jié)構(gòu)。火焰在初始加速階段的主要物理機理是前方未燃氣體受到前驅(qū)壓縮波作用而被加熱和壓縮的正反饋微分加速機制,之后在反射壓縮波影響下,火焰?zhèn)鞑ニ俣嚷杂薪档?。通過與實驗結(jié)果對比,數(shù)值模擬結(jié)果正確。

丙烷;火焰?zhèn)鞑ニ俣?k-ε湍流模型;4步反應(yīng)機制

在工業(yè)生產(chǎn)和日常生活中，可燃性氣體（如丙烷和甲烷）得到較為廣泛的應(yīng)用。但可燃性氣體在運輸、儲存、加工和使用過程中，在某些條件下極易發(fā)生意外爆炸。這些現(xiàn)象往往是因為可燃性氣體在被點燃后，在工業(yè)現(xiàn)場的機械設(shè)備、測試儀器、管道和建筑物等障礙物作用下，火焰的形狀、燃燒速率、傳播速度以及對周邊環(huán)境壓力有顯著變化。所有這些因素直接影響到可燃氣體災(zāi)害事故發(fā)生和演變的整個過程。在某些極端條件下火焰被不斷加速，甚至可能導(dǎo)致爆炸甚至爆轟的產(chǎn)生，造成重大的財產(chǎn)損失和人員傷亡[1-2]。因此，研究可燃氣體預(yù)混火焰?zhèn)鞑ヌ匦栽诳扇夹詺怏w災(zāi)害預(yù)防方面有著重要的實際意義。

目前，國內(nèi)外在可燃氣體預(yù)混火焰?zhèn)鞑ノ锢頇C理方面開展了一定的研究工作。在國外，F(xiàn)airweather[3]等針對甲烷/空氣預(yù)混氣體在半封閉管道內(nèi)開展燃燒爆炸實驗研究。發(fā)現(xiàn)火焰沿半封閉管的傳播大體上屬于層流，只有在爆炸實驗的最后階段才形成了快速湍流燃燒，并產(chǎn)生激波和超壓。Elicer-Cortes[4]等通過紋影法獲得氫氣/空氣預(yù)混火焰陣面結(jié)構(gòu)的變化過程，并建立了在實驗中測得的燃燒速度和混合氣體化學當量比的關(guān)系。Parsinejad[5]等采用陰影法建立了燃燒溫度和混合氣體化學當量比的函數(shù)關(guān)系。在國內(nèi)，盧捷[6]等研究了兩端密封管道中煤氣/空氣火焰?zhèn)鞑ヒ?guī)律及其加速機理，發(fā)現(xiàn)火焰加速是因前方未燃氣體被前驅(qū)壓縮波加熱和障礙物誘導(dǎo)的湍流區(qū)對燃燒過程的正反饋機制造成的。林柏泉[7]等開展了瓦斯爆炸過程中火焰?zhèn)鞑ヒ?guī)律及其加速機理的實驗研究。發(fā)現(xiàn)障礙物對瓦斯爆炸過程中火焰?zhèn)鞑ヒ?guī)律具有重要影響。障礙物能迅速提高瓦斯爆炸過程中火焰的傳播速度，這種火焰加速作用的機理可歸功于障礙物誘導(dǎo)的湍流區(qū)對燃燒過程的正反饋作用。周凱元[8]等開展了丙烷/空氣爆燃波的火焰面在直管道中加速運動的規(guī)律及其影響因素的實驗研究。探討了管道直徑和點火能量的變化以及當管道內(nèi)有障礙物時對火焰加速度的影響。

基于此，本文選用標準k-ε湍流模型，采用4步反應(yīng)機制?；榈难趸?，利用Fluent軟件對常溫常壓下丙烷/空氣預(yù)混氣體在光滑的方管道中火焰?zhèn)鞑ヌ匦缘臄?shù)值模擬研究，以期對丙烷在光滑管道中的加速機理做出一定的探究。

1 數(shù)值模擬

本文在光滑管道內(nèi)丙烷/空氣預(yù)混火焰?zhèn)鞑ヌ匦赃M行數(shù)值模擬過程中，做如下簡化假設(shè)：

1）管道內(nèi)的氣體為理想氣體，滿足理想氣體狀態(tài)方程；

2）點火前空氣和可燃氣充分混合，處于常溫常壓，并且處于靜止狀態(tài)；

3）混合氣體的狀態(tài)參數(shù)，如導(dǎo)熱系數(shù)、比熱容等，符合混合規(guī)則；

4）預(yù)混氣體燃燒為不可逆反應(yīng)；

5）把管道壁面視為剛性壁面，不考慮氣體流動與壁面的耦合作用。

1.1 控制方程

圖1 爆燃過程的雙波三區(qū)結(jié)構(gòu)Fig.1 Structure of double waves and three zones of flame front

眾所周知，爆燃波是一種帶有化學反應(yīng)的亞音速波，即燃燒波。燃燒波在傳播過程中形成了雙波三區(qū)結(jié)構(gòu)[1]，如圖1所示。其中，雙波是爆燃波和前驅(qū)沖擊波，三區(qū)分別是未燃區(qū)（0區(qū)）、燃燒區(qū)（1區(qū)）和已燃區(qū)（2區(qū)）。各區(qū)參數(shù)定義如下：e為比內(nèi)能，p為壓力，u為粒子速度，c為音速，T為溫度，γ為等熵指數(shù)。

跨過沖擊波陣面，兩側(cè)氣體狀態(tài)發(fā)生改變，但仍遵守質(zhì)量守恒定律、動量守恒定律和能量守恒定律[9]。

質(zhì)量守恒方程為

動量守恒方程為

能量守恒方程為

理想氣體狀態(tài)方程為

理想氣體本構(gòu)方程為

對多組分混合理想氣體而言，氣體常數(shù)為

式（6）中

氣體定壓比熱容為

本文利用Premix中的TWOPNT函數(shù)求解上述控制方程，采用有限差分法來離散微分方程[10]。

1.2 燃燒模型

丙烷-空氣預(yù)混燃燒總反應(yīng)方程式為

本文采用Hautman[11]提出的一個多步的準總包機理，用以下4步反應(yīng)來?；榈难趸?/p>

1.3 模型初始條件

本文只研究預(yù)混火焰?zhèn)鞑サ某跏茧A段（火焰?zhèn)鞑ニ俣鹊陀谖慈細怏w音速），湍流模型選用通用性最好的標準k-ε模型。采用的物理模型為長方形管道，長×寬分別為500 mm×80 mm，左側(cè)壁面中央點火，點火能量為1 mJ，管兩端均封閉。網(wǎng)格數(shù)為100×10；管道等效直徑為100 mm；當量比為1；丙烷-空氣混合濃度為6%；初始溫度為300 K；初始壓力為1 atm；湍流系數(shù)Ct=30；層流火焰拉伸系數(shù)Fs=0.15；k-ε模型常數(shù)為：Cμ=0.09，C1=1.44，C2=1.79，σk=1.0，σξ=1.3，渦流動能轉(zhuǎn)換系數(shù)Fk=0.1，選用丙烷的高位熱值46.357 mJ/kg，丙烷層流火焰速度為0.44 m/s。

2 數(shù)值模擬結(jié)果與分析

本文利用Fluent軟件對丙烷/空氣預(yù)混火焰在光滑管道中的傳播過程進行了數(shù)值模擬研究。圖2表示預(yù)混火焰沿管道軸線位置傳播過程中溫度變化云圖。由圖2可觀察到火焰加速過程中火焰形狀、結(jié)構(gòu)和表面積的變化情況。圖3表示火焰?zhèn)鞑ニ俣葀沿管道位置的變化情況，圖4表示火焰前鋒面到達管道位置、時刻變化。應(yīng)該指出的是，本文所選光滑管道長度為500 mm，測量范圍約為350 mm。通過記錄火焰到達管道中心軸線位置和時刻，可求得火焰?zhèn)鞑ニ俣葀，即

圖2 溫度云圖Fig.2 Temperature nephogram

其中：Xi-1、Xi分別為燃燒波到達管道軸線位置；ti-1、ti分別為燃燒波到達管道軸線不同位置的時刻。

圖3 火焰?zhèn)鞑ニ俣葀沿管道中心軸向位置X的變化圖Fig.3 Variations of flame propagation velocity v along central axial position X in channel

圖4 火焰前鋒面到達管道軸向位置X和時刻tFig.4 Variations of central axial position X and time-of-arrival t of flame front along channel

因此，由圖2～圖4可知：火焰?zhèn)鞑ニ俣认戎鸩皆黾?，在X=50 mm處到達最大值，之后逐步減小。預(yù)混氣體丙烷/空氣在點燃后，起初管道壁面對火焰阻礙作用不明顯，火焰呈現(xiàn)球形狀，如圖2（a）所示。此后，因管左端面的限制，火焰形狀逐漸變?yōu)榘肭蛐螤?，并以大約為10～15 m/s的速度向管右端傳播，如圖2（b）所示。之后因管道上下壁面的限制，半球形火焰逐漸演變?yōu)橹感位鹧?，如在t=3 ms時。在離點火端50 mm處火焰?zhèn)鞑ニ俣鹊竭_最大值，約為36 m/s，此后火焰?zhèn)鞑ニ俣乳_始降低。在t=1～3 ms這個階段（如圖2（a）～圖2（c）所示），火焰?zhèn)鞑ニ俣戎饾u增加，這主要是因為預(yù)混氣體丙烷/空氣燃燒產(chǎn)物密度下降，體積增加，燃燒產(chǎn)物不斷膨脹，火焰表面積不斷增加，從而提高了燃燒效率和火焰?zhèn)鞑ニ俣?。而火焰?zhèn)鞑ニ俣炔粩嘣黾?，在火焰前鋒面前形成微弱的壓縮波逐漸提高了未燃氣體的溫度和壓力，進一步提高了燃燒效率，火焰相應(yīng)地繼續(xù)加速并不斷向前發(fā)射壓縮波，此壓縮波的速度也逐漸提高，這是典型的火焰前未燃氣體受到前驅(qū)壓縮波作用而被加熱和壓縮的正反饋微分加速機制。

此后，因管道長度有限，火焰在管右端面反射的微弱壓縮波影響下逐步減速?；鹧媲颁h面變得越來越平，火焰表面積減小。在t=7 ms時（如圖2（d）所示），可觀察到火焰前鋒面表現(xiàn)為一個近似平面，此時火焰?zhèn)鞑ニ俣燃s為16 m/s。隨著位于火焰前鋒面的反射壓縮波強度不斷增強以及在管道上下壁面處的湍流作用影響下，火焰前鋒面開始凹向管左端面，呈現(xiàn)出反弓形火焰結(jié)構(gòu)，這就是典型的郁金香火焰，如圖2（e）和圖2（f）所示。由圖3和圖4可知，此階段火焰?zhèn)鞑ニ俣戎饾u從16 m/s減小到10 m/s左右。因此，燃燒波前的壓縮波在管右端反射波與火焰陣面?zhèn)鞑シ较蛳喾矗档土撕笳叩膫鞑ニ俣?。盡管，火焰表面積同時也不斷增加，但管右端反射的壓縮波在此階段是火焰?zhèn)鞑ニ俣冉档偷闹饕蛩亍?/p>

圖5表示文獻[12]通過陰影法觀察到丙烷/空氣在90°彎管中郁金香結(jié)構(gòu)火焰。因此，與本文圖2（e）和圖2（f）進行比較，數(shù)值模擬結(jié)果和實驗結(jié)果大體吻合。

圖5 郁金香結(jié)構(gòu)火焰示意圖Fig.5 Tulip structure of flame front

3 結(jié)語

本文通過選用標準k-ε湍流模型，采用4步反應(yīng)機制來?；榈难趸肍luent軟件對常溫常壓下丙烷/空氣預(yù)混氣體在光滑的方管道中火焰?zhèn)鞑ヌ匦缘臄?shù)值模擬研究，獲得了火焰?zhèn)鞑ニ俣?、火焰結(jié)構(gòu)、表面積、火焰到達位置和時刻沿管道的變化情況。結(jié)果表明：

1）火焰?zhèn)鞑ニ俣认仍黾雍笾饾u減小。在離點火端50 mm處到達最大值，之后逐漸降低。

2）在t=21 ms時，火焰開始逐漸呈現(xiàn)典型的郁金香結(jié)構(gòu)。

3）火焰在初始加速階段的主要物理機理是前方未燃氣體受到前驅(qū)壓縮波作用而被加熱和壓縮的正反饋微分加速機制，之后在反射壓縮波影響下，火焰?zhèn)鞑ニ俣嚷杂薪档汀?/p>

4）通過與文獻[12]實驗結(jié)果對比，數(shù)值模擬結(jié)果正確。

[1]關(guān)文玲，蔣軍成.我國化工企業(yè)火災(zāi)爆炸事故統(tǒng)計分析及事故表征物探討[J].中國安全科學學報，2008，18（3）：103-107.

[2]BJERKETVEDT D，BAKKE J R，WINGERDEN K V.Gas explosion handbook[J].Journal of Hazardous Materials，1997，52：1-150.

[3]FAIRWEATHER M，HARGRAVE G K，IBRAHIM S S，et al.Studies of premixed flame propagation in explosion tubes[J].Combustion and Flame，1999，116：504-518.

[4]ELICER-CORTES J C，RUZ C，HERNANDEZ R H，et al.Observation of preferred instability modes in a mechanically excited thermal plume using Schlieren visualizations[J].International Communication in Heat and Mass Transfer，2005，32：360-370.

[5]PARSINEJAD F，ARCARI C，METGHALCHI H.Flame structure and burning speed of JP-10 air mixtures[J].Combustion Science and Technology，2006，178：975-1000.

[6]盧 捷，寧建國，王 成，等.煤氣火焰?zhèn)鞑ヒ?guī)律及其加速機理研究[J].爆炸與沖擊，2004，24（4）：305-311.

[7]林柏泉，張仁貴，呂恒宏.瓦斯爆炸過程中火焰?zhèn)鞑ヒ?guī)律及其加速機理的研究[J].煤炭學報，1999，24（l）：56-59.

[8]周凱元，李宗芬.丙烷-空氣爆燃波的火焰面在直管道中的加速運動[J].爆炸與沖擊，2000，20（2）：137-142.

[9]李宇紅，祁海鷹，張宏武.甲烷預(yù)混燃燒火焰的詳細數(shù)值模擬[J].工程熱物理學報，2002，23（1）：119-122.

[10]劉 慧.多孔介質(zhì)內(nèi)預(yù)混氣體燃燒的實驗和數(shù)值模擬研究[D].沈陽：東北大學，2010.

[11]HAUTMAN D J，DRYER F，SCHUG K P，et al.A multiple-step overall kinetic mechanism for the oxidation of hydrocarbons[J].Combustion Science and Technology，1981，27：31-43.

[12]何學超.丙烷-空氣預(yù)混火焰在90°彎曲管道內(nèi)傳播特性的實驗和模擬研究[D].合肥：中國科學技術(shù)大學，2010.

（責任編輯：楊媛媛）

Numerical simulation study of propane/air premixed flame characteristics in square channel

LI Shu-ming，DAI Wen-yuan，CHENG Guan-bing
（College of Aeronautical Engineering，CAUC，Tianjin 300300，China）

A numerical simulation study on the characteristics of premixed flame of propane/air is carried out by software Fluent in a smooth square channel at normal pressure and temperature conditions.Both standardized k-ε turbulence model and four-step kinetic mechanisms of propane-oxidation are initiated.Flame propagation characteristics such as flame velocity，structure，surface area，position and time-of-arrival can be obtained along the channel.The main conclusions are：the flame propagation velocity first increases and then decreases. The flame velocity achieves maximum at the distance of 50 mm from the ignition point.When t=21 ms，a typical tulip structure flame can be observed.The positive-feedback differential loop resulted from the preheating and compressing unburnt mixtures by the precursor shock ahead of the flame front can be considered as the mainly mechanism of flame propagation in the channel.Subsequently，the flame propagation velocity decreases under the influence of compression waves reflected from the channel terminus.In contrast with the experimental results，the obtained numerical simulation results are reasonable.

propane；flame propagation velocity；k-ε turbulence model；four-step kinetic mechanisms

TK19

:A

:1674-5590(2015)02-0053-04

2013-09-24；

:2014-01-18

:中央高?；究蒲袠I(yè)務(wù)費專項（ZXH2012J001）

李書明（1964—），男，河北邢臺人，教授，博士，研究方向為航空發(fā)動機結(jié)構(gòu)與修理、動力燃燒.