李國慶，杜 揚，齊 圣,2，王世茂，李 蒙，李 潤

(1.陸軍勤務(wù)學院油料系，重慶 401311;2.62250部隊， 青海 格爾木 816099)

油氣爆炸是發(fā)生在化工領(lǐng)域和儲油場所最嚴重的災(zāi)害之一[1-2]。在實際的儲油場所，比如油料洞庫、覆土油罐巷道中不可避免地存在各種障礙物，如通風管道、各種閥門以及通道面積突然縮小的通道口(比如通道防護門)等。當火焰在傳播過程中與此類障礙物相互作用時，火焰形態(tài)、火焰?zhèn)鞑ニ俣?、流場結(jié)構(gòu)及爆炸超壓等與沒有障礙物的空間相比都會發(fā)生顯著的變化[3-5]。因此，研究火焰與障礙物的相互作用機制及由此引發(fā)的流場結(jié)構(gòu)變化對于有效預防油氣爆炸事故的發(fā)生及降低災(zāi)害損失意義重大。

對置障管道內(nèi)甲烷、氫氣和丙烷等可燃氣體的爆炸特性已有大量的研究[4,6-13]，研究重點主要為障礙物對火焰的湍流加速機理以及爆炸流場特性等。但受限于現(xiàn)有實驗方法和測試手段，很難捕捉到詳實的火焰?zhèn)鞑ミ^程和湍流流場。特別是在數(shù)值模擬方面，對于可燃氣體爆炸過程的研究大都采用傳統(tǒng)的雷諾平均模擬(RANS)方法[14-16]，其主要缺點在于它只能提供湍流的平均信息，這很難滿足研究置障管道內(nèi)油氣爆炸的復雜湍流過程的需求，無法有效預測油氣爆炸火焰-湍流耦合的行為。

大渦模擬將耗散尺度的脈動進行過濾，只對大尺度脈動進行求解，其空間分辨率比雷諾平均更高，可以有效捕捉油氣爆炸過程中的湍流特征。目前，關(guān)于置障管道內(nèi)油氣爆炸的大渦模擬研究偏少，油氣爆炸火焰與湍流的耦合特性還亟待進一步研究。本文中針對連續(xù)圓孔障礙物條件下油氣爆炸火焰與湍流耦合過程進行大渦模擬，并將數(shù)值模擬得到的火焰形態(tài)、火焰?zhèn)鞑ニ俣群捅ǔ瑝旱忍卣鲄?shù)與實驗結(jié)果進行比較，進而揭示油氣爆炸過程中火焰-湍流耦合機理。

1 數(shù)學模型

1.1 大渦模擬控制方程

由于在置障條件下的油氣爆炸是一個強湍流爆燃過程，因此本文中采用較能捕捉湍流特征的壁面自適應(yīng)局部渦粘模型(WALE)[17]的大渦模擬(LES)對實驗過程進行重現(xiàn)。大渦模擬對小尺度的脈動進行過濾，而只對大尺度的湍流脈動通過N-S方程進行計算求解，而對小尺度脈動產(chǎn)生的影響則通過亞網(wǎng)格模型進行模型假設(shè)來模擬，經(jīng)過過濾后大渦模擬的控制方程為:

(1)

(2)

(3)

(4)

建立亞網(wǎng)格模型使控制方程封閉,WALE模型[17]通過大渦速度場動態(tài)地求出亞網(wǎng)格模型的系數(shù)，能夠較為精確地捕捉層流到湍流的轉(zhuǎn)變，且不需要顯式過濾[18]，計算量較小。因此本文采用WALE模型作為計算的亞網(wǎng)格模型。其中渦粘模型方程為:

(5)

模型中

(6)

1.2 亞網(wǎng)格燃燒模型

在湍流預混燃燒的模擬中，火焰的厚度一般較小，若直接計算，對網(wǎng)格的要求一般會比較高，因此本文中采用的燃燒模型為Zimont模型[19]，在該模型中，通過加厚火焰前鋒來計算湍流火焰速度：

