王觀民王自興張凱博王文煜朱克敏欒素清解華華

1. 洛陽礦山機械工程設(shè)計研究院有限責(zé)任公司，河南 洛陽 471039

2. 礦山重型裝備國家重點實驗室（中信重工機械股份有限公司），河南 洛陽 471039

基于FLUENT的懸浮預(yù)熱器氣固耦合模擬與結(jié)構(gòu)優(yōu)化

王觀民1,2王自興1,2張凱博1,2王文煜1,2朱克敏1,2欒素清1,2解華華1,2

1. 洛陽礦山機械工程設(shè)計研究院有限責(zé)任公司，河南 洛陽 471039

2. 礦山重型裝備國家重點實驗室（中信重工機械股份有限公司），河南 洛陽 471039

根據(jù)某水泥廠應(yīng)用良好的懸浮預(yù)熱器C1筒，采用計算流體動力學(xué)軟件Fluent對其進行模擬，分析物料粒徑對分離效率的影響，得出進入預(yù)熱器物料的合理粒徑，再應(yīng)用到目前正在設(shè)計的懸浮預(yù)熱器的模擬優(yōu)化中，對其筒內(nèi)的煙氣和物料兩相流動進行三維湍流耦合模擬，得到筒內(nèi)煙氣的速度、壓力分布云圖、物料在旋風(fēng)筒內(nèi)的停留時間以及分離效率等。根據(jù)模擬分析結(jié)果確定合理的內(nèi)筒高度，降低了壓力損失并減少了耐熱鋼掛板的重量。

懸浮預(yù)熱器 Fluent 分離效率 優(yōu)化

0 引言

懸浮預(yù)熱器是預(yù)分解窯燒成系統(tǒng)的核心設(shè)備之一，承擔(dān)著氣固分散、物料加熱、氣固分離、物料輸送及部分物理、化學(xué)反應(yīng)等多項功能。分離效率和阻力損失是設(shè)計預(yù)熱器旋風(fēng)筒主要考慮的兩項性能指標。目前也有針對預(yù)熱器內(nèi)部流場和物料分離效率的研究，但是多局限于單一粒徑，沒有考慮物料的粒徑分布，水泥廠對進預(yù)熱器的粒度要求也僅是80μm篩余小于12%。故筆者首先針對某水泥廠應(yīng)用良好的懸浮預(yù)熱器C1筒結(jié)構(gòu)，應(yīng)用計算流體動力學(xué)軟件Fluent對其進行模擬，分析物料粒徑對分離效率的影響，得到滿足分離效率的物料粒徑臨界點，給出進入預(yù)熱器物料的合理粒徑，再應(yīng)用到某新型懸浮預(yù)熱器的模擬中，優(yōu)化其內(nèi)筒高度。

1 計算模型

1.1 物料相模型

為了計算的精確性，選用離散型模型(discrete phase model, DPM)，并運用其中的隨機軌道模型追蹤預(yù)熱器內(nèi)的顆粒相的運動。當(dāng)顆粒穿過流體運動時，顆粒的軌道以及傳熱量、傳質(zhì)量可通過當(dāng)?shù)亓黧w作用于顆粒上的各種平衡作用力、對流輻射引起的熱量質(zhì)量傳遞來進行計算[1]。離散相模型即可以通過在一個固定的流場中(非耦合方法)來預(yù)測離散型的分布，也可以在考慮離散型對連續(xù)相有影響的流場(相間耦合方法)中考察顆粒的分布。相間耦合計算中，離散相的存在影響了連續(xù)相的流場，而連續(xù)相的流場反過來又影響了離散相的分布[2]?？梢越惶嬗嬎氵B續(xù)相和離散相直到兩相的計算結(jié)果都達到收斂標準。

1.2 氣相模型

大量的研究已表明，預(yù)熱器單體內(nèi)部的流動狀態(tài)為不可壓縮湍流流動，而旋風(fēng)筒內(nèi)顆粒相的體積比率很低，滿足顆粒群軌道模型的基本條件。Reynolds應(yīng)力模型(RSM)是二階距模型以及在此基礎(chǔ)上經(jīng)簡化而得出的代數(shù)應(yīng)力模型。RSM采用二階關(guān)聯(lián)量來模擬未知的三階關(guān)聯(lián)量和其它關(guān)聯(lián)量，此模型對k方程和e方程中源相的模擬采用與k-e模型同樣的方法。在很多情況下，RSM模型能夠給出優(yōu)于標準k-e模型的結(jié)果[3-4]。針對懸浮預(yù)熱器內(nèi)三維強旋流的特點，以及懸浮預(yù)熱器不對稱的強旋流場的各向異性等特點，在模擬中氣相模型采用Reynolds應(yīng)力模型。

2 模擬分析與結(jié)構(gòu)優(yōu)化

2.1 不同粒徑顆粒的分離效率

根據(jù)某水泥廠應(yīng)用狀態(tài)良好的懸浮預(yù)熱器C1筒為依托來研究不同粒徑顆粒的分離效率，物料經(jīng)原料立磨后的粒度要求是80μm篩余小于12%，故研究顆粒粒徑的范圍在0～80μm。從表1可以看出，直徑越大的顆粒分離效應(yīng)越好，粒徑30μm是個臨界點，30μm～80μm的顆粒的分離效率能滿足工藝要求，即分離效率≥95%，所以在實際生產(chǎn)中需控制粒徑30μm以下的物料進入預(yù)熱器的含量，從而提高預(yù)熱器的分離效率。

