李 慧，王 飛，馬新輝

(1.黃河科技學院，河南 鄭州 450063；2.西安航天華威化工生物工程有限公司，陜西 西安 710100)

流化床是一種優(yōu)良的生物質(zhì)熱解轉(zhuǎn)化的反應器，而布風板作為流化床的關(guān)鍵部件，除了能支承物料以外，還能保證流化的穩(wěn)定，一定程度上降低壓降[1]。布風板布孔排列形式和布孔直徑一直是研究的重點問題之一。

單純通過實驗來研究布風板在流化床中的作用周期長、費用高，而精確的解析求解又極為困難。計算流體動力學(CFD)是一種有前景的方法，已經(jīng)有研究成功地使用CFD方法來預測流化床內(nèi)的流動狀態(tài)[2-5]。作者使用CFD商業(yè)軟件包Fluent軟件對循環(huán)流化床冷模布風板的開孔直徑和布孔形式做了數(shù)值模擬研究，以期確定布風板的最佳結(jié)構(gòu)形式。

1 布風板的結(jié)構(gòu)

1.1 流化床入口段平面示意圖

流化床入口段幾何形狀和尺寸和分布板開孔排列形式見圖1。其中，布風板開孔形式為正三角形和正方形[6-7]。

作者主要討論3種情況下，布風板開孔直徑d和孔間距L對流化床性能的影響。d=0.8 mm，L=2.5；d=1.5 mm，L=3.5 mm和d=3 mm，L=5.5 mm。

圖1 布板位置和結(jié)構(gòu)

1.2 流化床入口段三維簡化模型

為真實反映孔與孔之間氣流的相互作用與影響，從而更準確的模擬內(nèi)部流場，模擬實驗采用三維模型。布風板和流體的簡化縮小模型之一見圖2，其中，布風板開孔形式為正方形排列，開孔直徑為3 mm，開孔間距為5.5 mm。

圖2 布風板與流體三維模型

2 數(shù)學模型

在應用FLUENT軟件對圖1所示幾何模型進行模擬之前，利用相似原理和歐拉相似準則[7]對模型進行縮小。

相似原理：

(1)

歐拉相似準則：

(2)

操作工況相同，所以Pm=Pp且Tm=Tp，根據(jù)克拉柏倫方程式PM=ρRT，可得：ρm=ρp，

其中：Lp—實物總長；Lm—模型總長；dP—實物直徑；dm—模型直徑；CL—長度系數(shù)；Pp—實際氣體壓強；Pm—模擬氣體壓強；Tp—實際氣體溫度；Tm—模擬氣體溫度；Vp—實際氣體速度；Vm—模擬氣體速度；Eup—實際氣體歐拉系數(shù)；Eum—模擬氣體歐拉系數(shù)；ρP—實際氣體密度；ρm—模擬氣體密度；M—實際氣體分子量；R—摩爾氣體常數(shù)。

邊界條件為入口截面平均氣速為2.6 m/s，操作溫度為25 ℃，操作壓力為0.15 MPa，操作流體為氮氣。計算模型的網(wǎng)格劃分見圖3。

圖3 三維流體網(wǎng)格

3 結(jié)果與分析

3.1 正交實驗數(shù)值模擬

正交模擬實驗安排見表1和表2。

表1 布風板模擬水平與因子表

表2 布風板正交模擬表

3.2 模擬結(jié)果與分析

3.2.1 無布風板和有布風板(開孔正方形排列d=0.5 mm)比較

流化床內(nèi)沒有設置布風板時的模擬結(jié)果見圖4和圖6，流化床內(nèi)裝有布風板的情況見圖5和圖7。從圖4、圖6中可以看到流體的壓強沿流化床直管段持續(xù)下降，并且管內(nèi)呈現(xiàn)負壓狀態(tài)，壓強不穩(wěn)定。從圖5、圖7中可以清楚看到有了布風板之后，流體通過時，床內(nèi)壓強迅速穩(wěn)定。圖6中顯示無布風板時流體從入口到出口壓力損失約18.3 Pa；圖7顯示當氣體通過有布風板的床層時壓降達約3.2×104Pa。所以，盡管增加布風板可以實現(xiàn)穩(wěn)壓作用，形成等壓室，但是增加布風板所造成的壓力損失相當大，詳細模擬所得數(shù)據(jù)見表3。

