張賢，葛榮存，張守玉，劉青，張縵，楊海瑞，呂俊復(fù)

（1上海理工大學(xué)能源與動(dòng)力工程學(xué)院，上海 200093；2清華大學(xué)熱能工程系，北京 100084）

CFB密相區(qū)大顆粒橫向擴(kuò)散系數(shù)的CPFD模擬

張賢1,2，葛榮存1,2，張守玉1，劉青2，張縵2，楊海瑞2，呂俊復(fù)2

（1上海理工大學(xué)能源與動(dòng)力工程學(xué)院，上海 200093；2清華大學(xué)熱能工程系，北京 100084）

運(yùn)用一種離散單元法（DEM）計(jì)算顆粒流體力學(xué)（CPFD）對尺寸為900 mm×100 mm×1200 mm的準(zhǔn)三維流化床的密相區(qū)大顆粒擴(kuò)散行為進(jìn)行研究。模擬之前，依照前人實(shí)驗(yàn)研究對CPFD方法進(jìn)行驗(yàn)證，模擬結(jié)果與實(shí)驗(yàn)結(jié)果符合較好，證明了CPFD方法模擬的有效性。模擬中通過注入示蹤粒子的方法來研究大顆粒在密相區(qū)中的橫向擴(kuò)散系數(shù)，研究了流化風(fēng)速、顆粒直徑對顆粒橫向擴(kuò)散系數(shù)的影響。模擬結(jié)果顯示，氣泡是引起密相區(qū)內(nèi)顆粒混合的主要因素；隨著流化風(fēng)速增加，顆粒橫向擴(kuò)散系數(shù)變大；隨著顆粒直徑增大，顆粒橫向擴(kuò)散系數(shù)減小。

流化床；密相區(qū)；CPFD；氣泡；擴(kuò)散；流化風(fēng)速

Abstract:Computational particle fluid dynamics (CPFD),a discrete element method (DEM),was utilized to simulate the lateral diffusion behavior of large solid particles in a 900 mm×100 mm×1200 mm quasi-three dimensional fluidized-bed.Before the simulating,verification had been performed in accordance with previous experimental studies to ensure the reliability of CPFD.The effects of the tracer particle diameter and the fluidizing air velocity on lateral dispersion coefficient of large solid particles in dense zone were investigated by the method of tracer particle.The results showed that the bubble behavior was the major factor of solid mixing in dense zone,which had been confirmed in many literatures.An increased gas velocity or a reduced tracer particle diameter resulted in an increased lateral dispersion coefficient.

Key words:fluidized-bed; dense zone; CPFD; bubble; diffusion; fluidization air velocity

引 言

循環(huán)流化床廣泛應(yīng)用于石油、化工、能源、冶金行業(yè)，其氣固混合好、傳質(zhì)強(qiáng)烈的特點(diǎn)，為反應(yīng)的進(jìn)行提供了優(yōu)越的條件。隨著裝備的容量放大，床內(nèi)的反應(yīng)物濃度和溫度的均勻性問題逐漸顯示出來[1-3]。通常，對于循環(huán)流化床，沿提升管高度方向上自下而上可以分為密相區(qū)、過渡區(qū)和稀相區(qū)3個(gè)區(qū)域。底部密相區(qū)是氣固兩相接觸的初始區(qū)域，顆粒在密相區(qū)的混合和擴(kuò)散情況對整個(gè)床層的流動(dòng)、傳質(zhì)和傳熱產(chǎn)生重大影響。由于密相區(qū)流動(dòng)十分復(fù)雜，固含率大，很多測試手段的使用受到限制，測量研究結(jié)果存在很大差異[4]。在此條件下，數(shù)值模擬可以突破實(shí)驗(yàn)條件的限制，獲取流化床內(nèi)局部流動(dòng)和顆粒運(yùn)動(dòng)等詳細(xì)信息[5-7]。目前，流化床內(nèi)的多相氣固流動(dòng)模擬方法，依據(jù)對顆粒相處理方式的不同，大致可以分為兩類：Euler-Euler方法和Euler-Lagrange方法。通常流化床密相區(qū)的模擬方法采用前者，其缺點(diǎn)在于不能對顆粒的尺寸進(jìn)行定義，不能得到不同顆粒在系統(tǒng)中的具體運(yùn)動(dòng)行為，也不能得到粒徑分布對于流化床的流體力學(xué)性能的影響[7-9]。了解密相區(qū)中顆粒行為應(yīng)該采用Lagrange方法。計(jì)算顆粒流體力學(xué)（CPFD）中，對于顆粒相采用的是Lagrange模型，模型中顆粒的尺寸可以定義[10-11]，能研究不同粒徑顆粒擴(kuò)散行為，也能研究粒徑分布對于流化床的流體力學(xué)性能的影響[12]。相比于傳統(tǒng) Lagrange方法，離散單元法（DEM）CPFD方法中真實(shí)顆粒被打包成計(jì)算顆粒，顆粒間的接觸應(yīng)力用顆粒相應(yīng)力模型來描述，這樣可以顯著降低計(jì)算量[13]。本文利用CPFD模型，研究密相區(qū)中顆粒擴(kuò)散行為，分析流化風(fēng)速和示蹤顆粒粒徑等的影響。

