1、第 41卷第 3期 2010年 5月 鍋 爐 技 術(shù)BOIL ER TECHNOLO GYVol. 41, No. 3May. ,2010 收稿日期 :2009205221作者簡介 :王建軍 (19712 , 男 , 博士 , 副教授 , 主要從事流態(tài)化、 多相流分離的研究。 文章編號 : CN3121508(2010 0520021206循環(huán)流化床鍋爐爐膛內(nèi)氣固兩相流的數(shù)值模擬王建軍 1, 李東芳 2, 姬廣勤 1, 金有海 1(1. 中國石油大學(xué) (華東 機(jī)電工程學(xué)院 , 山東 東營 257061; 2. 海洋石油工程股份有限公司 ,河北 塘沽 300451關(guān)鍵詞 : 循環(huán)流化床鍋爐 ;

2、雙流體模型 ; 氣固兩相流 ; 數(shù)值模擬摘 要 : 利用 CFD 軟件 Fluent , ( 流的宏觀流動特性進(jìn)行了數(shù)值模擬。 準(zhǔn)確性。 通過定性與定量分析 , , 核” 流動結(jié)構(gòu)及顆粒軸向速度 中心處向上 , , 沿軸向爐膛中下部區(qū)域及沿 同時(shí) , 操作條件對顆粒軸向速度的影響都表現(xiàn)為中心區(qū)域顆粒向 邊壁處的氣固兩相流動規(guī)律還有待于進(jìn)一步研究。 中圖分類號 : T K 227. 1 文獻(xiàn)標(biāo)識碼 : A0 前 言 目前 , 對于循環(huán)流化床內(nèi)的氣固兩相流主要 集中在對循環(huán)流化床反應(yīng)器 1-2及鼓泡床 3-4的 研究 。 循環(huán)流化床鍋爐爐膛內(nèi)和循環(huán)流化床反 應(yīng)器內(nèi)的氣固兩相流動特性有一定的差別 ,

3、 不僅 體現(xiàn)在燃燒室的高徑比 , 循環(huán)系統(tǒng)中采用的顆粒 循環(huán)流率 , 床料的特性 , 而且循環(huán)流化床鍋爐有 二次風(fēng)的加入 , 對循環(huán)流化床鍋爐內(nèi)氣固兩相流 的研究并不多 5-6。 本文以歐拉雙流體模型和顆 粒動力學(xué)理論為基礎(chǔ)采用 CFD 軟件 Fluent 研究 對循環(huán)流化床鍋爐爐膛內(nèi)氣固兩相流動特性的 影響進(jìn)行數(shù)值模擬 。1 計(jì)算模型及數(shù)值方法1. 1幾何模型及計(jì)算條件 圖 1為整個(gè)循環(huán)流化床鍋爐循環(huán)系統(tǒng)幾何 模型及網(wǎng)格模型 , 模型按照工業(yè)裝置 12 1縮小 得到 。 沿流化床的徑向?yàn)?X 方向 , 沿爐膛的高度 方向?yàn)?Y 方向 。 1. 2計(jì)算方法 本文中選用歐拉雙流體模型 7, 固相

4、應(yīng)用顆 粒動力學(xué)理論 , 氣相采用 -s 湍流模型 。對于壓 力 -速度耦合 , 采用 Phase Coupled SIM PL E算圖 1 循環(huán)流化床鍋爐二維模型及爐膛區(qū)域網(wǎng)格法 。 速度的求解被相耦合 , 但是用分離的方式 , 采用多重網(wǎng)格方案同時(shí)地求解各相速度分量形 成的矢量方程 。空隙率選用 QU IC K 格式 , 以提 高計(jì)算精度 , 使得相之間界面清晰 。其余項(xiàng)選用 二階迎風(fēng)格式 。計(jì)算中氣相介質(zhì)為常溫空氣 , 固相介質(zhì)為石 英砂顆粒 , 顆粒平均直徑為 0. 245mm , 顆粒密度 =2300kg/m 3, 堆積密度 =1367kg/m 3。(1 入口邊界條件 :氣固兩相的入

5、口在計(jì)算中均采用均勻的速度入口邊界條件 , 一次風(fēng)進(jìn)口u 0=2. 0m/s , 二次風(fēng)入口處 u sec =2. 86m/s , 采用-s 模型 , 模型參數(shù)設(shè)置見文獻(xiàn) 8。(2 出口邊界條件 :出口處采用壓力出口邊 鍋 爐 技 術(shù) 第 41卷界條件 , 出口處設(shè)定為標(biāo)準(zhǔn)大氣壓 。顆粒的出口邊界條件設(shè)置為連續(xù)性邊界條件 。(3 壁面邊界條件 :顆粒采用 Johnson 和J ackson 部分滑移邊界條件 , 通過引入不同的鏡面反射系數(shù)及顆粒與壁面的恢復(fù)系數(shù) , 對顆粒在邊壁附近的運(yùn)動進(jìn)行描述 。氣相在壁面處為無滑移邊界 , 該處速度為零 。(4 初始邊界條件 :初始計(jì)算時(shí) , 設(shè)立管處有顆

6、粒堆積 , 以便于形成循環(huán)系統(tǒng) , 顆粒堆積靜止高度 h =1. 2m , 堆積密度為 0. 63; 在計(jì)算區(qū)域的其余部分均為氣體所充滿 。22. 1料量和氣固之間的動量交換。 床層壓降的計(jì)算能用于判斷數(shù)值模型計(jì)算的準(zhǔn)確性 9。 張瑞英 10等提出循環(huán)流化床的床層壓降的計(jì)算公式 , 發(fā)現(xiàn)床層壓降與顆粒固含率基本成線性關(guān)系。 該計(jì)算公式?jīng)]有考慮床層邊壁顆粒向下流動的影響。p =1. 1335p g (1- h +0. 0047(1本文也只選取了核心區(qū)域的數(shù)值來計(jì)算床層截面平均空隙率 。圖 2所示為數(shù)值計(jì)算值與理論分析值的對比結(jié)果 。如圖所示 , 沿爐膛從底部到頂部 , 床層壓降先急劇下降 , 隨

7、后逐漸減小至恒定值 。 從圖中可以看出兩者吻合較好 , 證實(shí)了本次數(shù)值計(jì)算的可靠性 。圖 2 軸向壓降的模擬值與文獻(xiàn) 11值的對比2. 2氣固兩相瞬時(shí)流場分布 圖 3為不同特征時(shí)刻循環(huán)流化床鍋爐爐膛內(nèi)顆粒固含率的分布 。由圖可直觀地看到 , 在計(jì)算的初始階段 , 顆粒在床層底部顆粒并不是均勻彌散分布的 , 而是部分顆粒會形成不規(guī)則形狀的團(tuán)聚物 , 在床層內(nèi)不斷的聚合與解體 , 同時(shí)與氣體發(fā)生強(qiáng)烈的混合 。未形成團(tuán)聚物的顆粒很容易被氣流攜帶懸浮上升到床層的上部 ,部分顆粒團(tuán)聚物由于其形態(tài)發(fā)生變化 , 使其迎 風(fēng)面積增加 , 也會被氣流帶著上升 。其中一部 分顆粒會被氣流攜帶出爐膛 , 進(jìn)入旋風(fēng)分

8、離器 進(jìn)行外循環(huán) 。隨著時(shí)間的推移 , 顆粒會懸浮于 整個(gè)爐膛內(nèi) 。在爐膛上部稀相段 , 沿床層徑向 方向 , 顆粒也存在團(tuán)聚現(xiàn)象 , 且容易發(fā)生在邊壁 處 , 因此顆粒固含率在邊壁處較高 , 中心處低 , 呈現(xiàn)典型 的環(huán) -核結(jié) 構(gòu) 沿床 層高度 方 向 , 圖 3 不同時(shí)刻爐膛內(nèi)顆粒固含率的瞬時(shí)分布 圖 4為爐膛內(nèi)顆粒速度場全貌圖及沿軸向 不同區(qū)域的局部放大圖 。由圖可知 , 顆粒在爐膛 內(nèi)中心氣速比較高 , 邊壁區(qū)域氣速較低 , 并且由 局部圖可看出 , 顆粒呈現(xiàn)在中心區(qū)域向上流動 , 在邊壁處向下流動的內(nèi)循環(huán)流動結(jié)構(gòu) 。這種內(nèi) 循環(huán)結(jié)構(gòu)沿爐膛高度方向的中部懸浮段比較明 顯 。 在爐膛頂

9、部 , 由于出口結(jié)構(gòu)的約束作用 , 顆 粒在該處做折轉(zhuǎn)運(yùn)動 , 部分顆粒與爐膛頂部相撞 后 , 被反彈沿邊壁向下運(yùn)動 , 形成內(nèi)循環(huán)物料 。 在爐膛底部 , 可看到沿邊壁下落的顆粒在到達(dá)床 層底部之后 , 會被一次風(fēng)加速 , 折轉(zhuǎn)向上運(yùn)動 , 形 成爐膛內(nèi)循環(huán)的流動特性 。同時(shí)在二次風(fēng)口區(qū)圖 4 爐膛內(nèi)顆粒速度場分布22 第 3期 王建軍 , 等 :循環(huán)流化床鍋爐爐膛內(nèi)氣固兩相流的數(shù)值模擬 域 , 氣固兩相的混合比較強(qiáng)烈 , 二次風(fēng)的擾動使 得沿二次風(fēng)口一側(cè)邊壁下落的顆粒運(yùn)動發(fā)生改 變 , 因此改善了氣固之間沿床層軸向與徑向上的 混合 。 因此 , 這部分顆粒很難到達(dá)床層底部 , 會 在一次風(fēng)

10、和二次風(fēng)的作用下 , 改變運(yùn)動方向折轉(zhuǎn) 向上運(yùn)動 , 可以看出二次風(fēng)對顆粒沿軸向向上運(yùn) 動會有加速的作用 , 同時(shí)會阻止顆粒到達(dá)下部密 相區(qū) , 減小底部密相區(qū)的燃燒負(fù)荷 。2. 3顆粒時(shí)均流動規(guī)律 對氣固兩相流場進(jìn)行非穩(wěn)態(tài)計(jì)算 , 計(jì)算時(shí)間 持續(xù)到 35s 以獲得穩(wěn)定的結(jié)果 , 并取最后 20s 為時(shí)間樣本 。圖 5固含率分布 , , 在近 壁處固含率高的 -核” 流動結(jié)構(gòu) :在爐膛中心 區(qū)域 , 顆粒固含率較小且分布比較均勻 , 并沿軸 向變化較小 ; 而在爐膛邊壁區(qū) , 顆粒固含率則明 顯高于核心區(qū)且隨著半徑的增大而增大 , 在壁面 處達(dá)到最大值 。同時(shí) , 隨著爐膛高度的增加 , 顆

11、粒固含率無論是在中心區(qū)域還是邊壁區(qū)域都是 逐漸減小的 , 但邊壁處顆粒固含率降低較為明 顯 , 而且中心區(qū)域與邊壁區(qū)域的顆粒固含率差也 會減小 , 即顆粒固含率沿徑向分布趨于均勻 。也 正如圖 3中看到的爐膛內(nèi)呈明顯的上稀下濃的 分布狀況 。 在二次風(fēng)進(jìn)口附近 , 由于二次風(fēng)的橫 向作用致使進(jìn)風(fēng)口對面的壁面上顆粒固含率比 進(jìn)風(fēng)口側(cè)壁面上顆粒固含率高 , 但隨著高度的增 加 , 氣固混合愈加均勻 , 顆粒固含率分布 趨于 平緩 。圖 5 不同爐膛高度處顆粒時(shí)均固含率分布圖 圖 6是顆粒固含率沿床層高的軸向分布 。 從圖中可知 , 顆粒固含率總體呈現(xiàn)上稀下濃的分 布規(guī)律 。 顆粒固含率的最大值出現(xiàn)

12、在返料口附 近 (h =0. 3m 處 。這主要是由于二次風(fēng)的加入在該處外循環(huán)物料的返回固含率比較高 。隨著 床層高度的增加 , 顆粒固含率在 h =0. 3m 0. 7m 的范圍內(nèi)又迅速下降 , 即存在一個(gè)快速轉(zhuǎn)變 區(qū) , 主要是因?yàn)樵谠搮^(qū)域由于二次風(fēng)的加入 , 使 得床層孔隙率增加 。隨后床層中上部顆粒固含 率逐漸減小且趨于均勻 , 占據(jù)大部分床層 。在二 次風(fēng)上部區(qū)域的稀相段 , 顆粒固含率沿軸向的分 布是呈單調(diào)指數(shù)分布的 5中的結(jié)論 是一致的 圖 6 顆粒固含率沿床高的軸向分布圖 圖 7為不同爐膛高度處 , 顆粒時(shí)均軸向速度分布 。 選取軸向 3個(gè)截面 y =1m ,2m ,3m ,

13、主 要是考慮二次風(fēng)位置 y =0. 9m 以上到爐膛出口 y =3. 75m 以下爐膛高度范圍內(nèi) , 顆粒的流動特 征 。 從圖中可以看出 , 在各個(gè)位置上 , 顆粒的軸 向速度分布相似 , 均為在中心區(qū)域?yàn)檎?, 在邊壁 處為負(fù) , 這意味著顆粒的運(yùn)動形式為中心處向 上 、 近壁面處向下的內(nèi)循環(huán)流動結(jié)構(gòu) 。這與文獻(xiàn) 5研究的規(guī)律相符 。在該操作條件下 , 沿爐膛 高度方向 , 中心區(qū)域顆粒的上行速度基本隨著爐 膛高度的增加而增加 , 同時(shí)曲線由分布平緩變成 拋物線分布趨勢 。在邊壁處 , 顆粒下行速度隨著 爐膛高度的增加而減小 。圖 7 不同爐膛高度處顆粒時(shí)均軸向速度分布32 鍋 爐 技 術(shù)

14、 第 41卷 圖 8為不同爐膛高度截面處顆粒徑向速度分布 , 在爐膛稀相段 , 顆粒徑向速度分布呈 “倒 S ”型分布 , 表明顆粒在環(huán)形區(qū)域和核心區(qū)域之間的運(yùn)動方向 , 顆粒由核心區(qū)域向環(huán)形區(qū)域運(yùn)動 , 而后形成絮狀物向下運(yùn)動 。如圖所示 , 顆粒在徑向大部分區(qū)域內(nèi) , 是向左側(cè)邊壁 (即 X 軸負(fù)方向 運(yùn)動的 。 在 y =1. 0m 處 , 由于二次風(fēng)的影響 , 以及在 y =3. 0m 處 , 由于爐膛上部出口結(jié)構(gòu)的影響 ,都增加了顆粒沿徑向的運(yùn)動趨勢 ; 說明沿爐膛高度中部懸浮段 , 顆粒沿徑向的運(yùn)動趨勢不大 , “環(huán)-核” 流動結(jié)構(gòu)是比較穩(wěn)定的 。圖 8 不同爐膛高度處顆粒時(shí)均徑向

15、速度分布2. 4操作條件對顆粒流動特性的影響 通過改變操作氣速以及顆粒循環(huán)流率 , 研究顆粒宏觀流動特性 , 主要從顆粒沿床層軸向與徑向的非均勻流動特性來分析 。圖 9為不同操作氣速及質(zhì)量循環(huán)流率下的 顆粒時(shí)均固含率徑向分布 。由圖 9(a 為一定顆 粒循環(huán)流率下 , 改變操作氣速 , 爐膛軸向截面上 顆粒固含率徑向分布規(guī)律 。由圖可知 :隨著操作 氣速的增大 , 床層截面平均固體顆粒固含率下 降 , 顆粒固含率的徑向分布相對趨于均勻 , 即邊 ; 氣速增 大 ,;, 基本處于稀相氣力輸送 的狀態(tài) , 同時(shí)顆粒固含率沿床層高度方向基本保 持不變 。圖 9(b 為一定操作氣速下 , 改變顆粒循環(huán)

16、流 率 , 爐膛軸向截面上顆粒固含率徑向分布規(guī)律 。 顆粒循環(huán)流率越大 , 整個(gè)截面上顆粒固含率越 大 。 在爐膛中下部區(qū)域 , 顆粒循環(huán)流率對邊壁處 的顆粒固含率影響要顯著于對中心處的影響 。 在爐膛上部區(qū)域 , 顆粒循環(huán)流率增大使得整個(gè)截 面沿徑向顆粒固含率都呈現(xiàn)明顯增大趨勢 , 這主 要是因?yàn)槌隹诮Y(jié)構(gòu)的影響 , 使得到達(dá)爐膛頂部 的顆粒回流量增加 , 出口結(jié)構(gòu)的影響區(qū)域延伸至 y =3m 處 。圖 9 不同操作條件下 , 顆粒時(shí)均固含率徑向分布 由圖 9可見 , 在不同的操作條件下 , 顆粒固含率沿徑向都呈現(xiàn)中心低 , 邊壁高的的分布規(guī)律 。 雖然爐膛各個(gè)截面顆粒固含率均隨表觀氣速 U

17、g 的減小或顆粒循環(huán)速率 G s 的增大而增加 ,但在爐膛不同高度處 , 不同徑向區(qū)域的顆粒固含率對操作條件變化的敏感程度不同 。操作條件 改變 , 主要體現(xiàn)在對爐膛中下部區(qū)域的影響比較 大 , 同時(shí)邊壁區(qū)的顆粒固含率受 U g 和 G s 的影響 相對比較顯著 。圖 10為在不同操作氣速及質(zhì)量循環(huán)流率條 42 第 3期 王建軍 , 等 :循環(huán)流化床鍋爐爐膛內(nèi)氣固兩相流的數(shù)值模擬 件下 , 不同爐膛高度截面上的顆粒軸向速度沿徑 向的分布規(guī)律。圖 10(a 表明隨著操作氣速增 大 , 爐膛內(nèi)任一截面的顆粒速度增加。中心區(qū)域 顆粒向上運(yùn)動的速度增大的效果要顯著于邊壁處 向下運(yùn) 動 速 度 的 變

18、化。在 爐 膛 下 部 區(qū) 域 y =1. 02. 0m 高度范圍內(nèi) , 顆粒加速能力增強(qiáng) , 中 心區(qū)域向上運(yùn)動的速度增大 , 顆粒的加速效果顯 著。 這主要是由于在保持一次風(fēng)率不變的條件 下 , 增大一次風(fēng)速 , 二次風(fēng)速也會增加 , 因此提高 了二次風(fēng)對顆粒的再次加速作用。 隨著高度的增 加 , 二次風(fēng)作用減弱 , 小。 , 分區(qū)域 , 隨著氣速的增加而增加 , 而在頂部區(qū)域 y =3m 處 , 隨著氣速的增加而減小。這與黎小 波 11的結(jié)論不完全相符。分析原因可能是由于 :氣速提高 , 顆粒的橫向擴(kuò)散率增加 , 顆粒向邊壁處 運(yùn)動的趨勢增強(qiáng) ; 在床層下部區(qū)域 , 顆粒固含率比 較大

19、, 操作氣速使顆粒在邊壁區(qū)域的聚集效應(yīng)更 加明顯 , 顆粒所受重力和曳力都增加 , 但對于大的 , 所以 , 而在爐 , 邊壁處顆 , 顆粒速度 圖 10 不同操作條件下 , 時(shí)均顆粒軸向速度徑向分布 由圖 10(b 可知 , 隨著顆粒循環(huán)流率的增加 ,不同徑向區(qū)內(nèi)顆粒速度受影響的程度不同 。同 樣表現(xiàn)為中心區(qū)域顆粒的向上速度增大的效果 要顯著于邊壁處向下運(yùn)動速度的變化 。在中心 區(qū)域 , 顆粒向上的運(yùn)動速度隨著顆粒循環(huán)流率的 增加而增大 , 這是由于顆粒循環(huán)流率的增加 , 爐 膛內(nèi)截面顆粒固含率增加 , 在床層中心區(qū)域局部 氣速增大 , 因此氣固曳力增加 , 顆粒向上的速度 增加 。 而邊壁

20、處向下的顆粒速度 , 在床層大部分 區(qū)域 , 隨著顆粒循環(huán)流率的增加而增加 , 在頂部 區(qū)域 , 影響不顯著 。這是因?yàn)?:邊壁區(qū)局部氣速 減小 , 氣固曳力減小 , 同時(shí)由于顆粒固含率增加 , 顆粒更容易團(tuán)聚 , 重力作用給顆粒的影響也會使 顆粒向下的速度增加 , 在爐膛頂部區(qū)域顆粒充分 發(fā)展 , 顆粒聚集效應(yīng)減弱 , 因此變化不明顯 。同 時(shí) , 顆粒質(zhì)量循環(huán)流率越大 , 在爐膛內(nèi)更容易形 成邊壁 下 行 顆 粒 流 , 即 邊 壁 流 的 延 伸 長 度 也越長 。在不同的操作條件下 , 爐膛內(nèi)顆粒軸向速 度沿徑 向 基 本 呈 現(xiàn) 中 心 區(qū) 域 向 上 , 且 速 度 值 大 , 邊

21、壁 區(qū) 域 向 下 , 且 速 度 值 小 的 分 布 趨 勢 。 同時(shí)操作條件對顆粒軸向速度的影響都 表現(xiàn) 為中心區(qū)域顆粒的向上速度變化的效果 要顯 著于邊壁處向下運(yùn)動速度 的變 化 。邊壁 處的 下行顆粒流由于該區(qū)域內(nèi)復(fù)雜的氣固相 互作 用及顆粒團(tuán)聚的影響使得其規(guī)律還有待 于進(jìn) 一步研究 。3 結(jié) 論 結(jié)合顆 粒動力 學(xué)理 論 , 采 用雙 流體模 型 , 對爐膛內(nèi)氣固兩相流的宏觀流動特性的 影響 進(jìn)行數(shù)值模擬 。將爐膛壓降 沿軸向 的分 布規(guī) 律與理論值對比 , 驗(yàn)證了本次數(shù)值模擬計(jì)算的 準(zhǔn)確性 。5226 鍋 爐 技 術(shù) 第 41 卷 ( 1 通過定性與定量分析 , 得到顆粒固含率 2

22、 吳東垠 ,李吉武 , 許暉 . 水煤漿燃燒技術(shù)在勝利油田的應(yīng)用 J . 工業(yè)鍋爐 ,2005 (3 :23 - 26. 3 D. Gidaspow , R. Bezburuah , J . Ding. Hydrodynamics of cir2 culating fluidized beds C . Fluidization V II proceedings of t he 7t h engineering Foundation conference on fluidization , 1992 :75 - 82. 4 白丁榮 ,金涌 ,俞芷青 ,等 . 快速流化床截面平均空隙率軸向 沿爐膛

23、軸向呈典型的上稀下濃及沿徑向方向呈 環(huán) - 核流動分布 。顆粒速度在中心區(qū)域較高且 向上 ,邊壁區(qū)速度值較低且向下的內(nèi)循環(huán)流動 結(jié)構(gòu) 。 ( 2 爐膛各個(gè)截面處 , 顆粒固含率均隨表觀 氣速 U g 的減小或顆粒循環(huán)速率 Gs 的增大而增 加 ,但在爐膛不同高度處 , 不同徑向區(qū)域的顆粒 固含率對操作條件變化的敏感程度不同 。操作 條件改變 ,主要體現(xiàn)在對爐膛中下部區(qū)域的影響 比較大 ,同時(shí)邊壁區(qū)的顆粒固含率受 U g 和 Gs 的 影響相對比較顯著 。 ( 3 操作條件對顆粒軸向速度的影響都表 現(xiàn)為中心區(qū)域顆粒的向上速度變化的效果要顯 著于邊壁處向下運(yùn)動速度的變化 。邊壁處下降 流的變化還有

24、待于進(jìn)一步研究 。 參考文獻(xiàn) : 1 姜秀民 ,馬玉峰 ,崔志剛 ,等 . 水煤漿流化懸浮高效潔凈燃燒 分布及其影響因素 J . 化學(xué)反應(yīng)工程與工藝 , 1990 ,6 (1 :63 - 71. 技術(shù)研究與應(yīng)用 J . 化學(xué)工程 ,2006 ,34 (1 :62 - 65. (1. College of Mechanical and Elect rical Engineering ,China U niversity of Pet roleum ( East China ,Dongying 257061 ,China ; Off shore Oil Engineering Co , L td.

25、 , Tanggu 300451 , China 2. field ; numerical simulatio n field was performed fo r a circulating fluidized bed boiler . Co mpariso n of t he simulated p ressure drop wit h t he t heoretical value test s t he validatio n of t his numerical simulatio n. Thro ugh t he qualitative and quantitative analy

26、sis , it is analyzed t hat t he particle pat tern t hat t he particle movement is upward in t he core zo ne and down in t he annulus zo ne. In additio n , under different operating co nditio n , particle volume f ractio n dist ributio n is affected p ro minently at t he lower2middle zo ne alo ng t h

27、e axial directio n and at t he wall area alo ng t he radial directio n. Furt hermore , t he effect of operating co nditio n o n t he particle velocit y reflect s t hat t he changing of upward particle velocit y at t he core area is f urt her investigatio n. more o bvio us t han at t he wall area . T

28、he gas2solid flow pat tern at t he wall area needs Abs t ra c t : U sing t he CFD sof t ware Fluent , based o n t he t wo2fluid model co mbining wit h co ncent ratio n dist ributio n show s a co re2annular flow f rame where t he solid is denser in t he annular regio n t han core regio n and t he particle ax