氣體分布板對(duì)制粒流化床流化質(zhì)量的影響

劉訓(xùn)濤, 劉福燦, 沙宇峰, 劉金輝

(1.黑龍江科技大學(xué) 機(jī)械工程學(xué)院, 哈爾濱 150022; 2.哈爾濱納諾機(jī)械設(shè)備有限公司, 哈爾濱 150078)

0 引 言

在國(guó)內(nèi)能源轉(zhuǎn)型和碳中和背景下,流化床在能源、化工、制藥、污水處理等領(lǐng)域的應(yīng)用得到進(jìn)一步拓展,通過流化床可以讓固體顆粒呈現(xiàn)流化狀態(tài),促進(jìn)顆粒間的有效混合[1-4]。尤其在工業(yè)制粒過程中,流化床發(fā)揮著越來越重要的作用。制粒過程中流化床內(nèi)部的氣固兩相返混劇烈,傳熱、傳質(zhì)效率較高,流化后混合顆粒的顆粒粒度均勻,流動(dòng)性、壓縮成形性好。氣體分布板作為制粒流化床的重要組成部件之一,既起著支撐物料的作用,又承擔(dān)著均勻分布?xì)怏w,調(diào)整氣流方向、大小的作用,能夠改善流化床內(nèi)氣體分布,創(chuàng)造一個(gè)良好的起始流化條件,形成一個(gè)好的氣固接觸環(huán)境。朱沈瑾等[5]針對(duì)錐形非均勻分布板流化床進(jìn)行數(shù)值模擬計(jì)算。Li等[6]設(shè)計(jì)了一個(gè)帶有一個(gè)中央進(jìn)氣口和六個(gè)輔助進(jìn)氣口的錐形分配器,即使在較低的表觀氣速下,死區(qū)現(xiàn)象也得到了很好地抑制,有效地改善了流態(tài)化質(zhì)量。魏新利等[7]利用Fluent數(shù)值模擬軟件對(duì)流化床進(jìn)行流場(chǎng)模擬,重點(diǎn)討論布風(fēng)板正方形和正三角形布孔排列形式和三種開孔直徑下對(duì)流態(tài)的影響。Li等[8]采用減速法研究了空氣分布器結(jié)構(gòu)對(duì)片狀顆粒流化特性的綜合影響。李少華等[9]應(yīng)用CFD軟件Fluent對(duì)流化床爐內(nèi)布風(fēng)板進(jìn)行流場(chǎng)模擬,研究在一定的開孔率下不同的布風(fēng)板開孔布置與不同開孔直徑下對(duì)流化床冷態(tài)流場(chǎng)與布風(fēng)板壓降的影響。

雖然專家學(xué)者對(duì)流化床分布板進(jìn)行了大量的研究,但是實(shí)際工作中,仍然存在一定量由氣體分布板引起的顆粒循環(huán)流動(dòng)的均勻性差,產(chǎn)生偏流等問題。在壁面處,由于物料黏性和氣體分布不均,還會(huì)存在一定的渦旋和循環(huán)死區(qū),導(dǎo)致顆粒流化效果差,制粒效率低,筆者應(yīng)用CFD-DEM方法探究分布板均勻分布和非均勻分布排列形式和三種開孔直徑下對(duì)流化床內(nèi)部流場(chǎng)的影響,得到流化過程顆粒和氣體的動(dòng)態(tài)特性。

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

由于流化床內(nèi)氣固兩相流動(dòng)的復(fù)雜性,數(shù)值模擬中通常采用顆粒軌道模型或雙流體模型進(jìn)行計(jì)算[10-11]。其中,雙流體模型又稱歐拉-歐拉法(Euler-Euler),著眼于空間,以流場(chǎng)中每一空間位置為描述對(duì)象,描述這些位置上流體的物理屬性隨時(shí)間的變化狀態(tài)。因此,被廣泛應(yīng)用于研究流化床內(nèi)氣相和固相流動(dòng)的復(fù)雜的流體力學(xué)行為。采用雙流體模型仿真模擬三種氣體分布板下流化床內(nèi)部的氣固兩相流動(dòng)狀態(tài),其質(zhì)量守恒方程和動(dòng)量守恒方程[12]如下。

當(dāng)氣相時(shí),質(zhì)量守恒方程為

(1)

當(dāng)固相時(shí),質(zhì)量守恒方程為

(2)

式中:εg——?dú)怏w體積分?jǐn)?shù),%;