1 數(shù)學(xué)模型

CPFD基于 MP-PIC（multi-phase particle-incell）方法與流體耦合的微分方程來描述顆粒的Lagrange運(yùn)動(dòng)[14]。其控制方程如下。

氣體的連續(xù)性方程和動(dòng)量方程為

式中，ug為氣體速度，ρg為氣體密度，εg為氣體體積分?jǐn)?shù)，p為氣體壓力，g為重力加速度，τg為氣體應(yīng)力張量，F(xiàn)為單位體積內(nèi)氣固之間的動(dòng)量交換率，μg為黏性系數(shù)，δij為Kronecker符號(hào)。

顆粒的加速度為

式中，up為顆粒速度，εp為顆粒體積分?jǐn)?shù)，τp為顆粒的法向應(yīng)力。

顆粒的屬性通過內(nèi)插算子被內(nèi)插到Euler網(wǎng)格中，將xp點(diǎn)的顆粒屬性內(nèi)插到i點(diǎn)的網(wǎng)格單元中的內(nèi)插算子為

對于z方向的內(nèi)插算子采用相同的形式，因此任意單元內(nèi)的顆粒體積分?jǐn)?shù)εi,j,k可以表示為

式中，Np為計(jì)算顆粒數(shù)，np,m為計(jì)算顆粒內(nèi)包含的實(shí)際顆粒數(shù)，Vi,j,k為計(jì)算單元體積，Vp,m為顆粒體積。

顆粒速度的隱式方程為

式中，upn+1為顆粒位置處的插值流體速度，?ppn+1為顆粒位置處的插值壓力梯度，?τpn+1為顆粒位置處的插值顆粒應(yīng)力梯度。

一個(gè)時(shí)間步長后新的顆粒位置為

氣固之間的動(dòng)量耦合是通過相間動(dòng)量交換，動(dòng)量單元ξ的相間動(dòng)量交換為

顆粒間相互作用力用顆粒的正應(yīng)力來描述

式中，ps為常數(shù)，1 Pa；ε為10?7數(shù)量級的一個(gè)很小的數(shù)來消除奇點(diǎn)的影響；β的推薦值為2～5；εcp為堆積密度/顆粒密度。

2 橫向擴(kuò)散系數(shù)計(jì)算

橫向擴(kuò)散系數(shù)可基于擴(kuò)散模型[3,15-16]求解。根據(jù)Fick定律和整體質(zhì)量守恒，可以確立擴(kuò)散的模型方程為

式中，Dx、Dy、Dz分別為x、y、z方向的擴(kuò)散系數(shù)，m2·s?1；wx、wy、wz分別為x、y、z方向的顆粒整體對流平均速度，m·s?1；γs為反應(yīng)速率，kg·m?3。

依據(jù)特定條件對模型進(jìn)行簡化，根據(jù)初始條件和邊界條件求解擴(kuò)散方程，可以獲得橫向擴(kuò)散系數(shù)。

橫向擴(kuò)散系數(shù)也可以基于單顆粒軌跡求取[8,17]。假如當(dāng)t=0時(shí)，在（x0,y0,z0）處注入N個(gè)粒子，每個(gè)粒子的瞬時(shí)位移可表示為

示蹤顆粒在t時(shí)刻的橫向平均位移為

t時(shí)刻顆粒橫向位移的方差為

顆粒的橫向擴(kuò)散系數(shù)可以表示為

本文模擬中，顆粒相采用Lagrange方法，因此密相區(qū)顆粒橫向擴(kuò)散系數(shù)采用基于顆粒軌跡的方法求解。

3 模型及模擬條件

