外部空間與初始溫度對(duì)氫氣與空氣混合氣體爆炸過(guò)程的影響?

趙懿明 楊振欣 張 欣 張傳彪 許張歸 周莊紅 張 云 趙亮亮 曹衛(wèi)國(guó)

中北大學(xué)環(huán)境與安全工程學(xué)院(山西太原，030051)

引言

氫氣因能量密度大、燃燒產(chǎn)物無(wú)污染等優(yōu)點(diǎn)而得到應(yīng)用和推廣[1]。但從化學(xué)本質(zhì)上講，氫氣屬于第2類危險(xiǎn)化學(xué)品，在儲(chǔ)存、運(yùn)輸和使用中易發(fā)生爆炸事故[2-3]。因此，安全使用氫能也是在氫氣推廣使用中亟需解決的問(wèn)題。

為了更好地認(rèn)識(shí)氫氣的危險(xiǎn)性，大量研究人員對(duì)氫氣爆炸的壓力、火焰等燃爆特性及影響因素開展了相關(guān)研究。Kim等[4]對(duì)圓形管道內(nèi)氫氣的自燃機(jī)理進(jìn)行流動(dòng)可視化研究，發(fā)現(xiàn)自燃首先發(fā)生在氫氣與空氣混合區(qū)前鋒面后管壁處的高度混合點(diǎn)處。Shen等[5-7]探究了方形管道中氫氣與空氣預(yù)混火焰?zhèn)鞑バ螒B(tài)的變化，發(fā)現(xiàn)郁金香變形總是伴隨著火焰尖端速度的脈動(dòng)。Wang等[8]在具有單個(gè)孔板的圓柱形容器中進(jìn)行氫氣與空氣混合氣體的爆炸特性實(shí)驗(yàn)，發(fā)現(xiàn)孔板阻塞對(duì)最大爆炸壓力上升速率的影響顯著。進(jìn)一步增大孔板阻塞率發(fā)現(xiàn)，隨著阻塞率的增加，火焰?zhèn)鞑ニ俣群捅▔毫τ猩呲厔?shì)[9-10]。

雖然從事相關(guān)研究的專業(yè)人員對(duì)氫氣爆炸的風(fēng)險(xiǎn)具有較為豐富的認(rèn)識(shí)，但是由于對(duì)氫氣爆炸發(fā)生機(jī)制還未完全了解，導(dǎo)致開展氫氣爆炸實(shí)驗(yàn)存在一定的安全風(fēng)險(xiǎn)。而隨著計(jì)算機(jī)技術(shù)的發(fā)展，數(shù)值模擬可以對(duì)現(xiàn)階段實(shí)驗(yàn)無(wú)法捕捉到的爆炸參數(shù)進(jìn)行仿真再現(xiàn)，輔助研究人員對(duì)爆炸過(guò)程中微觀動(dòng)力學(xué)過(guò)程及爆炸機(jī)理進(jìn)行認(rèn)識(shí)，因此被廣泛地應(yīng)用在氣體爆炸的研究中。周寧等[11]采用大渦模擬(LES)模型研究了大長(zhǎng)徑比空間對(duì)氫氣與空氣爆炸過(guò)程的影響，發(fā)現(xiàn)火焰陣面結(jié)構(gòu)發(fā)展受管道截面尺寸影響較大。任家帆等[12]基于Navier-Stokes方程組以及k-ε湍流模型，在球形爆炸容器內(nèi)對(duì)不同比例的氫氣與空氣混合氣體的燃爆壓力及壓力上升速率進(jìn)行數(shù)值模擬研究，并通過(guò)實(shí)驗(yàn)驗(yàn)證了該模型的準(zhǔn)確性。Ugarte等[13]對(duì)氫氣與空氣混合氣體的泄爆場(chǎng)景建立數(shù)值計(jì)算模型，并改變?nèi)萜鞔笮?、形狀和點(diǎn)火位置進(jìn)行實(shí)驗(yàn)，發(fā)現(xiàn)容器內(nèi)壓力傳播規(guī)律受容器大小和形狀影響，而泄爆口火焰形狀、大小受點(diǎn)火位置影響較大。通過(guò)對(duì)以上研究結(jié)果進(jìn)行分析可以看出:氫氣爆炸的燃爆特性受外部環(huán)境影響，尤其是受限制空間影響；同時(shí)，初始溫度作為重要的環(huán)境因素，對(duì)氫氣的燃爆特性也有重要影響。

