許曉元，孫金華，劉晅亞

（1. 中國(guó)科學(xué)技術(shù)大學(xué)，安徽 合肥 230026；2. 應(yīng)急管理部天津消防研究所，天津 300381）

含有可燃?xì)怏w的連通型裝置在現(xiàn)實(shí)中普遍存在?？扇?xì)怏w通常在密封的封閉空間中，如反應(yīng)器、反應(yīng)釜、儲(chǔ)罐等各類化工容器或生產(chǎn)裝置，這些裝置之間并不孤立，而是經(jīng)由管道連接形成連通型裝置?？扇?xì)怏w在生產(chǎn)、儲(chǔ)存、運(yùn)輸和使用過(guò)程中由于人為失誤、容器和管道的缺陷、外部環(huán)境等因素，導(dǎo)致容器管道封閉空間內(nèi)混入空氣等其他助燃性氣體，極易發(fā)生爆炸。經(jīng)過(guò)前人研究，連通裝置爆炸強(qiáng)度比單一容器更高，這是因?yàn)楸úê腿紵鹧嫱ㄟ^(guò)管道傳播，會(huì)導(dǎo)致氣體的壓縮并使得氣體湍流程度更高。研究表明，當(dāng)連通型裝置間的管道較長(zhǎng)時(shí)，長(zhǎng)、徑比超過(guò)40，爆燃就有可能在管道內(nèi)演變?yōu)楸Z，此時(shí)爆炸傳播速度可達(dá)2 000 m/s，壓力達(dá)到甚至超過(guò)3 MPa。除連通型工藝裝置外，實(shí)際生產(chǎn)活動(dòng)中也存在多種形式的連通型結(jié)構(gòu)。如建筑物中暢通的兩個(gè)或多個(gè)房間，礦井中的空間及通道，核電保護(hù)裝置與冷卻循環(huán)裝置等。這些連通型裝置內(nèi)部都有可能發(fā)生危害性極大的可燃?xì)怏w爆炸事故。在實(shí)際生產(chǎn)裝置爆炸事故中，由于管道較長(zhǎng)、裝置結(jié)構(gòu)復(fù)雜，因此，發(fā)生可燃?xì)怏w與主燃?xì)怏w混合過(guò)程中，可燃?xì)怏w混合物在裝置內(nèi)分布不均勻，存在體積分?jǐn)?shù)梯度。另外，煤礦井下通風(fēng)不暢、密閉的隧道中瓦斯也存在體積分?jǐn)?shù)梯度。因此，研究具有體積分?jǐn)?shù)梯度的連通裝置的可燃?xì)怏w爆炸特性對(duì)指導(dǎo)爆炸防控工作更具有實(shí)際意義。

針對(duì)連通裝置內(nèi)可燃?xì)怏w爆炸特性和泄爆技術(shù)的研究，國(guó)內(nèi)外學(xué)者做了大量的研究，前人針對(duì)連通裝置爆炸的絕大多數(shù)研究都基于連通裝置內(nèi)的氣體為均一體積分?jǐn)?shù)，前人的研究主要集中在爆炸危害后果、爆炸影響因素和泄爆技術(shù)幾個(gè)方面。連通裝置可燃?xì)怏w爆炸的影響因素主要為尺寸效應(yīng)（容器容積比例、管道的長(zhǎng)度、管徑的大小等）、點(diǎn)火位置、火焰?zhèn)鞑シ较?、障礙物阻塞率和位置等。在連通裝置爆炸危害后果研究方面，Bartkneckt[1]指出，當(dāng)連通裝置兩個(gè)容器的容積比為1∶1 時(shí)，與獨(dú)立容器相比壓力上升速率提高了10 倍，壓力上升速率增高的原因是氣體湍流和噴射火焰；Phylakton 等[2]研究發(fā)現(xiàn)，連通裝置的爆炸強(qiáng)度及壓力上升速率與燃燒速度關(guān)系密切，而燃燒速度又與湍流程度相關(guān)，通過(guò)測(cè)定，連通裝置內(nèi)最大燃燒速度可達(dá)370 m/s，最大壓力上升速率可達(dá)2.086×108Pa/s。在連通裝置爆炸影響因素研究方面，Lunn 等[3]通過(guò)粉塵爆炸實(shí)驗(yàn)研究了連通裝置容積比及連接管徑的影響，研究表明，管道體積與容積體積之間的比例也是影響燃燒爆炸的因素之一；Holbrow 等[4]針對(duì)連通裝置的粉塵爆炸做了大量試驗(yàn)研究，可燃物質(zhì)采用煤粉、石墨粉等，實(shí)驗(yàn)發(fā)現(xiàn)，管道長(zhǎng)度、直徑、容器容積比、粉塵的燃爆性能和泄放面積均對(duì)爆炸壓力有很大影響；Holbrow 等[5]開(kāi)展了實(shí)驗(yàn)研究，實(shí)驗(yàn)裝置為容積為2～20 m3的圓柱形容器，管道長(zhǎng)15 m、半徑分別為0.15、0.25、0.50 m 的管道容器組合的連通裝置，研究了容積比例、管徑、泄爆面積等因素對(duì)連通裝置內(nèi)粉塵爆炸和泄爆的影響；嚴(yán)建駿等[6]、尤明偉等[7]、王志榮等[8]采用實(shí)驗(yàn)方法研究了連通裝置內(nèi)氣體爆炸過(guò)程，主要從初始條件、點(diǎn)火位置、火焰?zhèn)鞑シ较騺?lái)分析連通容器內(nèi)壓力變化；尤明偉等[9]分析了障礙物阻塞率和位置、連接管道直徑和長(zhǎng)度、火焰?zhèn)鞑シ较?、點(diǎn)火位置等因素對(duì)火焰?zhèn)鞑ゼ氨◤?qiáng)度的影響。在連通裝置泄爆方面，王志榮等[10]研究了連通裝置內(nèi)氣體的泄爆，結(jié)果表明一個(gè)泄爆口不能明顯降低容器內(nèi)的超壓，而兩個(gè)泄壓口同時(shí)泄爆才能更好地實(shí)現(xiàn)泄爆。

