基于欠膨脹噴射火模型的含氫天然氣管道失效后果分析

，，，

(浙江大學(xué) 化工機(jī)械研究所，杭州 310027)

0 引言

含氫天然氣管道輸送是氫能利用的重要途徑之一，特別是我國(guó)獨(dú)特的地域特點(diǎn)，西部地區(qū)風(fēng)電充足，無(wú)法并網(wǎng)的棄風(fēng)、棄電可制氫后摻入現(xiàn)有西氣東輸管網(wǎng)，進(jìn)而輸送到用戶(hù)終端，具有廣闊的市場(chǎng)前景。隨著含氫天然氣輸送管道的廣泛應(yīng)用，其安全性問(wèn)題已引起全球范圍內(nèi)的廣泛關(guān)注。由于甲烷和氫氣性質(zhì)的不同[1]，氫氣的摻入往往會(huì)增大管道失效引起的后果，其中氣體的泄漏燃燒行為一直是其安全性問(wèn)題研究的重點(diǎn)[2]。當(dāng)含氫天然氣的輸送管道氣體發(fā)生泄漏時(shí)，通常會(huì)在泄漏噴口處發(fā)展成為欠膨脹噴射，若被點(diǎn)燃，由其引起的噴射火會(huì)延伸數(shù)米甚至超過(guò)十米[3]。因此，為了保證含氫天然氣管道運(yùn)行過(guò)程中的安全性，亟待研究含氫天然氣的欠膨脹噴射火行為。

關(guān)于氫氣或含氫混合氣體欠膨脹噴射火行為的研究主要有兩種手段：試驗(yàn)研究和數(shù)值模擬。國(guó)際上已開(kāi)展了一些與氫氣相關(guān)的欠膨脹噴射火試驗(yàn)研究[4-6]，基于這些試驗(yàn)結(jié)果，Birch等[7]首次提出了偽直徑計(jì)算公式，該公式根據(jù)質(zhì)量流率守恒定義了偽直徑。之后,Birch等[8]又通過(guò)膨脹過(guò)程中的質(zhì)量守恒和動(dòng)量守恒重新定義了偽直徑，同樣假設(shè)偽直徑處的壓力降至環(huán)境壓力。然而，由于在世界范圍內(nèi)研究欠膨脹噴射火的試驗(yàn)裝置數(shù)量有限，以及大規(guī)模的試驗(yàn)研究通常耗資巨大且每次試驗(yàn)結(jié)果都是在特定的條件下得到的，數(shù)據(jù)適用條件的可拓展性差，因此，在實(shí)際工程中很難進(jìn)行推廣。

數(shù)值模擬手段由于具有良好的經(jīng)濟(jì)性、可操作性及可重復(fù)性等優(yōu)點(diǎn)，在工程計(jì)算和科學(xué)研究中應(yīng)用較多[9]。目前，國(guó)內(nèi)外已針對(duì)欠膨脹噴射火開(kāi)展了一定的數(shù)值模擬研究。國(guó)際上， Ekoto等[10]通過(guò)采用Birch模型以及修正的地面反射率，建立了較為恰當(dāng)?shù)难芯看蟪叽鐨錃鈬娚浠鹧孑椛涮卣鞯臄?shù)值模型。在國(guó)內(nèi)，陳國(guó)華等[11]利用Fluent軟件對(duì)大尺寸甲烷欠膨脹噴射火進(jìn)行了模擬，研究發(fā)現(xiàn)，火焰內(nèi)部溫度分布和火焰周?chē)鸁彷椛浞植嫉臄?shù)值模擬結(jié)果與試驗(yàn)數(shù)據(jù)一致性較好，火焰尺寸的數(shù)值模擬結(jié)果與經(jīng)驗(yàn)公式也能較好地吻合，從而驗(yàn)證了Fluent軟件模擬大尺寸噴射火的適用性。但是，目前國(guó)內(nèi)外針對(duì)含氫混合氣體欠膨脹噴射火行為及相應(yīng)失效后果的研究較為缺乏。

