基于CFD-DEM 耦合的條縫出流自由下落流場特性研究

伍吉富 孫宏發(fā) 張靜

西安建筑科技大學(xué)建筑設(shè)備科學(xué)與工程學(xué)院

在工業(yè)生產(chǎn)過程中，散料在重力作用下自由下落至料堆或傳送帶時(shí)，會散發(fā)大量的粉塵，使得地下環(huán)境受到污染，不僅會影響設(shè)備的運(yùn)行，更會危及工作人員的健康。上世紀(jì)60 年代以來，各國學(xué)者對散料顆粒自由下落過程進(jìn)行了大量的研究，Hemeon[1]基于單顆粒在靜止空氣中下落最先提出了預(yù)測自由下落微粒流的卷吸空氣量模型。Cooper[2]對卷吸空氣量進(jìn)行了更深一步的研究，從理論分析提出了不同顆粒粒徑的卷吸空氣量模型。Ogata[3]通過實(shí)驗(yàn)研究了顆粒流的速度分布特性，發(fā)現(xiàn)顆粒流軸心速度大于單顆粒的沉降速度，并給出了不同顆粒雷諾數(shù)下的卷吸空氣量模型。隨著計(jì)算能力及求解精度的快速提升，采用CFD模擬進(jìn)行顆粒下落的研究已較為普遍[4-6]。現(xiàn)有研究主要針對圓孔出流的點(diǎn)狀塵源，鮮有對條縫出流線狀塵源的研究，且相關(guān)模擬研究大都忽略了顆粒碰撞的影響，文獻(xiàn)[7]指出顆粒體積分?jǐn)?shù)即便只有0.04%，顆粒碰撞的影響依舊不能被忽略。基于計(jì)算流體力學(xué)（CFD）與離散元素法[8]（DEM）的CFD-DEM 耦合模型[9]能更好地考慮顆粒之間的碰撞，因此本文采用CFD-DEM耦合模型對線狀塵源的條縫出流不同出口長度及粒徑的自由下落過程進(jìn)行了研究，從微觀尺度分析微粒羽流運(yùn)動及卷吸空氣的流場特性，以期為通風(fēng)除塵設(shè)計(jì)提供相關(guān)指導(dǎo)。

1 CFD-DEM 耦合模型

近年來，CFD-DEM 耦合模型作為一種全新的仿真方法已經(jīng)得到廣泛應(yīng)用。該模型能夠更好的模擬顆粒與顆粒，顆粒與壁面之間的碰撞，獲得顆粒和流場的微觀信息，保證對顆粒運(yùn)動的精確模擬，被譽(yù)為最有前途的氣固兩相流模擬方法之一[10]。

1.1 流體相控制方程

流體相基于連續(xù)介質(zhì)模型，考慮顆粒相作用的流體運(yùn)動連續(xù)性方程和動量方程分別為：

式中：ρf，uf，P，τ 分別代表顆粒的密度，速度，壓力以及應(yīng)力張量；g 為重力加速度；ε 為空隙度，是網(wǎng)格單元中流體所占的體積分?jǐn)?shù)；S 為動量匯，是作用在網(wǎng)格單元內(nèi)流體阻力的總和。

1.2 顆粒相控制方程

顆粒相的運(yùn)動基于離散元素法，根據(jù)牛頓第二定律求解顆粒運(yùn)動方程：

式中：mi，Ii，up,i，ωp,i分別為顆粒i 的質(zhì)量，轉(zhuǎn)動慣量，平移和旋轉(zhuǎn)速度；Fn,ij，F(xiàn)t,ij是顆粒i，j 之間的法向和切向接觸阻力；Mij是顆粒之間的接觸扭矩；Ff-p，Mf-p分別是流體對顆粒的作用力和作用力矩。文獻(xiàn)[4]指出在自由下落的過程中Saffman 升力對顆粒流的擴(kuò)散有著重要的影響，因此在本文研究中，流體對顆粒的作用力主要考慮阻力和Saffman 升力。

2 模型建立及驗(yàn)證

2.1 物理模型建立

在實(shí)際工程中，條縫出流更為普遍。建立如圖1 所示顆粒條縫出流自由下落的物理模型和三維坐標(biāo)系。計(jì)算區(qū)域?yàn)橐粋€長方體，模型尺寸為600 mm×600 mm×1000 mm（長×寬×高）。在模型頂部設(shè)置顆粒流出口，高為5 mm，寬度b 為2 mm，長度l 隨模擬工況變化。采用ICEM 軟件對模型進(jìn)行計(jì)結(jié)構(gòu)化網(wǎng)格劃分，并對顆粒流出口及下落部分進(jìn)行局部網(wǎng)格加密。

