發(fā)送罐內(nèi)煤粉流動(dòng)特性和機(jī)理的數(shù)值研究

王燕豐 范贏(yíng) 孫姣

摘要 為研究氣力輸送發(fā)送罐中煤粉出料穩(wěn)定性及兩相流動(dòng)機(jī)理，采用基于Realizable k-ε雙方程湍流模型并結(jié)合顆粒動(dòng)力學(xué)理論的Two Fluid Model（TFM）對(duì)其進(jìn)行數(shù)值模擬。根據(jù)實(shí)驗(yàn)數(shù)據(jù)對(duì)Gidaspow模型中的曳力系數(shù)進(jìn)行修訂，得到自定義的曳力模型，使之符合實(shí)際工況。對(duì)煤粉在出料口的質(zhì)量流量、體積分?jǐn)?shù)的波動(dòng)特性進(jìn)行分析，發(fā)現(xiàn)隨著入口氣速增大，煤粉出料穩(wěn)定性呈現(xiàn)出先增加后降低的規(guī)律。最后對(duì)發(fā)送罐中空氣-煤粉兩相耦合運(yùn)動(dòng)流動(dòng)機(jī)理進(jìn)行分析，發(fā)現(xiàn)罐中空氣主要以氣流渦的形式進(jìn)行運(yùn)動(dòng)，具體表現(xiàn)為多種尺度、相互碰撞、疊加的渦驅(qū)動(dòng)煤粉沿流向擴(kuò)散和展向遷移。

關(guān) 鍵 詞 氣力輸送；TFM；曳力系數(shù)；數(shù)值模擬；流動(dòng)機(jī)理；渦流

中圖分類(lèi)號(hào) TQ022.4 文獻(xiàn)標(biāo)志碼 A

0 引言

氣力輸送煤粉廣泛應(yīng)用于電力、冶金和環(huán)保等行業(yè)中，而發(fā)送罐是實(shí)現(xiàn)此類(lèi)粉體輸送的核心裝置之一[1]。氣體輸送煤粉涉及到復(fù)雜的氣-固耦合流動(dòng)，國(guó)內(nèi)外關(guān)于煤粉氣力輸送發(fā)送罐的研究針對(duì)不同操作條件下的輸送規(guī)律較多[2]，對(duì)發(fā)送罐內(nèi)氣固兩相流動(dòng)機(jī)理進(jìn)行分析還比較少，尚未有較為全面深入的研究。崔越等[3]對(duì)栓塞式氣力輸送灰罐進(jìn)行了數(shù)值模擬，得到了物料發(fā)送過(guò)程中可分為起始、穩(wěn)定和下降3個(gè)階段，發(fā)現(xiàn)穩(wěn)定輸送階段基本不受輸送壓力的影響，但是發(fā)送罐出料的穩(wěn)定性會(huì)直接影響氣力輸送過(guò)程中管路內(nèi)的流型[4-5]。徐貴玲，陳曉平等[6，7]在不同壓力和物性條件的上出料式發(fā)送罐氣力輸送試驗(yàn)臺(tái)上，以氮?dú)鉃檩斔徒橘|(zhì)，研究煤粉氣力輸送特性，與高壓輸送相比，常壓輸送的固氣質(zhì)量比較高，氣耗率及輸送能耗較低。賀春輝，沈湘林等[8]以?xún)贾轃熋好悍圩鬏斔臀锪?，發(fā)現(xiàn)上引式出料氣力輸送系統(tǒng)的煤粉質(zhì)量流量、表觀(guān)氣速和固氣質(zhì)量比均略高于下引式出料系統(tǒng)，但二者均具有良好的輸送穩(wěn)定性。可見(jiàn)，不同工況下的研究會(huì)得出不同的結(jié)果，研究成果存在局限性。因此，對(duì)于罐內(nèi)兩相流動(dòng)的機(jī)理研究顯得尤為重要。