本文針對(duì)含氫天然氣的欠膨脹噴射火行為和失效后果，利用Birch模型計(jì)算得到的偽直徑處的狀態(tài)參數(shù)來(lái)描述噴射出口處的工況，再利用計(jì)算流體力學(xué)軟件Fluent建立含氫天然氣欠膨脹噴射火的模擬模型。首先，通過(guò)對(duì)比模擬結(jié)果與試驗(yàn)結(jié)果，驗(yàn)證模型的有效性；然后，基于該模型，分析氫氣含量、管道壓力和噴口尺寸對(duì)含氫天然氣欠膨脹噴射火火焰特性的影響規(guī)律。相關(guān)研究可為含氫天然氣高壓管道輸送涉及的氫氣含量選取、管道布置、管道監(jiān)測(cè)檢查、風(fēng)險(xiǎn)管理等提供參考。

1 控制方程及數(shù)值計(jì)算模型

本文利用氣相流體的質(zhì)量、動(dòng)量、能量、化學(xué)組分等方程的解，描述氫氣/天然氣(甲烷)混合氣體欠膨脹噴射火在空氣中的擴(kuò)散燃燒?？刂品匠贪ㄟB續(xù)性方程、動(dòng)量守恒方程、能量守恒方程、組分質(zhì)量守恒方程以及氣體狀態(tài)方程，具體形式分別如下。

(1)連續(xù)性方程。

(1)

式中ρ——流體密度,kg/m3;

u,v,w——x，y，z方向上的速度分量,m/s。

(2)動(dòng)量守恒方程。

(2)

(3)

(4)

式中μ——流體的動(dòng)力黏度,Pa·s;

p——流體微元體上的壓力,MPa;

Su，Sv，Sw——?jiǎng)恿渴睾惴匠痰膹V義源項(xiàng)。

(3)能量守恒方程。

(5)

式中T——溫度,K;

k——流體的傳熱系數(shù),W/(m2·K);

cp——比熱容,J/(kg·K);

ST——流體的內(nèi)源熱及由于黏性作用流體機(jī)械能轉(zhuǎn)換為熱能的部分,J。

(4)組分質(zhì)量守恒方程。

(6)

式中cs——組分s的質(zhì)量分?jǐn)?shù)；

Ds——組分s的擴(kuò)散系數(shù),m2/s;

Ss——組分s的化學(xué)反應(yīng)速率,kg/(m3·s)。

(5)氣體狀態(tài)方程。

當(dāng)氫氣壓力低于17.2 MPa時(shí)，可作為理想氣體處理[12]，本文的研究?jī)H針對(duì)10 MPa以下的混合氣體，故該模型的計(jì)算采用理想氣體狀態(tài)方程。

p=ρRT

(7)

式中R——理想氣體常數(shù)。

對(duì)于湍流模型的選擇，大量計(jì)算結(jié)果表明k-ε模型對(duì)模擬射流的擴(kuò)散和衰減有較高的準(zhǔn)確性[13]，故本文采用k-ε湍流模型研究氫氣/甲烷混合氣體的欠膨脹噴射火。

Fluent軟件中提供了5種計(jì)算輻射換熱的模型，根據(jù)各種模型的使用條件，本文采用使用條件為光學(xué)厚度大于1的P-1輻射模型進(jìn)行計(jì)算。

關(guān)于燃燒模型，F(xiàn)luent軟件中提供了5種方法模擬燃燒化學(xué)反應(yīng)。本文采用其中最常用的有限速率模型進(jìn)行模擬，該模型是組分質(zhì)量分?jǐn)?shù)的輸運(yùn)方程，反應(yīng)速率以源項(xiàng)的形式出現(xiàn)在組分輸運(yùn)方程中。利用渦耗散模型求解流場(chǎng)噴射火燃燒反應(yīng)速率，反應(yīng)r中物質(zhì)i的產(chǎn)生速率Ri,r由以下表達(dá)式中較小的一個(gè)給出：

