白建平，范健強，王　越，李向富，王　雨

(重慶科技學院 安全工程學院，重慶 401331)

0　引言

近年來，國內(nèi)外粉塵爆炸事故屢屢頻發(fā)，造成了巨大的人員傷亡和財產(chǎn)損失，給相關(guān)行業(yè)的安全生產(chǎn)和社會安穩(wěn)大局帶來了重大影響。據(jù)統(tǒng)計，2005—2015年期間我國大陸發(fā)生72起粉塵爆炸事故[1]?；谶@種形勢，粉塵爆炸的研究已成為熱點，據(jù)統(tǒng)計1986—2016年期間國家自然基金委資助162項粉塵災(zāi)害領(lǐng)域的項目，資助金額高達6 128萬，其中，2007—2016年 10 a時間資助項目數(shù)和資助金額分別占據(jù)70%和86.4%[2]。

粉塵最大爆炸壓力、最大爆炸壓力上升速率、爆炸極限濃度等爆炸特性參數(shù)常用來用來評估粉塵爆炸危險性，但是爆炸特性參數(shù)不是粉塵固有屬性，它不僅取決于粉塵粒徑、濃度、濕度、揮發(fā)性等物理性質(zhì)[3-6]，還與實驗設(shè)備密切相關(guān)[7-8]。20 L球形爆炸裝置是目前用于測試粉塵爆炸特性參數(shù)最廣泛的設(shè)備之一。粉塵顆粒在20 L球形爆炸裝置球罐(以下簡稱“球罐”)內(nèi)分散的過程受到重力、曳力、顆粒之間及顆粒與壁面之間的碰撞等力的綜合作用，使得粉塵在球罐內(nèi)的分布狀態(tài)處于不斷變化中。粉塵在分散過程中存在一個最佳的分布狀態(tài)，此時點火，熱量可以傳遞到較多粉塵顆粒上，粉塵體系的活化分子的密度將達到最大，第一時間參加爆炸反應(yīng)的粉塵將會更多，粉塵體系釋放出熱量也會更多更迅速，實驗結(jié)果也更能客觀真實全面反映粉塵爆炸危險性，因此，合理的點火延遲時間對客觀真實全面評估粉塵爆炸危險性具有重要意義。

近些年來，針對點火延遲時間對粉塵爆炸特性的影響規(guī)律，諸多學者[9-11]進行了大量的實驗研究，結(jié)果表明，在粉塵粒徑和濃度不同時，最大爆炸壓力和最大爆炸壓力上升速率對應(yīng)的最佳點火時間也不一樣。針對噴粉后球罐內(nèi)流場的演化規(guī)律，陳嘉琛[12]研究了粒徑對鋁粉擴散的影響規(guī)律，結(jié)果表明，點火延遲時間段隨粉塵顆粒粒度的增大而減??；Benedetto[13]和Murillo[14]的研究結(jié)果表明，由于粉塵顆粒的作用，3D模型球罐內(nèi)的各流場參量并不嚴格對稱征，而2D模型則無法反映這一特征；Sarli[15-17]研究了粉塵粒徑、濃度及噴嘴形狀對球罐內(nèi)流場的影響，結(jié)果表明，隨著粉塵粒徑的增大，氣相和顆粒相的運動軌跡發(fā)生剝離，粉塵濃度越大沉降越明顯，安裝穿孔環(huán)形噴嘴的球罐湍流動能分布相當均勻，但其值顯著低于回彈噴嘴。但目前為止，尚未有研究粉塵密度對球罐內(nèi)流場參數(shù)的影響規(guī)律。不同密度的粉塵顆粒，因單位體積所受的重力不同，粉塵沉降和分散、流場速度及湍流動能會呈現(xiàn)不同的規(guī)律，筆者通過CFD軟件對不同密度粉塵顆粒在罐內(nèi)分散規(guī)律進行研究，揭示粉塵密度對球罐內(nèi)流場參量和點火延遲時間的影響規(guī)律，為后期研究20 L球型爆炸裝置的粉塵爆炸提供指導。

