非貯壓手持式滅火器氣固兩相流數(shù)值模擬仿真研究

曹書華，莊弘煒，余鵬飛

?

非貯壓手持式滅火器氣固兩相流數(shù)值模擬仿真研究

曹書華，莊弘煒，余鵬飛

（武警工程大學，陜西西安，710086）

對非貯壓手持式滅火器的氣固兩相流進行研究，運用Flunet軟件進行數(shù)值模擬仿真，通過氣體速度矢量圖和顆粒速度流線圖，分析了兩相流體在筒體內(nèi)和外的運動狀態(tài)，并通過試驗加以驗證。結(jié)果表明仿真與試驗符合程度較高，本研究為滅火器結(jié)構(gòu)優(yōu)化提供了理論指導。

滅火器；氣固兩相流；仿真；試驗

非貯壓手持式滅火器是一種新型的反自焚滅火裝備，其原理為產(chǎn)氣藥劑燃燒生成的煙火型氣溶膠氣體將筒體內(nèi)部的超細干粉全部穩(wěn)定、持續(xù)、迅速地噴發(fā)出去，作用于人體或其它物體的火焰燃燒表面。由于超細干粉的噴射是一個氣固兩相流運動，其噴射流場的特性是影響滅火效果的關(guān)鍵。因此，需要對非貯壓手持式滅火器的氣固兩相流進行研究，為滅火器性能改進和實戰(zhàn)應用提供技術(shù)指導。

1 氣固兩相流數(shù)值模型分析

非貯壓手持式滅火器使用時，其滅火劑噴射過程是一種具有復雜的顆粒運動規(guī)律的稠密氣固兩相流動，采用離散單元法（Discrete Element Method，DEM）對氣固兩相流流場進行描述[1]，即對氣相流體采用宏觀尺度上的描述，在歐拉坐標系下考察流體的運動情況，對固相顆粒采用顆粒尺度上的描述，在拉格朗日坐標系下考察顆粒的運動情況[2-5]。

（1）氣相控制方程

根據(jù)質(zhì)量守恒定律和牛頓第二定律，在歐拉坐標下描述氣相性質(zhì)。

（2）固相控制方程

滅火劑顆粒在流場中的運動可以分作平動運動和轉(zhuǎn)動運動，并進行非勻速運動。為了更準確描述滅火劑顆粒的運動狀態(tài)，把整個滅火劑顆粒的非勻速運動看作無數(shù)個勻速運動組成，即用微分的概念，將滅火劑顆粒在每一個非常小的時間間隔△內(nèi)看作勻速運動[6-7]。可得，滅火劑顆粒在三維坐標中的平動和轉(zhuǎn)動方程：

式（2）~（8）中：為滅火劑顆粒位移；f為滅火劑顆粒的碰撞力；f為氣體對滅火劑顆粒的總作用力；為流場的壓力梯度；V為單個滅火劑顆粒的體積；為滅火劑顆粒轉(zhuǎn)動的角速度；T為滅火劑顆粒受到的轉(zhuǎn)動合力矩；為滅火劑顆粒的轉(zhuǎn)動慣量。

2 氣固兩相流數(shù)值模擬仿真過程

（1）幾何模型建立：根據(jù)筒體內(nèi)腔的尺寸，建立幾何模型的主體，長162mm，直徑47mm，噴口處內(nèi)錐角為130°，噴口直徑7mm。對產(chǎn)氣發(fā)生器結(jié)構(gòu)進行忽略，改為在模型底部設(shè)置一個進氣口。已知產(chǎn)氣發(fā)生器上有8個產(chǎn)氣孔，每個孔直徑5mm，因此進氣口大小近似為8個產(chǎn)氣孔大小之和，其直徑為17mm。器體模型如圖1所示。

圖1 筒體模型

由試驗觀察，將流場區(qū)域大小設(shè)置為4m×3m ×3m，滅火器噴口方向為水平方向。整個計算區(qū)域模型如圖2所示。

圖2 計算區(qū)域模型

（2）網(wǎng)格劃分：整個非貯壓手持式滅火器氣固兩相流場的幾何模型具有復雜的外形，因此用T-grid方法對其進行四面體網(wǎng)格劃分。為了更精確地研究滅火器的氣固兩相流動，在劃分計算區(qū)域網(wǎng)格時，最大體積約8.5×10-3m3，最小體積約6.3×10-11m3，總數(shù)約38萬左右。