(8)

(9)

Mw,i——第i種物質(zhì)的分子量;

A——經(jīng)驗(yàn)常數(shù),A=4.0;

YR——反應(yīng)物組分的質(zhì)量分?jǐn)?shù);

Mw,R——反應(yīng)物組分的分子量;

B——經(jīng)驗(yàn)常數(shù),B=0.5;

YP——生成物組分的質(zhì)量分?jǐn)?shù);

Mw,j——生成物組分的分子量。

目前,關(guān)于氫氣燃燒化學(xué)反應(yīng)的模型主要采用單步、雙步和多步反應(yīng)，但針對(duì)大尺度的燃燒、爆炸，考慮計(jì)算效率與計(jì)算結(jié)果的準(zhǔn)確性[14]，本文采用最常見(jiàn)的單步總包反應(yīng),研究氫氣/甲烷混合氣體與空氣中氧氣的反應(yīng)。

采用軸對(duì)稱(chēng)二維模型研究氫氣/甲烷混合氣體自管道泄漏噴口垂直于地面射出產(chǎn)生的噴射火，計(jì)算區(qū)域選取對(duì)稱(chēng)區(qū)域的一半，形狀為矩形，如圖1所示(圖中Y方向與地面平行)。氫氣/甲烷混合氣體由矩形的左下方的速度進(jìn)口進(jìn)入計(jì)算區(qū)域；與速度進(jìn)口相鄰的管壁設(shè)置為固壁；矩形的上邊界和右邊界面為空氣環(huán)境，設(shè)置為壓力出口邊界條件；矩形的下邊界面設(shè)置為對(duì)稱(chēng)軸。

圖1 氫氣/甲烷混合氣體噴射火計(jì)算區(qū)域示意

速度入口處的湍流強(qiáng)度I按照下式計(jì)算：

I=0.16Re1/8

(10)

其中:

(11)

式中Re——雷諾數(shù);

ρ——噴口處的氣體密度,kg/m3;

u——噴口處的流速,m/s;

D——噴口直徑,m。

常壓狀態(tài)下,混合氣體的動(dòng)力黏度按照下式計(jì)算：

(12)

式中yi——組分i的摩爾分?jǐn)?shù)；

μi——組分i的黏度,Pa·s;

Mi——組分i的相對(duì)分子質(zhì)量。

2 Birch模型

本文利用Birch模型來(lái)描述所有真實(shí)噴射出口處的工況。具體模型如圖2所示。氣體從噴口噴射出到整個(gè)膨脹過(guò)程皆遵守質(zhì)量守恒和動(dòng)量守恒。

狀態(tài)0-壓力容器或管道中氣體的狀態(tài)；狀態(tài)1-氣體剛從噴口噴出的狀態(tài)；狀態(tài)2-噴射氣流經(jīng)膨脹之后，壓力降為環(huán)境壓力時(shí)的狀態(tài)；D2-需要求取的偽直徑

圖2 Birch模型

假設(shè)整個(gè)膨脹面上的黏性力可以忽略不計(jì)，且無(wú)環(huán)境流體流入，則質(zhì)量和動(dòng)量守恒方程分別為：

ρ2u2A2=ρ1u1A1CD

(13)

(14)

式中ρ1,ρ2——噴口、偽噴口處的氣體密度,kg/m3;

u1，u2——噴口、偽噴口處的氣體速度,m/s;

A1，A2——噴口、偽噴口的面積,m2;

CD——?dú)怏w泄漏系數(shù)，當(dāng)噴口形狀為圓形時(shí)，取值1.0;

p1，p2——噴口、偽噴口處的氣體壓力,MPa。

分別求取u2和A2，此處假設(shè)p2等于環(huán)境壓力p∞，得到：

(15)

(16)