1　 數(shù)值模型建立

1.1　幾何模型

20 L球形爆炸裝置的噴粉過程：儲粉罐(體積為0.6 L)里的粉塵在壓力2 MPa的壓縮空氣的作用下，噴入初壓為-0.06 MPa球罐(體積為2 L)內(nèi)，在回彈噴嘴和罐壁的反彈作用下，形成均勻彌散的粉塵云。為客觀模擬粉塵的分散情況，本文建立的3D模型最大限度地保留了20 L球形爆炸裝置的主要特征，如圖1所示。并對回彈噴嘴處網(wǎng)格進行加密，網(wǎng)格總數(shù)為986 865，最小網(wǎng)格體積為1.2×10-10m3，球罐壁面及噴嘴的邊界條件均設(shè)為wall。

圖1　20 L球形爆炸裝置幾何模型Fig.1　Geometric model of 20 L spherical explosive device

1.2　模型假設(shè)與控制方程

罐內(nèi)粉塵分散過程為非定常氣固兩相流問題，鑒于粉塵動力源為壓縮空氣，故假設(shè)空氣為理想氣體，湍流流動模型采用標準 k-ε模型。氣相流動采用雷諾平均N-S(RANS)方程進行求解，雷諾平均N-S控制方程：

連續(xù)方程：

(1)

動量守恒方程：

(2)

顆粒選取3種常見爆炸性粉塵，假設(shè)粉塵尺寸為圓形顆粒，粒徑為50 μm，球罐內(nèi)粉塵表觀濃度為250 g/m3，粉塵密度和體積分數(shù)如表1所示。由于顆粒體積分數(shù)小于10%，故選用離散相模型。

表1　粉塵密度及體積分數(shù)

顆粒分散過程中考慮重力、壁面的反射作用及顆粒相和氣相之間的相互作用，忽略虛擬質(zhì)量力、熱泳力、布朗力、薩夫曼力等的影響。顆粒所受作用力平衡方程為：

(3)

式中：F為顆粒所受重力，N；up和u分別為顆粒和氣體速度，m·s-1；ρp和ρ分別為顆粒和氣體密度，kg·m-3；FD(u-up)為顆粒單位質(zhì)量的曳力，N。可由式(4)得到：

(4)

式中：μ為流體動力粘度，Pa·s；dp為顆粒直徑，m；Re為顆粒雷諾數(shù)，可由式(5)得到；CD為拖曳力系數(shù)，可由式(6)得到：

(5)

(6)

1.3　參數(shù)設(shè)置