表1 C1筒對不同粒徑顆粒物料的分離效率

2.2 模擬分析及結(jié)構(gòu)優(yōu)化

為了得到更為滿意的結(jié)果，進行了多次內(nèi)筒高度的模擬優(yōu)化，限于篇幅大小，只選取三種內(nèi)筒高度的模型進行論述。圖1、圖2和圖3中按(a)、(b)、(c)的順序，內(nèi)筒高度依次降低，其中(a)為經(jīng)驗公式計算的內(nèi)筒高度。

圖1 剖面1速度云圖

圖2 剖面2速度云圖

圖3 剖面1壓力云圖

氣體進入旋風(fēng)筒后沿邊壁旋轉(zhuǎn)向下運動，在底部發(fā)生折返，從中心螺旋流出，速度呈現(xiàn)出雙峰分布的特點。從圖1可以看到，模擬計算準確地捕捉到了這一特性，在平均入口風(fēng)速18 m/s的情況下，這三種內(nèi)筒高度的速度分布類似，速度大小接近，但是隨著內(nèi)筒高度的降低，高速區(qū)域向上移動。從圖2可以看到內(nèi)筒內(nèi)的速度分布呈環(huán)形分布，中間速度最低，向周圍遞增；隨著內(nèi)筒高度的降低，環(huán)形高速區(qū)的速度升高。另外可知被煙氣帶走的物料主要集中在內(nèi)筒內(nèi)壁附件，會加劇對內(nèi)筒內(nèi)壁附件的磨損。

從圖3可以看出，靜壓的分布在此剖面沿徑向由外向內(nèi)遞減，中心軸線附近靜壓較低，且為負壓區(qū)。該負壓區(qū)的出現(xiàn)，說明了在懸浮預(yù)熱器的中心位置存在真空區(qū)，這主要是由于在懸浮預(yù)熱器軸線附近存在著強度很高的強制渦。不同之處在于隨著內(nèi)筒高度的降低，壓力最高處和最低處的壓差不斷減小，由(a)結(jié)構(gòu)中的1 072 Pa減小到(c)結(jié)構(gòu)的896 Pa；預(yù)熱器進口和出口的壓損隨著內(nèi)筒高度的降低也不斷減小，由(a)結(jié)構(gòu)中的669 Pa減小到(c)結(jié)構(gòu)的603 Pa。

從表2可知，(a)、(b)兩種內(nèi)筒高度的旋風(fēng)筒分離效率都能滿足工藝要求，即≥90%，而(c)結(jié)構(gòu)的旋風(fēng)筒分離效率不能滿足工藝要求，三種結(jié)構(gòu)的物料停留時間相差不大。綜合以上因素確定(b)結(jié)構(gòu)的旋風(fēng)筒更合理，壓損比經(jīng)驗公式計算出的(a)結(jié)構(gòu)降低41 Pa，并且降低了耐熱鋼掛板的使用量，節(jié)省了成本。

表2 不同結(jié)構(gòu)的分離效率和停留時間

各級懸浮預(yù)熱器采用以上方法進行模擬分析及結(jié)構(gòu)優(yōu)化前后的性能參數(shù)對比如表3所示，優(yōu)化改進后的預(yù)熱器系統(tǒng)總體壓損降低321 Pa。經(jīng)過優(yōu)化后的預(yù)熱器成功應(yīng)用到唐山某水泥廠，反饋的數(shù)據(jù)表明預(yù)熱器系統(tǒng)總體壓損降低約300 Pa，這與本研究結(jié)果較吻合。

表3 各級懸浮預(yù)熱器優(yōu)化前后參數(shù)的對比

3 結(jié)束語

水泥生料顆粒直徑在30μm以上時，懸浮預(yù)熱器具備較高的分離效率。

優(yōu)化改進后的預(yù)熱器系統(tǒng)總體壓損降低321 Pa，降幅為12%，旋風(fēng)筒的總分離效率僅降低1.5%，說明優(yōu)化后的預(yù)熱器為系統(tǒng)節(jié)能減排要求的實現(xiàn)創(chuàng)造了條件。

[1] Raoufi A, Shams M, Kanani H. CFD analysis of flow field in square cyclones[J]. Powder Technol, 2009, 191(3): 349-357．

[2] Cortes C, Gil A. Modeling the gas and particle flow inside cyclone separators[J]. Prog Energy Combust Sci, 2007, 33(5): 409-452.

[3] 王博. 旋風(fēng)式預(yù)熱器內(nèi)氣固兩相運動的數(shù)值仿真研究[D].西安: 西安建筑科技大學(xué).2003.

[4] 陳全德、陳晶、崔素萍,等.水泥預(yù)分解技術(shù)與熱工系統(tǒng)工程[M]. 北京:中國建材工業(yè)出版社, 2005.

TQ172.622.24

A

1008-0473(2015)02-0035-03

10.16008/j.cnki.1008-0473.2015.02.009

2014-12-28）