鑒于此，結(jié)合實(shí)驗(yàn)和數(shù)值模擬，使用20 L球形爆炸實(shí)驗(yàn)裝置和20 L圓柱形爆炸實(shí)驗(yàn)裝置對(duì)不同氫氣濃度時(shí)氫氣與空氣混合氣體的燃爆壓力特性進(jìn)行了研究。之后，對(duì)理論當(dāng)量濃度氫氣在兩種容器內(nèi)的燃爆過(guò)程進(jìn)行數(shù)值模擬仿真，并通過(guò)實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證仿真結(jié)果的正確性。在此基礎(chǔ)上，對(duì)比分析了容器形狀和初始溫度對(duì)氫氣與空氣混合氣體爆炸壓力及火焰?zhèn)鞑ヌ匦缘挠绊?，為氫氣的安全使用提供基礎(chǔ)數(shù)據(jù)。

1 實(shí)驗(yàn)儀器

球形爆炸容器和圓柱形爆炸容器如圖1所示。球形容器的主體是一個(gè)內(nèi)徑為0.34 m的球體；圓柱形容器是一個(gè)高為0.30 m、內(nèi)徑為0.29 m的圓柱體。兩者都為不銹鋼材質(zhì)，體積均為20 L。容器內(nèi)設(shè)有抽真空、排氣、進(jìn)氣、清洗等接口，抽真空接口附近安裝真空表。容器都包括配氣系統(tǒng)、點(diǎn)火系統(tǒng)、加熱系統(tǒng)、同步控制系統(tǒng)等。在實(shí)驗(yàn)開始之前，容器內(nèi)部抽真空，氫氣和空氣通過(guò)氣體質(zhì)量流量計(jì)控制進(jìn)行預(yù)混合，然后將預(yù)混氣瓶中的氣體注入到爆炸容器中，在爆炸容器中靜止10 min，以消除氣體湍流的影響。同時(shí)，為了保證點(diǎn)火成功，采用100 mJ的電火花點(diǎn)火。球形容器和圓柱形容器的壓力傳感器均設(shè)置在容器側(cè)壁中心處。球形爆炸容器在進(jìn)行壓力測(cè)試的同時(shí)，采用高速攝影拍攝球體內(nèi)部的火焰?zhèn)鞑ミ^(guò)程。

圖1 不同形狀的氣體爆炸實(shí)驗(yàn)裝置示意圖Fig.1 Schematic diagram of experimental apparatus of gas explosion with different shapes

2 數(shù)值模擬

2.1 計(jì)算模型

建立了符合實(shí)驗(yàn)情況的三維幾何模型，兩種不同形狀的容器中都包括點(diǎn)火電極。為反映兩種不同規(guī)格容器內(nèi)的動(dòng)態(tài)反應(yīng)壓力，對(duì)爆炸容器中不同位置處的壓力進(jìn)行了采集，計(jì)算模型及具體壓力采集位置見圖2。圖2中，A點(diǎn)代表點(diǎn)火中心位置；B點(diǎn)代表容器的右側(cè)壁面中心位置；C點(diǎn)代表上壁面中心位置。

圖2三維幾何模型及壓力采集位置Fig.2 3D geometric model and pressure acquisition location

使用ICEM軟件對(duì)幾何模型進(jìn)行網(wǎng)格劃分。整個(gè)流場(chǎng)采用四面體網(wǎng)格，為了保證網(wǎng)格的獨(dú)立性和計(jì)算結(jié)果的收斂性，在計(jì)算前對(duì)網(wǎng)格的無(wú)關(guān)性進(jìn)行了驗(yàn)證。計(jì)算最終采用的網(wǎng)格數(shù)量:球形密閉容器中為366 888個(gè)，圓柱形密閉容器中為448 282個(gè)。網(wǎng)格分布如圖3所示。