圖4 無布風板壓強云圖(X=0截面)

圖5 有布風板壓強云圖(X=0截面)

軸向位置/mm圖6 無布風板軸線壓強分布

軸向位置/mm圖7 有布風板軸線壓強分布

表3 模擬結(jié)果

3.2.2 有布風板時的速度和壓強分布(正方形和三角形排列)

依據(jù)表2正交實驗表進行模擬，結(jié)果見圖8、圖9、圖10、圖11、圖12、圖13。

圖8、圖10和圖12顯示了流體沿軸線的速度分布曲線。從圖8、圖10、圖12都可以看到，在相同的布孔直徑的情況下，正三角形布孔形式與正方形布孔形式相比，正三角形排列的速度流場更容易達到穩(wěn)定，且都有峰值速度出現(xiàn)。布風板開孔0.8 mm和1.5 mm時出現(xiàn)了超過60 m/s的最大速度峰值，這主要是因為布風板中心有空的原因。

軸向位置/mm圖8 開孔直徑d=0.8 mm軸向速度圖

軸向位置/mm圖9 開孔直徑d=0.8 mm軸向壓強圖

軸向位置/mm圖10 開孔直徑d=1.5 mm軸向速度圖

軸向位置/mm圖11 開孔直徑d=1.5 mm軸向壓強圖

軸向位置/mm圖12 開孔直徑d=3 mm軸向速度圖

從壓強圖9、圖11、圖13可看出，氣體通過布風板的壓力很快就達到穩(wěn)定，且正三角形排列形式下床內(nèi)的壓降比正方形排列形式下的壓力損失小，壓降低。因此，從減少壓降、控制能量損失這方面來講，選擇正三角形排列的布風板更為合適。就曲線發(fā)展趨勢看，相同的布孔形式下，開孔越小對應位置上的速度將越小，而壓降則越大。

分析圖8～圖13，可明顯看出，在布風板中心開孔對流化床內(nèi)的速度場影響很大，但對壓力場的影響相對較小。因此，從流場穩(wěn)定角度來看，盡量避免在布風板中心開孔。

軸向位置/mm圖13 開孔直徑d=3 mm軸向壓強圖

4 結(jié) 論

(1) 在流化床內(nèi)設置布風板有利于床內(nèi)速度場和壓力場的穩(wěn)定；

(2) 布風板開孔直徑越小，壓降越大，但流化床內(nèi)的速度場越容易達到穩(wěn)定；

(3) 布風板開孔直徑相同時，排列方式為等邊三角形比正方形時壓降更??；

(4) 為了得到較穩(wěn)定的速度場和壓力場，應避免在布風板中心處開孔。

[ 參 考 文 獻 ]

[1] A IBRAHIM.A generalized fludized bed reactor model across the flow regimes[D].America：The University of British Columbia，2001：14-33.

[2] Y Tsuji,T Ifawaguchi，et al.Discrete particle simulation of two-dimensional fluidized bed[J].Powder Technology，1993,79(7):77-79.

[3] P A Langston,U Tuzun，D M.Heyes.Discrete element simulation of granular flow in 2D and 3D hoppers：dependence of discharge rate and wall stress on particle interactions[J].Chemical Engineering Science，1995，50(6)：967-987.

[4] B P B Hoomans，J A M Kuipers.Discrete particle simulation of bubble formation in a two-dimensional fast-fluidized bed：a hard-sphere approach[J].Powder Technology，1993,79(10):99-118.

[5] Masayuki Horio.Numerical simulation of fluidized bed processes-review of latest achievements of numerical simulation and discussion on the future[C].Japan：Conferences Papers Seventh China-Japan Symposium,2000：1126-1128.

[6] 林宗虎，魏敦崧，安恩科.循環(huán)流化床鍋爐[M].北京：化學工業(yè)出版社,2004:124-128.

[7] 姚運金,張學才,陳明強.噴動流化床內(nèi)木屑流動特性的研究(Ⅲ)— 床錐底角的影響[J].淮南工業(yè)學院學報,2002,22(4)：57-60.