本文主要是研究大顆粒在密相區(qū)擴(kuò)散行為，所關(guān)心區(qū)域主要是提升管底部，而CFB鍋爐內(nèi)顆粒的粒徑分布范圍較大，較大的顆粒在密相區(qū)中無法被攜帶到爐膛上部，與上部快速流態(tài)化不同，爐膛底部密相區(qū)為鼓泡流化狀態(tài)，因此鼓泡床內(nèi)的顆粒擴(kuò)散行為可以定性地反映循環(huán)流床內(nèi)顆粒的混合[18-19]，因此，本文的模擬對象實(shí)為鼓泡流化床。文中采用準(zhǔn)三維實(shí)體模型，模擬對象的計(jì)算域?yàn)?900 mm(長)×100 mm(寬)×1200 mm(高)。劃分網(wǎng)格時(shí)，考慮到當(dāng)網(wǎng)格尺寸較小時(shí)，所需計(jì)算時(shí)間較長；當(dāng)網(wǎng)格尺寸較大時(shí)，模擬精度達(dá)不到要求。綜合考慮計(jì)算效率和精度，采用網(wǎng)格尺寸為10 mm中等精度結(jié)構(gòu)化網(wǎng)格。

計(jì)算中，床料顆粒粒徑采用600 μm。顆粒均按球形度0.66、密度2600 kg·m?3考慮。示蹤顆粒直徑dp分別采用0.6、1、1.5、3、5 mm，除粒徑外示蹤顆粒與床料其他性質(zhì)相同。流化風(fēng)速依據(jù)流化特點(diǎn)選擇為 0.75、1.00、1.25、1.60 和 1.75 m·s?1。流體黏性系數(shù)為 1.8×10?5kg·(m?s)?1，密度為 1.17 kg·m?3。顆粒間碰撞恢復(fù)系數(shù)為0.4，正向碰撞恢復(fù)系數(shù)為0.3，壁面切向碰撞恢復(fù)系數(shù)為0.99。靜止床高為 0.4 m，初始狀態(tài)由重力作用下顆粒自由堆積得到，堆積孔隙率約為0.42，對于不同粒徑顆粒的接觸應(yīng)力采用混合加速度模型。氣固兩相間作用，當(dāng)顆粒體積分?jǐn)?shù)εp＞0.85εcp時(shí)采用Ergun公式，當(dāng)顆粒體積分?jǐn)?shù)εp＜0.75εcp時(shí)采用Wen-Yu公式，當(dāng)顆粒體積分?jǐn)?shù) 0.75εcp＜εp＜0.85εcp時(shí)，按二者線性計(jì)算。氣體對壁面的條件為無滑移，顆粒對壁面的條件為部分滑移。入口采用均勻布風(fēng)方式。模擬時(shí)間為50 s，時(shí)間步長由軟件中CFL（0.8～1.5）模塊自動(dòng)確定，約為0.00025 s。

示蹤法是大多數(shù)學(xué)者研究密相區(qū)中顆粒擴(kuò)散行為所廣泛采用的一種方法[3,16-17,20-21]，該方法簡單且易于實(shí)現(xiàn)，能夠很好地捕捉顆粒具體運(yùn)動(dòng)行為。本文模擬中采用示蹤粒子法，注入方式為脈沖注入，注入時(shí)間為模擬開始5 s后。

4 模擬結(jié)果與討論

4.1 模型驗(yàn)證

為了驗(yàn)證模擬的可信度，利用CPFD方法對文獻(xiàn)[3,22]中實(shí)驗(yàn)結(jié)果分別進(jìn)行了預(yù)測，并將實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了比較，見圖1。

圖1 模擬結(jié)果與實(shí)驗(yàn)結(jié)果對比Fig.1 Comparison of simulation results with experimental data

文獻(xiàn)[3]實(shí)驗(yàn)條件為長900 mm、寬100 mm、高5200 mm的冷態(tài)流化床，床料粒徑分別為600、1000μm，靜止床高17.5 cm，流化風(fēng)速2.2、2.5、3.0、3.7、4.1 m·s?1。模擬結(jié)果與實(shí)驗(yàn)結(jié)果對比見圖1(a)。文獻(xiàn)[22]實(shí)驗(yàn)條件為長600 mm、寬50 mm、高350 mm 的冷態(tài)流化床，床料粒徑 491 μm，靜止床高3.05、4.356、5.23 cm，流化風(fēng)速 0.87 m·s?1[18]。模擬結(jié)果與實(shí)驗(yàn)結(jié)果對比見圖1(b)。從圖1中可以看出，CPFD方法模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù)的變化趨勢完全一致，量級相同，模擬結(jié)果與實(shí)驗(yàn)結(jié)果符合較好。表明用CPFD預(yù)測密相區(qū)顆粒橫向擴(kuò)散行為是可信的。