基于歐拉-歐拉框架的TFM模型已經(jīng)廣泛的應(yīng)用于流化床[9]、噴動(dòng)床[10]的數(shù)值模擬，其在模擬氣固兩相流的可行性也已經(jīng)被驗(yàn)證[11-12]。顆粒動(dòng)力學(xué)理論（KTGF）是基于稠密的氣體分子動(dòng)力學(xué)理論而建立的描述顆粒相運(yùn)動(dòng)的方程。該方程引入了顆粒溫度和徑向分布函數(shù)來(lái)描述顆粒相的速度脈動(dòng)與顆粒碰撞等。該理論中[13]，顆粒相被處理成擬流體，將顆粒運(yùn)動(dòng)類(lèi)比于氣體分子的熱運(yùn)動(dòng)，顆粒的隨機(jī)運(yùn)動(dòng)動(dòng)能被定義為“顆粒擬溫度”，其正比于顆粒速度脈動(dòng)的均方。本文主要分析發(fā)送罐內(nèi)氣固兩相的流動(dòng)特性，屬于介觀(guān)層次[14]的考察，選用TFM模型和KTGF相耦合的方法可滿(mǎn)足模擬需求。

1 模型描述

1.1 TFM

在歐拉框架下構(gòu)建的雙流體模型屬于比較復(fù)雜的多相流模型。它把每一相看做具有不同運(yùn)動(dòng)特性的連續(xù)性介質(zhì)，而且各相之間存在相互滲透、耦合的現(xiàn)象，從而建立了一套完整的包含有各相的連續(xù)性方程與動(dòng)量方程的方程組進(jìn)行求解。方程把各界面間的交換系數(shù)與壓力項(xiàng)耦合到一起，而耦合的方式則根據(jù)所含相的情況進(jìn)行相應(yīng)的選擇，譬如顆粒流與非顆粒流的處理方式就有所不同。針對(duì)顆粒流，可以采用顆粒流體力學(xué)中的分子運(yùn)動(dòng)理論來(lái)得到流動(dòng)特性。跟單流體模型比較而言，TFM優(yōu)勢(shì)在于考慮了氣固兩相間由于存在相互滑移而引起的阻力以及固相的湍流輸運(yùn)。因此在大多數(shù)情形下，其計(jì)算結(jié)果與實(shí)際情況更加吻合。由于每一相都有其單獨(dú)的守恒方程，各相可以進(jìn)行單獨(dú)的計(jì)算，因此該模型適應(yīng)性很強(qiáng)。當(dāng)然，該方法也有其局限性。顯然，方程數(shù)的增加必然導(dǎo)致計(jì)算量的增加。此外，顆粒相的導(dǎo)熱系數(shù)，擴(kuò)散系數(shù)等物性參數(shù)的確定依然存在一些問(wèn)題，雙流體模型在計(jì)算過(guò)程中也可能會(huì)發(fā)生偽擴(kuò)散現(xiàn)象。

連續(xù)性方程（質(zhì)量守恒方程）

[??tαqρq+?·αqρqυq=p=1nmpq]， （1）

式中：[αq]為q相的體積分?jǐn)?shù)；[ρq]為q相的密度； [υq]為q相的速度； [mpq]為從第p相到q相的質(zhì)量傳遞。

動(dòng)量守恒方程

[??tαqρqυq+?·αqρqυqυq=-αq?p+?·τq+p=1nRpq+mpqυpq+αqρqFq+Flift，q+FVm，q]， （2）

式中：[τq]是第q相的壓力應(yīng)變張量；[Fq]是外部體積力；[Flift，q]是升力；[FVm，q]是虛擬質(zhì)量力；[Rpq]是相間作用力；p為所有相共享的壓力；[υpq]為相間速度。

1.2 顆粒動(dòng)力學(xué)模型