接著使用等熵流關(guān)系式關(guān)聯(lián)狀態(tài)1和狀態(tài)0的氣體狀態(tài)參數(shù)：

(17)

(18)

(19)

式中p0——管道中氣體的壓力,MPa;

γ——絕熱指數(shù);

R——理想氣體常數(shù);

T0——管道中氣體的溫度,K;

M——?dú)怏w的相對(duì)分子質(zhì)量,kg/mol。

對(duì)于具有s種組分的混合氣體的γ和M，計(jì)算公式[15-16]如下：

(20)

(21)

式中aVi——混合氣體中組分i的體積分?jǐn)?shù)；

γi——組分i的絕熱指數(shù)。

對(duì)于氫氣和甲烷氣體，絕熱指數(shù)分別取1.410,1.314，相對(duì)分子質(zhì)量分別取0.002,0.016 kg/mol，由此可以得到不同氫氣含量的氫氣/甲烷混合氣體的絕熱指數(shù)和相對(duì)分子質(zhì)量。

此外，試驗(yàn)研究發(fā)現(xiàn)，偽噴口處氣體的溫度T2等于管道中的溫度T0，故有：

(22)

(23)

式中D2——偽直徑,mm;

D1——實(shí)際噴口直徑,mm。

結(jié)合計(jì)算式(22)和(23),可以得到偽噴口處的氣體速度和偽直徑大小。

3 模擬分析及討論

3.1 Birch模型的驗(yàn)證

在模擬不同氫氣含量的氫氣/甲烷混合氣體的噴射火之前，首先需要對(duì)所采用模型的正確性進(jìn)行驗(yàn)證，本文根據(jù)Schefer等[12,17]的試驗(yàn)進(jìn)行模型的驗(yàn)證，兩個(gè)案例的試驗(yàn)參數(shù)及對(duì)應(yīng)的根據(jù)Birch模型求得的模擬參數(shù)(偽噴口參數(shù))如表1所示。

表1 兩個(gè)案例的試驗(yàn)參數(shù)及對(duì)應(yīng)的偽噴口參數(shù)

火焰長(zhǎng)度是研究大尺寸噴射火的重要特征參數(shù)，故本文將用Birch模型得到的模擬火焰長(zhǎng)度與試驗(yàn)值和理論值進(jìn)行比較，對(duì)Birch模型進(jìn)行驗(yàn)證。本文采用溫度閾值的方法對(duì)火焰長(zhǎng)度的模擬結(jié)果進(jìn)行計(jì)算，Schefer等[12，18]在試驗(yàn)及模擬中以燃燒溫度1 300～1 500 K作為噴射火火焰長(zhǎng)度的判據(jù)，因此參考這一溫度范圍，設(shè)定溫度閾值為1 500 K。根據(jù)表1中的兩種工況得到的火焰云圖如圖3所示，在溫度云圖中，可以通過(guò)找到1 500 K的等溫線(xiàn)的方式獲得火焰長(zhǎng)度的模擬值。該參數(shù)的理論值采用Delichatsios[16]提出的基于無(wú)量綱弗勞德數(shù)Frf的火焰長(zhǎng)度計(jì)算公式進(jìn)行計(jì)算。

圖3 火焰溫度云圖

利用上述方法得到火焰長(zhǎng)度的模擬值、計(jì)算值及試驗(yàn)值如表2所示?？梢钥闯?，工況1中的模擬火焰長(zhǎng)度為7.51 m，試驗(yàn)得到的火焰長(zhǎng)度6.7 m，模擬火焰長(zhǎng)度比試驗(yàn)火焰長(zhǎng)度高出12.1%，與理論火焰長(zhǎng)度7.58 m相比，模擬火焰長(zhǎng)度減小了0.9%；對(duì)于工況2關(guān)于氫氣/甲烷混合氣體噴射火的模擬，得到的模擬火焰長(zhǎng)度(7.42 m)比試驗(yàn)火焰長(zhǎng)度(6.84 m)高出8.5%，比理論值高出7.4%。對(duì)于目前的大尺寸噴射火模擬研究來(lái)說(shuō)，這種誤差是可以接受的，因此，本文采用的噴射火模型能夠較好地模擬高速氫氣以及氫氣/甲烷混合氣體的欠膨脹噴射火。