4.2 流動(dòng)形態(tài)分析

為了便于觀察氣泡行為，選取ug=0.75 m·s?1的工況，來觀察床內(nèi)流動(dòng)結(jié)構(gòu)的發(fā)展?fàn)顩r，見圖2。整個(gè)床層由顆粒濃度很稀甚至不含顆粒的氣泡相和顆粒濃度較高的乳化相組成。當(dāng)風(fēng)速超過最小流化風(fēng)速umf時(shí)，多余的氣體將會(huì)以氣泡的形式穿過床層，氣泡由底部布風(fēng)板處開始產(chǎn)生，床層也隨之膨脹，剛產(chǎn)生時(shí)氣泡體積較小，隨著氣泡的上升，氣泡周圍環(huán)境壓力逐漸變小，體積逐漸變大，到達(dá)床面后便爆裂開，這是由于到達(dá)床面后氣泡失穩(wěn)。氣泡上升過程中，伴有合并、破碎等現(xiàn)象。氣泡穿過床層、氣泡爆裂以及氣泡之間的合并、破碎都會(huì)引起顆粒的混合，與文獻(xiàn)[15,23-24]中觀察到的現(xiàn)象一致。

圖2 床內(nèi)不同時(shí)刻流動(dòng)形態(tài)（ug=0.75 m·s?1）Fig.2 Flow patterns at different time in bed

圖3(a)為ug=0.75 m·s?1時(shí)，最初 5 s內(nèi)床層膨脹比與時(shí)間的關(guān)系。床層膨脹比為不同條件下床層高度與靜止床高之間的比值。觀察圖 3(a)發(fā)現(xiàn)，模擬開始1 s內(nèi)，床層膨脹比隨時(shí)間增加而增大，兩者之間接近于線性關(guān)系，超過1 s后，有所減小，達(dá)到1.5 s后基本趨于穩(wěn)定。表明當(dāng)t=1.5 s時(shí)，床內(nèi)氣固兩相的流動(dòng)結(jié)構(gòu)已發(fā)展完全。

圖3 床層膨脹比Fig.3 Bed expansion ratio

圖4是不同流化風(fēng)速下，t=5 s時(shí)單一粒徑床內(nèi)流動(dòng)形態(tài)。當(dāng)流化風(fēng)速較小時(shí)，床內(nèi)氣泡體積小，合并現(xiàn)象少，到達(dá)床面爆裂后的拋射能力弱；隨著流化風(fēng)速的增加，氣泡的尺度變大，合并現(xiàn)象變多，床層膨脹度也增大，氣泡在床面爆裂后的拋射能力更強(qiáng)；當(dāng)流化風(fēng)速較大時(shí)，床層膨脹度變化不大，氣泡尺寸和數(shù)量基本不變，局部位置上的氣泡相連，部分氣體短路穿過床層。圖3(b)為床層膨脹比與流化風(fēng)速的關(guān)系。由圖3(b)發(fā)現(xiàn)，床層膨脹比先隨著流化風(fēng)速的增加變大，超過一定風(fēng)速后，床層膨脹比增加幅度變小，到最后基本不變。與圖4中觀察到的現(xiàn)象一致。

圖4 床內(nèi)不同風(fēng)速下流動(dòng)形態(tài)（t=5 s）Fig.4 Flow patterns under different fluidization velocities in bed

4.3 橫向擴(kuò)散系數(shù)求解

圖5為示蹤顆粒橫向位移方差X2與時(shí)間t的典型關(guān)系。在示蹤顆粒5 s注入后至12 s之前，其斜率基本不變，當(dāng)12 s之后，因?yàn)槟Ｐ统叽缬邢蓿聚櫫Ｗ訖M向位移受到限制，不再正確反映顆粒擴(kuò)散行為。本文顆粒橫向擴(kuò)散系數(shù)Dx的計(jì)算范圍選取-t曲線前面斜率相近段，即示蹤粒子注入之后的7 s范圍內(nèi)的擴(kuò)散行為?？紤]到流態(tài)化的床內(nèi)波動(dòng)性，本文針對同一工況選取不同時(shí)刻，分別求橫向擴(kuò)散系數(shù)，最后求得平均值。

圖5 橫向位移方差與時(shí)間的關(guān)系Fig.5 Averaged square of lateral displacement with time