物性參數(shù)的設(shè)置如前文所述，初始時刻，儲粉罐和連接管初始壓力為2 MPa，球罐(初始壓力為-0.06 MPa。采用SIMPLE算法，迭代時間步長為0.000 1 s，迭代時間步數(shù)為2 000步，每隔1個時間步長迭代20次，計算時間為200 ms。

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

2.1　模型驗證

對無塵空氣數(shù)值模擬的結(jié)果如圖2和3所示，40 ms之后儲粉罐和球罐壓力均相對穩(wěn)定，震蕩幅度越來越小，儲粉罐和球罐壓力變化曲線與文獻[13]的結(jié)果基本一致。經(jīng)測量，本文的幾何模型高壓區(qū)域(圖1儲粉罐部分)的實際體積為0.601 L，低壓區(qū)域(圖1球罐部分)的體積為20.019 L，根據(jù)理想氣體方程計算2個區(qū)域最終的壓力為1 411.9 Pa。40 ms之后球罐的壓力與理論值之間的相對誤差不斷減小，圖3為132 ms時計算域的壓力云圖，該時刻相對誤差為10.29%?？梢灶A見，最終儲粉罐和球罐壓力會無限逼近理論值。因此，本文建立的數(shù)值模型是可靠的。

圖2　儲粉罐和球罐內(nèi)的壓力變化趨勢Fig.2　Variation trend of pressure in the storage tank and the ball tank

圖3　儲粉罐和球罐內(nèi)的壓力分布(132 ms)Fig.3　Pressure distribution in the storage tank and the ball tank ( 132 ms)

2.2　密度對球罐內(nèi)流場參量的影響規(guī)律

本文以鋁粉、鋯粉、鋅粉為對象，對球罐球心(點火位置)處的湍流動能、流場速度、粉塵濃度的變化趨勢進行實時監(jiān)測，結(jié)果如圖4～6所示。

圖4　球罐內(nèi)球心處湍流動能的變化規(guī)律Fig.4　Variation of turbulent kinetic energy at the center of the spherical tank

由圖4可知，球罐內(nèi)的湍流動能在噴粉期間迅速達到峰值，3種粉塵對應(yīng)的峰值分布為3 324.06，3 968.88，4 038.57 m2/s2，這說明粉塵密度越小，峰值越小，但密度相差不大的粉塵，湍流動能的變化曲線和峰值差異很小。而且3種粉塵對應(yīng)的湍流動能在14.5 ms同時上升至峰值，這說明粉塵密度對湍流動能的增值速率沒有影響。由于噴粉壓力的減小，14.5 ms湍流動能開始急劇衰減，衰減速率逐漸變緩，在110 ms之后逐漸趨于穩(wěn)定。

圖5　球罐內(nèi)球心處流場速度的變化規(guī)律Fig.5　Variation of flow velocity at the center of the spherical tank

由圖5可知，球罐內(nèi)的流場速度在噴粉期間迅速達到峰值，3種粉塵對應(yīng)的峰值分布為102.01，92.91，92.44 m/s，這說明粉塵密度越小，流場速度峰值越大，且抵達峰值的時間越短；密度接近的粉塵，流場速度峰值和變化曲線很接近。隨后因噴粉壓力的減小，流場速度開始急劇減小，減小速率逐漸也變小，在100 ms之后也逐漸趨于穩(wěn)定。

圖6　球罐內(nèi)球心處粉塵云濃度的變化規(guī)律Fig.6　Variation of dust cloud concentration at the center of the spherical tank

由圖6可知，球罐內(nèi)的粉塵云濃度變化曲線呈“鋸齒”狀，這是由于粉塵顆粒反復被罐壁反彈和氣流的沖擊。粉塵濃度在噴粉結(jié)束后才達到峰值，然后又急劇下降。3種粉塵對應(yīng)的粉塵濃度的峰值分布為2.36，1.81，1.72 g/m2，這說明粉塵密度越小，粉塵濃度增值速率越快，濃度峰值越大，也越先降低至穩(wěn)定值。

2.3　球罐內(nèi)流場參量分布及演化規(guī)律

以鋁粉為例，分析球罐內(nèi)湍流動能、流場速度、粉塵云濃度及顆粒軌跡4種流場參量的演化規(guī)律和內(nèi)在聯(lián)系。截取10，60，100，200 ms 4個具有代表性的時間段的流場變量云圖，如圖7所示。

圖7　球罐內(nèi)流場參量的分布及變化過程Fig.7　Distribution and variation of flow field parameters in spherical tank

由圖7可知，在噴粉階段，湍流動能和速度的最大值集中在噴嘴處；噴粉結(jié)束，峰值逐漸向上移動。噴粉初始階段，球罐內(nèi)湍流動能、流場速度、粉塵云濃度和顆粒軌跡云圖呈現(xiàn)對稱結(jié)構(gòu)。流場存在2個對稱旋渦，而隨后這2個旋渦順時針轉(zhuǎn)動，并逐漸合二為一，如圖7(b)，流場速度和湍流動能的隨即也逐漸變?yōu)椴粚ΨQ，如圖7(a)和(b)，同時粉塵顆粒在旋渦的推擠下向罐壁堆積，如圖7(c)和(d)，球心處的粉塵濃度急劇下降，如圖6，而罐壁處的粉塵濃度遠大于標稱濃度(0.25 kg/m3)。旋渦變化的原因，一方面，粉塵顆粒因重力和噴粉速度的減小引發(fā)沉降現(xiàn)象；另一方面，球罐內(nèi)粉塵的增多，致使密度梯度增大。因此，用20 L球形爆炸裝置測試粉塵爆炸特性時，應(yīng)該在粉塵尚未大量堆積在罐壁、流場速度和湍流動能足夠穩(wěn)定時點燃粉塵。粉塵密度不同，球心處的湍流濃度、流場速度、粉塵濃度3個變量的峰值也不同，粉塵濃度衰減至穩(wěn)定值的時間也不同。結(jié)合圖4～7分析得知，對于鋁粉，點火延遲時間建議在50～60 ms之間，而對于鋯粉和鋅粉，點火延遲時間建議在60～80 ms之間。