ρg——?dú)怏w密度,kg/m3;

vg——?dú)怏w速度,m/s;

εs——顆粒體積分?jǐn)?shù);

ρs——顆粒密度,kg/m3;

vs——顆粒速度,m/s。

動(dòng)量方程基于動(dòng)量守恒原理,氣相與固相的動(dòng)量守恒方程進(jìn)行構(gòu)建。

當(dāng)氣相時(shí)

β(vs-vg)+?εgτg+ρgεgg。

(3)

當(dāng)固相時(shí)

εs?τs+β(vg-vs)+ρsεsg,

(4)

式中:vg——?dú)怏w速度向量;

vs——固體速度向量;

τg——?dú)怏w剪應(yīng)力向量;

τs——固體剪應(yīng)力向量;

g——重力加速度向量;

β——相間摩擦系數(shù)(氣固相間曳力系數(shù));

?ps——固體壓力梯度;

?pg——?dú)怏w壓力梯度。

顆粒接觸力模型選用Hertz-Mindlin無滑移接觸模型[13],該模型可精確得出顆粒在彈性碰撞下所受的剪切應(yīng)力,采用軟球線性碰撞模型,將碰撞過程簡(jiǎn)化,用一個(gè)彈簧、滑移器和阻尼器的組合模型來描述,如圖1所示。

圖1 流化顆粒線彈性碰撞模型

當(dāng)兩個(gè)顆粒碰撞時(shí),其法向疊合量為

?=R1+R2-|r1-r2|,

(5)

式中:R1、R2——球形顆粒半徑,mm;

r1、r2——球心位置矢量。

合力和合力矩通過顆粒碰撞過程中切向接觸力和法向接觸力的疊加進(jìn)行計(jì)算。顆粒與壁面的接觸與顆粒間的碰撞過程相同,因此在計(jì)算顆粒與壁面之間的作用力時(shí)將壁面假設(shè)置成顆粒。顆粒間的法向接觸力為Fn,其表達(dá)式為

式中:E*——等效彈性模量,MPa;

E1、E2——顆粒彈性模量,MPa;

R*——等效顆粒半徑,mm;

μ1、μ2——顆粒泊松比。

顆粒碰撞時(shí)的切向接觸力Ft,其表達(dá)為

Ft=-Stδ,

式中:St——切向剛度,N/m;

δ——切向疊合量,mm;

G1、G2——顆粒的彈性剪切模量,MPa。

2 氣體分布板結(jié)構(gòu)對(duì)流化效果的影響

2.1 仿真初始參數(shù)的設(shè)置

運(yùn)用Solidworks軟件建立流化床三維模型,采用六面體網(wǎng)格對(duì)模型進(jìn)行劃分,如圖2所示。為了避免其他因素對(duì)氣固流化床流場(chǎng)計(jì)算結(jié)果造成影響,保證生成網(wǎng)格的質(zhì)量,將網(wǎng)格增長(zhǎng)率設(shè)置為1.1,網(wǎng)格尺寸為71.42 mm。劃分后的網(wǎng)格的單元數(shù)量為1 392 892,節(jié)點(diǎn)數(shù)為258 178。

圖2 流化床模型尺寸

在工程實(shí)際中,流化床內(nèi)部填充的固體顆粒往往是不規(guī)則形狀,而且顆粒的直徑大小也不相同。為了簡(jiǎn)化計(jì)算,假設(shè)模擬過程中采用的固體顆粒為理想的單顆粒球體,粒徑均勻一致。流化氣體為常溫、常壓下空氣,氣流沿進(jìn)氣面均勻分布,密度保持一致,不會(huì)隨著壓強(qiáng)的變化而改變,為不可壓縮流體。模擬的過程中氣體速度呈線性增長(zhǎng),初始時(shí)刻氣體的速度為0 m/s,0.5 s后氣體速度趨于穩(wěn)定,穩(wěn)定狀態(tài)下氣體速度為10 m/s。氣體速度的函數(shù)表達(dá)為

IF(t≤0.5[s],t/1[s]*-10[m/s],-5[m/s])。

流化床初始床層填充高度為16.5 mm,填充顆粒的直徑為1.4 mm,顆粒初始速度為0 m/s,床層空隙率為46%。設(shè)定求解器的步長(zhǎng)為0.001,步數(shù)為3 000,每一步最大迭代次數(shù)20。參數(shù)設(shè)置見表1。

表1 流化床參數(shù)