（3）求解器選擇：由試驗已知滅火器的噴射持續(xù)時間約為1.13s，結(jié)合其噴射距離，可知兩相流動狀態(tài)為低速流動。所以，模擬采用基于壓力的求解器。該求解器適用于中等可壓縮流動，適合用于低速流動的CFD模型。FLUENT的求解流程圖如圖3所示。

（4）模擬步驟：模擬滅火劑從滅火器筒體中噴出分為2個步驟：第1步是進行穩(wěn)態(tài)模擬，獲得氣流場；第2步是在已經(jīng)獲得的氣流場中加入離散相，模擬追蹤離散相顆粒的運動軌跡。

圖3 FLUENT求解流程圖

穩(wěn)態(tài)模擬計算時，采用標準-湍流模型[3]，壁面的流動也采用標準壁面方程。氣體入口給定特定的氣體組分以及速度大小，氣體組分按照“12.1%K2O、26.6%N2、24.7%H2O、25.2%CO、2.7%CO2”在“Species”中給出。速度大小為11.5m/s，這是能將滅火劑顆粒在1.13s左右完全從滅火器筒體中吹完的速度大小。初始時刻滅火器筒體和計算區(qū)域都是空氣。邊界條件設(shè)置時，將計算區(qū)域的地面設(shè)成“wall”，代表地面；計算區(qū)域其它面則都設(shè)成0壓力的壓力出口，也即自由開口，模擬空氣邊界；滅火器筒體設(shè)成“wall”，表示固體壁面。松弛因子為：壓力0.3，動量0.7，其它保持默認大小。

穩(wěn)態(tài)模擬收斂后，加入離散相進行瞬態(tài)模擬。離散相以“Injections”加入，選定面噴射源，再選擇整個滅火器器體作為噴射源。滅火劑顆粒初速度為零，平均直徑大小為15μm，總質(zhì)量150g。顆粒追蹤選擇離散隨機游動模型。拽力模型選用“Spherical”，離散相的物理模型選用“Saffman Lift Force”和DEM碰撞模型，這是由于滅火劑顆粒較輕，會受到空氣的浮升力作用，同時顆粒之間也會有碰撞，模擬考慮器性碰撞。此外，將重力項打開，加速度大小設(shè)為9.81m/s2，方向豎直向下。瞬態(tài)模擬的時間步長為0.01s，也就是每隔0.01s對顆粒進行一次追蹤，得出運動軌跡及速度大小等情況。

（5）求解控制壓力速度耦合采用“SIMPLE”算法。空間離散方法：梯度為“Least Squares Cell Based”，壓力采用“Standard”，動量采用二階迎風離散，其它都采用一階迎風離散。

3 模擬結(jié)果及分析

以往對非貯壓手持式滅火器的氣固兩相流研究，局限于檢測條件的不足，主要通過試驗觀測的方法，對滅火劑的射流形態(tài)與滅火器的滅火效能進行分析。通過氣固兩相流的數(shù)值模擬，可以對滅火劑顆粒的運動狀態(tài)有一個更深的認識。

從圖4器筒豎直截面上的氣體速度矢量圖可以看出，動力源氣體進入器筒后，與超細干粉滅火劑顆?；旌希跍缁饎╊w粒間流動，給予了滅火劑顆粒一定的動量。氣固兩相流在中軸線上的運動最為激烈，大體向著噴口方向運動。除此之外，流體在器筒壁面上發(fā)生碰撞，反彈后加劇了顆粒間碰撞運動，在局部形成更無序的狀態(tài)，破壞了兩相流體向噴口處運動的軌跡[8]。由圖5流場中的氣體速度矢量圖可以看出，氣體裹挾了滅火劑粉體從噴口噴出后，在空氣中做擴散運動，當擴散到最大狀態(tài)后，氣體原本具有的動量耗散干凈，開始做分子熱運動。

圖4 筒體豎直截面上氣體速度矢量圖

圖5 流場中的氣體速度矢量圖

圖6 滅火劑顆粒速度流線圖

由圖6滅火劑顆粒的速度流線圖可以看出，滅火劑顆粒在噴射過程中，射流中心的集束性非常強。這樣不利于超細干粉滅火劑的散布，覆蓋不到足夠的面積，使滅火劑在某一小區(qū)域過于集中，降低了其他區(qū)域的滅火劑顆粒濃度，以至于達不到預期的滅火效能。整個模擬過程中滅火劑顆粒的噴射持續(xù)時間為1.07s，與實際噴射持續(xù)時間1.13s非常接近，兩者相差4.4%，說明仿真結(jié)果與試驗很相近，在可接受誤差范圍之內(nèi)。