學(xué)者們?cè)谘芯恐袑Ⅲw積分?jǐn)?shù)梯度方向與爆炸波傳播方向垂直的情形定義為垂直體積分?jǐn)?shù)梯度，將體積分?jǐn)?shù)梯度方向與爆炸波傳播方向平行的情形定義為平行體積分?jǐn)?shù)梯度，見(jiàn)圖1。前人的體積分?jǐn)?shù)梯度氣體爆炸研究主要用的方法為實(shí)驗(yàn)研究和數(shù)值模擬研究，研究的對(duì)象大多為單一容器或管道。在受限空間內(nèi)平行體積分?jǐn)?shù)梯度爆炸研究方面，Thomas 等[11]在22 mm×10 mm×500 mm 的受限空間及直徑為50 mm 的爆轟管內(nèi)研究了平行體積分?jǐn)?shù)梯度對(duì)爆炸傳播速度的影響；Kuznetsov 等[12]在內(nèi)徑為174 mm 長(zhǎng)度為6.0～11.2 m 的爆轟管研究了氫氣-空氣平行體積分?jǐn)?shù)梯度對(duì)爆轟波的影響，探討爆轟波熄滅、二次起爆、爆燃轉(zhuǎn)爆轟等動(dòng)力學(xué)行為。在受限空間內(nèi)垂直體積分?jǐn)?shù)梯度爆炸研究方面，Vollmer 等[13-14]研究了截面60 mm×300 mm長(zhǎng)5.4 m 的管道內(nèi)氫氣-空氣垂直體積分?jǐn)?shù)梯度、障礙物等對(duì)氫氣-空氣火焰加速、爆燃轉(zhuǎn)爆轟距離的影響；Kessler 等[15]研究了垂直體積分?jǐn)?shù)梯度中爆炸波的穩(wěn)定性和淬火概率問(wèn)題；Ishii 等[16]在截面40 mm×20 mm 長(zhǎng)0.5 m 的管道內(nèi)研究了垂直體積分?jǐn)?shù)梯度如何影響爆轟波陣面形狀和爆轟不穩(wěn)定性；Wang 等[17]對(duì)前人的實(shí)驗(yàn)開(kāi)展了數(shù)值模擬，對(duì)氫氣爆炸和丙烷爆炸流場(chǎng)進(jìn)行研究；Han 等[18]利用數(shù)值模擬的方法研究了垂直體積分?jǐn)?shù)梯度對(duì)爆轟波前沿結(jié)構(gòu)穩(wěn)定性的影響。

綜上所述，封閉空間具有體積分?jǐn)?shù)梯度的可燃?xì)怏w爆炸研究，由于其與實(shí)際工程問(wèn)題最接近，且通過(guò)研究能解釋實(shí)際工程問(wèn)題所面臨的難點(diǎn)，更為重要的是非均勻氣體爆炸所涉及的火焰加速、爆燃轉(zhuǎn)爆轟、爆轟波傳播機(jī)理與均勻氣體存在較大的差別，近年來(lái)逐漸成為國(guó)外學(xué)者的研究熱點(diǎn)，而針對(duì)連通裝置內(nèi)體積分?jǐn)?shù)梯度的爆炸研究還鮮有報(bào)道。本文中將通過(guò)數(shù)值模擬的方法研究連通裝置內(nèi)均一體積分?jǐn)?shù)和具有體積分?jǐn)?shù)梯度的CH4爆炸過(guò)程中的爆炸參數(shù)、爆炸波傳播等特點(diǎn)，以期為連通裝置內(nèi)可燃?xì)怏w爆炸事故防控、泄爆技術(shù)提供理論指導(dǎo)。

圖1 體積分?jǐn)?shù)梯度分類Fig.1 Volume fraction gradient classification

1 數(shù)值模擬研究

1.1 物理模型模擬軟件及工況設(shè)置介紹

選取容器管道連通裝置作為研究對(duì)象，研究具有體積分?jǐn)?shù)梯度的容器管道連通裝置內(nèi)甲烷氣體爆炸特性。該裝置由兩個(gè)圓柱形爆炸容器通過(guò)方形管道連通，水平布置在地面上，見(jiàn)圖2。大容器的容積為60 L，小容器的容積為20 L，兩個(gè)容器直徑與容器高度之比為1∶1，在大小容器底部由一個(gè)截面積為0.035 m×0.035 m，長(zhǎng)為3 m 的方形管道連接。為了研究體積分?jǐn)?shù)梯度，將連通裝置劃分為5 個(gè)區(qū)域，區(qū)域1 為大容器，區(qū)域2 為連接大容器的長(zhǎng)為1 m 的管道1，區(qū)域3 為連接管道1 的長(zhǎng)為1 m 的管道2，區(qū)域4 為連接管道2 和小容器的長(zhǎng)為1 m 的管道3，區(qū)域5 為小容器。