圖3三維幾何模型的網(wǎng)格分布Fig.3 Grid distribution of 3D geometric model

在模擬計(jì)算中，每一個(gè)時(shí)間步都進(jìn)行20步迭代，且能量殘差、動(dòng)量殘差和質(zhì)量殘差都小于10-6，從而保證了計(jì)算結(jié)果的收斂性。在容器中心區(qū)域點(diǎn)火，以直徑15 mm的高溫區(qū)域模擬點(diǎn)火源，點(diǎn)火源初始溫度1 200 K，其他參數(shù)與未反應(yīng)區(qū)域相同，近似模擬電火花引發(fā)的氣體強(qiáng)迫點(diǎn)火。

2.2 數(shù)值計(jì)算方法

利用計(jì)算流體力學(xué)軟件ANSYS Fluent求解基本控制方程。主要包括質(zhì)量方程、動(dòng)量方程和能量方程[14]。因本文中氣體爆炸模擬是一個(gè)化學(xué)反應(yīng)過(guò)程，因此還需考慮組分輸運(yùn)模型[15]。

模擬時(shí)采用壓力基瞬態(tài)求解器；壓力-速度耦合方程采用半隱式Simple算法求解；梯度差分采用Green-Gauss節(jié)點(diǎn)基格式；壓力、密度、動(dòng)量、湍流變量以及燃燒組分等均采用二階迎風(fēng)格式。組分輸運(yùn)模型在湍流流動(dòng)中包括了詳細(xì)的化學(xué)反應(yīng)機(jī)理，選用氫氣燃燒的19步詳細(xì)化學(xué)反應(yīng)機(jī)理進(jìn)行數(shù)值計(jì)算[16]，假定化學(xué)反應(yīng)都發(fā)生在小渦當(dāng)中，反應(yīng)時(shí)間由小渦的生存時(shí)間和化學(xué)反應(yīng)本身所需要的時(shí)間共同控制。

3 結(jié)果與討論

3.1 不同氫氣濃度爆炸時(shí)容器內(nèi)典型的壓力特性

不同容器中，氣體爆炸的火焰?zhèn)鞑ニ俣燃捌漤憫?yīng)壓力等參數(shù)并不完全一致。對(duì)氫氣體積分?jǐn)?shù)為

10.0%、20.0%、28.6%、30.0%、40.0%、50.0%、60.0%和70.0%的混合氣體在20 L球形容器內(nèi)的爆炸火焰形貌傳播特征及其壓力變化規(guī)律進(jìn)行了研究。其中，體積分?jǐn)?shù)28.6%為氫氣的理論當(dāng)量濃度。利用高速攝影裝置和壓力傳感器分別對(duì)密閉容器內(nèi)火焰發(fā)展?fàn)顩r及爆炸壓力變化情況進(jìn)行了測(cè)試，分析火焰?zhèn)鞑グl(fā)展過(guò)程中的爆炸壓力響應(yīng)規(guī)律。氫氣體積分?jǐn)?shù)為10.0%、30.0%、50.0%時(shí)，典型的爆炸壓力和火焰形貌組合圖如圖4所示。

通過(guò)圖4可知，在點(diǎn)火初期，火焰從點(diǎn)燃位置向球體四周自由擴(kuò)散。由于點(diǎn)火初期容器內(nèi)未燃?xì)怏w溫度較低，因此火焰面緩慢發(fā)展，火焰前鋒陣面近似球形。隨著爆炸的發(fā)展，未燃?xì)怏w逐漸被加熱。當(dāng)氫氣體積分?jǐn)?shù)為30.0%時(shí)，火焰面發(fā)展最快。點(diǎn)火后3 ms，火焰面迅速靠近觀察視窗并于5 ms后充滿整個(gè)球體；在7 ms時(shí)，通過(guò)視窗拍攝到的火焰尺度及亮度均達(dá)到最大，在此范圍內(nèi)，容器壁面壓力迅速增加到最大；氫氣爆炸反應(yīng)結(jié)束后，因球體壁面向外部散熱，容器內(nèi)溫度和壓力不斷減小，火焰亮度減弱并逐漸熄滅。