4.3.1 示蹤顆粒數(shù)目對橫向擴(kuò)散系數(shù)的影響 考慮到示蹤顆粒數(shù)目不能太多，太多會(huì)影響床內(nèi)氣固流動(dòng)行為；由于顆粒運(yùn)動(dòng)的隨機(jī)性，示蹤顆粒不能太少，否則模擬結(jié)果不具有重復(fù)性和代表性。模擬中示蹤粒子數(shù)目是通過質(zhì)量流量及注入時(shí)間來控制，本文示蹤顆粒數(shù)目在102～103數(shù)量級上，各工況下示蹤粒子質(zhì)量不超過床層總質(zhì)量的 2%，示蹤粒子的加入對床層影響很小。

為了解示蹤顆粒數(shù)目對橫向擴(kuò)散系數(shù)的影響，模擬中分析了ug=1.5 m·s?1及3種不同示蹤顆粒粒徑條件下示蹤顆粒數(shù)目對顆粒橫向擴(kuò)散系數(shù)Dx的影響,模擬結(jié)果如表1所示。相同示蹤顆粒粒徑dp條件下，由不同顆粒數(shù)量模擬得到的顆粒橫向擴(kuò)散系數(shù)Dx誤差在3%以內(nèi)。表明在本文所采用的示蹤顆粒數(shù)目數(shù)量級范圍內(nèi)，示蹤顆粒數(shù)目對顆粒橫向擴(kuò)散系數(shù)Dx影響較小，基本可以忽略不計(jì)。

表1 示蹤顆粒數(shù)目相關(guān)性分析Table 1 Correlation analysis of tracer particles number

4.3.2 橫向擴(kuò)散系數(shù)與風(fēng)速的關(guān)系 橫向擴(kuò)散系數(shù)Dx與流化風(fēng)速ug的關(guān)系如圖6所示。

由圖6可見，隨著流化風(fēng)速ug的增加，顆粒的橫向擴(kuò)散系數(shù)Dx增大，曲線的斜率先變大后變小，不同粒徑間橫向擴(kuò)散系數(shù)Dx差距變大。依據(jù)對床內(nèi)流動(dòng)形態(tài)分析知道，當(dāng)流化風(fēng)速增加時(shí)，氣泡體積變大，氣泡在上升、合并和爆裂時(shí)所引起的顆粒混合更加強(qiáng)烈，增強(qiáng)了顆粒橫向位移，從而增大了顆粒橫向擴(kuò)散系數(shù)Dx。受限于床體尺寸，最大氣泡尺寸存在上限，部分氣體以短路的方式穿過床層，導(dǎo)致橫向擴(kuò)散系數(shù)Dx增長變慢。

4.3.3 橫向擴(kuò)散系數(shù)與粒徑的關(guān)系 顆粒橫向擴(kuò)散系數(shù)Dx與顆粒粒徑dp之間的關(guān)系見圖7。如圖7所示，隨著示蹤顆粒粒徑dp的增大，顆粒橫向擴(kuò)散系數(shù)Dx減小，這是由于當(dāng)粒徑dp增大時(shí)，體積增大，顆粒向前運(yùn)動(dòng)將要克服更大的阻力，而且當(dāng)密度相同時(shí)，體積增大后，單個(gè)顆粒的質(zhì)量增加，使其運(yùn)動(dòng)將需要更大的動(dòng)力。綜合來說，當(dāng)粒徑dp增大其擴(kuò)散將會(huì)受到抑制，其橫向擴(kuò)散系數(shù)減小。在圖7的曲線族中，從下往上流化風(fēng)速依次增加。當(dāng)流化風(fēng)速ug較小時(shí)，不同粒徑顆粒得到的橫向擴(kuò)散系數(shù)相近，低流化風(fēng)速下，顆粒粒徑dp對顆粒橫向擴(kuò)散系數(shù)Dx影響較??；當(dāng)流化風(fēng)速ug增加時(shí)，顆粒橫向擴(kuò)散系數(shù)Dx曲線先下降較快，dp超過1500 μm后，變得較為平緩。表明流化風(fēng)速ug增大時(shí)，顆粒粒徑dp對顆粒橫向擴(kuò)散系數(shù)Dx影響增強(qiáng)，超過一定粒徑后，顆粒橫向擴(kuò)散系數(shù)Dx基本不變。

圖7 橫向擴(kuò)散系數(shù)與粒徑的關(guān)系Fig.7 Relation between lateral dispersion coefficient and tracer particle diameter