為了研究氣體分布板開孔參數(shù)對(duì)制粒流化效果的影響,在開孔率10.15%的條件下給出了3-5-9(A型)、5-5-9(B型)、5-3-9(C型)三種氣體分布板,氣體分布板的開孔尺寸,A、B、C型分布板的內(nèi)圈直徑分別為3、5、5 mm,外圈直徑分別為5、5、3 mm,孔間距均為9 mm。不同分布板開孔示意如圖3所示。采用單一變量方法來進(jìn)行模擬分析,通過保持入口進(jìn)氣速度、顆粒屬性和顆粒初始填充高度不變,研究不同氣體分布板流化床的流場(chǎng)特征和顆粒體積分?jǐn)?shù)分布規(guī)律。

圖3 不同分布板開孔方式示意

2.2 流化床氣固兩相流瞬態(tài)值模擬

采用EDEM軟件與Fluent軟件進(jìn)行聯(lián)合仿真,在Fluent中通過求解流體力學(xué)方程,可以獲得流化床中氣體速度和壓力等參數(shù),在EDEM中通過離散單元法對(duì)顆粒的運(yùn)動(dòng)與碰撞進(jìn)行計(jì)算,可以獲得顆粒間的相互作用力、顆粒的位置和速度等參數(shù)。

2.2.1 不同分布板下顆粒分布變化情況

在相同送風(fēng)條件、相同初始床層高度、不同型號(hào)的氣體分布板下,流化床中心截面(y=0)固體顆粒由啟動(dòng)到穩(wěn)定工況時(shí)的分布情況,如圖4所示。

圖4 不同分布板下顆粒隨時(shí)間變化情況

由圖4a、b可見,在氣體速度尚未達(dá)到穩(wěn)定狀態(tài)時(shí),t=0.5 s時(shí),分布板A與分布板B中的顆粒在氣體的作用下逐步進(jìn)入湍流狀態(tài),流化床層高度較為均勻,床層表面跳動(dòng)幅度不大,顆粒均勻分布啟動(dòng)階段分布板A中流化床層表面高度高于分布板B。分布板C中心部位的顆粒率先涌出床層,靠近內(nèi)壁附近的顆粒響應(yīng)較慢。

在流化過程中,2.0 s以后分布板A中的顆粒開始周期性的運(yùn)動(dòng),分布板B與分布板C在1.5 s時(shí)顆粒達(dá)到穩(wěn)定狀態(tài),相比于分布板B與分布板C,分布板A中顆粒到達(dá)穩(wěn)定狀態(tài)的時(shí)間較晚。在流化床穩(wěn)定運(yùn)行時(shí),分布板A與分布板B中的顆粒沿著流化床中心軸圍繞內(nèi)壁做周期性旋轉(zhuǎn)運(yùn)動(dòng),分布板C中的顆粒由中心部位流化至最大高度,然后沿壁面部位回落至床層中,并不斷循環(huán)此過程。

分布板A中的顆粒大多集中在z=0.4 m以下的區(qū)域,與分布板C相比,床層內(nèi)部顆粒比較集中,流化過程中床層高度無法達(dá)到一個(gè)穩(wěn)定值,床層最大高度存在周期性的變化。在一個(gè)變化周期中,床層膨脹比的最大值為9.6,最小值為23.7,可以認(rèn)為此時(shí)顆粒在床層中并未能夠進(jìn)行充分流化,在靠近流化床內(nèi)壁的位置出現(xiàn)明顯的顆粒堆積現(xiàn)象。布設(shè)分布板B時(shí),顆粒集中分布在z=0.4 m以下的區(qū)域,與分布板A相比,在床層膨脹階段,顆粒的膨脹高度較分布板A增加,進(jìn)入穩(wěn)定階段,顆粒同樣存在周期性運(yùn)動(dòng),顆粒的周期性運(yùn)動(dòng)的形式改變,呈現(xiàn)振蕩運(yùn)動(dòng)。布設(shè)分布板C時(shí),顆粒均勻分布在z=0.6 m以下的區(qū)域,顆粒的軸向濃度變化較小,床內(nèi)顆粒無明顯的堆積現(xiàn)象,與分布板A、B相比,床層內(nèi)部顆粒運(yùn)動(dòng)穩(wěn)定,沒有明顯的震蕩現(xiàn)象。

2.2.2 不同氣體分布板下流化床內(nèi)部流場(chǎng)分布

