郭涵予　陶 剛　張禮敬　嚴景藝

南京工業(yè)大學安全科學與工程學院江蘇省城市與工業(yè)安全重點實驗室

球形容器內(nèi)H2-Air爆燃特性的數(shù)值模擬*

郭涵予陶 剛張禮敬嚴景藝

南京工業(yè)大學安全科學與工程學院江蘇省城市與工業(yè)安全重點實驗室

為了研究氫氣在密閉容器中燃燒的發(fā)展過程，基于20L球形容器中H2-air爆炸實驗（Crowl and Jo，2009），采用數(shù)值模擬的方法對密閉容器中Φ=1的混合氣體的爆炸壓力和火焰鋒面位置進行了詳細分析。通過對比實驗和模擬的結(jié)果可知：燃燒過程中受重力影響使得不同方向的火焰?zhèn)鞑ニ俣炔煌?，導致實際燃燒過程與模擬相比較為緩慢；壁面的阻礙作用導致火焰在密閉容器內(nèi)傳播速度先增加后降低；實驗過程中的熱耗散和浮力影響是導致誤差的主要原因。

預混火焰；爆炸壓力；數(shù) 值模擬；k-ε湍流模型

在國外，Liu等[1]研究了初始壓力對點火延滯期和點燃溫度的影響，探討氫氧均勻混合體系中基元反應對第三爆炸極限的影響。Young-Do Jo[2]利用火焰增長模型研究了密閉容器內(nèi)H2-air混合體系的爆炸壓力-時間曲線以及最大爆炸壓力、爆燃指數(shù)和燃燒速率等隨濃度的變化關系，結(jié)果表明最大爆炸壓力并非出現(xiàn)在化學當量比處。Tao等[3-4]采用化學平衡計算的方法研究了固定氮氣濃度下氫氣最大爆炸壓力和爆燃指數(shù)的變化規(guī)律，并擬合了相關的經(jīng)驗公式。在國內(nèi)，董剛[5]等采用二維Euler方程求解了H2-air混合氣體的瞬態(tài)爆轟過程，分析了爆轟波陣面的變化特征以及爆轟波后的組分濃度及溫度分布。李書明等[6]采用Fluent模擬了常溫常壓下H2-air預混火焰在光滑管道中的傳播特性，分析了火焰?zhèn)鞑ニ俣?、火焰結(jié)果、表面積等的變化情況，分析了郁金香火焰形成的原因。孫伯剛等[7]采用CFD仿真軟件對氫氣層流燃燒進行了數(shù)值模擬，采用TFSC模型對定容燃燒彈進行了模擬，結(jié)果在常壓下、溫度范圍300～900K，當量比的范圍0.5～1.2時較為準確。付佳佳等[8]采用氫的單步化學反應模型對不同的噴射速度和噴口直徑的噴射火進行了大渦模擬，結(jié)果表明噴射速度和噴口直徑的增加會導致火焰高度增加。

目前國內(nèi)外的研究大多是圍繞氫氣燃燒過程中的燃燒速度及爆炸壓力隨濃度的變化特點展開研究，而物理模型也多為管道模型，對于流動模型則大多采用層流燃燒模型進行分析，因此本文基于20L球形容器中的H2-air混合氣體爆炸實驗數(shù)據(jù)（Crowl and Young-Do Jo）進行研究，運用CFD軟件進行模擬，基于標準湍流k-ε模型采用組分概率密度輸運燃燒模型（composition PDF transport）對球形容器中的二維爆炸過程進行數(shù)值分析，探討爆炸過程，通過對比實驗結(jié)果分析爆炸特性參數(shù)的變化特點。

1 實驗數(shù)據(jù)分析

實驗裝置采用20L的球形爆炸容器，初始壓力為1atm，初始溫度為300k，點火能量大約為10J，點火位置為裝置中心點火，壓力傳感器位于球形容器點火位置右側(cè)壁面處[9]。在實際燃燒過程中，點火以后火焰開始由中心位置向四周傳播，火焰內(nèi)核受到浮力的影響導致上浮，因此在燃燒過程中火焰向上的傳播速度最大，向下的傳播速度最小，位于容器上方的預混氣體優(yōu)先反應，因此火焰形狀并非規(guī)則的球形[10]。圖1為化學當量比為1處的壓力隨時間變化曲線。壓力波的傳播滯后于火焰?zhèn)鞑ィ捌趬毫Σ蛔?，隨著燃燒的不斷加快，火焰鋒面對未燃預混氣體產(chǎn)生擠壓，大約在10ms處壓力波傳至壓力傳感器出，壓力開始上升。燃燒過程中壓力逐漸增大，由于壁面的阻礙作用，導致火焰?zhèn)鞑ニ俣冉档?，但隨著壓力的上升化學反應速率加快，在17ms達到最大爆炸壓力，此時預混氣體幾乎全部參與反應，壓力波峰值達到最大。17ms后由于火焰?zhèn)鞑ミ^程中容器壁面的反射作用，導致壓力波疊加，以及爆炸過程中的聲波影響，導致最大壓力出現(xiàn)振蕩，隨著燃燒的結(jié)束，容器內(nèi)燃燒壓力逐漸趨向穩(wěn)定，此時最大爆炸壓力約為7.8atm，隨著燃燒的結(jié)束，熱損耗導致燃燒產(chǎn)物溫度逐漸降低，最大爆炸壓力逐漸減小。