(7)

2 數(shù)值模擬

2.1 物理模型

本文中大渦模擬以文獻[11]中實驗結(jié)果為基礎(chǔ)進行，實驗裝置和實驗系統(tǒng)如圖1所示。采用的物理模型和實驗管道尺寸相同，如圖2所示，管道內(nèi)部尺寸為100 mm×100 mm×1 000 mm，實驗管道材質(zhì)為有機玻璃，厚度為20 mm。在管道內(nèi)安裝了4個連續(xù)圓孔障礙物，第一個障礙物距離點火端100 mm，每個障礙物之間的間距為100 mm，圓孔半徑為40 mm，障礙物阻塞率為49.8%。管道右側(cè)端部采用盲板密封，盲板中心安裝有點火頭，距離點火頭20 mm位置安裝有一支壓力傳感器。左側(cè)端部采用薄膜密封，薄膜的破口壓力很小，對流場影響可以忽略[6]。采用實驗室自主研制的油氣霧化裝置進行配氣，實驗中初始壓力和溫度分別為實驗室環(huán)境壓力(約96.4～97.3 kPa)和溫度(約298～303 K)。采用初始體積分數(shù)φ=1.7%的油氣進行實驗，并用電火花引爆油氣。利用壓力傳感器記錄管道內(nèi)超壓-時間曲線，高速攝影儀拍攝油氣爆炸火焰發(fā)生和發(fā)展過程。為保證實驗的可靠性，每組實驗在相同工況下重復進行至少3次。

在數(shù)值計算時，為了減小管道出口氣體回流對管道內(nèi)部壓力的影響，物理模型在管道外沿x、y、z方向分別延伸500、500、1 000 mm，即設(shè)置了一個尺寸為500 mm×500 mm×1 000 mm的外部空間區(qū)域。這個外部空間區(qū)域和管道出口相連接，可以模擬管道內(nèi)火焰沖出管道的反應(yīng)進程，也可以模擬管道泄壓過程。

2.2 網(wǎng)格劃分

在大渦模擬中，亞網(wǎng)格尺度的物理擴散是隨著網(wǎng)格的加密而減小的，因此，對于大渦模擬而言，并不存在通常數(shù)值計算中所謂的“網(wǎng)格無關(guān)性”概念，較好的做法是采用十分精細的網(wǎng)格和很小的時間步長，以確保其數(shù)值擴散遠小于雷諾平均方法。但是，在實際計算中，由于受到計算機計算能力的限制，通常仍然會進行網(wǎng)格無關(guān)性檢驗，兼顧計算精度和計算成本。鑒于此，本文中在網(wǎng)格無關(guān)性檢驗的基礎(chǔ)上，采用六面體網(wǎng)格對計算模型進行網(wǎng)格劃分。整個計算區(qū)域網(wǎng)格總數(shù)為1 200 000個，其中管道內(nèi)部大約有950 000個網(wǎng)格，網(wǎng)格尺寸大小為2～2.5 mm，外部空間網(wǎng)格總數(shù)為250 000左右，網(wǎng)格尺寸為2～10 mm。

2.3 邊界條件和初始條件

管道右側(cè)封閉端、管道內(nèi)壁面和障礙物表面都設(shè)置成絕熱無滑移壁面邊界。另外，由于實驗中管道左側(cè)薄膜能夠在很低的壓力條件下破裂，對油氣泄壓爆炸過程影響較小，因此在大渦模擬計算時，忽略薄膜的影響[7]。管道外部擴展區(qū)域的壁面設(shè)置為壓力出口，壓力出口表壓設(shè)置為0 Pa。

未燃氣體和已燃氣體的比熱容按照近似溫度的五階分段多項式函數(shù)來計算。氣體分子黏度通過Sutherland法則來計算。初始反應(yīng)物選用φ(CH)=1.7%的汽油蒸汽[1-2]，層流火焰速度假設(shè)為隨溫度和壓力變化保持常數(shù)不變，為0.43 m/s[21]。