本文中將容積總計(jì)83.7 L 的連通裝置劃分為156244 個(gè)非結(jié)構(gòu)化網(wǎng)格，網(wǎng)格劃分情況見(jiàn)圖3?？扇?xì)怏w初始條件設(shè)置：壓力為105Pa，溫度為300 K，爆炸容器的壁面設(shè)置為粗糙的絕熱壁面。本研究中做了網(wǎng)格無(wú)關(guān)化分析，當(dāng)網(wǎng)格劃分為18 348 個(gè)非結(jié)構(gòu)化網(wǎng)格時(shí)，某一工況條件下連通裝置爆炸的最大壓力為665.5 kPa，當(dāng)網(wǎng)格劃分為85 648 個(gè)非結(jié)構(gòu)化網(wǎng)格時(shí)，相同工況條件下連通裝置爆炸的最大壓力為680.2 kPa，當(dāng)網(wǎng)格劃分為156 244 個(gè)非結(jié)構(gòu)化網(wǎng)格時(shí)，相同工況條件下連通裝置爆炸的最大壓力為676.9 kPa，粗糙網(wǎng)格、中等網(wǎng)格與精密網(wǎng)格壓力值的誤差分別為?1.68%和0.488%，誤差較小，可知三種網(wǎng)格劃分情況對(duì)結(jié)果的影響很小，因此，本研究所有工況的網(wǎng)格設(shè)置為156 244 個(gè)非結(jié)構(gòu)化網(wǎng)格。為了驗(yàn)證Fluidyn 軟件對(duì)爆炸模擬的可靠性，本文利用本軟件模擬了文獻(xiàn)[19]中同工況條件9.5%甲烷爆炸的實(shí)驗(yàn)，物理模型的對(duì)比圖見(jiàn)圖4，根據(jù)實(shí)驗(yàn)工況將數(shù)值模擬中爆炸容器壁面設(shè)置為傳熱時(shí)，通過(guò)分析結(jié)果可知，數(shù)值模擬的結(jié)果比實(shí)驗(yàn)值高7.10%，偏差可以接受。由于當(dāng)壁面假設(shè)為絕熱時(shí)，爆炸危害后果最嚴(yán)重，因此，本文所有的研究中均將壁面設(shè)置為絕熱條件。

圖2 管道容器連通裝置Fig.2 A connection device with pipes and vessels

數(shù)值模擬軟件Fluidyn 求解Navier-Stokes 方程以及描述理想氣體混合物的物種體積分?jǐn)?shù)、質(zhì)量和能量守恒的方程。求解紊流方程的雷諾平均，雷諾應(yīng)力采用線性渦黏模型（LEVM）進(jìn)行建模?？刂品匠倘缦拢?/p>

式中： ρ 為密度，U 為速度矢量，ym為物質(zhì)m 的質(zhì)量分?jǐn)?shù)，t 為溫度，Dm為物質(zhì)m 的有效擴(kuò)散系數(shù)，Sm為物質(zhì)m 方程的源項(xiàng) m =1,···,n。

Fluidyn 中線性渦黏性模型為：

圖3 網(wǎng)格劃分Fig.3 Model mesh generation

圖4 實(shí)驗(yàn)裝置和物理模型的比較Fig.4 Comparison between the experimental apparatus and the physical model

式中：Γ 為輸運(yùn)系數(shù)，下標(biāo)l 表示依賴于流體的層流（分子）成分；下標(biāo)t 表示湍流成分。

數(shù)值模擬軟件Fluidyn 中湍流模型運(yùn)用的為k-ε 模型，k-ε 模型是一個(gè)雙方程線性渦黏性模型。Fluidyn CFD 求解器在三維空間和時(shí)間上求解的控制方程可以表示為一般對(duì)流擴(kuò)散方程，其形式如下：

數(shù)值模擬軟件Fluidyn 中爆炸模型的化學(xué)反應(yīng)流模型包含Arrhenius 化學(xué)動(dòng)力學(xué)模型和渦流耗散模型。在Arrhenius 模型中，反應(yīng)速率用 Kj表示，渦流耗散模型中，反應(yīng)速率用 K 表示。Kj完全由化學(xué)動(dòng)力學(xué)決定。Arrhenius 形式如下：

式中：Aj為指數(shù)前因子，T 為熱力學(xué)溫度，ζj為溫度指數(shù)，Ej為活化能，R 為通用氣體常數(shù)， xi為物質(zhì) i 的摩爾分?jǐn)?shù)，vi為物種 i 速率指數(shù)。

Arrhenius 模型只考慮了反應(yīng)的化學(xué)動(dòng)力學(xué)部分，但沒(méi)有考慮湍流對(duì)反應(yīng)速率的影響。由于湍流微混合，瞬時(shí)體積分?jǐn)?shù)變化很大。數(shù)值模擬軟件Fluidyn 建立了以燃料和氧化劑為反應(yīng)物的燃燒反應(yīng)的渦流耗散模型。在該模型中，反應(yīng)速率計(jì)算如下：