表2 火焰長(zhǎng)度模擬值、試驗(yàn)值和理論值的對(duì)比

3.2 網(wǎng)格敏感性分析

在流體力學(xué)模型計(jì)算過(guò)程中，網(wǎng)格尺寸的選取直接影響到計(jì)算精度和計(jì)算規(guī)模。通常情況下，在計(jì)算之前都需要對(duì)網(wǎng)格進(jìn)行敏感性分析，保證獲得計(jì)算工況收斂的解的同時(shí)，得到計(jì)算工況下最小計(jì)算量對(duì)應(yīng)的最佳網(wǎng)格尺度。本文的網(wǎng)格敏感性分析以管道壓力10 MPa、噴口直徑5 mm、噴口溫度293.15 K、環(huán)境壓力0.101 MPa、純氫氣介質(zhì)的工況為例，計(jì)算區(qū)域?yàn)?0 m×2 m。計(jì)算采用三類(lèi)不同的網(wǎng)格尺度(a類(lèi)，b類(lèi)，c類(lèi))，分析網(wǎng)格敏感性對(duì)氫氣/甲烷混合氣體噴射火火焰軸線(xiàn)溫度和火焰長(zhǎng)度的影響。三種情況隨遠(yuǎn)離噴口方向沿軸向和徑向均采用遞增的網(wǎng)格尺寸劃分方式，其中，b類(lèi)網(wǎng)格尺度對(duì)a類(lèi)的軸向和徑向進(jìn)行加密，c類(lèi)網(wǎng)格尺度對(duì)b類(lèi)的軸向和徑向進(jìn)行加密，使得b類(lèi)的單元總數(shù)對(duì)a類(lèi)增加一倍，c類(lèi)的單元總數(shù)對(duì)b類(lèi)增加一倍，噴口處的最小網(wǎng)格尺寸接近于1 mm(c類(lèi))。

表3和圖4中給出了不同網(wǎng)格尺度下氫噴射火火焰長(zhǎng)度和火焰軸線(xiàn)溫度的計(jì)算結(jié)果。可以看出，對(duì)于軸線(xiàn)的平均溫度，三種網(wǎng)格尺度的計(jì)算結(jié)果偏差較小，但就火焰長(zhǎng)度來(lái)說(shuō)，由a類(lèi)網(wǎng)格尺度得到的結(jié)果明顯大于b類(lèi)和c類(lèi)，b類(lèi)和c類(lèi)的計(jì)算結(jié)果基本相同。通過(guò)上述對(duì)比，考慮到計(jì)算精度和計(jì)算規(guī)模的最佳配合，本文將選用網(wǎng)格尺度b類(lèi)來(lái)計(jì)算分析氫氣/甲烷混合氣體噴射火的火焰特征。

表3 不同網(wǎng)格尺度下氫噴射火的火焰長(zhǎng)度

圖4 不同網(wǎng)格尺度下氫噴射火軸線(xiàn)溫度 分布對(duì)比

3.3 失效后果影響規(guī)律研究

對(duì)由管道噴射口射出的氫氣/甲烷混合氣體的欠膨脹噴射火的研究，會(huì)有諸多因素影響其火焰特征，本文主要研究氫氣含量、管道壓力以及噴口尺寸對(duì)氫氣/甲烷混合氣體欠膨脹噴射火火焰軸線(xiàn)溫度分布和火焰長(zhǎng)度的影響規(guī)律。