初始條件包括壓力、流場初速度、初始溫度和反應(yīng)進程變量。初始溫度設(shè)為300 K，其余初始參數(shù)設(shè)置為0。采用Patch功能實現(xiàn)點火引爆，具體操作方法為在管道右側(cè)封閉端中心位置設(shè)置一個半徑為5 mm的半球形區(qū)域，并且設(shè)置該區(qū)域的反應(yīng)進程為1，以此來實現(xiàn)模擬點火功能[7]。

2.4 計算方法

大渦模擬計算通過ANSYS公司的Fluent計算平臺來實現(xiàn)。采用SIMPLE算法對壓力和速度進行耦合計算，對流項和擴散項分別采用二階迎風格式和二階中心差分格式來進行離散。采用聯(lián)想ThinkServer TD350服務(wù)器(64位、Xeon E5-2609 v4 CPU (8 processes)、8G RAM)實現(xiàn)模擬計算。迭代2計算步長設(shè)置為1×10-5s，每一個時間步長要求迭代20次，保證能量方程殘差小于1×10-6，反應(yīng)進程變量方程殘差小于 1×10-4，動量方程殘差小于1×10-5。完成一次模擬計算的時間大約為52 h。

3 模型驗證

以往大多數(shù)模擬研究，采用將模擬計算結(jié)果和實驗結(jié)果進行對比來說明數(shù)值模擬模型的可靠性[22-23]。本文中采用對比數(shù)值計算和實驗所得火焰?zhèn)鞑バ螒B(tài)、壓力-時序曲線、火焰鋒面位置和火焰?zhèn)鞑ニ俣葋砼袛鄶?shù)值模擬計算結(jié)果的精度。為了證明大渦模擬的精確性，在數(shù)值模擬計算時還采用了RNGk-ε湍流模型來進行計算，并將計算結(jié)果和大渦模擬計算結(jié)果進行對比，相關(guān)的湍流控制方程[14]在此不做贅述。

圖3為選取的模擬計算和實驗結(jié)果中8個不同時刻火焰?zhèn)鞑バ螒B(tài)對比圖。從模擬結(jié)果可見，大渦模擬和RNGk-ε湍流模型都能比較直觀地反映出火焰在管道內(nèi)外的傳播形態(tài)變化過程。但是，大渦模擬相比RNGk-ε湍流模型能更好地表現(xiàn)出火焰結(jié)構(gòu)的褶皺、彎曲等精細化演變過程。圖4所示是數(shù)值模擬和實驗所得火焰鋒面位置隨時間的變化關(guān)系曲線，從圖4中可見，大渦模擬計算結(jié)果和實驗所得火焰鋒面位置基本一致，體現(xiàn)出大渦模擬對火焰鋒面位置較好的預測性。但是，采用RNGk-ε湍流模型計算所得的火焰鋒面位置與實驗結(jié)果相比較誤差較大，尤其是預測的火焰鋒面的移動速度明顯偏高。

圖5為數(shù)值模擬和實驗相同測點監(jiān)測所得壓力-時間曲線對比關(guān)系。從圖中可見，實驗、大渦模擬和RNGk-ε湍流模型計算結(jié)果中取得最大超壓峰值的時間分別為27.4、25.5、16 ms，后兩者與實驗結(jié)果相比較誤差分別為6.93%和41.6%。三者最大爆炸超壓峰值分別為194.77、165.35、83 kPa，后兩者與實驗結(jié)果相比較誤差分別誤差為15.1%和57.4%。由此可見，大渦模擬對爆炸超壓的預測較RNGk-ε湍流模型精度更高。值得注意的是，本文中實驗和數(shù)值模擬得到的最大爆炸超壓峰值相比其余大多數(shù)文獻中獲得的最大超壓峰值更大，原因主要為：一是本文中采用的實驗管道的長徑比相較其余大多數(shù)文獻中的管道長徑比更大，因此油氣爆炸火焰在管道內(nèi)加速距離更長，導致爆炸發(fā)展更充分，引起更劇烈的爆炸；二是本文中采用了4個連續(xù)障礙物，顯著增強了流場湍流度，提高了已燃氣體和未燃氣體的熱交換速率。由此，增強了障礙物對火焰的加速效應(yīng)，并能誘導更高的火焰?zhèn)鞑ニ俣龋鶕?jù)文獻[11]中關(guān)于油氣爆炸火焰?zhèn)鞑ニ俣群捅ǔ瑝貉葑冞^程的耦合關(guān)系可知，最大爆炸超壓峰值會變得更大。