式中：k 為湍流動(dòng)能，ε 為湍流耗散率， ρ 為密度，Mf為燃料分子量，Aedc、Bedc為常數(shù)，yf為燃料質(zhì)量分?jǐn)?shù)，yo為氧化劑質(zhì)量分?jǐn)?shù)， yp為產(chǎn)物質(zhì)量分?jǐn)?shù)， φo為氧化劑燃料比，φp產(chǎn)物燃料比。

利用Fluidyn 軟件來(lái)對(duì)具有體積分?jǐn)?shù)梯度的容器管道連通裝置進(jìn)行研究，設(shè)置的工況見(jiàn)表1。CH4的可爆體積分?jǐn)?shù)區(qū)間為5.0%～15.0%，本文研究了整個(gè)連通裝置內(nèi)CH4體積分?jǐn)?shù)均在爆炸區(qū)間內(nèi)的情形，為保證CH4的爆炸反應(yīng)，將CH4的體積分?jǐn)?shù)取為6.0%～14.0%。區(qū)域1 的體積分?jǐn)?shù)保持6.0%不變，其他各區(qū)域體積分?jǐn)?shù)按一定數(shù)值依次增高。將體積分?jǐn)?shù)梯度定義為區(qū)域5 的體積分?jǐn)?shù)值減去區(qū)域1的體積分?jǐn)?shù)值。本文中研究了體積分?jǐn)?shù)梯度分別為2.0%、4.0%、6.0%、6.5%、7.0%、7.5%、8.0%的工況。將具有體積分?jǐn)?shù)梯度的工況時(shí)，連通裝置內(nèi)平均CH4體積分?jǐn)?shù)作為此體積分?jǐn)?shù)梯度工況的均一體積分?jǐn)?shù)對(duì)比工況，通過(guò)計(jì)算，體積分?jǐn)?shù)梯度為2.0%、4.0%、6.0%、6.5%、7.0%、7.5%、8.0%對(duì)應(yīng)的均一體積分?jǐn)?shù)分別為6.517%、7.034%、7.551%、7.697%、7.827%、7.958%、8.067%。所有工況中，點(diǎn)火位置位于區(qū)域1 即大容器的中心。為研究不同工況下的爆炸壓力等參數(shù)，在每個(gè)區(qū)域的中心位置設(shè)置測(cè)點(diǎn)，該測(cè)點(diǎn)可測(cè)量壓力、溫度、速度。

表1 模擬工況Table 1 Simulated working conditions

1.2 體積分?jǐn)?shù)梯度對(duì)連通裝置內(nèi)最大壓力的影響

對(duì)不同工況條件下最大爆炸壓力、最大爆炸壓力所處的位置、達(dá)到最大爆炸壓力時(shí)刻進(jìn)行分析，見(jiàn)表2 和表3，表2 為連通裝置不同體積分?jǐn)?shù)時(shí)最大爆炸壓力情況，表3 為連通裝置不同體積分?jǐn)?shù)梯度時(shí)最大爆炸壓力情況。通過(guò)分析表2 和表3 可知，連通裝置內(nèi)均一體積分?jǐn)?shù)為6.517%～8.067%及區(qū)域1 體積分?jǐn)?shù)6.0%體積分?jǐn)?shù)梯度為2.0%～8.0%，且區(qū)域1 中心點(diǎn)火時(shí)，最大爆炸壓力位于區(qū)域5。由于大容器作為起爆容器，小容器作為傳爆容器，爆炸激波從大容器經(jīng)管道的持續(xù)加速傳播至小容器，小容器內(nèi)的氣體受到擠壓并產(chǎn)生強(qiáng)烈化學(xué)反應(yīng)，因此導(dǎo)致小容器的爆炸壓力最大。