引入顆粒流體力學(xué)來(lái)處理發(fā)送罐中顆粒相的相關(guān)問(wèn)題。顆粒動(dòng)力學(xué)理論源自于氣體分子運(yùn)動(dòng)理論，借助速度分布函數(shù)的Boltzmann方程來(lái)推導(dǎo)顆粒相的物性參數(shù)及守恒方程，與熱力學(xué)溫度表征微觀(guān)尺度分子的能量波動(dòng)類(lèi)似，顆粒相擬溫度表達(dá)了宏觀(guān)尺度顆粒隨機(jī)運(yùn)動(dòng)的動(dòng)能，通過(guò)顆粒相擬溫度Θ來(lái)建立顆粒相宏觀(guān)的黏度、壓力、擴(kuò)散系數(shù)等參數(shù)的本構(gòu)關(guān)系：

[τs=[-ps+αsξs??us]I-2αsμsSs]， （3）

[Ss=12?us+?usT-13??usI]， （4）

[μs=10ρsdsΘπ0.596αs1+eg01+45g0αs1+e2+45αsρsdsg01+eΘπ0.5]， （5）

[ξs=43α2sρsdsg01+eΘπ0.5]， （6）

[ps=αsρsΘ+2ρs（1+e）α2sg0Θ]， （7）

式中：[ξs]是固相體積黏度，源于顆粒之間碰撞造成的動(dòng)量交換；[μs]是顆粒相剪切黏度，由動(dòng)力部分和碰撞部分組成；[ps]是顆粒相壓力。

顆粒相擬溫度[Θ]和顆粒隨機(jī)運(yùn)動(dòng)的動(dòng)能成正比，顆粒相脈動(dòng)能守恒方程

[32??tαsρsΘ+??αsρsusΘ=τs：?us+??q-γ-3βΘ]， （8）

其中：

[q=k?Q ，k=2rsa2sds1+eg0Qp0.5+150rsdsQp0.53841+eg01+65asg01+e2 ，g=31-e2as2rsg0Q4dsQp0.5 ，g0=1-asas，max13-1 ，] （9）