圖6所示是實驗和數(shù)值模擬所得火焰?zhèn)鞑ニ俣葘Ρ?，火焰?zhèn)鞑ニ俣仁峭ㄟ^計算火焰鋒面在管道軸向移動速度得到，計算公式為：

Sf=(xn+1-xn)/Δtn

(8)

式中：Sf表示火焰?zhèn)鞑ニ俣?xn+1-xn表示相鄰的兩幅高速攝影照片中(在大渦模擬中指反應(yīng)進程變量c=0.5的相鄰兩幅圖像)火焰鋒面的真實距離之差, Δtn表示選取的兩幅相鄰圖像的時間差(本文中Δtn=0.001 s)。

從圖6中可見，實驗和大渦模擬所得火焰?zhèn)鞑ニ俣鹊淖兓?guī)律基本一致，在火焰鋒面沒有接觸障礙物之前，火焰保持較低速度傳播，大致為9 m/s；當火焰鋒面經(jīng)過障礙物之后，火焰速度開始急劇增大，并且在27 ms取得最大值。實驗和大渦模擬所得最大火焰速度分別為410.71、390 m/s，大渦模擬和實驗結(jié)果相對誤差為5%。當火焰速度達到最大值之后，由于可燃氣體量的減少和熱量的大量損耗，火焰?zhèn)鞑ニ俣燃眲∠陆?。對于RNGk-ε湍流模型計算結(jié)果而言，其火焰?zhèn)鞑ニ俣仍谧兓厔萆吓c實驗結(jié)果和大渦模擬計算結(jié)果相比誤差較大，不能體現(xiàn)出障礙物對火焰的加速過程，并且最大火焰速度僅為140 m/s，相比實驗結(jié)果誤差為65.9%。

由上述對比分析可見，本文中采用的大渦模擬計算模型能夠較好地預測半受限空間汽油-空氣混合物泄壓爆炸過程火焰和爆炸超壓特征。

4 置障管道油氣泄壓爆炸火焰特性分析

4.1 障礙物對火焰?zhèn)鞑バ螒B(tài)影響

圖7所示為大渦模擬所得7個不同時刻火焰?zhèn)鞑ソ?jīng)過障礙物過程中的3維立體火焰圖像，火焰結(jié)構(gòu)采用反應(yīng)進程變量為c=0.5的等值面渲染得到。從圖7可見，t=14 ms時，火焰還沒受到管道側(cè)壁和障礙物的影響，火焰鋒面保持半球形向前傳播；t=18 ms時，火焰鋒面?zhèn)鞑ソ?jīng)過第1個障礙物，火焰形狀由半球形轉(zhuǎn)變成蘑菇形，并且火焰鋒面面積顯著增大；t=20～23 ms之間，火焰穿越剩余3個連續(xù)障礙物，在這一過程中，火焰鋒面形態(tài)變得更加不規(guī)則，出現(xiàn)明顯的褶皺和卷曲變形現(xiàn)象，并且當火焰鋒面受到障礙物壁面阻礙后，還出現(xiàn)火焰向點火端回傳現(xiàn)象(如圖7中20、21 ms火焰圖像所示)。當火焰鋒面穿越4個連續(xù)障礙物之后，火焰出現(xiàn)明顯的破碎現(xiàn)象，如t=25，26 ms所示。此外，還可以觀察到火焰鋒面在傳播經(jīng)過障礙物時，在障礙物和管道內(nèi)壁形成的腔體內(nèi)還有一部分未燃氣體滯留，文獻[23-24]中的研究表明這一部分未燃氣體與爆炸超壓后期振蕩有密切的關(guān)系。