圖5 為不同工況條件下最大爆炸壓力隨體積分?jǐn)?shù)/體積分?jǐn)?shù)梯度變化曲線圖。通過(guò)圖5 可知，連通裝置最大爆炸壓力隨均一體積分?jǐn)?shù)的增大而增大，約成線性關(guān)系；連通裝置最大爆炸壓力隨體積分?jǐn)?shù)梯度的變化不是線性關(guān)系，最大爆炸壓力與體積分?jǐn)?shù)梯度的變化曲線出現(xiàn)2 個(gè)波峰，一個(gè)是4.0%一個(gè)是7.0%。這是因?yàn)镃H4與O2的化學(xué)反應(yīng)當(dāng)量比時(shí)C H4的空氣體積分?jǐn)?shù)為9.5%[20]，均一體積分?jǐn)?shù)越接近當(dāng)量體積分?jǐn)?shù)，爆炸壓力越大。當(dāng)體積分?jǐn)?shù)梯度為2.0%～7.0%時(shí)，所對(duì)應(yīng)的工況比相應(yīng)工況最大爆炸壓力大，即具有體積分?jǐn)?shù)梯度所發(fā)生爆炸危害后果更嚴(yán)重。當(dāng)體積分?jǐn)?shù)梯度為7.5%～8.0%時(shí)，所對(duì)應(yīng)的工況比相應(yīng)工況的最大爆炸壓力小，即均一體積分?jǐn)?shù)工況發(fā)生爆炸危害后果更嚴(yán)重。從圖5可以推斷，在體積分?jǐn)?shù)梯度7.0%～7.5%之間，存在一個(gè)體積分?jǐn)?shù)梯度值，在此體積分?jǐn)?shù)梯度情況下，所對(duì)應(yīng)的工況與相應(yīng)工況所產(chǎn)生的效果一致。由于工況11 和工況13 的體積分?jǐn)?shù)與化學(xué)當(dāng)量比較接近，化學(xué)反應(yīng)完全，釋放的能量大，因此這種情況下，均一體積分?jǐn)?shù)比具有體積分?jǐn)?shù)梯度的工況產(chǎn)生的最大爆炸壓力大；工況1、3、5、7、9 的體積分?jǐn)?shù)與化學(xué)當(dāng)量比差距較大，而工況2、4、6、8、10 中，連通裝置的部分區(qū)域的體積分?jǐn)?shù)與CH4化學(xué)當(dāng)量比體積分?jǐn)?shù)較為接近，因此這種情況下，體積分?jǐn)?shù)梯度比均一體積分?jǐn)?shù)工況所產(chǎn)生的最大爆炸壓力大。圖6 所示為不同工況條件下最大壓力時(shí)刻隨體積分?jǐn)?shù)/體積分?jǐn)?shù)梯度的變化。由圖6 可知，最大壓力時(shí)刻隨體積分?jǐn)?shù)的增大而提前，并約成線性關(guān)系；而最大壓力時(shí)刻值隨體積分?jǐn)?shù)梯度的增大先提前再滯后。這是因?yàn)?，?.517%～8.067%范圍內(nèi)，隨著體積分?jǐn)?shù)的增大，越來(lái)越接近CH4的化學(xué)當(dāng)量體積分?jǐn)?shù)，越接近化學(xué)當(dāng)量體積分?jǐn)?shù)，爆炸反應(yīng)越快，因此，達(dá)到最大爆炸壓力最大值的時(shí)刻越提前。

表2 連通裝置不同體積分?jǐn)?shù)時(shí)最大爆炸壓力相關(guān)參數(shù)Table 2 Related parameters of the maximum explosion pressures for different volume fractions of connected devices

表3 連通裝置不同體積分?jǐn)?shù)梯度時(shí)最大爆炸壓力相關(guān)參數(shù)Table 3 Related parameters of the maximum explosion pressures for different volume fraction gradients of connected devices

圖5 不同工況條件下最大爆炸壓力隨體積分?jǐn)?shù)和體積分?jǐn)?shù)梯度變化曲線Fig.5 Changes of the maximum explosion pressure with volume fraction and volume fraction gradient under different working conditions

1.3 體積分?jǐn)?shù)梯度對(duì)連通裝置內(nèi)最大壓力上升速率的影響

對(duì)不同工況條件下最大壓力上升速率、最大壓力速率所處的位置、達(dá)到最大壓力上升速率時(shí)刻進(jìn)行了分析，見(jiàn)表4 和表5，表4 為連通裝置不同體積分?jǐn)?shù)時(shí)最大壓力上升速率情況，表5 為連通裝置不同體積分?jǐn)?shù)梯度時(shí)最大壓力上升速率情況。通過(guò)分析表4 和表5 可知，連通裝置內(nèi)均一體積分?jǐn)?shù)為6.517%～8.067%及區(qū)域1 體積分?jǐn)?shù)6.0%體積分?jǐn)?shù)梯度為2.0%～8.0%，且區(qū)域1 中心點(diǎn)火時(shí)，最大壓力上升速率位于區(qū)域2 或區(qū)域3。這是因?yàn)?，壓力上升速率代表著爆炸氣流的湍流程度，爆炸壓力波從區(qū)域1 傳播至區(qū)域2，空間驟然變小，導(dǎo)致爆炸氣流湍流程度增強(qiáng)，因此，區(qū)域2 附近的管道位置壓力上升速率最大。

圖7 所示為不同工況條件下最大壓力上升速率隨體積分?jǐn)?shù)/體積分?jǐn)?shù)梯度變化曲線。通過(guò)圖7 可知，連通裝置最大壓力上升速率隨均一體積分?jǐn)?shù)的增大而增大，約成線性關(guān)系；連通裝置最大壓力上升速率隨體積分?jǐn)?shù)梯度的變化不是線性關(guān)系，最大壓力上升速率與體積分?jǐn)?shù)梯度的變化曲線出現(xiàn)2 個(gè)波峰，一個(gè)是4.0% 一個(gè)是7.5%。這是因?yàn)榫惑w積分?jǐn)?shù)越接近當(dāng)量比體積分?jǐn)?shù)，爆炸越劇烈，最大壓力上升速率也越大。當(dāng)體積分?jǐn)?shù)梯度為2.0%～4.0%時(shí)，體積分?jǐn)?shù)梯度工況比對(duì)應(yīng)的均一體積分?jǐn)?shù)工況最大壓力上升速率大，當(dāng)體積分?jǐn)?shù)梯度為6.0%～8.0%時(shí)，體積分?jǐn)?shù)梯度工況比對(duì)應(yīng)的均一體積分?jǐn)?shù)工況最大壓力上升速率小。從圖7可以推斷，在體積分?jǐn)?shù)梯度4.0%～6.0%之間，存在一個(gè)體積分?jǐn)?shù)梯度值，在此體積分?jǐn)?shù)梯度情況下，爆炸的最大湍流程度與對(duì)應(yīng)的均一體積分?jǐn)?shù)工況一致。圖8所示為不同工況條件下最大壓力上升速率時(shí)刻隨體積分?jǐn)?shù)/體積分?jǐn)?shù)梯度的變化。由圖8 可知，最大壓力上升速率時(shí)刻隨體積分?jǐn)?shù)的增大而提前，并約成線性關(guān)系；而最大壓力上升速率時(shí)刻隨體積分?jǐn)?shù)梯度的增大先提前再滯后。這是因?yàn)?，最大壓力上升速率存在于壓力快速上升即爆炸的過(guò)程中，在6.517%～8.067%范圍內(nèi)，隨著體積分?jǐn)?shù)的增大，越來(lái)越接近CH4的化學(xué)當(dāng)量體積分?jǐn)?shù)，越接近化學(xué)當(dāng)量體積分?jǐn)?shù)，爆炸反應(yīng)越快，因此，達(dá)到最大壓力上升速率最大值的時(shí)刻越提前。