式（8）中：[Θ]是顆粒相擬溫度，定義為[Θ=13u'su's]；[u's]為顆粒的脈動(dòng)速度，表征顆粒之間因碰撞而產(chǎn)生的小尺度脈動(dòng)強(qiáng)度；式右邊第1項(xiàng)代表由于變形作用導(dǎo)致的能量生成；第2項(xiàng)是能量擴(kuò)散；第3項(xiàng)是由于非彈性碰撞導(dǎo)致的能量耗散；最后1項(xiàng)是在2個(gè)相之間脈動(dòng)能量的凈轉(zhuǎn)移率；[γ]是碰撞能量耗散；q是脈動(dòng)能量通量；k是脈動(dòng)能量傳導(dǎo)率；e是顆粒碰撞的恢復(fù)系數(shù)；g0是徑向分布函數(shù)；[αs，max]是對(duì)于任意填充的最大固體體積分?jǐn)?shù)。

1.3 曳力模型

當(dāng)黏性流體繞流物體時(shí)，物體會(huì)受到壓力和切向應(yīng)力的作用。在沿著物體橫截面的流動(dòng)平面中，物體所受的合力F可以分解為與來(lái)流方向同向的作用力FD和垂直于來(lái)流方向的升力FL。FD與流體和物體之間的相對(duì)速度有關(guān)，當(dāng)物體運(yùn)動(dòng)速度大于流體流速時(shí)，F(xiàn)D即為物體在流體中受到的阻力。反之，當(dāng)流體流速大于物體運(yùn)動(dòng)速度時(shí)，F(xiàn)D受到的就是流體的曳力。

本文對(duì)氣固兩相間的曳力系數(shù)β選用Gidaspow模型，該模型適用于密相氣固流。當(dāng)αg<0.8時(shí)系數(shù)[β]采用Ergun提出的公式，而當(dāng)αg≥0.8時(shí)系數(shù)β采用Wen、Yu的公式：

[β=150α2sμgαgd2s+1.75ρgαsug-usds ，αg<0.8，β=34Cdαsαgρgug-usdsα-2.65g ，αg≥0.8，] （10）

[Cd=24Res1+0.15Res0.687 ，Res<1 000 ，Cd=0.44 ，Res≥1 000 ，] （11）

[Res=αsρgug-usdsμg]。 （12）

2 數(shù)值計(jì)算說(shuō)明

2.1 幾何模型及網(wǎng)格劃分

以上引式氣力輸送發(fā)送罐作為研究對(duì)象，幾何模型尺寸如圖1所示。模型整體高度為700 mm；填料高度為550 mm；上端直徑為200 mm；底部為空氣入口，直徑為50 mm；出料管位于模型中心處，豎直向上，直徑為10 mm。由圖1所示可知，發(fā)送罐模型為軸對(duì)稱(chēng)模型，內(nèi)部流場(chǎng)流動(dòng)情況亦呈對(duì)稱(chēng)分布，因此取整體模型四分之一進(jìn)行研究。利用ICEM軟件對(duì)幾何模型進(jìn)行結(jié)構(gòu)網(wǎng)格劃分，如圖2所示。

2.2 網(wǎng)格無(wú)關(guān)性檢驗(yàn)

對(duì)于氣固兩相流流動(dòng)，為使網(wǎng)格捕捉到所需流場(chǎng)的細(xì)節(jié)，要求網(wǎng)格的最大尺寸應(yīng)低于最小的細(xì)節(jié)特征，即網(wǎng)格尺寸要小于顆粒直徑的1/10[15]。如圖3所示，分別劃分10萬(wàn)、40萬(wàn)、80萬(wàn)、100萬(wàn)網(wǎng)格，發(fā)現(xiàn)80萬(wàn)網(wǎng)格滿(mǎn)足求解需要，繼續(xù)增加網(wǎng)格數(shù)對(duì)模擬結(jié)果影響較小。

2.3 數(shù)值模擬條件設(shè)置

容器底部進(jìn)氣口類(lèi)型設(shè)置為速度入口，湍流強(qiáng)度和水力直徑分別為5%，50 mm；出料口類(lèi)型設(shè)置為出流條件；壁面條件選用無(wú)滑移壁面；初始物料填充高度設(shè)置為550 mm；采用SIMPLE算法即求解壓力耦合方程組的半隱式方法進(jìn)行求解，空間離散格式選用具有較高精度的二階迎風(fēng)格式和QUICK格式，防止一階精度導(dǎo)致的數(shù)值上的擴(kuò)散誤差。湍流模型選用Realizable k-ε雙方程模型，但該湍流模型適用于充分發(fā)展的湍流，不能應(yīng)用于近壁面區(qū)域，因此要把第1層網(wǎng)格的計(jì)算節(jié)點(diǎn)放置在湍流核心區(qū)內(nèi)。而過(guò)渡層與粘性底層區(qū)域中的物理量由壁面函數(shù)來(lái)計(jì)算完成。[Y+]的計(jì)算公式：[Y+=u*y/ν]；式中u*是近壁面摩擦速度；y為第1層網(wǎng)格節(jié)點(diǎn)與壁面的間距；[ν]為流體的運(yùn)動(dòng)粘度。在計(jì)算過(guò)程中設(shè)置邊界層，第1層的高度為0.1 mm，增長(zhǎng)率為1.1，計(jì)算得到[Y+]在30～300之間，滿(mǎn)足壁面函數(shù)法的要求。

3 模擬的結(jié)果與討論

3.1 曳力系數(shù)的修正和誤差分析

首先將模擬數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行校正，由圖4可以發(fā)現(xiàn)Gidaspow模型在速度較低時(shí)與實(shí)驗(yàn)數(shù)據(jù)的趨勢(shì)一致，在數(shù)值上低于實(shí)驗(yàn)數(shù)據(jù)；但當(dāng)速度增大到0.113 m/s時(shí)，需要對(duì)Gidaspow模型中的曳力系數(shù)進(jìn)行修正。校正之后如圖5所示，模擬數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)基本吻合。