4.3.4 橫向擴(kuò)散系數(shù)對比 研究結(jié)果表明[8,25]，顆粒橫向擴(kuò)散系數(shù)Dx實(shí)驗(yàn)值大致在 0.0001～0.1 m2·s?1范圍內(nèi)，存在 3個(gè)數(shù)量級以上的差距，這是由于不同學(xué)者實(shí)驗(yàn)條件和測試手段差別很大。當(dāng)流化風(fēng)速較低或床體尺寸較小時(shí)[26-28]，顆粒橫向擴(kuò)散系數(shù)Dx處于 0.0001～0.001 m2·s?1數(shù)量級；當(dāng)流化風(fēng)速較高或床體尺寸較大時(shí)[16,29-30]，顆粒橫向擴(kuò)散系數(shù)Dx處于 0.01～0.1 m2·s?1數(shù)量級。本文模擬床體橫截面尺寸為0.9 m×0.1 m，尺寸較小，流化風(fēng)速為0.75～2 m·s?1，流化風(fēng)速較低，模擬所得顆粒橫向擴(kuò)散系數(shù)Dx為 0.0002～0.003 m2·s?1，與前人實(shí)驗(yàn)結(jié)果數(shù)量級上符合較好。進(jìn)一步說明運(yùn)用CPFD方法來模擬實(shí)驗(yàn)室規(guī)模的循環(huán)流化床密相區(qū)中氣固流動(dòng)、顆粒的擴(kuò)散行為是可行的，這種方法預(yù)測密相區(qū)中顆粒運(yùn)動(dòng)行為具有明顯的優(yōu)勢，其在高濃度氣固兩相流數(shù)值模擬的研究中具有廣闊的應(yīng)用前景。

5 結(jié) 論

本文運(yùn)用CPFD方法對循環(huán)流化床密相區(qū)大顆粒的擴(kuò)散行為進(jìn)行研究，采用示蹤顆粒方法得到大顆粒的運(yùn)動(dòng)軌跡，并根據(jù)單顆粒位移計(jì)算得到顆粒平均橫向擴(kuò)散系數(shù)Dx。主要結(jié)論如下。

（1）利用CPFD預(yù)測顆粒橫向擴(kuò)散系數(shù)與實(shí)驗(yàn)結(jié)果比較吻合，表明CPFD對循環(huán)流化床密相區(qū)模擬，可以反映床內(nèi)氣固流動(dòng)結(jié)構(gòu)的演化過程，以及氣泡的產(chǎn)生、上升以及爆裂對密相區(qū)中顆?；旌系淖饔?。

（2）文中分析了不同條件下床層膨脹比的變化規(guī)律，隨著時(shí)間的增加，床層膨脹比先增加后減小，最后趨于不變；隨著流化風(fēng)速的增加，床層膨脹比變大，超過一定風(fēng)速后，其增加幅度變小，到最后基本不變。

（3）通過流化風(fēng)速ug、顆粒粒徑dp對顆粒橫向擴(kuò)散系數(shù)Dx的影響的CPFD模擬發(fā)現(xiàn)，隨流化風(fēng)速ug增加，顆粒橫向擴(kuò)散系數(shù)Dx增大；隨顆粒粒徑dp增加，顆粒橫向擴(kuò)散系數(shù)Dx減小，當(dāng)流化風(fēng)速ug較小時(shí)，顆粒粒徑dp變化對顆粒橫向擴(kuò)散系數(shù)Dx影響較小，當(dāng)流化風(fēng)速ug較大時(shí)，顆粒粒徑dp對顆粒橫向擴(kuò)散系數(shù)Dx影響變大。

（4）根據(jù)顆粒軌跡計(jì)算，密相區(qū)顆粒橫向擴(kuò)散系數(shù)Dx為 0.0002～0.003 m2·s?1，不同學(xué)者通過實(shí)驗(yàn)得到的橫向擴(kuò)散系數(shù)Dx分布在 0.0001～0.1 m2·s?1，模擬結(jié)果落在實(shí)驗(yàn)范圍之內(nèi)，表明CPFD在高濃度氣固兩相流數(shù)值模擬的研究中具有廣闊的應(yīng)用前景。

[1]劉道銀,陳曉平,唐智,等.側(cè)面進(jìn)料在循環(huán)流化床密相區(qū)混合特性的試驗(yàn)研究[J].工程熱物理學(xué)報(bào),2009,30(3): 529-532.LIU D Y,CHEN X P,TANG Z,et al.Experimental study on the mixing of particles feeding into the bottom zone of a CFB[J].Journal of Engineering Thermophysics,2009,30(3): 529-532.