圖4 球形容器中不同時(shí)間的爆炸壓力和火焰形貌Fig.4 Explosion pressure and flame morphology in spherical container at different times

不同氫氣濃度時(shí)容器內(nèi)的最大爆炸壓力p見表1。從表1可以看出，在測(cè)試濃度范圍內(nèi)，兩種容器測(cè)試得出的最大爆炸壓力的變化規(guī)律一致。這是因?yàn)橛绊懽畲蟊▔毫Φ闹饕蛩厥菤錃獾臐舛?。氫氣濃度過(guò)低時(shí)，氧氣供應(yīng)充分。從傳熱方面分析，氫氣濃度越高，分子數(shù)目越多，參與化學(xué)燃燒的分子數(shù)增多，釋放的總熱量越大，整個(gè)燃燒程度加強(qiáng)，爆炸壓力隨之增加并達(dá)到最大。此外，實(shí)驗(yàn)測(cè)得的氫氣的最佳爆炸濃度(體積分?jǐn)?shù)30.0%)略高于當(dāng)量比濃度(體積分?jǐn)?shù)28.6%)，造成這一現(xiàn)象的主要原因是實(shí)驗(yàn)過(guò)程中氣體混合不均勻以及測(cè)試的不確定性等。進(jìn)一步增加氫氣濃度，單位體積內(nèi)的氫氣分子數(shù)增多，有限密閉空間內(nèi)氧氣供應(yīng)不足，氫氣不能完全燃燒，參與反應(yīng)的有效氫分子數(shù)減少，最終導(dǎo)致爆炸壓力隨氫氣濃度的增加而減小。

表1 不同氫氣濃度時(shí)容器內(nèi)的最大爆炸壓力Tab.1 Maximum explosion pressure in containers at different hydrogen concentrations

造成兩種不同形狀容器爆炸壓力差異的原因?qū)⑼ㄟ^(guò)數(shù)值模擬進(jìn)行探討。

3.2 計(jì)算模型的驗(yàn)證

圖5為293 K模擬條件下兩種容器中理論當(dāng)量濃度(體積分?jǐn)?shù)28.6%)時(shí)右壁面爆炸壓力隨時(shí)間的變化規(guī)律，并與相同條件下的實(shí)驗(yàn)值進(jìn)行了對(duì)比。從圖5中可以看出:20 L球形容器中，爆炸壓力的實(shí)驗(yàn)結(jié)果在8 ms左右達(dá)到最大值0.82 MPa，模擬結(jié)果在9 ms左右達(dá)到最大值0.86 MPa；20 L圓柱形容器中，爆炸壓力的實(shí)驗(yàn)結(jié)果在9 ms附近達(dá)到最大值0.78 MPa，模擬結(jié)果也在9 ms左右達(dá)到最大值0.82 MPa。模擬結(jié)果與實(shí)驗(yàn)結(jié)果具有較高的重復(fù)性。

圖5 兩種容器內(nèi)右壁面爆炸壓力的模擬與實(shí)驗(yàn)對(duì)比Fig.5 Comparison of explosion pressure on the right wall of two kinds of containers between simulation and experiment

3.3 容器形狀對(duì)氫氣爆炸過(guò)程的影響

為了直觀重現(xiàn)不同形狀的密閉容器中氫氣多維瞬態(tài)的爆炸過(guò)程，深入分析動(dòng)態(tài)爆炸過(guò)程中容器內(nèi)爆炸壓力、火焰溫度等的變化規(guī)律，通過(guò)ANSYS Fluent軟件對(duì)理論當(dāng)量濃度下氫氣在20 L球形容器和20 L圓柱形容器的爆炸過(guò)程進(jìn)行3D模擬。圖6~圖7為模擬得出的理論當(dāng)量濃度時(shí)氫氣在不同容器中的動(dòng)態(tài)爆炸過(guò)程。

圖6 不同容器中爆炸壓力的空間分布Fig.6 Space distribution of explosion pressure in different containers