圖1 顆粒條縫出流自由下落物理模型

2.2 邊界條件及模擬工況設(shè)置

對流體相，顆粒流出口設(shè)為速度入口邊界（Velocity inlet），取值0.001 m/s[6]，出流口側(cè)面以及底面為墻邊界（Wall），計(jì)算區(qū)域其他部分設(shè)置為自由出流邊界（Outflow）。對顆粒相，設(shè)置顆粒初速為0 m/s，顆粒的材料力學(xué)性質(zhì)詳見表1[3，11-13]，接觸參數(shù)詳見表2[14-15]。

表1 材料力學(xué)性質(zhì)

表2 接觸參數(shù)

為研究不同出口尺寸和粒徑對顆粒流場特性的影響，進(jìn)行了多工況模擬，不同工況質(zhì)量流量的計(jì)算公式[16]如下：

式中：mp是顆粒流的質(zhì)量流量；ρb是顆粒的堆積密度；A0是孔口面積；b 是顆粒出流寬度；dp則是顆粒的平均粒徑。

選取工業(yè)散料氧化鋁顆粒進(jìn)行分析，設(shè)置不同的條縫出流口尺寸，并參照文獻(xiàn)[5]設(shè)置顆粒粒徑分別為300 μm，400 μm，500 μm，具體模擬工況見表3。

表3 模擬計(jì)算工況

2.3 湍流模型選擇及驗(yàn)證

前人對圓孔出流自由下落過程進(jìn)行了許多研究，為了驗(yàn)證CFD-DEM 耦合模型的在自由下落研究的可行性，對圓孔出流過程進(jìn)行驗(yàn)證。

在顆粒流下落研究中常用的湍流模型[17]有：標(biāo)準(zhǔn)k-ε 模型（SKE）、重整化k-ε 模型（RNG）、可實(shí)現(xiàn)的k-ε 模型（RKE），雷諾應(yīng)力模型（RSM）以及剪切應(yīng)力輸運(yùn)模型(SST)。Ogata 等人通過實(shí)驗(yàn)研究了不同孔徑自由下落顆粒流的速度分布[3-4]，參照實(shí)驗(yàn)設(shè)置出流圓孔直徑D 為4 mm，采用平均粒徑為454 μm（幾何標(biāo)準(zhǔn)差為0.116）玻璃微珠進(jìn)行模擬，玻璃微珠的初速度為0.77 m/s，堆積密度ρb為1460 kg/m3，質(zhì)量流量為0.002 kg/s。為選擇合適的湍流模型，采用這五種湍流模型與Ogata 的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對比，如圖2（a）所示。從圖中可以看出，不同湍流模型與實(shí)驗(yàn)值的變化規(guī)律基本一致，其平均預(yù)測偏差分別為6.7%、5.7%、5.8%、4.6%、7.2%，RSM 雷諾應(yīng)力湍流模型的偏差最小，因此后續(xù)數(shù)值模擬選用RSM 湍流模型進(jìn)行分析。

針對圓孔出流自由下落過程Hemeon[1]，Cooper[2]及Ogata[3，5]等先后提出了不同的卷吸空氣量模型。采用氧化鋁顆粒進(jìn)行分析，粒徑取為300 μm，質(zhì)量流量取為0.0015 kg/s[18]，以對流體相卷吸空氣進(jìn)行驗(yàn)證。在某一高度截面上，對顆粒流卷吸直徑內(nèi)的氣流速度進(jìn)行積分計(jì)算模擬所得的卷吸空氣量，并將其與理論計(jì)算值進(jìn)行對比，詳見圖2（b）。可以看出，Ogata 提出的卷吸空氣模型計(jì)算值相比于Hemeon 和Cooper 的計(jì)算值要大，這可能是因?yàn)镠emeon 低估了顆粒群中單顆粒的速度，模擬值與Ogata 提出的卷吸空氣模型擬合的較好，為CFD-DEM 耦合模型對自由下落流體相的模擬提供了理論支持。

圖2 模型驗(yàn)證

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

3.1 不同下落高度流場特性分析

圖3 和圖4 分別給出了工況1 不同下落高度顆粒流的速度分布圖。圖3 所示，由于顆粒流核心區(qū)大，在z=100 mm 處Y 方向核心區(qū)顆粒速度大小基本一致，隨著下落高度的增加，顆粒流逸散寬度增大且短邊方向增長更為明顯，速度也逐漸增大，但增長趨緩。相同高度不同方向上顆粒的速度分布見圖4，在z=500 mm處不同方向上的顆粒速度分布差異很大，而在z=900 mm 處不同方向上速度分布趨于一致，這和三維射流在軸對稱衰減區(qū)空氣速度剖面相似的結(jié)論相仿[19]。