從入口速度與出口質(zhì)量流量的關(guān)系可以看出，隨著入口氣速的增大，出口的質(zhì)量流量先增大后減小，在u=0.118 m/s時(shí)，出口質(zhì)量流量達(dá)到最大值為0.021 kg/s，再增大入口氣速，出口的質(zhì)量流量逐漸降低。此外，衡量輸送能力強(qiáng)弱還與輸送總耗時(shí)有關(guān)。從圖6可以看出，輸送總耗時(shí)在速度較低的時(shí)候隨著入口氣速的增大而減小，當(dāng)超過(guò)某一值（u=0.12 m/s）時(shí)，輸送總耗時(shí)基本保持不變，最終維持在85 s左右。

這對(duì)于工程中實(shí)際的意義是獲取最具性?xún)r(jià)比的操作條件。所以，以增大入口氣體速度來(lái)提高輸送量的方式前期作用效果明顯，但當(dāng)速度達(dá)到一定值之時(shí)，改善效果不再明顯，只會(huì)造成成本顯著提高。

3.2 出口體積分?jǐn)?shù)和質(zhì)量流量的波動(dòng)

工程運(yùn)用上，輸送的穩(wěn)定性一直是關(guān)注的熱點(diǎn)。本節(jié)主要對(duì)輸送過(guò)程中出口煤粉的質(zhì)量流量和體積分?jǐn)?shù)的波動(dòng)進(jìn)行分析。圖7是不同入口氣速下出口質(zhì)量流量隨著時(shí)間的變化，圖8是不同入口氣速下出口固體體積分?jǐn)?shù)隨時(shí)間的變化。從圖中可以看出質(zhì)量流量和體積分?jǐn)?shù)都呈上下波動(dòng)狀，其波動(dòng)程度可用標(biāo)準(zhǔn)差來(lái)表征。標(biāo)準(zhǔn)差值越大表示波動(dòng)程度越大，輸送過(guò)程越不穩(wěn)定。波動(dòng)程度劇烈會(huì)增大對(duì)管路的沖擊，造成管路顫動(dòng)，甚至是引發(fā)安全事故。

從圖7可以看出，在入口氣速較低的時(shí)候，出口的質(zhì)量流量波動(dòng)的整體振幅較小，隨著入口氣速的增大，波動(dòng)的振幅也逐漸增大，其波動(dòng)情況與入口氣體速度不是單純的線(xiàn)性關(guān)系。從圖8中可以看出出口固體體積分?jǐn)?shù)隨氣速變化并不明顯，在速度為0.1 m/s時(shí)，波動(dòng)程度相對(duì)較小。

從圖9可以看出，出口固體體積分?jǐn)?shù)的標(biāo)準(zhǔn)差隨著入口氣速的增大，先減小后增大；而出口固體體積分?jǐn)?shù)的均方根則隨著入口氣速的增大，先增大后減小。通過(guò)對(duì)出口固體體積分?jǐn)?shù)的均方根（αRMS）和標(biāo)準(zhǔn)差（αS）進(jìn)行曲線(xiàn)擬合可得出方程

[αRMS=0.207+u-4.10u2]， （13）

[αS=0.15-0.50u+2.19u2。] （14）

通過(guò)計(jì)算可以得出以下結(jié)論：在入口氣速u(mài)=0.10～0.12 m/s時(shí)，出口的質(zhì)量流量最大、出口的固體體積分?jǐn)?shù)的均方根最大并且出口固體體積分?jǐn)?shù)的標(biāo)準(zhǔn)差最小。

3.3 發(fā)送罐內(nèi)氣固兩相流流動(dòng)機(jī)理分析