從圖6爆炸壓力云圖中可以看出:在點(diǎn)火瞬間，球形容器和圓柱形容器中點(diǎn)火電極處的壓力都達(dá)到0.1 MPa；隨后，壓力波向壁面處傳播，使壁面處的壓力在10 ms時(shí)達(dá)到0.8 MPa左右。從圖6中還可以看出，兩種容器內(nèi)隨后的壓力傳播規(guī)律并不一致。由于壁面的反射作用，球形容器中的壓力波在球體中心處疊加，使得球體中心處的壓力在20 ms時(shí)達(dá)到最大；隨后，壓力再均勻地向壁面處傳播，直至球體中各處壓力達(dá)到平衡。而圓柱形容器的壓力傳播較球形容器不規(guī)則，使得四周反射的壓力波在圓柱體中心疊加的時(shí)間和強(qiáng)度并不一致；圓柱中心壓力也在20 ms達(dá)到最大，但是強(qiáng)度較球形容器中的??；隨后，壓力向上、下壁面和左、右壁面?zhèn)鞑サ膹?qiáng)度也不一致。

從圖7爆炸溫度云圖中可以看出，球形容器和圓柱形容器中溫度的傳播規(guī)律較為一致。點(diǎn)火瞬間，溫度達(dá)到2 000 K，火焰以球形向四周傳播；5 ms時(shí)，火焰?zhèn)鬟f到容器的一半；10 ms時(shí)，火焰已充滿整個(gè)容器，此時(shí)容器內(nèi)溫度達(dá)到最高；32 ms時(shí)，由于反應(yīng)后壁面的冷卻作用，壁面處的溫度開始下降。

3.4 容器內(nèi)不同位置處的壓力特性

在實(shí)驗(yàn)測(cè)試過(guò)程中，壓力傳感器一般設(shè)置在容器壁面。數(shù)值模擬能得到任一時(shí)刻、任何位置的壓力分布。以氫氣體積分?jǐn)?shù)28.6%時(shí)的爆炸壓力為例，A、B、C3個(gè)壓力采集點(diǎn)分別監(jiān)測(cè)的壓力變化情況如圖8所示。

圖8 不同位置處的爆炸壓力特性Fig.8 Explosion pressure characteristics at different locations

通過(guò)圖8可以看出:兩種容器中，點(diǎn)火電極處的壓力都在5 ms時(shí)達(dá)到最大值0.4 MPa；15 ms前，兩種形狀的密閉容器中，電極中心處的壓力傳播規(guī)律基本一致，結(jié)合圖6~圖7可知，此時(shí)的爆炸火焰已傳播到壁面，還未反射到球體中心處；20 ms時(shí)，由于壓力波均勻的反射和疊加作用，球形容器中的點(diǎn)火電極處的壓力達(dá)到極大值0.96 MPa，而由于不均勻的反射和疊加，圓柱形容器中點(diǎn)火電極處的壓力達(dá)到極大值0.81 MPa；隨后，圓柱形容器中電極處壓力下降的幅度不如球形容器中明顯。

兩種容器中，上、右壁面處的壓力都是在9 ms時(shí)就達(dá)到最大。球形容器上、右壁面處的壓力分別為0.79 MPa和0.86 MPa。造成球形容器內(nèi)部壓力不一致的原因主要是由于點(diǎn)火電極存在一定的阻礙和冷卻效應(yīng)，造成爆炸傳播過(guò)程存在差異。圓柱形容器中，上、右壁面處的壓力為0.75 MPa和0.82 MPa，主要是由于圓柱形容器側(cè)壁面距離中心的距離小于上壁面，壓力在側(cè)壁面首先形成集聚效應(yīng)。球形容器中，上、右壁面壓力的傳播趨勢(shì)基本一致，振蕩幅度也基本一致；而圓柱形容器中，上、右壁面壓力在第1個(gè)峰值前的傳播規(guī)律基本一致，到達(dá)壁面后，由于向容器內(nèi)四周不均勻地傳播，使得第1個(gè)峰值后的壓力振蕩傳播規(guī)律并不完全同步。隨著時(shí)間的推移，壓力振蕩周期逐漸一致，容器內(nèi)最大壓力穩(wěn)定在0.5 MPa左右，這也是造成兩種容器測(cè)試得出的最大爆炸壓力不一致的原因。