對(duì)于含氫天然氣的噴射火研究，掌握氫氣含量，即氫氣體積分?jǐn)?shù)對(duì)火焰特征的影響是首要任務(wù)。本節(jié)以工作壓力10 MPa的管道為例，假設(shè)泄漏噴口直徑5 mm，管道內(nèi)部氣體溫度293.15 K，研究常溫、常壓環(huán)境下(p∞=0.101 MPa，T∞=293.15 K)不同氫氣含量，即氫氣含量分別為0,20%,40%,60%,80%及100%，對(duì)氫氣/甲烷噴射火的影響規(guī)律。模擬中采用Birch模型確定的混合氣體在不同氫氣含量介質(zhì)對(duì)應(yīng)的偽噴口狀態(tài)，再進(jìn)行模擬計(jì)算。

圖5示出了不同氫氣含量所對(duì)應(yīng)噴射火火焰軸線(xiàn)溫度分布。可以看出，火焰軸線(xiàn)上的溫度在遠(yuǎn)離噴口方向上先升高后降低，隨著氫氣含量不斷增加，火焰軸線(xiàn)上的最高溫度依次升高，純氫氣火焰軸線(xiàn)最高溫度相對(duì)于純甲烷提高了8.0%。

圖5 不同氫氣含量下火焰軸線(xiàn)溫度分布

圖6示出了不同氫氣含量下對(duì)應(yīng)的噴射火火焰長(zhǎng)度及計(jì)算得到的混合氣體燃燒放出的凈功率。可以看出，混合氣體燃燒放出的功率變化趨勢(shì)與火焰長(zhǎng)度變化趨勢(shì)一致，由此可知，噴射火火焰長(zhǎng)度主要與氣體燃燒釋放的功率有關(guān)。

工作壓力是管道操作過(guò)程中的重要參數(shù)，同時(shí)也是影響氫氣/甲烷混合氣體噴射火火焰特征的重要參數(shù)之一。本文假設(shè)管道上的泄漏噴口直徑為5 mm，管道內(nèi)部氣體溫度為293.15 K，研究常溫、常壓環(huán)境下不同管道壓力，即泄漏壓力分別為2，4，6，8，10 MPa時(shí)，對(duì)氫氣/甲烷噴射火的影響規(guī)律。

圖7示出不同管道壓力下火焰軸線(xiàn)溫度分布，其中圖7(a)，(b)，(c)分別為甲烷、含50vol%氫氣的氫氣/甲烷混合氣體(50%H2)及純氫氣的模擬結(jié)果，經(jīng)對(duì)比可以發(fā)現(xiàn)，圖7(a)，(b),(c)火焰軸線(xiàn)上的溫度隨壓力變化的趨勢(shì)一致，均為在遠(yuǎn)離噴口方向先升高后降低，且最高溫度基本不變，分別在2 270，2 330，2 460 K附近。由此可知，火焰線(xiàn)軸上最高溫度與管道壓力無(wú)關(guān)，但與混合氣體中氫氣的含量有關(guān)，隨著氫氣含量的增加、最高溫度升高，這主要是因?yàn)闅錃獾膿饺爰觿×嘶旌蠚怏w的燃燒反應(yīng)。

圖8示出了不同管道壓力下甲烷、含50vol%氫氣的氫氣/甲烷混合氣體以及純氫氣噴射火的火焰長(zhǎng)度?？梢钥闯觯瑢?duì)于不同氫氣含量的混合氣體，噴射火火焰長(zhǎng)度隨著泄漏壓力的升高而增加，但火焰長(zhǎng)度增速隨著混合氣體中氫氣含量的升高逐漸減緩，這主要與壓力升高引起的質(zhì)量泄漏率增大有關(guān)。

欠膨脹噴射火的試驗(yàn)研究發(fā)現(xiàn)，噴射火火焰特征受?chē)娍谛螤詈统叽绲挠绊懨黠@。本部分研究具有不同直徑的圓形噴口對(duì)氫氣/甲烷噴射火火焰特征的影響規(guī)律，噴口直徑分別取5，7，10，15 mm。仍以工作壓力10 MPa的管道為例，管道溫度為293.15 K，周?chē)h(huán)境為常溫、常壓環(huán)境。