取三維模型1/2中心對(duì)稱(chēng)面作為觀(guān)測(cè)面。圖中X代表氣力輸送發(fā)送罐軸向尺寸，Y代表氣力輸送發(fā)送罐徑向尺寸。如圖10、圖11所示，進(jìn)入發(fā)送罐的空氣主要沿著阻力較小的罐體壁面進(jìn)行運(yùn)動(dòng)，從流線(xiàn)運(yùn)動(dòng)規(guī)律可以看出，在壁面處形成強(qiáng)烈漩渦，這些強(qiáng)烈漩渦對(duì)罐內(nèi)平行于Z軸截面上的煤粉運(yùn)動(dòng)產(chǎn)生了兩個(gè)方向的作用力：一個(gè)是展向的運(yùn)動(dòng)，即Y方向，該方向作用力使煤粉沿著展向進(jìn)行移動(dòng)，增強(qiáng)煤粉與空氣的混合程度；另一個(gè)是流向的運(yùn)動(dòng)，即X方向，該方向的作用力是煤粉運(yùn)動(dòng)的主動(dòng)力，力越大煤粉運(yùn)動(dòng)的越快，進(jìn)而進(jìn)入輸送管的流量越大。

從空氣流線(xiàn)圖可以看出，由入口進(jìn)入的空氣可分為兩部分，一部分穿過(guò)罐中物料沉積層，加速物料整體的流動(dòng)，另一部分沿著壁面運(yùn)動(dòng)，在壁面附近與煤粉進(jìn)行能量傳遞，渦流現(xiàn)象非常明顯。從圖中可以發(fā)現(xiàn)，在接近發(fā)送罐頂部的空氣域存在兩個(gè)旋轉(zhuǎn)方向相反的空氣漩渦，因?yàn)檠刂苈愤\(yùn)動(dòng)的空氣流首先達(dá)到罐體頂部，而沿著罐體壁面的空氣流隨后到達(dá)罐體頂部區(qū)域，這兩股空氣流在罐體頂部區(qū)域相遇，形成了外層逆時(shí)針?lè)较蚝蛢?nèi)層順時(shí)針?lè)较虻目諝鉁u。從煤粉流線(xiàn)圖可以看出，首先煤粉受重力的作用整體往下運(yùn)動(dòng)，其次受到氣流的曳力作用，煤粉的密度分布發(fā)生變化。隨著輸送時(shí)間的增加，罐中物料不斷減少，頂部的空氣渦尺度也隨之越來(lái)越大。從圖中可以看出煤粉在流向和展向上的運(yùn)動(dòng)速度都受到了空氣渦強(qiáng)度和方向的影響。

對(duì)于流動(dòng)產(chǎn)生大量的漩渦可以從兩方面進(jìn)行考濾。一方面，大量的漩渦會(huì)造成劇烈的湍流耗散，增大能量的消耗；另一方面，漩渦可以增大湍流強(qiáng)度，增強(qiáng)脈動(dòng)程度，使氣固之間的混合更加均勻，有利于兩相之間的動(dòng)量和能量傳遞。如圖12所示，在軸向Z = 550 mm的截面觀(guān)測(cè)煤粉沿徑向的運(yùn)動(dòng)情況，可以看出煤粉的運(yùn)動(dòng)受到空氣流的影響，沿Z軸方向翻滾，呈沸騰狀，加速了氣固兩相間的動(dòng)量和能量的交換。

4 結(jié)論

本文采用結(jié)合了顆粒動(dòng)力學(xué)理論和Realizable k-ε的歐拉雙流體模型，對(duì)上引式氣力輸送發(fā)送罐的流動(dòng)機(jī)理和出料過(guò)程流動(dòng)特性進(jìn)行了數(shù)值模擬研究。對(duì)Gidaspow曳力模型進(jìn)行了修正，使模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù)相吻合。對(duì)出料口固體質(zhì)量流量和固體體積分?jǐn)?shù)的波動(dòng)特性進(jìn)行分析，考察了不同入口氣速下對(duì)出料穩(wěn)定性的影響，在此基礎(chǔ)上，對(duì)罐內(nèi)氣固兩相的流動(dòng)機(jī)理進(jìn)行分析，得出以下結(jié)論：

1）隨入口氣速的增加，出料口煤粉的質(zhì)量流量和體積分?jǐn)?shù)都呈現(xiàn)先增大后減小的變化趨勢(shì)，其穩(wěn)定性呈現(xiàn)先增加后降低的變化趨勢(shì)，入口氣速在v=0.10～0.12 m/s范圍內(nèi)存在一個(gè)最優(yōu)值，可使波動(dòng)最小、煤粉體積分?jǐn)?shù)和出料量最大。