在流化床內(nèi)部顆粒進(jìn)入穩(wěn)定運(yùn)行階段時(shí),選取顆粒一個(gè)循環(huán)運(yùn)動(dòng)周期,分析不同氣體分布板的氣體跡線分布規(guī)律,一個(gè)周期內(nèi)流化床中心截面(y=0)氣體跡線如圖5所示。

圖5 一個(gè)周期內(nèi)流化床中心截面氣體跡線

由圖5可知,布設(shè)氣體分布板A時(shí),氣體軌跡在底部分布比較均勻,隨著高速氣體穿過分布板不斷向上運(yùn)動(dòng),氣體的方向發(fā)生偏轉(zhuǎn),逐漸偏向于流化床右側(cè)氣體出口的位置,氣體的軌跡變得混亂,形成環(huán)流。在z=0.2～0.4 m的位置,氣體軌跡變得復(fù)雜,局部形成環(huán)流區(qū)域,環(huán)流區(qū)域的位置和面積不固定,隨著時(shí)間不斷發(fā)生變化,流場(chǎng)的穩(wěn)定性較差。

氣體分布板的開孔形式為B時(shí),氣體生成的環(huán)流區(qū)域面積較大,環(huán)流區(qū)域出現(xiàn)在較低的位置。布設(shè)分布板C時(shí),氣體在z=0.3～0.1 m的位置,氣體跡線對(duì)稱性較好,軸向速度較大。隨著高度增加,氣體逐漸向四周擴(kuò)散,床層內(nèi)部未出現(xiàn)明顯的環(huán)流,在z=0.7 m的位置,氣體軌跡開始均勻分布,氣體速度降低。氣體分布板的內(nèi)圈直徑擴(kuò)大時(shí),氣體產(chǎn)生的環(huán)流區(qū)域的位置逐漸降低,當(dāng)外圈直徑縮小時(shí),氣體向中間部位集中,有利于改善流化床內(nèi)部的氣體的流動(dòng),有利于加強(qiáng)顆粒在流化床內(nèi)的均勻分布。

2.2.3 不同氣體分布板下固體顆粒動(dòng)力學(xué)特性

流化床內(nèi)中心截面(y=0)載體顆粒軸向速度分布云圖見圖6。規(guī)定沿z軸正向?yàn)檩S向速度的正值。布設(shè)分布板A時(shí),床層內(nèi)部形成偏心環(huán)流,顆粒最大速度為0.825 m/s,最小速度為-0.825 m/s。速度分布在床層內(nèi)部不均勻,右側(cè)大部分區(qū)域?yàn)檎蛩俣?左下角小部分區(qū)域?yàn)樨?fù)值,床層水平速度梯度較大,加劇了顆粒之間的碰撞與摩擦,容易導(dǎo)致產(chǎn)生過粉碎現(xiàn)象,可能會(huì)堵塞流化床頂部濾袋。

圖6 中心截面顆粒軸向速度云圖

不同分布板z軸高度下床層顆粒的徑向速度分布如圖7所示。

圖7 不同床層下氣體分布板對(duì)顆粒徑向速度的影響

布設(shè)分布板B時(shí),床層內(nèi)部環(huán)流的幅值較高,顆粒速度較高的區(qū)域集中在流化床左下角區(qū)域,在右下角的位置顆粒速度較低,且顆粒多向下運(yùn)動(dòng),容易在此處產(chǎn)生顆粒堆積現(xiàn)象。床層內(nèi)部顆粒的最大速度為0.793 m/s,最小速度為-1.598 m/s,顆粒速度沿軸向與徑向變化較小,床層內(nèi)部無過大的速度產(chǎn)生。布設(shè)分布板C時(shí),顆粒在流化床內(nèi)形成明顯的對(duì)稱環(huán)流,形成向上的主流區(qū)和壁面附近的顆粒下落區(qū),主流區(qū)顆粒的最大速度為1.987 m/s,下落區(qū)顆粒的速度為-1.419 m/s。兩側(cè)區(qū)域的流速差有利于主流區(qū)與環(huán)流區(qū)顆粒的相互摻雜和交換,環(huán)流區(qū)內(nèi)的固體顆粒能夠得到及時(shí)地補(bǔ)充,有利于氣體和顆粒之間傳質(zhì)作用的發(fā)生。