圖1　實驗中壓力隨時間變化曲線

2 控制方程和數(shù)值方法

2.1控制方程

為了簡化計算過程，根據(jù)球形對稱性的幾何關系，在進行數(shù)值模擬時將物理模型簡化為2D模型，控制方程如下[11]。

質(zhì)量方程，又稱為連續(xù)性方程：

動量守恒方程：

能量守恒方程：

式中：

t—時間，s；

ρ—密度，g/m3；

E—流體微團的總能，J；

p—靜壓，Pa；

keff—有效熱傳導系數(shù)，W/(m2·K)；

Jj—組分j的擴散通量，mol/(m2·s)；

Sh—包含化學反應熱以及其他用戶定義的體積熱源項，J。

湍流方程，k-ε模型是目前使用最廣泛的湍流模型，其中標準k-ε模型只適合完全湍流的流動過程模擬，其輸運方程為：

式中：

G—平均速度梯度引起的湍動能，J；C1、C2—經(jīng)驗常數(shù)。

2.2燃燒模型與初始條件

可燃氣體的燃燒過程伴隨著化學反應的發(fā)生，所以在模擬過程中應當選擇合適的化學反應來表征燃燒過程，而在火焰?zhèn)鞑ルA段，由于火焰的拉伸效果導致燃燒并非穩(wěn)定的層流傳播，湍流模型更有助于對燃燒過程的模擬。因此本文采用標準k-ε模型作為湍流模型，化學反應模型選擇組分概率密度輸運燃燒模型(composition PDF transport)進行模擬，該模型能夠根據(jù)導入的反應機理模擬湍流火焰中的有限速率反應，能夠根據(jù)導入的反應機理對容器中混合氣體的組分概率密度進行計算，可以求解預混燃燒，但對計算量的要求較高，適合在2D條件下進行模擬。模擬選用Chemkin中的氫氧反應機理，該機理包含了9種物質(zhì)的19個基元反應，能夠很好的表征氫氧反應中各組分反應的變化特征，見表1。物理模型基于Jo和Crowl的實驗裝置[9]，即選用20L球形爆炸容器。模擬選取化學當量比Φ=1的理想H2-air混合氣體，容器壁面均做絕熱處理，初始條件為常溫常壓，中心位置賦予2000K的高溫區(qū)域作為點火源，監(jiān)測點位于容器右側(cè)壁面位置，與實驗裝置保持一致。

表1　氫氣-氧氣反應動力學機理

3 模擬結(jié)果分析

根據(jù)上述條件，本文采用SIMPLE算法對燃燒過程進行模擬，燃燒過程中的壓力變化如圖2所示。從圖中可以得知，壓力波在前3ms尚未傳播至監(jiān)測點。在3ms時壓力開始上升，隨著反應的進行，燃燒受到湍流的影響導致火焰并非層層傳播，而是處于擴散傳播，因此前期壓力變化并非一個穩(wěn)定上升過程。模擬過程為高溫引燃，因此壓力上升從2ms開始，經(jīng)過6ms達到最大爆炸壓力，并最終穩(wěn)定在7.9atm，與實驗得到的最大爆炸壓力相比，誤差僅為1.2%。由于模擬過程中未考慮壁面反射波的疊加效果，因此最終壓力與實驗結(jié)果相比，沒有明顯的震蕩效果。

圖3為不同時間點火焰的傳播過程。模擬初始時刻在中心位置規(guī)定一個2000K的高溫區(qū)域，0.5ms時火焰開始向周圍傳播，周圍預混氣體的溫度開始上升，此時火焰?zhèn)鞑ニ俣容^慢，火焰形成初期，火焰反應區(qū)和預熱區(qū)均較為狹窄，隨著火焰的傳播，燃燒速度逐漸增加，火焰的反應區(qū)和預熱區(qū)逐漸增厚，4ms時火焰的預熱區(qū)相比于初始階段已經(jīng)十分明顯，此時的火焰?zhèn)鞑ニ俣容^快，火焰前鋒即將到達容器壁面位置。隨后由于壁面反射波的作用，導致火焰?zhèn)鞑ニ俣戎饾u降低，5.85ms時容器內(nèi)最低溫度為1290K，此時容器內(nèi)預混氣體已經(jīng)全部點燃。

圖2　數(shù)值模擬中壓力隨時間的變化

圖3　不同時刻的火焰?zhèn)鞑ミ^程

圖4為不同時刻火焰前鋒的位置，由圖知，初始時刻火焰?zhèn)鞑ポ^慢，隨著火焰的傳播，火焰鋒面的預熱區(qū)逐漸增大，且溫度逐漸增高，而隨著未燃區(qū)域的溫度逐漸增加，導致未燃區(qū)域與已燃區(qū)域的溫差逐漸減小，火焰?zhèn)鞑ニ俣戎饾u增加[12]，但當火焰?zhèn)鞑ミ_到4ms時，由于容器壁對壓力波的阻礙作用，火焰?zhèn)鞑ニ俣戎饾u降低，壁面引起的反射波導致火焰?zhèn)鞑ニ俣冉档停瑫r湍流效果增強，火焰緩慢向壁面?zhèn)鞑ィ也⒎菍訉觽鞑?，最終引燃整個區(qū)域的預混氣體。