[2]肖顯斌,楊海瑞,呂俊復(fù),等.CFB密相區(qū)內(nèi)顆粒橫向擴(kuò)散對燃燒的影響[J].煤炭轉(zhuǎn)化,2003,26(4): 55-59.XIAO X B,YANG H R,Lü J F,et al.Effect of the solid lateral dispersion on combustion in the dense zone of CFB[J].Coal Conversion,2003,26(4): 55-59.

[3]楊海瑞,呂俊復(fù),劉青,等.循環(huán)流化床鍋爐密相區(qū)內(nèi)顆粒的橫向擴(kuò)散研究[J].熱能動(dòng)力工程,2001,16(4): 395-399.YANG H R,Lü J F,LIU Q,et al.Study on the solid lateral dispersion in the dense zone of circulating fluidized bed[J].Thermal and Power Engineering,2001,16(4): 395-399.

[4]劉寶勇,魏緒玲,張斌,等.氣-固循環(huán)流化床底部密相區(qū)研究進(jìn)展[J].化工技術(shù)與開發(fā),2014,43(3): 38-42.LIU B Y,WEI X L,ZHANG B,et al.Research development of bottom dense phase zone of gas-solid circulating fluidized bed[J].Technology & Development of Chemical Industry,2014,43(3):38-42.

[5]王翀.流化床氣固兩相流數(shù)值模擬及實(shí)驗(yàn)研究[D].重慶: 重慶大學(xué),2011.WANG C.Numerical simulation and experimental for flow characteristics in fluidized-bed[D].Chongqing: Chongqing University,2011.

[6]汪琦.氣固流化床兩相流動(dòng)的 CFD模型研究和實(shí)驗(yàn)驗(yàn)證[D].武漢: 華中科技大學(xué),2012.WANG Q.Comparative analysis of CFD models of gas-solid fluidized bed and experimental verification[D].Wuhan: Huazhong University of Science and Technology,2012.

[7]閆盛楠.鼓泡流化床不規(guī)則形狀顆粒氣固兩相流動(dòng)特性研究[D].哈爾濱: 哈爾濱工業(yè)大學(xué),2014.YAN S N.Investigation on gas-solid two-phase flow characteristics of irregularly shaped particle in a bubbling fluidized bed[D].Harbin:Harbin Institute of Technology,2014.

[8]劉道銀,陳曉平,陸利燁,等.流化床密相區(qū)顆粒擴(kuò)散系數(shù)的CFD數(shù)值預(yù)測[J].化工學(xué)報(bào),2009,60(9): 2183-2190.LIU D Y,CHEN X P,LU L Y,et al.Prediction of solids dispersion coefficient in fluidized bed dense zone using CFD simulation[J].CIESC Journal,2009,60(9): 2183-2190.

[9]OKE O,LETTIERE P,SALATINO P,et al.Eulerian modeling of lateral solid mixing in gas-fluidized suspensions[J].Procedia Engineering,2015,102: 1491-1499.

[10]SNIDER D M,O'ROURKE P J,ANDREWS M J.Sediment flow in inclined vessels calculated using a multiphase particle-in-cell model for dense particle flows[J].International Journal of Multiphase Flow,1998,24(8): 1359-1382.

[11]SNIDER D M.An incompressible three-dimensional multiphase particle-in-cell model for dense particle flows[J].Journal of Computational Physics,2001,170(2): 523-549.

[12]SHI X,LAN X,LIU F,et al.Effect of particle size distribution on hydrodynamics and solids back-mixing in CFB risers using CPFD simulation[J].Powder Technology,2014,266: 135-143.

[13]CHEN C,WERTHER J,HEINRICH S,et al.CPFD simulation of circulating fluidized bed risers[J].Powder Technology,2013,235(2):238-247.

[14]XU B H,YU A B.Numerical simulation of the gas-solid flow in a fluidized bed by combining discrete particle method with computational fluid dynamics[J].Chemical Engineering Science,1997,52(16): 2785-2809.

[15]沈來宏,章名耀.鼓泡流化床內(nèi)顆?；旌系膶α?擴(kuò)散模型[J].中國電機(jī)工程學(xué)報(bào),1995,15(1): 45-53.SHEN L H,ZHANG M Y.Modeling for solids mixing in bubbling fluidized beds[J].Proceeding of The Chinese Society for Electrical Engineering,1995,15(1): 45-53.

[16]李炳順,孫運(yùn)凱,雍玉梅,等.二維循環(huán)流化床密相區(qū)物料擴(kuò)散冷態(tài)試驗(yàn)與模型建立[J].鍋爐技術(shù),2007,38(2): 28-31.LI B S,SUN Y K,YONG Y M,et al.Coal-feeding diffusion cold experiment and model-building in dense phase of two-dimensional circulating fluidized bed[J].Boiler Technology,2007,38(2): 28-31.