表4 連通裝置不同體積分?jǐn)?shù)時(shí)最大壓力上升速率相關(guān)參數(shù)Table 4 Related parameters of the maximum pressure rise rates for different volume fractions of connected devices

表5 連通裝置不同體積分?jǐn)?shù)梯度時(shí)最大壓力上升速率相關(guān)參數(shù)Table 5 Parameters related to the maximum pressure rise rates in the connected devices with different volume fraction gradients

圖7 不同工況下最大壓力上升速率隨體積分?jǐn)?shù)和體積分?jǐn)?shù)梯度的變化Fig.7 Changes of the maximum pressure rise rate with volume fraction and volume fraction gradient under different working conditions

圖8 不同工況下最大壓力上升速率時(shí)刻隨體積分?jǐn)?shù)和體積分?jǐn)?shù)梯度的變化Fig.8 Changes of the arrival time of the maximum pressure rise rate with volume fraction and volume fraction gradient under different working conditions

1.4 體積分?jǐn)?shù)梯度對(duì)連通裝置內(nèi)最高溫度的影響

對(duì)不同工況條件下最高溫度、最高溫度所處的位置、達(dá)到最高溫度時(shí)刻進(jìn)行了分析，見(jiàn)表6 和表7，表6 為連通裝置不同體積分?jǐn)?shù)時(shí)最高溫度情況，表7 為連通裝置不同體積分?jǐn)?shù)梯度時(shí)最高溫度情況。通過(guò)分析表6 和表7 可知，連通裝置內(nèi)均一體積分?jǐn)?shù)為6.517%～8.067%，且大容器區(qū)域1 中心點(diǎn)火時(shí)，最高溫度位于區(qū)域1；連通裝置區(qū)域1 體積分?jǐn)?shù)6.0%體積分?jǐn)?shù)梯度為2.0%～8.0%，且區(qū)域1 中心點(diǎn)火時(shí)，最高溫度位于區(qū)域5。

表6 連通裝置不同均一體積分?jǐn)?shù)下最高溫度相關(guān)參數(shù)Table 6 Related parameters of the maximum temperatures in connected devices with different volume fractions

表7 連通裝置不同體積分?jǐn)?shù)梯度時(shí)最高溫度相關(guān)參數(shù)Table 7 Rrelated parameters of the maximum temperatures in connected devices with different volume fraction gradients

圖9 所示為不同工況條件下最高溫度隨體積分?jǐn)?shù)/體積分?jǐn)?shù)梯度變化曲線圖。通過(guò)圖9 可知，連通裝置最高溫度隨均一體積分?jǐn)?shù)值增大而升高，約成線性關(guān)系；連通裝置最高溫度隨體積分?jǐn)?shù)梯度的變化不是線性關(guān)系，最高溫度與體積分?jǐn)?shù)梯度的變化曲線出現(xiàn)先升高后降低的趨勢(shì)，體積分?jǐn)?shù)梯度為4.0%時(shí)，達(dá)到最高溫度。這是因?yàn)楫?dāng)均一體積分?jǐn)?shù)值越接近當(dāng)量比體積分?jǐn)?shù)時(shí)，爆炸越激烈，因此最高溫度也越高。研究中不同體積分?jǐn)?shù)梯度工況下比對(duì)應(yīng)的均一體積分?jǐn)?shù)工況的最高溫度偏高，說(shuō)明連通裝置內(nèi)相同總量的CH4，區(qū)域1 體積分?jǐn)?shù)低，區(qū)域2 至5 體積分?jǐn)?shù)依次增高，區(qū)域1 中心點(diǎn)火時(shí)，具有體積分?jǐn)?shù)梯度的工況比對(duì)應(yīng)的均一體積分?jǐn)?shù)工況釋放的能量多。圖10 所示為不同工況條件下最高溫度時(shí)刻值隨體積分?jǐn)?shù)/體積分?jǐn)?shù)梯度的變化。由圖10 可知，最高溫度時(shí)刻隨體積分?jǐn)?shù)的增大而提前，并約成線性關(guān)系；而最高溫度時(shí)刻隨體積分?jǐn)?shù)梯度的增大先提前再滯后。這是因?yàn)?，?.517%～8.067%范圍內(nèi)，隨著體積分?jǐn)?shù)的增大，越來(lái)越接近CH4的化學(xué)當(dāng)量體積分?jǐn)?shù)，越接近化學(xué)當(dāng)量體積分?jǐn)?shù)，爆炸反應(yīng)越劇烈，釋放的能量越大，因此，達(dá)到最高溫度的時(shí)刻越提前。