圖5為實驗數(shù)據(jù)與模擬數(shù)據(jù)的對比分析圖，由圖可知，模擬的爆炸過程中壓力上升速率高于實驗記錄的壓力上升速率，這主要是由于兩方面原因?qū)е碌?。燃燒初期，火焰受到空氣浮力的作用導致火焰?nèi)核上升，此時火焰?zhèn)鞑シ较蛑饕譃橄蛏蟼鞑?、向下傳播和四周傳播，其中燃燒速度差別較大的是向上和向下傳播，不同的燃燒方向?qū)е禄鹧鎯?nèi)核下方的未燃氣體并不能及時被引燃，因此在燃燒前期，實驗過程中的已燃區(qū)域要小于模擬得到的已燃區(qū)域，因此壓力上升速率也會低于模擬得到的數(shù)據(jù)。在模擬過程中，選取的假設條件為不可壓縮的理想氣體，且反應容器為絕熱反應容器，因此在模擬過程中容器內(nèi)外沒有發(fā)生熱交換，但在實驗過程中會有一部分熱損失，這會導致已燃氣體產(chǎn)生的能量不能全部用來引燃未燃氣體，因此使得容器壁面附近的預混氣體不能及時點燃。

圖4　不同時刻火焰前鋒的位置

圖5　實驗值與模擬值對比

4 結(jié)論

本文通過采用CFD軟件模擬H2-air預混氣體爆炸過程，對比分析實驗數(shù)據(jù)與模擬結(jié)果，主要結(jié)論：

（1）實驗過程中預混氣體受到重力作用導致火焰內(nèi)核上浮，致使球形容器中的燃燒方式為非球形燃燒，導致燃燒速度降低。同時由于壓力波的疊加以及聲波的影響，使得最大爆炸壓力出現(xiàn)振動。

（2）模擬結(jié)果表明，忽略重力影響下火焰呈球形燃燒，在燃燒過程中火焰?zhèn)鞑ニ俣认纫蚧鹧婕铀贆C理而增加，后因壁面阻礙壓力波的傳播而降低。

（3）與實驗相比，模擬最終達到的爆炸壓力與實驗結(jié)果誤差僅為1.2%，但由于模擬不能考慮壁面反射效果以及熱耗散，因此實驗達到最大爆炸壓力所需時間比模擬達到最大爆炸壓力所需時間長。

[1] Liu Y,Pei P.Asymptotic Analysis on Autoignition an Explosion Limits of Hydrogen–oxygen Mixtures in Homogeneous Systems[J].International Journal of Hydrogen Energy,2006,31(5):639-647

[2] Jo Y D,Crowl D A.Explosion Characteristics of Hydrogenair Mixtures in a Spherical Vessel[J].Process Safety Progress,2009,29(3):216-223

[3] Tao G,Crowl D A.Comparison of the Maximum Gas Combustion Pressure of Hydrogen/Oxygen/Nitrogen-between Chemical Equilibrium Calculations and Experimental Data [J]. Procedia Engineering,2013,62:786-790

[4] Tao G,Crowl D A,Guo HY.Equations for Flammability Parameters,Pmax and KG[C].American Institute of Chemical Engineers/Global Congress on Process Safety,New Orleans,2014

[5] 董剛,范寶春,謝波.氫氣-空氣混合物中瞬態(tài)爆轟過程的二維數(shù)值模擬[J].高壓物理學報,2004,18(1):40-46

[6] 李書明,張旭龍,程關兵.氫氣-空氣預混火焰?zhèn)鞑ヌ匦缘臄?shù)值模擬研究[J].裝備制造技術,2013(11):23-26

[7] 孫柏剛,趙建輝,劉福水.預混氫氣層流燃燒速度的數(shù)值模擬[J].燃燒科學與技術,2010,16(5):430-435

[8] 付佳佳,王昌建,秦俊,等.氫氣噴射火的大渦模擬[J].燃燒科學與技術,2013,19(5):473-477

[9] Crowl D A,Jo Y.A Method for Determining the Flammable Limits of Gases in a Spherical Vessel[J]．Process Safety Progress,2009,28(3):227-236

[10] Kuznetsov M,Kobelt S,Grune J,et al.Flammability Limits and Laminar Flame Speed of Hydrogen–air Mixtures at Subatmospheric Pressures [J].International Journal of Hydrogen Energy,2012,37(22):17580–17588

[11] 朱紅鈞.FLUENT流體分析及仿真實用教程[M].北京:人民郵電出版社,2010

[12] 岑可法.燃燒理論與污染控制[M].北京:機械工業(yè)出版社, 2004

江蘇省2014年度普通高校研究生科研創(chuàng)新計劃項目(KYLX_0783)