黃濤 范進(jìn) 丁建國
(南京理工大學(xué),江蘇南京 210094)
封閉管道內(nèi)氣體燃燒火焰?zhèn)鞑ヌ匦匝芯?/p>
黃濤 范進(jìn) 丁建國
(南京理工大學(xué),江蘇南京 210094)
為了研究可燃?xì)怏w在管道內(nèi)的燃燒規(guī)律,選取甲烷與空氣預(yù)混氣體為可燃?xì)怏w,建立內(nèi)置障礙物情況不同的封閉管道模型,運(yùn)用CFD計算流體動力學(xué)方法,對管道內(nèi)氣體燃燒過程進(jìn)行了數(shù)值模擬。監(jiān)測了火焰面的形狀、位置、速度的變化;分析了擋板位置、數(shù)量以及阻塞率對火焰面?zhèn)鞑ミ^程的影響。研究表明:火焰面?zhèn)鞑ニ俣入S擋板數(shù)量、阻塞率的增加而增大,擋板和點(diǎn)火區(qū)的距離與管徑的比值在0.5~3范圍內(nèi)時,火焰面?zhèn)鞑ニ俣冗_(dá)到最大值。研究為工程結(jié)構(gòu)的防爆設(shè)計提供依據(jù)。
封閉管道 氣體燃燒 障礙物 計算流體動力學(xué) 火焰?zhèn)鞑?/p>
2014年6月3日,重慶市萬圣經(jīng)濟(jì)開發(fā)區(qū)硯石臺煤礦發(fā)生重大瓦斯爆炸事故,造成22人死亡;8月1日,臺灣高雄前鎮(zhèn)區(qū)發(fā)生丙烯氣體泄漏并爆炸,造成28人死亡,302人受傷。無論是工業(yè)生產(chǎn),還是日常生活中,可燃?xì)怏w安全問題已引起人們的高度關(guān)注。
氣體爆炸問題的研究方法,主要有解析法,實(shí)驗(yàn)方法和數(shù)值方法。解析法中,常用的解析模型有等溫模型、絕熱模型和一般模型[1]。模型對實(shí)際情況進(jìn)行了較大的簡化,分析存在一定的誤差。相比實(shí)驗(yàn)方法,數(shù)值方法成本低,安全性高,已經(jīng)被許多研究人員采用。采用數(shù)值方法主要有集中參數(shù)方法(LMP)和計算流體動力學(xué)方法(CFD)。CFD方法能夠反映變量在空間的分布,可以適應(yīng)形狀復(fù)雜的情況,得到了廣泛的應(yīng)用。運(yùn)用CFD方法,王志榮等[2]模擬了連通容器內(nèi)氣體的爆炸,結(jié)果表明連通管后的容器內(nèi)產(chǎn)生預(yù)壓縮,湍流導(dǎo)致爆炸壓力增加;李帆等[3]模擬了光滑管道內(nèi)氣體的爆炸,得到了不同時刻管內(nèi)壓力、湍流動能的分布;畢明樹等[4]對管道內(nèi)可燃?xì)獗ㄟM(jìn)行了一維數(shù)值模擬。
實(shí)際上,輸氣管道在使用過程中,可燃?xì)庵泻械碾s質(zhì)將在管道內(nèi)沉積、聚集,導(dǎo)致管道流通面積減小。在長徑比L/D較大(L/D>>5)且內(nèi)部有較密集的隔板和設(shè)備等障礙物的密封體內(nèi),氣體燃燒過程更加復(fù)雜[1、5-6]。為更好的對管道及其圍護(hù)結(jié)構(gòu)進(jìn)行抗爆設(shè)計,有必要對不同阻塞率下管道內(nèi)火焰?zhèn)鞑ヌ匦赃M(jìn)行研究。
本文利用計算流體動力學(xué)軟件Fluent對內(nèi)徑為20mm,長300mm(L/D=15)的管道內(nèi)甲烷與空氣預(yù)混氣體燃燒過程進(jìn)行數(shù)值模擬。
建立內(nèi)徑為20mm,長為300mm的管道二維模型,網(wǎng)格大小為1mm×1mm,共約6000個網(wǎng)格。四周邊界條件為固定壁面,無滑移。由于氣體爆炸發(fā)生時間較短,故在模擬過程中忽略熱損失,假設(shè)壁面為絕熱。管道內(nèi)充滿甲烷與空氣的混合氣體,甲烷含量為10%,混合氣體溫度為300K。化學(xué)反應(yīng)方程式為:
圖1 管截面示意圖
表1 管道中擋板設(shè)置位置及阻塞率
CH4+2O2-CO2+2H2O
以管道左端面中心為圓點(diǎn),設(shè)置半徑為3mm的半圓形點(diǎn)火區(qū),點(diǎn)火區(qū)溫度為1800K。如圖1中,紅色區(qū)域所示。
本文共建立7個管道模型。0號為內(nèi)壁光滑管道模型,管內(nèi)無擋板。其他模型,分別在管道內(nèi)距離點(diǎn)火區(qū)l處,設(shè)置1道或2道擋板,對管道形成一定的阻塞。管道截面如圖1所示。具體情況見表1。
3.1 內(nèi)壁光滑管道內(nèi)氣云燃燒模擬
圖2 火焰面位置及形狀變化
圖3 火焰面位置、速度隨時間的變化關(guān)系
圖4 0、1、2號模型火焰面速度隨時間的變化關(guān)系
內(nèi)壁光滑管道內(nèi)氣體模擬燃燒過程中,火焰面位置及形狀變化如圖2所示。點(diǎn)火初期,火焰面近似成半圓形向外擴(kuò)散,由于側(cè)壁溫度與環(huán)境溫度相同,為300K,火焰面沿軸向延伸,距離壁面較近處還未燃燒。0.5ms時刻,已燃區(qū)域逐漸飽滿,火焰面成平面,化學(xué)反應(yīng)主要集中于已燃區(qū)域右側(cè)。
由于湍流作用,火焰面有一定的厚度,繼續(xù)向前傳播。2ms時刻,火焰面形成指向已燃區(qū)域的尖端,形成“郁金香”火焰面,火焰發(fā)展為明顯的湍流,火焰面厚度進(jìn)一步增加。此后,火焰面大致維持該形狀,逐漸向右傳播,至所有氣體燃燒完畢。
每隔0.1ms在FLUENT中監(jiān)測管道中心線上火焰面中心位置,得到火焰面位置與時間的關(guān)系曲線。對該曲線進(jìn)行擬合并求導(dǎo),可得火焰面中心的速度與時間關(guān)系。火焰面的位置、速度與時間關(guān)系如圖3所示。
點(diǎn)火初期,點(diǎn)火溫度較高,由阿累尼烏斯方程[1]可知,化學(xué)反應(yīng)速率較大,火焰面以較快速度向前傳播,隨著燃燒的進(jìn)行,火焰面速度逐漸減小。當(dāng)火焰面以穩(wěn)定的“郁金香”形向前傳播時,一方面,因?yàn)楸O(jiān)測點(diǎn)的局限性以及湍流的存在,造成監(jiān)測數(shù)據(jù)不能準(zhǔn)確反映整個火焰面的位置,引起一定的誤差;另一方面,由于受到側(cè)壁及管道右端面反射的超壓影響,傳播速度呈現(xiàn)一定的波動性。
圖5 燃燒腔內(nèi)已燃?xì)怏w形狀(上)與燃燒腔超壓(下)
圖6 火焰斷層示意圖
圖7 3、4、5號管火焰面速度隨時間的變化關(guān)系
圖8 4、6號管火焰面速度隨時間的變化關(guān)系
3.2 阻塞率對火焰?zhèn)鞑ニ俣鹊挠绊?/p>
分別監(jiān)測0、1、2號模型火焰面中心位置。記錄火焰面中心速度與時間關(guān)系,如圖4所示。在有擋板的管道中,火焰面?zhèn)鞑ニ俣仍诒诿娣瓷涑瑝旱淖饔孟乱廊挥胁▌樱畲笏俣扰c平均速度均遠(yuǎn)大于光滑管道。擋板對火焰面?zhèn)鞑ビ泻艽笥绊憽?/p>
由于擋板的存在,擋板左側(cè)形成燃燒腔,腔內(nèi)超壓值較大,擋板處形成超壓出口,造成氣體噴涌,如圖5所示。燃燒氣體在超壓作用下通過擋板,火焰?zhèn)鞑ニ俣仍黾樱纬赏黄馉罨鹧婷妗醢鍖饬饔忻黠@的擾動,產(chǎn)生大尺寸湍流,燃燒劇烈,進(jìn)一步增加火焰?zhèn)鞑ニ俣取T趦烧吖餐饔孟禄鹧媲瓣嚸娓咚偻ㄟ^擋板,同時,擋板與管壁形成的角落處氣體燃燒速度較慢,形成火焰斷層,如圖6。隨著角落可燃?xì)馊紵M,火焰斷層消失。
3.3 擋板位置對火焰?zhèn)鞑サ挠绊?/p>
分別監(jiān)測3、4、5號模型火焰面中心位置。記錄火焰面速度與時間關(guān)系,如圖7所示。
在點(diǎn)火初期,5號模型中擋板距離點(diǎn)火區(qū)較遠(yuǎn),對火焰面?zhèn)鞑ビ绊懖淮螅?、4號模型中,管道內(nèi)火焰面?zhèn)鞑ニ俣仍鲩L較快,4號模型速度峰值大于其他管道。火焰面?zhèn)鞑ブ衅冢?號模型受擋板影響,火焰面?zhèn)鞑ニ俣确逯蹈哂谄渌艿馈_@意味著,距離點(diǎn)火區(qū)遠(yuǎn)的擋板對火焰面?zhèn)鞑コ跗谟绊戄^小,對達(dá)到擋板前后的火焰面?zhèn)鞑ゲ庞休^明顯的影響。當(dāng)擋板處于合適位置上時,將產(chǎn)生最大火焰面?zhèn)鞑ニ俣龋撐恢镁嚯x點(diǎn)火區(qū)0.5~3倍管徑。
3.4 擋板數(shù)量對火焰?zhèn)鞑サ挠绊?/p>
監(jiān)測6號模型火焰面中心位置變化。4、6號模型火焰面速度隨時間變化曲線如圖8所示。6號模型火焰面速度峰值大于其他管道。相對于4號模型,6號模型中火焰面速度峰值增加近30%。由此可知,兩塊擋板之間形成的燃燒腔對火焰面?zhèn)鞑ト匀挥屑铀僮饔谩S^測顯示,50mm擋板處火焰面?zhèn)鞑ネ瑯映霈F(xiàn)斷層現(xiàn)象。
本文對7個充滿甲烷-空氣預(yù)混氣體,內(nèi)部擋板設(shè)置情況不同的封閉管道模型進(jìn)行燃燒模擬,得出以下結(jié)論:
(1)在L/D較大的密閉容器預(yù)混氣體燃燒過程中,最大火焰面速度在燃燒初期形成。傳播過程中火焰面速度受邊界條件影響出現(xiàn)波動。
(2)當(dāng)管道內(nèi)存在障礙物時,火焰面?zhèn)鞑ニ俣却蟠笤黾樱枞试礁撸俣仍黾釉酱螅⒃趽醢逄幊霈F(xiàn)火焰斷層現(xiàn)象。
(3)存在最不利距離,當(dāng)障礙物和點(diǎn)火區(qū)的距離與管徑的比值在0.5~3范圍內(nèi)時,出現(xiàn)火焰面?zhèn)鞑ニ俣茸畲笾怠U系K物離點(diǎn)火區(qū)太遠(yuǎn)或者太近時,其對火焰面?zhèn)鞑ニ俣鹊挠绊懗潭冉档汀?/p>
(4)增加障礙物數(shù)量會形成多個燃燒腔,燃燒腔對火焰面?zhèn)鞑ビ屑铀僮饔谩?/p>
[1]趙衡陽.氣體和粉塵爆炸原理[M].北京:北京理工大學(xué)出版社,1996.
Zhao Hengyang. Principle of gas and dust explosion [M]. Beijing:Beijing Institute of Technology Press, 1996(in Chinese)
[2]王志榮,蔣軍成,鄭楊艷.連通容器氣體爆炸的CFD模擬[J].化工學(xué)報,2007,58(4):854-861.
Wang Zhirong,Jiang Juncheng,Zheng Yangyan. CFD simulation on gas explosion field in linked vessels[J].Journal of Chemical Industry and Engineering,2007,58(4):854-861(in Chinese).
[3]李帆,譚迎新.管道內(nèi)甲烷-空氣預(yù)混爆炸燃燒的數(shù)值模擬[J].工業(yè)安全與環(huán)保,2009,35(4):28-30.
Li Fan,Tan Yingxin. Numerical Simulation of Detonation and Combustion of Methane-air Premix in Pipeline[J]. Industrial Safety and Environmental Protection,2009,35(4):28-30(in Chinese).[4]畢明樹,尹旺華,丁信偉.管道內(nèi)可燃?xì)怏w爆炸的一維數(shù)值模擬[J].天然氣工業(yè),2003,23(4):89-92.
Bi Mingshu,Yin Wanghua,Ding Xinwei,et al.One-dimensional numerical simulation of gas explosion in pipe[J]. natural gas industry,2003,23(4):89-92(in Chinese).
[5]Phylakton H,Andrews G E. Gas explosions in long closed vessels [J].Combustion Science and Technology,1991,77(1):27-39
[6]Bjerketvedt D,Bakke J R,Wingerder K V. Gas explosion handbook [J].Journal of Hazardous Materials,1997,52:1-150