2）初始階段隨著入口氣速的增加，輸送物料所需時(shí)間逐漸減小，但當(dāng)入口氣速超過(guò)0.12 m/s時(shí)，耗時(shí)基本保持恒定。

3）發(fā)送罐中的氣體運(yùn)動(dòng)以各種類(lèi)型的渦流運(yùn)動(dòng)為主，作用方式是各種尺度渦的疊加、碰撞和干涉；運(yùn)動(dòng)形式分解為流向和展向，驅(qū)動(dòng)煤粉的展向擴(kuò)散和流向移動(dòng)，展向擴(kuò)散加速煤粉和空氣的混合程度，氣流曳力驅(qū)動(dòng)煤粉離開(kāi)輸送罐。

參考文獻(xiàn)：

[1] 楊倫，謝一華. 氣力輸送工程 [M]. 北京：機(jī)械工業(yè)出版社，2006：172-179.

[2] 陸慧林，劉文鐵，孫永立，等. 稠密氣固兩相湍流流動(dòng)的實(shí)驗(yàn)和數(shù)值模擬 [J]. 力學(xué)學(xué)報(bào)，2000，32（4）：385-391.

[3] YUE C，HONG G，SUN J，et al. Numerical simulation of gas-solid flow in a conveying vessel [J]. Powder Technology， 2012，226（8）：34-42.

[4] JONES M G，WILLIAMS K C. Predicting the mode of flow in pneumatic conveying systems [J]. Particuology， 2008，6（5）：289-300.

[5] Setia G，Mallick S S. Modelling fluidized dense-phase pneumatic conveying of fly ash[J]. Powder Technology， 2015，270：39-45.

[6] 徐貴玲，陳曉平，梁財(cái)，等. 煤粉上出料式發(fā)送罐氣力輸送特性[J]. 化工學(xué)報(bào)，2012，63（6）：1709-1716.

[7] 徐貴玲，盧平，許盼，等. 煤粉外水含量對(duì)上出料式發(fā)送罐供料特性的影響[J]. 化工學(xué)報(bào)，2016，67（7）：2767-2776.

[8] 賀春輝，陳先梅，倪紅亮，等. 發(fā)送罐出料方式對(duì)煤粉高壓密相氣力輸送特性的影響[J]. 中國(guó)電機(jī)工程學(xué)報(bào)，2012，32（23）：32-39.

[9] ZHANG Y，LEI F，WANG S，et al. A numerical study of gas-solid flow hydrodynamics in a riser under dense suspension upflow regime[J]. Powder Technology，2015，280 （20）：227-238.

[10] FATTAHI M，HOSSEINI S H，AHMADI G. CFD simulation of transient gas to particle heat transfer for fluidized and spouted regimes[J]. Applied Thermal Engineering，2015，105（6）：385-396.

[11] 沈志恒，李巍，陸慧林，等. 不同曳力模型對(duì)提升管內(nèi)氣固流動(dòng)特性的影響[J]. 石油機(jī)械，2011，39 （10）：32-38.

[12] CLOETE S，JOHANSEN S T，AMINI S. Investigation into the effect of simulating a 3D cylindrical fluidized bed reactor on a 2D plane [J]. Powder Technology，2013，239（17）：21-35.

[13] CHAN C K，GUO Y C，LAU K S. Numerical modeling of gas-particle flow using a comprehensive kinetic theory with turbulence modulation [J]. Powder Technology，2005，150（1）：42-55.

[14] TSUJI Y. Multi-scale modeling of dense phase gas-particle flow[J]. Chemical Engineering Science，2007，62（13）：3410-3418.

[15] AGRAWAL K，LOEZOS P N，SYAMLAL M，et al. The role of meso-scale structures in rapid gas-solid flows[J]. Journal of Fluid Mechanics，2001，445：151-185.

[責(zé)任編輯 田 豐]