3　結(jié)論

1)在其他條件一致的情況下，粉塵密度越小，球罐內(nèi)的湍流動能的峰值越小，而粉塵云濃度和流場速度的峰值則越大，但密度相近的粉塵，峰值相差不大。

2)在其他條件一致的情況下，粉塵密度對湍流動能的增值速率沒有影響，而粉塵密度越小，流場速度和粉塵濃度增值速率越快，粉塵濃度衰減至穩(wěn)定值的時間也越短。

3)根據(jù)本文研究結(jié)果，鋁粉最佳點火延遲時間在50～60 ms之間，鋯粉和鋅粉在60～80 ms之間，這說明粉塵密度越小，點火延遲時間越小。

4)噴粉前期，球罐內(nèi)流場是對稱分布的，而隨著球罐內(nèi)粉塵顆粒沉降和密度梯度增大，逐漸演化為非對稱分布。

[1]鐘圣俊. 我國粉塵爆炸防護問題與建議[J]. 勞動保護,2016(9):14-16.

ZHONG Shengjun. Problems and suggestions on dust explosion protection in China[J]. Labour Protection,2016(9):14-16.

[2]韓方偉,李世虎,王天暘. 粉塵災(zāi)害領(lǐng)域國家自然科學基金資助情況分析[J]. 中國安全生產(chǎn)科學技術(shù),2016,12(12):166-170.

HAN Fangwei, LI Shihu, WANG Tianyang. Investigation and analysis on projects in the field of dust disaster supported by National Natural Science Foundation of China[J]. Journal of Safety Science and Technology, 2016, 12(12):166-170.

[3]王林元,呂瑞琪,鄧洪波. 不同粒徑鎂鋁合金粉塵爆炸與抑爆特性研究[J]. 中國安全生產(chǎn)科學技術(shù),2017, 13(1):34-38.

WANG Linyuan, LYU Ruiqi, DENG Hongbo. Study on characteristics of explosion and explosion suppression for magnesium aluminum alloy dust with different particle size[J]. Journal of Safety Science and Technology, 2017, 13(1):34-38

[4]馬小明,周可. 不飽和聚酯樹脂鈕扣粉塵靜電爆炸敏感性研究[J]. 中國安全生產(chǎn)科學技術(shù), 2016, 12(3):149-154.

MA Xiaoming, ZHOU Ke. Study on sensitivity of electrostatic explosion for dust of unsaturated polyester resin button[J]. Journal of Safety Science and Technology, 2016, 12(3):149-154.

[5]代濠源,樊建春,孫莉,等. 初始溫度對濕法成型硫磺燃燒爆炸特性影響的試驗研究[J].中國安全生產(chǎn)科學技術(shù), 2015, 11(3):24-28.

DAI Haoyuan, FAN Jianchun, SUN Li, et al. Experimental study on influence by initial temperature to combustion and explosion characteristics of wet -granulation sulfur [J].Journal of Safety Science and Technology, 2015, 11(3):24-28.

[6]蒯念生,黃衛(wèi)星,袁旌杰,等. 點火能量對粉塵爆炸行為的影響[J]. 爆炸與沖擊,2012,32(4):432-438.