[17]MOSTOUFI N,CHAOUKI J.Local solid mixing in gas-solid fluidized beds[J].Powder Technology,2001,114(1/2/3): 23-31.

[18]YUE G X,LU J F,ZHANG H,et al.Design theory of circulating fluidized bed boilers[C]//18th International Conference on Fluidized Bed Combustion.Fairfield,NJ: American Society of Mechanical Engineers,2005: 135-146.

[19]楊石,楊海瑞,呂俊復(fù),等.新一代節(jié)能型循環(huán)流化床鍋爐燃燒技術(shù)[J].動(dòng)力工程學(xué)報(bào),2009,29(8): 728-732.YANG S,YANG H R,Lü J F,et al.The new generation combustion technology for energy saving circulating fluidized bed boilers[J].Journal of Power Engineering,2009,29(8): 728-732.

[20]RHODES M J,ZHOU S,HIRAMA T,et al.Effects of operating conditions on longitudinal solids mixing in a circulating fluidized bed riser[J].AIChE Journal,1991,37(10): 1450-1458.

[21]KOJIMA T,ISHIHARA K I,GUILIN Y,et al.Measurement of solids behaviour in a fast fluidized bed[J].Journal of Chemical Engineering of Japan,1989,22(4): 341-346.

[22]SHI Y F,FAN L T.Lateral mixing of solids in batch gas-solids fluidized beds[J].Industrial & Engineering Chemistry Process Design and Development,1984,23(2): 337-341.

[23]金涌.流態(tài)化工程原理[M].北京: 清華大學(xué)出版社,2002.JIN Y.Theory of Fluidized Engineering[M].Beijing: Tsinghua University Press,2002.

[24]徐旭常.燃燒理論與燃燒設(shè)備[M].北京: 機(jī)械工業(yè)出版社,1990.XU X C.Combustion Theory and Combustion Equipment[M].Beijing: China Machine Press,1990.

[25]胡南,郭兆君,楊海瑞,等.CFB鍋爐爐膛內(nèi)顆粒橫向擴(kuò)散系數(shù)研究[J].動(dòng)力工程學(xué)報(bào),2016,36(3): 168-171.HU N,GUO Z J,YANG H R,et al.Experimental study on lateral dispersion coefficient of solid particles in a CFB boiler[J].Chinese Journal of Power Engineering,2016,36(3): 168-171.

[26]SALAM T F,REN Y,GIBBS B M.Lateral solid and thermal dispersion in fluidized bed combustors[C]//9th InternationalConference on Fluidized Bed Combustion.New York: American Society of Mechanical Engineers,1987: 541-545.

[27]WINAYA I N S,SHIMIZU T,YAMADA D.A new method to evaluate horizontal solid dispersion in a bubbling fluidized bed[J].Powder Technology,2007,178(3): 173-178.

[28]BERRUTI F,SCOTT D S,RHODES E.Measuring and modelling lateral solid mixing in a three-dimensional batch gas-solid fluidized bed reactor[J].Canadian Journal of Chemical Engineering,2010,64(1): 48-56.

[29]XIAO P,YAN G,WANG D.Investigation on horizontal mixing of particles in dense bed in circulating fluidized bed (CFB)[J].Journal of Thermal Science,1998,7(2): 78-84.

[30]SCHLICHTHAERLE P,WERTHER J.Solids mixing in the bottom zone of a circulating fluidized bed[J].Powder Technology,2001,120(1): 21-33.

Lateral dispersion coefficient of large solid particles in CFB dense zone using CPFD method

ZHANG Xian1,2,GE Rongcun1,2,ZHANG Shouyu1,LIU Qing2,ZHANG Man2,YANG Hairui2,Lü Junfu2
(1School of Energy and Power Engineering,University of Shanghai for Science and Technology,Shanghai200093,China;2Department of Thermal Engineering,Tsinghua University,Beijing100084,China)

TK 224.1

A

0438—1157（2017）10—3725—08

10.11949/j.issn.0438-1157.20170394

2017-04-13 收到初稿，2017-05-19收到修改稿。

聯(lián)系人：呂俊復(fù)。

張賢（1990—），男，碩士研究生。

國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目（2016YFB0600201）。

Received date:2017-04-13.

Corresponding author:Lü Junfu,lvjf@mail.tsinghua.edu.cn

Foundation item:supported by the National Key Research and Development Program (2016YFB0600201).