從圖7a可以看出,于分布板A上方z=0.05 m處,流化床內(nèi)徑向顆粒速度分布呈現(xiàn)左側(cè)高右側(cè)低的趨勢(shì)。對(duì)于分布板B,流化床內(nèi)徑向顆粒速度分布呈現(xiàn)右側(cè)高左側(cè)低的趨勢(shì)中心區(qū)域速度最低。對(duì)于分布板C,在x=-0.025～0.025 m處,流化床內(nèi)顆粒速度存在一定波動(dòng),在x=-0.025～-0.075 m與x=0.025～0.075 m處,顆粒的速度分布呈現(xiàn)先下降后上升的趨勢(shì),在靠近邊壁位置,顆粒速度為0,整體流化速度的對(duì)稱性較好。

從圖7b可以看出,在z=0.1 m的高度上分布板B顆粒速度顯著降低,顆粒速度在0 ～0.5 m/s之間波動(dòng),整體速度分布呈現(xiàn)出兩側(cè)高中間低的變化趨勢(shì),顆粒處于主流化區(qū)域,顆粒速度的對(duì)稱性增強(qiáng)。對(duì)于分布板C,流化床中間區(qū)域的顆粒速度增加,在三種分布板中,中間區(qū)域達(dá)到最大值。在兩側(cè)位置,顆粒速度仍呈現(xiàn)先下降后上升的趨勢(shì)。對(duì)于分布板A,隨著軸向高度增加至z=0.1 m,床層內(nèi)的顆粒速度分布趨勢(shì)與z=0.05 m相比沒有太大的改變。

從圖7c可以看出,在z=0.15 m高度時(shí),對(duì)于分布板A,顆粒速度與z=0.1 m時(shí)具有相似的變化趨勢(shì),顆粒速度左側(cè)高右側(cè)低,波動(dòng)范圍較小,且顆粒的速度顯著減小。布設(shè)分布板B的流化床內(nèi)的顆粒流化高度已經(jīng)達(dá)到極限,僅有極少數(shù)的顆粒存在,顆粒的速度趨近于0。對(duì)于分布板C,顆粒的速度仍呈現(xiàn)中間高兩邊低的趨勢(shì)。從圖8a～c可知,分布板開孔形式的不同,對(duì)顆粒的流化效果影響不同。對(duì)于均勻孔徑的開孔方式,顆粒呈現(xiàn)非對(duì)稱性流化狀態(tài),床層內(nèi)部的顆粒速度集中分布于壁面附近,顆粒能達(dá)到的最大高度較低。對(duì)于分布板A與C非均勻開孔的方式,隨著軸向高度的增加,床層內(nèi)部的顆粒速度分布變化比較小,顆粒能達(dá)到較大的高度。

3 結(jié) 論

不同的氣體分布板結(jié)構(gòu)會(huì)影響流化床內(nèi)氣固兩相的流動(dòng)特性。改變開孔直徑與開孔方式,將會(huì)直接影響床層內(nèi)部的氣體的流向,同時(shí)顆粒的分布狀態(tài)與動(dòng)力學(xué)特性也會(huì)隨之改變,進(jìn)而影響流化床的流化效果。

(1)布設(shè)分布板C時(shí),流化床內(nèi)的顆粒運(yùn)行狀態(tài)穩(wěn)定,顆粒與氣體的摻雜更加充分,有利于增強(qiáng)氣體對(duì)顆粒的流化作用。布設(shè)分布板A與B時(shí),流化床內(nèi)的顆粒呈現(xiàn)周期性的運(yùn)動(dòng),顆粒多聚集在壁面附近,不利于氣固兩相充分混合。

(2)布設(shè)分布板A與分布板B的流化床內(nèi)部產(chǎn)生位置不穩(wěn)定氣體環(huán)流,顆粒不能達(dá)到均勻分布的狀態(tài),分布板C內(nèi)部氣體分布較為集中,有利于顆粒的快速達(dá)到流化狀態(tài)。

(3)布設(shè)分布板C時(shí),顆粒在流化床內(nèi)形成明顯的對(duì)稱環(huán)流,有利于氣體和顆粒之間傳質(zhì)作用的發(fā)生。分布板A與分布板B的流化床內(nèi)部產(chǎn)生偏心環(huán)流區(qū),顆粒之間劇烈摩擦,容易產(chǎn)生過粉碎現(xiàn)象。

(4)均勻孔徑開孔的分布板,床層內(nèi)部的顆粒呈現(xiàn)非對(duì)稱的流化狀態(tài),顆粒多集中在內(nèi)壁附近。非均勻開孔方式的分布板,床層內(nèi)部的顆粒速度變化較小,能達(dá)到較高的位置。