圖9 不同工況下最高溫度隨體積分?jǐn)?shù)和體積分?jǐn)?shù)梯度的變化Fig.9 Changes of the maximum temperature with volume fraction and volume fraction gradient under different working conditions

圖10 不同工況下最高溫度時(shí)刻值隨體積分?jǐn)?shù)和體積分?jǐn)?shù)梯度的變化Fig.10 Changes of the arrival time of the maximum temperature with volume fraction and volume fraction gradient under different working conditions

1.5 體積分?jǐn)?shù)梯度對(duì)連通裝置內(nèi)最大速度的影響

對(duì)不同工況條件下最大速度、最大速度所處的位置、達(dá)到最大速度時(shí)刻進(jìn)行了分析，見(jiàn)表8和表9，表8 為連通裝置不同體積分?jǐn)?shù)時(shí)最大速度情況，表9 為連通裝置不同體積分?jǐn)?shù)梯度時(shí)最大速度情況。通過(guò)分析表8 和表9 可知，連通裝置內(nèi)均一體積分?jǐn)?shù)為6.517%～8.067% 及區(qū)域1 體積分?jǐn)?shù)6.0%體積分?jǐn)?shù)梯度為2.0%～8.0%，且區(qū)域1 中心點(diǎn)火時(shí)，最大速度位于區(qū)域4。這是因?yàn)楫?dāng)連通裝置大容器發(fā)生爆炸時(shí)，細(xì)長(zhǎng)管道內(nèi)從大容器方向向小容器方向氣流速度持續(xù)加大，呈不斷加速過(guò)程，因此，與小容器相連的管道內(nèi)速度最大，可達(dá)400～600 m/s。

連通裝置最大速度隨均一體積分?jǐn)?shù)值增大而增大，約成線性關(guān)系；連通裝置最大速度隨體積分?jǐn)?shù)梯度的變化不是線性關(guān)系，最大速度隨體積分?jǐn)?shù)梯度增大先增大后減小。這是因?yàn)?，均一體積分?jǐn)?shù)值越接近當(dāng)量體積分?jǐn)?shù)，爆炸反應(yīng)越激烈，連通裝置內(nèi)最大速度越大。最大速度時(shí)刻值隨體積分?jǐn)?shù)的增大而降低，并約成線性關(guān)系；而最大速度時(shí)刻隨體積分?jǐn)?shù)梯度的增大先提前再滯后。這是因?yàn)椋?.517%～8.067%范圍內(nèi)，隨著體積分?jǐn)?shù)值的增大，越來(lái)越接近CH4的化學(xué)當(dāng)量體積分?jǐn)?shù)，越接近化學(xué)當(dāng)量體積分?jǐn)?shù)，爆炸反應(yīng)越激烈，因此，達(dá)到最大速度的時(shí)刻越提前。

表8 連通裝置不同體積分?jǐn)?shù)時(shí)最大速度相關(guān)參數(shù)Table 8 Related parameters of the maximum velocities in connected devices with different volume fractions

表9 連通裝置存在不同體積分?jǐn)?shù)梯度時(shí)最大速度相關(guān)參數(shù)Table 9 Related parameters of the maximum velocities in connected devices with different volume fraction gradients

1.6 結(jié)果分析

在區(qū)域1 中點(diǎn)火爆炸后，未燃?xì)怏w迅速膨脹，產(chǎn)生火焰鋒面前后的密度差導(dǎo)致流體力學(xué)不穩(wěn)定，膨脹比越大，流體動(dòng)力學(xué)不穩(wěn)定易使火焰失穩(wěn)，熱膨脹比可用如下公式計(jì)算[21]：

式中：σ 是熱膨脹比，ρu是未燃?xì)怏w密度，ρb是燃燒產(chǎn)物體積分?jǐn)?shù)，ρu和 ρb可通過(guò)GASEQ 軟件計(jì)算得到。

通過(guò)計(jì)算，區(qū)域1 中6.0%體積分?jǐn)?shù)的甲烷燃燒熱膨脹比為5.598，當(dāng)區(qū)域1 氣體完全燃燒，體積完全膨脹后為335.88 L，該數(shù)值遠(yuǎn)大于其他區(qū)域的體積總和，這意味著區(qū)域1 的未燃?xì)怏w會(huì)將區(qū)域2～4 的氣體全部壓縮至區(qū)域5 中。

從整個(gè)爆炸過(guò)程來(lái)看，各區(qū)域內(nèi)的氣體混合物溫度隨容器中壓力的升高而增加，因此可近似的看作絕熱壓縮的過(guò)程，絕熱壓縮使得未燃?xì)怏w溫度升高，對(duì)火焰燃燒速率產(chǎn)生了影響[22]，即：

式中：Su為火焰速度，Kr為室溫下的燃燒速速率，Tu為未燃?xì)鉁囟?，T0為初始溫度， p0為初始?jí)毫?，p 為爆炸壓力，β是壓力對(duì)燃速的影響因子。