4.2 火焰和流場耦合關(guān)系分析

4.2.1 內(nèi)場火焰和流場耦合關(guān)系

圖8為管道內(nèi)部置障條件下油氣混合物預混火焰?zhèn)鞑ズ土鲌鼋Y(jié)構(gòu)圖，圖中綠色代表火焰鋒面(c=0.5)，彩色箭頭代表速度矢量。在t=14 ms之前，管道內(nèi)氣體流動是層流狀態(tài)，流場強度較小，流線形態(tài)較規(guī)則，流場方向基本保持從點火端向管道開口端發(fā)展；t=18 ms時，火焰?zhèn)鞑ソ?jīng)過第1個障礙物，流場形態(tài)開始變得不規(guī)則，并且在相鄰障礙物之間的流場區(qū)域出現(xiàn)渦旋結(jié)構(gòu)，在渦旋結(jié)構(gòu)作用下，火焰鋒面形狀從半球形向蘑菇狀轉(zhuǎn)變，增大了已燃氣體和未燃氣體的接觸面積；t=20 ms時，火焰?zhèn)鞑ソ?jīng)過第2個障礙物，此時流場強度進一步增大，渦旋結(jié)構(gòu)更加明顯，并且在渦旋結(jié)構(gòu)的卷曲作用下，火焰在第1～2個障礙物之間出現(xiàn)向點火端回傳現(xiàn)象，同時可見此時火焰鋒面前部的流場強度更加劇烈；t=21～23 ms之間，火焰?zhèn)鞑ソ?jīng)過第3～4個障礙物，在這個過程中，火焰結(jié)構(gòu)變得更加褶皺、彎曲，增大了已燃氣體與未燃氣體的接觸面積，提高了燃燒速率，增大了火焰?zhèn)鞑ニ俣龋虼藢е禄鹧驿h面前方流場渦旋更加明顯，流場強度急劇增大。同時，隨著流場強度的增大，渦旋更加明顯，對氣體的卷曲作用明顯增強，因此提高了已燃氣體和未燃氣體之間的交換速率，反過來提高燃燒速率和火焰?zhèn)鞑ニ俣?。因此，在火焰?zhèn)鞑ズ土鲌錾鲜稣钭饔玫恼T導下，火焰形態(tài)產(chǎn)生嚴重變形，火焰鋒面前流場強度急劇增大。

4.2.2 外場火焰和流場耦合關(guān)系

圖9所示為外場火焰和流場結(jié)構(gòu)圖。t=26 ms以前，火焰還沒傳播至管道開口端，但是從t=21 ms起，管道外流場逐漸形成“雙渦旋”結(jié)構(gòu)，且隨著時間的推移，渦旋結(jié)構(gòu)越來越顯著，并且以管道中軸線為分界線，中軸線上部渦旋方向為順時針，中軸線下部渦旋方向為逆時針方向。t=26.5 ms時，火焰以射流火焰形式從管道內(nèi)部傳播到管道外部，并進入管道外部渦旋區(qū)域，在渦旋作用下，火焰鋒面在t=27 ms時從柱狀火焰向蘑菇狀火焰轉(zhuǎn)變，并于t=28 ms時形成形狀較完整的蘑菇狀火焰。通過圖中流場矢量圖顏色可以判斷得到，在蘑菇狀火焰區(qū)域，軸向流場速度明顯大于徑向流場速度，可以推斷，蘑菇狀火焰的形成和流場在軸向和徑向的速度差存在密切的內(nèi)在聯(lián)系。此后，火焰鋒面的蘑菇狀火焰逐漸向啞鈴狀火焰轉(zhuǎn)變(t=30 ms所示)，并且火焰鋒面向前傳播的同時，渦旋結(jié)構(gòu)也向前移動，并且始終和火焰鋒面重合。