(a)CH4

(b)50%H2

(c)H2

圖9示出不同氫氣含量所對(duì)應(yīng)噴射火火焰軸線(xiàn)溫度分布，其中圖9(a)，(b)，(c)分別為甲烷、含50vol%氫氣的氫氣/甲烷混合氣體及純氫氣的模擬結(jié)果，可以看出，火焰軸線(xiàn)上的溫度在遠(yuǎn)離噴口方向上先升高后降低，不同噴口尺寸對(duì)應(yīng)火焰軸線(xiàn)上的最高火焰溫度基本不變，分別在2 270,2 330,2 460 K附近。由此可知，火焰軸線(xiàn)上最高溫度與噴口尺寸無(wú)關(guān)，但與混合氣體中氫氣的含量有關(guān)，隨著氫氣含量的增加，最高溫度升高，這主要是因?yàn)闅錃獾膿饺爰觿×嘶旌蠚怏w的燃燒反應(yīng)。

圖8 不同管道壓力下的火焰長(zhǎng)度

(a)CH4

(b)50%H2

(c)H2

圖10示出了不同噴口尺寸下甲烷、含50vol%氫氣的氫氣/甲烷混合氣體以及純氫氣噴射火的火焰長(zhǎng)度?？梢钥闯觯瑢?duì)于不同氫氣含量的混合氣體，噴射火火焰長(zhǎng)度隨著噴口直徑的增大而增加，且火焰長(zhǎng)度與噴口直徑基本呈線(xiàn)性關(guān)系增大，這主要是由于噴口直徑增大、使得質(zhì)量泄漏率增大引起的。

圖10 不同噴口尺寸下的火焰長(zhǎng)度

4 結(jié)論

本文基于計(jì)算流體軟件建立了含氫天然氣(氫氣/甲烷混合氣體)欠膨脹噴射火模擬模型，利用該模型分析研究了含氫天然氣的欠膨脹噴射火行為和失效后果，研究要點(diǎn)及結(jié)論如下。

(1)計(jì)算模型中采用Birch模型描述噴口處的工況，通過(guò)與Schefer和Studer等學(xué)者開(kāi)展的氫氣及氫氣/甲烷混合氣體的噴射火試驗(yàn)的對(duì)比，表明模擬結(jié)果和試驗(yàn)結(jié)果具有較好的一致性，驗(yàn)證了采用Birch模型的欠膨脹噴射火模型能夠較好地模擬氫氣/甲烷混合氣體的欠膨脹噴射火。

(2)氫氣/甲烷混合氣體噴射火火焰軸線(xiàn)溫度分布的模擬結(jié)果表明，火焰軸線(xiàn)上的溫度分布在遠(yuǎn)離噴口方向上呈先升高、后降低的變化趨勢(shì)，軸線(xiàn)上的最高溫度大小僅與混合氣體中氫氣的含量有關(guān)；氫氣/甲烷混合氣體噴射火火焰長(zhǎng)度的模擬結(jié)果表明，火焰長(zhǎng)度與混合氣體中的氫氣含量、管道壓力以及噴口尺寸的大小有關(guān)。

(3)氫氣含量的增加使得氫氣/甲烷混合氣體噴射火火焰軸線(xiàn)上的最高溫度升高，這主要是因?yàn)闅錃庀鄬?duì)于甲烷氣體具有更高的擴(kuò)散系數(shù)和更快的火焰?zhèn)鞑ニ俣龋瑲錃夂康纳呒觿×嘶旌蠚怏w燃燒反應(yīng)的劇烈程度，故導(dǎo)致溫度升高；噴射火火焰長(zhǎng)度則主要取決于泄漏氣體燃燒放出的功率大小，泄漏氣體燃燒放出的功率越大，則火焰長(zhǎng)度越長(zhǎng)，因此，火焰長(zhǎng)度隨管道壓力、噴口尺寸的增大而增加。