文中涉及的7 種體積分?jǐn)?shù)梯度工況，剛好經(jīng)歷貧燃至富燃的燃燒過(guò)程，而由貧燃向富燃轉(zhuǎn)變過(guò)程中，最大爆炸壓力、最大爆炸壓力上升速率、火焰溫度及最大速度總體均呈現(xiàn)先升后降的趨勢(shì)。

區(qū)域1 內(nèi)點(diǎn)火后發(fā)生層流燃燒，通過(guò)區(qū)域2、區(qū)域3、區(qū)域4 傳播后，由于容器形狀改變和管道剛性壁面的作用，火焰受到擾動(dòng)，由層流燃燒轉(zhuǎn)變?yōu)橥牧魅紵?，加速向區(qū)域5 傳遞，并對(duì)區(qū)域5 進(jìn)行了預(yù)壓縮，在壓力累積與噴射火焰的共同作用下，區(qū)域5 直接發(fā)生了湍流燃燒，因此，壓力高，上升速率快。

區(qū)域2～4 均為管道內(nèi)，管內(nèi)流團(tuán)燃燒時(shí)，因體積膨脹而對(duì)周圍介質(zhì)作功，從而在火焰前方形成壓縮波，誘導(dǎo)和加速當(dāng)?shù)刭|(zhì)點(diǎn)的運(yùn)動(dòng)?；鹧嬉虼硕冃危紵娣e增大，燃燒速率加快，加速火焰的進(jìn)一步推動(dòng)又導(dǎo)致火焰陣面的更大變形。由于壁面附近存在粘性邊界層，流速的增大會(huì)使管內(nèi)出現(xiàn)湍流。大渦湍流使火焰皺褶，增加燃燒面積，小渦湍流可提高火焰內(nèi)部的有效輸運(yùn)能力，這些都有助于提高燃燒速率。湍流強(qiáng)度隨軸向流速的增加而增加。當(dāng)湍流足夠強(qiáng)時(shí)，未燃流團(tuán)將被火焰所吞食形成燃燒的相干流塊。這些流塊在軸心附近溫度高、流速大，壁面附近溫度低流速小[23]。推測(cè)由于可燃?xì)怏w在軸向體積分?jǐn)?shù)呈不均勻分布，導(dǎo)致區(qū)域2～4 內(nèi)湍流程度呈現(xiàn)較為復(fù)雜的變化，在總體先升后降的趨勢(shì)下發(fā)生波動(dòng)，根據(jù)式（7）可發(fā)現(xiàn)最大爆炸壓力、最大爆炸壓力上升速率、最大溫度及最大速度均會(huì)受到一定影響。

2 結(jié) 論

（1）連通裝置均一甲烷體積分?jǐn)?shù)6.517%至8.067%，且大容器中心點(diǎn)火工況時(shí)，最大爆炸壓力、最大爆炸壓力上升速率、最高溫度和最大速度，以及這些爆炸參數(shù)達(dá)到最大值時(shí)的時(shí)刻值隨體積分?jǐn)?shù)值的變化約呈線性關(guān)系；（2）連通裝置大容器甲烷體積分?jǐn)?shù)6.0%，體積分?jǐn)?shù)梯度為2.0%～8.0%，且大容器中心點(diǎn)火時(shí)，最大爆炸壓力、最大爆炸壓力上升速率、最高溫度和最大速度隨體積分?jǐn)?shù)梯度總體上先增大后減小，但由于管道內(nèi)湍流程度高呈現(xiàn)較為復(fù)雜的變化，在體積分?jǐn)?shù)梯度6.0%～7.5%之間可能出現(xiàn)波動(dòng)；（3）連通裝置內(nèi)均一甲烷體積分?jǐn)?shù)在6.517%～8.067%范圍內(nèi)，且大容器內(nèi)甲烷的體積分?jǐn)?shù)6.0%，體積分?jǐn)?shù)梯度為2.0%～8.0%，大容器中心點(diǎn)火時(shí)，最大爆炸壓力位于小容器內(nèi)；（4）連通裝置內(nèi)均一甲烷體積分?jǐn)?shù)在6.517%～8.067%范圍內(nèi)，且大容器內(nèi)甲烷的體積分?jǐn)?shù)為6.0%，體積分?jǐn)?shù)梯度為2.0%～8.0%，大容器中心點(diǎn)火時(shí)，最大壓力上升速率位于管道1 或管道2 中；（5）連通裝置內(nèi)均一甲烷體積分?jǐn)?shù)為6.517%～8.067%范圍內(nèi)，大容器中心點(diǎn)火時(shí)，最大溫度位置在大容器內(nèi)；連通裝置大容器內(nèi)甲烷的體積分?jǐn)?shù)為6.0%，體積分?jǐn)?shù)梯度為2.0%～8.0%時(shí)，大容器中心點(diǎn)火時(shí)，最高溫度位置在小容器內(nèi)；（6）由于火焰?zhèn)鞑ピ诠艿纼?nèi)呈現(xiàn)持續(xù)加速現(xiàn)象，因此連通裝置內(nèi)均一甲烷體積分?jǐn)?shù)為6.517%～8.067%范圍內(nèi)，大容器中甲烷體積分?jǐn)?shù)6.0%，體積分?jǐn)?shù)梯度為2.0%～8.0%，且大容器中心點(diǎn)火時(shí)，最大速度位于管道3，速度值可達(dá)400～600 m/s。