KUAI Niansheng, HUANG Weixing, YUAN Jingjie, et al. Influence of ignition energy on dust explosion behavior[J]. Explosion and Shock Waves, 2012,32(4):432-438.

[7]丁小勇,譚迎新,秦澗,等. 垂直哈特曼管與水平管道中鋁粉爆炸特性[J].消防科學與技術(shù), 2012,31(6):561-563.

DING Xiaoyong, TAN Yingxin, QIN Jian, et al. Study on alumimun dust explosion characteristics in vertical Hartmann tube and horizontal pipeline[J]. Fire Science and Technology, 2012,31(6):561-563.

[8]苑春苗,王麗茸,陳寶智,等. 1.2 L Harttman管式與20 L球型爆炸測試裝置爆炸猛度實驗研究[J]. 中國安全生產(chǎn)科學技術(shù),2008,4(1):108-111.

YUAN Chunmiao, WANG Lirong, CHEN Baozhi, et al. Experimental study on the violence of dust explosion with the 20 litre sphere and with the 1.2 L Harttman tube [J]. Journal of Safety Science and Technology,2008,4(1):108-111.

[9]陳春燕,龍思華,肖國清,等. 點火延遲時間對甘薯粉塵爆炸的影響研究[J]. 中國安全生產(chǎn)科學技術(shù),2016,12(12):171-175.

CHEN Chunyan, LONG Sihua, XIAO Guoqing, et al. Study on influence of ignition delay time on explosion of sweet potato dust[J]. Journal of Safety Science and Technology, 2016,12(12):171-175.

[10]衛(wèi)少潔,譚迎新. 不同濃度下點火延長時間對煤粉點火能的影響[J]. 中國粉體技術(shù),2016,12(6):32-34.

WEI Shaojie, TAN Yingxin. Effect of ignition delay time on pulverized coal ignition energy at different concentrations[J]. China Powder Science and Technology,2016,12(6):32-34.

[11]張博,張奇,譚汝媚. 噴粉壓力及點火延遲時間對粉塵爆炸參數(shù)的影響[J]. 高壓物理學報,2014,28(2):183-190.

ZHANG Bo, ZHANU Qi, TAN Rumei. Influence of dispersion pressure and Ignition delay time on the dust explosion parameters[J]. Chinese Journal of High Pressure Physics, 2014,28(2):183-190.

[12]陳嘉琛,張奇,馬秋菊,等. 20 L球型罐內(nèi)不同粒徑鋁粉擴散的數(shù)值模擬[J]. 高壓物理學報,2014,28(2):202-208.

CHEN Jiachen, ZHANG Qi, MA Qiuju, et al. Numerical simulation of aluminite dust dispersion for different particle sizes in 20 L vessel[J]. Chinese Journal of High Pressure Physics, 2014,28(2):202-208.

[13]Benedetto A D, Russo P, Sanchirico R, et al. CFD simulations of turbulent fluid flow and dust dispersion in the 20 liter explosion vessel[J]. Aiche Journal, 2013, 59(7):2485-2496.

[14]Murillo C., Bardin-monnier N., Blanchard C., et al. CFD to improve the repeatability and accuracy of dust explosion tests in the 20-liters sphere[J]. Chemical Engineering Transactions, 2016(48)：115-120.

[15]Sarli V D, Sanchirico R, Russo P, et al. CFD modeling and simulation of turbulent fluid flow and dust dispersion in the 20 L explosion vessel equipped with the perforated annular nozzle[J]. Journal of Loss Prevention in the Process Industries, 2015(38):204-213.

[16]Sarli V D, Russo P, Sanchirico R, et al. CFD simulations of the effect of dust diameter on the dispersion in the 20 L bomb[J]. Chemical Engineering Transactions, 2013(31):727-732.

[17]Sarli V D, Russo P, Sanchirico R, et al. CFD simulations of dust dispersion in the 20 L vessel: effect of nominal dust concentration[J]. Journal of Loss Prevention in the Process Industries, 2014, 27(1):8-12.