圖3 不同高度顆粒速度

圖4 不同方向顆粒速度

圖5 和圖6 分別給出了工況1 不同下落高度顆粒流卷吸空氣速度分布圖。由圖5 可知，卷吸空氣寬度遠(yuǎn)大于顆粒流逸散寬度，且隨著下落高度增加而增大，但卷吸空氣軸心速度先增大后減小。這是因?yàn)樵谙侣涑跏茧A段，顆粒流向兩側(cè)的逸散寬度相對較小，顆粒對單位流體的動量匯隨著顆粒速度的增大而增大，而隨著顆粒逸散寬度的增大，顆粒速度增長趨緩，作用在單位流體上的動量匯S 逐漸減少，因此卷吸空氣速度減小。圖6 所示，在遠(yuǎn)離出口的z=500 mm 和z=900 mm 處不同方向上卷吸空氣速度分布基本趨于一致。

圖5 不同高度卷吸空氣速度分布

圖6 不同方向卷吸空氣速度分布

3.2 不同出口長度流場特性分析

對于不同出口長度20 mm，30 mm，40 mm，圖7給出了核心區(qū)顆粒流和卷吸空氣速度的對比。考慮到不同方向顆粒和卷吸空氣軸心速度一致，后續(xù)對比只從X 方向進(jìn)行分析。由圖7 可知，出口長度越大顆粒和卷吸空氣速度越大，并在l＞30 mm 后趨于定值。其原因是，出口長度越大顆粒流核心區(qū)越大，更不易向外擴(kuò)散，受空氣阻力的影響越小，但隨著出口長度不斷增大，顆粒流核心區(qū)受空氣影響不明顯因而顆粒和卷吸空氣速度逐漸趨于定值。

圖7 不同出口長度顆粒和卷吸空氣速度對比

參考Wypych 等人[20]提出的方法，結(jié)合單位質(zhì)量流量的卷吸空氣量進(jìn)行分析。由圖8 可知卷吸空氣量隨出口長度增加而增大，但單位質(zhì)量流量卷吸空氣量隨出口長度增加而減小。這可能是因?yàn)轭w粒流核心區(qū)越大，孔隙度越小，顆粒與空氣接觸面積越小，所以單位質(zhì)量流量卷吸空氣量越少。

圖8 不同出流長度卷吸空氣量對比

3.3 不同顆粒直徑流場特性分析

針對粒徑分別為300 μm，400 μm，500 μm 的氧化鋁顆粒，核心區(qū)顆粒和卷吸空氣速度的對比分析見圖9。由圖9 可知，顆粒粒徑越大顆粒速度越大，卷吸空氣速度越小。這是由于顆粒粒徑越大質(zhì)量流量越小，顆粒數(shù)越少，所受空氣阻力越小，因此顆粒速度越大，同時(shí)流體內(nèi)顆粒數(shù)越少，卷吸空氣速度越小。

不同粒徑卷吸空氣量對比如圖10 所示，可以看出相同高度處粒徑越小顆粒的卷吸空氣量越大，且單位質(zhì)量流量卷吸空氣量也隨粒徑減小而增大。這是因?yàn)樵趩挝毁|(zhì)量流量下，粒徑越小顆粒數(shù)越多，顆粒傳遞給單位體積流體的動量匯S 越多，所以單位質(zhì)量流量卷吸空氣量也越大。

圖9 不同粒徑顆粒和卷吸空氣速度對比

圖10 不同粒徑卷吸空氣量對比

4 結(jié)論

本文采用基于計(jì)算流體力學(xué)和離散元素法的CFD-DEM 耦合方法，研究了條縫出流不同下落高度，出口長度及粒徑對流場的影響，得到了以下結(jié)論：

1）條縫出流的逸散寬度在短邊方向增長速度相比長邊更快。在z=500 mm 處不同方向上顆粒速度分布差異較大，而在z=900 mm 處不同方向的速度分布基本趨于一致。遠(yuǎn)離出口處相同高度不同方向上的卷吸空氣速度分布趨于一致。

2）隨著出口長度度的增加，顆粒和卷吸空氣的速度逐漸增大并在l＞30 mm 后趨于定值。粒徑越大，顆粒速度越大，而卷吸空氣速度越小。

3）對于實(shí)際工程設(shè)計(jì)，卷吸空氣量隨出口長度的增大而增大，但單位質(zhì)量流量的卷吸空氣量隨之減小。粒徑越小，卷吸空氣量和單位質(zhì)量流量的卷吸空氣量越大。