從圖8中還可以看出，部分點(diǎn)(例如A點(diǎn))最大爆炸壓力大于最大壁面壓力，說(shuō)明此時(shí)間段為整個(gè)體系的爆炸猛烈程度最大的時(shí)刻。整個(gè)反應(yīng)體系爆炸壓力最猛烈的位置并不在壁面處，從側(cè)面說(shuō)明了僅通過(guò)壁面壓力傳感器來(lái)測(cè)試最大爆炸壓力存在一定的局限性，壁面壓力傳感器并不能完全反映出整個(gè)爆炸體系的猛烈程度。

3.5 不同初始溫度下氫氣爆炸壓力特性

為進(jìn)一步研究氫氣爆炸壓力特性的影響因素，以20 L球形容器為例，對(duì)不同初始溫度(293、313、333、353 K)條件下、氫氣體積分?jǐn)?shù)28.6%時(shí)的爆炸壓力進(jìn)行了模擬，得到了最大爆炸壓力與時(shí)間的變化關(guān)系，如圖9所示。

圖9 不同初始溫度條件下的爆炸壓力特性Fig.9 Explosion pressure characteristics at different initial temperatures

從圖9可以看出:初始溫度從293 K上升到353 K時(shí)，達(dá)到最大爆炸壓力的時(shí)間由9 ms提前至7 ms；在體系初始?jí)毫Σ蛔兊那疤嵯?，密閉容器內(nèi)最大爆炸壓力由0.86 MPa降低到0.73 MPa。當(dāng)初始溫度提高20%，最大爆炸壓力降低15%左右。

溫度的升高使分子運(yùn)動(dòng)速度加快，單位面積內(nèi)的有效碰撞次數(shù)增加，因而反應(yīng)速率加快。而到達(dá)最大爆炸壓力的時(shí)間與反應(yīng)速率直接相關(guān)。因此，溫度升高，到達(dá)最大爆炸壓力的時(shí)間就越快。根據(jù)標(biāo)準(zhǔn)氣體狀態(tài)方程可知，當(dāng)初始溫度升高時(shí)，在體系壓力不變的前提下，單位體積的物質(zhì)的量減小，參與反應(yīng)的分子數(shù)降低。因此，密閉容器內(nèi)的最大爆炸壓力也隨之減小。

4 結(jié)論

利用20 L球形容器和20 L圓柱形容器對(duì)氫氣與空氣混合氣體爆炸過(guò)程進(jìn)行了實(shí)驗(yàn)和模擬研究，得到主要結(jié)論如下:

1)兩種容器得出的不同氫氣濃度下的爆炸壓力隨氫氣濃度的變化規(guī)律基本一致，最大爆炸壓力出現(xiàn)在氫氣體積分?jǐn)?shù)為30.0%的條件下；模擬結(jié)果與實(shí)驗(yàn)結(jié)果基本一致。

2)兩種容器中，氣體爆炸后的火焰均呈球面往外發(fā)展。容器中，上壁面的壓力均低于右壁面的壓力。由于壁面的反射作用，導(dǎo)致兩種容器中壓力振蕩趨勢(shì)并不一致，最終使得球形容器的最大爆炸壓力高于圓柱形容器。

3)溫度的升高使分子運(yùn)動(dòng)速度加快，單位面積內(nèi)的有效碰撞次數(shù)增加，因而反應(yīng)速率加快。初始溫度從293 K上升到353 K時(shí)，球形容器內(nèi)達(dá)到最大爆炸壓力的時(shí)間由9 ms提前至7 ms。當(dāng)初始溫度升高時(shí)，在體系壓力不變的前提下，單位體積的物質(zhì)的量減小，參與反應(yīng)的分子總數(shù)降低，最大爆炸壓力由0.86 MPa降低到0.73 MPa。