4 試驗驗證

試驗過程中，在噴口處設(shè)置一塊黑色幕布作為背景，幕布上每隔1m有一條黃色標示帶，用以提高射流分辨程度。將非貯壓手持式滅火器氣固兩相流模擬得到的滅火劑顆粒散布圖與高速攝影機拍攝到的滅火器試驗噴射圖進行比較，如圖7所示。

在非貯壓手持式滅火器啟動后，產(chǎn)氣發(fā)生器反應生成熱氣溶膠氣體并送入器筒內(nèi)，氣體積蓄并產(chǎn)生壓力，當壓力突破噴口薄膜所能承受的壓力值后，氣體裹挾著超細干粉滅火劑噴射而出。在外部大氣阻力和內(nèi)部氣流推動的作用下，滅火劑射流的邊界產(chǎn)生了湍流運動，邊界變得不穩(wěn)定，并形成“羽流現(xiàn)象”。當滅火劑粉體的動量越來越小，直至消耗完畢后，射流前端的滅火劑顆粒主要被空氣攜帶著隨之流動，受氣粉兩相流體的后續(xù)推動以及滅火劑顆粒的分子運動而繼續(xù)向前擴散。根據(jù)模擬射流與試驗射流之間的比較，兩者有近似相同的流體形態(tài)變化。因而，對非貯壓手持式滅火器的氣固兩相流仿真具有可信性，可以作為下一步研究的方法和手段。

5 結(jié)語

本文應用FLUNET軟件實現(xiàn)了非貯壓手持式滅火器的氣固兩相流場模擬，通過氣體速度矢量圖和顆粒速度流線圖，分析了兩相流體在筒體內(nèi)和外的運動狀態(tài)。結(jié)合高速攝影機采集的滅火劑顆粒射流圖像，說明了仿真與試驗的吻合程度較高，用仿真模擬代替試驗研究的方法是可行的，為下一步研究提出新的研究方法和思路。

[1] 吳昌寧,丁宇龍,金涌,程勇.氣固催化反應流動的CFD-DEM模擬研究[J].中國科技論文在線,2008,3（12）:897-904.

[2] 昌澤舟, A.Berlemont, G.Gouesbet.氣固兩相流中顆粒彌散的拉格朗日模擬[J].計算力學學報,2001,19(8):284-289.

[3] 袁竹林,朱立平,耿凡,彭正標.氣固兩相流動與數(shù)值模擬[M].南京:東南大學出版社,2013.

[4] 陳敏,袁竹林,郜時旺.移動床過濾除塵的直接數(shù)值模擬研究及試驗對比[J].中國電機工程學報,2003,26(11):199-203.

[5] 蔡桂英,袁竹林.用離散顆粒數(shù)值模擬對陶瓷過濾器過濾特性的研究[J].中國電機工程學報,2003,26(12):206-210.

[6] 熊源泉,袁竹林,章名耀.加壓條件下氣固噴射器輸送特性的三維數(shù)值模擬[J].化工學報,2004,35(10):1 638-1 643.

[7] 熊源泉.氣固噴射器輸送特性的試驗研究及其三維數(shù)值模擬[D].南京:東南大學,2003.

[8] 李希川.豎直窄矩形流道液泛特性研究[D].哈爾濱:哈爾濱工程大學,2013.

Numerical Simulation of Gas Solid Two Phase Flow in A Non Storage Pressure Hand Held Fire Extinguisher

CAO Shu-hua，ZHUANG Hong-wei，YU Peng-fei

（Engineering University of PAP, Xi’an, 710086）

Based on the study of gas-solid two-phase flow of the non pressure storage portable fire extinguisher, the gas velocity vector and particle velocity streamlines were applied to analyze the motion state of two-phase flow in and out of the cylinder body, by numerical simulation based on Flunet simulation software, as well as the correctness of the simulation results were verified by test. It showed that the results of simulation are accordance with that of test, which provide theoretical guidance for structure optimization of fire extinguisher.

Fire extinguisher；Gas-solid two-phase flow；Simulation；Test

1003-1480（2017）02-0049-04

TQ567.9

A

2017-02-28

曹書華（1992 -），男，在讀碩士研究生，主要從事非致命武器研究。