4.3 爆炸超壓與火焰?zhèn)鞑ヱ詈详P(guān)系

爆炸超壓是油氣爆炸過程中一個重要的爆炸參數(shù)，為了研究其演變規(guī)律與火焰?zhèn)鞑ブg的耦合關(guān)系，根據(jù)爆炸超壓曲線的變化特點將其劃分為圖10(a)中所示的4個變化階段。在分析時，將大渦模擬計算得到的火焰?zhèn)鞑ニ俣取⒒鹧婷娣e與爆炸超壓進行聯(lián)合分析，對應(yīng)的曲線如圖10所示?；鹧婷娣e根據(jù)下式計算：

(9)

式中：Ai代表等值面上微單元的面積；n代表等值面上微單元的數(shù)量(注：計算時等值面取c=0.5)。

由圖10可見，在15 ms以前(階段Ⅰ)，爆炸超壓平緩上升，對應(yīng)的火焰?zhèn)鞑ニ俣群突鹧婷娣e變化幅度也很小，此時火焰還沒受到障礙物的擾動，火焰處于層流傳播階段。從15～25 ms(階段Ⅱ)，爆炸超壓急劇上升，并在大約25 ms處取得最大值。在這一階段，火焰穿越障礙物，火焰發(fā)生卷曲變形，層流火焰向湍流火焰轉(zhuǎn)變，導致火焰?zhèn)鞑ニ俣群突鹧婷娣e急劇增大(如圖10(b)～(c)所示)，提高燃燒速率和熱釋放率，進而增大爆炸超壓。在第Ⅲ階段(大約25～30 ms)，爆炸超壓和火焰?zhèn)鞑ニ俣榷技眲∠陆?，并且都在大約30 ms左右取得極小值，同時火焰面積也小幅下降。在前三個階段，爆炸超壓、火焰?zhèn)鞑ニ俣群突鹧婷娣e具有“同升同降”的相似變化趨勢，耦合關(guān)系較好。在第Ⅳ階段，火焰鋒面?zhèn)鞑サ焦艿雷钭蠖?，火焰逐漸熄滅，火焰面積在短時間內(nèi)保持上升趨勢，隨后急劇下降。隨著火焰面積的改變，爆炸超壓呈現(xiàn)出振蕩下降的趨勢，并形成多個超壓峰值，其中振蕩期的第一個超壓峰值對應(yīng)火焰面積的最大峰值。由上述分析可見，油氣在置障管道內(nèi)泄壓爆炸過程的爆炸超壓、火焰?zhèn)鞑ニ俣群突鹧婷娣e之間耦合關(guān)系明顯，三者的變化趨勢存在較好的一致性。

5 結(jié) 論

(1)通過比較大渦模擬和RNGk-ε湍流模型計算結(jié)果發(fā)現(xiàn)，大渦模擬在預測油氣爆炸超壓、火焰?zhèn)鞑ニ俣燃盎鹧嫘螒B(tài)變化等方面比RNGk-ε湍流模型精確度更高，且能表現(xiàn)出更多流場的精細化結(jié)構(gòu)；

(2)障礙物誘導管道內(nèi)形成湍流度較高的流場區(qū)域，導致火焰產(chǎn)生褶皺彎曲變形，增大火焰面積，加速火焰?zhèn)鞑ィ?/p>

(3)通過比較分析爆炸超壓、火焰?zhèn)鞑ニ俣群突鹧婷娣e的演變規(guī)律，發(fā)現(xiàn)三者內(nèi)在聯(lián)系密切，具有顯著的耦合性，且隨時間的變化趨勢存在高度的一致性。