DD分解爐流場(chǎng)仿真分析

DD分解爐流場(chǎng)仿真分析

我院設(shè)計(jì)的某項(xiàng)目DD型分解爐，其結(jié)構(gòu)尺寸設(shè)計(jì)主要依據(jù)經(jīng)驗(yàn)值，現(xiàn)業(yè)主方要求對(duì)設(shè)計(jì)結(jié)構(gòu)合理性進(jìn)行驗(yàn)證。為降低投資成本及縮短驗(yàn)證周期，本文采用數(shù)值模擬仿真技術(shù)取代傳統(tǒng)試驗(yàn)方法，對(duì)所設(shè)計(jì)的分解爐建立相應(yīng)三維模化理論模型，通過(guò)模擬其內(nèi)部流場(chǎng)，得到相應(yīng)的壓力場(chǎng)、溫度場(chǎng)、顆粒運(yùn)動(dòng)流場(chǎng)等值。根據(jù)模擬結(jié)果驗(yàn)證設(shè)計(jì)結(jié)構(gòu)的合理性，同時(shí)找出優(yōu)化設(shè)計(jì)方案。

分解爐；三維；優(yōu)化；模擬仿真

1 前言

針對(duì)目前水泥工業(yè)蓬勃發(fā)展的預(yù)分解窯技術(shù)中的核心設(shè)備——分解爐的設(shè)計(jì)，設(shè)計(jì)人員大多依據(jù)經(jīng)驗(yàn)值及參考國(guó)外相應(yīng)數(shù)據(jù)進(jìn)行設(shè)計(jì)，在自主研發(fā)方面缺少創(chuàng)新性。目前我院對(duì)某項(xiàng)目設(shè)計(jì)的DD型分解爐，亦采用經(jīng)驗(yàn)進(jìn)行設(shè)計(jì)。該尺寸在傳統(tǒng)DD爐設(shè)計(jì)上進(jìn)行了改進(jìn)，因此業(yè)主方及設(shè)計(jì)人員均對(duì)設(shè)計(jì)結(jié)構(gòu)的合理性存在疑惑。若對(duì)設(shè)計(jì)結(jié)構(gòu)的準(zhǔn)確性采用傳統(tǒng)試驗(yàn)方法進(jìn)行驗(yàn)證，其周期長(zhǎng)，投資成本高，目前不可行。而國(guó)內(nèi)外很多研究者在實(shí)驗(yàn)和理論的指導(dǎo)下采用數(shù)值模擬手段對(duì)分解爐進(jìn)行研究，通過(guò)分析和統(tǒng)計(jì)各種物理量來(lái)分析所設(shè)計(jì)的分解爐內(nèi)復(fù)雜的物理化學(xué)過(guò)程和運(yùn)動(dòng)規(guī)律，得到分解爐內(nèi)各種物理化學(xué)場(chǎng)的分布情況和規(guī)律，以助于指導(dǎo)分解爐的優(yōu)化設(shè)計(jì)，達(dá)到分解爐內(nèi)復(fù)雜過(guò)程的可視化，同時(shí)為技術(shù)人員提供詳盡可靠的操作信息，并且可以對(duì)整個(gè)水泥生產(chǎn)過(guò)程的優(yōu)化起到促進(jìn)作用[1,2]。作為實(shí)驗(yàn)手段的一種延伸，數(shù)值模擬技術(shù)有著更加優(yōu)越的研究效率[2,3]。

因此，本文采用數(shù)值模擬技術(shù)，運(yùn)用Fluent軟件，對(duì)該項(xiàng)目中設(shè)計(jì)的分解爐內(nèi)三維流場(chǎng)合理性進(jìn)行分析，并將結(jié)構(gòu)設(shè)計(jì)進(jìn)行優(yōu)化。本文主要對(duì)設(shè)計(jì)結(jié)構(gòu)合理性分析中數(shù)值模擬部分，以及優(yōu)化結(jié)構(gòu)的一種方案進(jìn)行闡述。

2 解決思路

為驗(yàn)證該結(jié)構(gòu)合理性，按下面思路進(jìn)行模擬：

（1）先求單相流場(chǎng)運(yùn)動(dòng)，來(lái)求得穩(wěn)定的內(nèi)流場(chǎng)，目的是為得到分解爐內(nèi)空氣流場(chǎng)的速度分布及渦耗散分布情況，從而分析結(jié)構(gòu)合理性。此時(shí)湍流模型采用RNG k-ε模型（Re-Normalization Group k-ε turbu?lence model的簡(jiǎn)稱）。

（2）然后加入氣固兩相流模型，得出CaCO3的停留時(shí)間及分布，進(jìn)一步分析結(jié)構(gòu)合理性。采用DPM（Dis?crete partical model）模型的歐拉-拉格朗日（Eulerian-Lagrangian）模擬。

（3）對(duì)設(shè)計(jì)結(jié)構(gòu)進(jìn)行比對(duì)，優(yōu)化結(jié)構(gòu)設(shè)計(jì)。

3 計(jì)算模型

3.1 建立模型

根據(jù)流量選擇，設(shè)計(jì)分解爐尺寸如圖1a。模擬分析主要針對(duì)內(nèi)流場(chǎng)，因此只保留影響內(nèi)流場(chǎng)的結(jié)構(gòu)尺寸，將模型進(jìn)行簡(jiǎn)化，不考慮澆筑層厚度，尺寸見(jiàn)圖1b。設(shè)計(jì)中物料從C4進(jìn)入分解爐，采用導(dǎo)流板將物料導(dǎo)入分解爐。此處模擬，將導(dǎo)流板結(jié)構(gòu)簡(jiǎn)化，同等位置、同等角度采用同截面積的方形接口取代，按照?qǐng)D1尺寸，建立三維模型，如圖2。

表1 邊界條件指定

3.2 網(wǎng)格劃分（圖3）

要進(jìn)行流場(chǎng)的數(shù)值模擬計(jì)算，首先需要將所要計(jì)算的空間或區(qū)域離散化，區(qū)域的網(wǎng)格化是把微分方程離散化并奠定數(shù)值求解的基礎(chǔ)，模擬仿真中70%～80%的工作量都在網(wǎng)格劃分上，網(wǎng)格劃分的方式將直接影響到方程離散化的難易，影響計(jì)算速度和所需的存儲(chǔ)量，并影響數(shù)值解的收斂性和準(zhǔn)確性。在劃分網(wǎng)格時(shí)，要求網(wǎng)格能夠適應(yīng)計(jì)算區(qū)域中物理量的變化情形，即在物理量變化比較劇烈的地方適當(dāng)?shù)丶用芫W(wǎng)格數(shù)量，以保證計(jì)算精度，在物理量變化平緩的地方適當(dāng)減少網(wǎng)格數(shù)量，以減少計(jì)算量[4]。

采用ANSYS Workbench的Meshing網(wǎng)格劃分工具來(lái)劃分分解爐的CFD網(wǎng)格，在規(guī)則形狀上，采用垂直于流速方向的六面體網(wǎng)格來(lái)劃分。網(wǎng)格尺寸給定100mm，核心六面體網(wǎng)格的劃分方法，總數(shù)量為164354，采用斜率（skewness）評(píng)判標(biāo)準(zhǔn)，網(wǎng)格質(zhì)量，平均值為0.44。

3.3 數(shù)值計(jì)算方法

將通用輸運(yùn)方程在網(wǎng)格上利用控制容積方法進(jìn)行離散，標(biāo)量和矢量的控制容積采用交錯(cuò)網(wǎng)格。在求解動(dòng)量方程時(shí)，對(duì)于對(duì)流項(xiàng)和擴(kuò)散項(xiàng)采用混合差分格式進(jìn)行離散，求解差分方程用壓力－速度修正的SIMPLE算法。差分方程組在切向方向上采用逐線掃描的CTDMA迭代求解，在徑向和軸向方向上采用逐線掃描的TDMA迭代求解。計(jì)算時(shí)采用亞松弛技術(shù)保持計(jì)算的穩(wěn)定性[5,6,7]。

3.4 邊界條件

該結(jié)構(gòu)設(shè)計(jì)為日產(chǎn)150t磁化焙燒系統(tǒng)使用，分解爐下端接熱風(fēng)爐。為保證系統(tǒng)足夠的風(fēng)量及還原氣氛，根據(jù)熱工計(jì)算，要求分解爐入口風(fēng)速為35m/s。在Fluent軟件中給定如圖4及表1所示的邊界及相應(yīng)參數(shù)。

4 DD分解爐模擬結(jié)果合理性分析

4.1 單相計(jì)算結(jié)果分析

4.1.1 壓力場(chǎng)

根據(jù)質(zhì)量加權(quán)平均，計(jì)算得出該分解爐結(jié)構(gòu)的壓降為386.6Pa。圖5為分解爐內(nèi)壓力分布云圖和等值圖。

可以看出，分解爐內(nèi)區(qū)域靜壓值均為負(fù)值，壓力分布較均勻。分解爐下半段三分之二處至縮口處，壓力值最大。物料入口附近壓力大于窯尾煙氣入口段，這種壓力分布規(guī)律對(duì)噴入的生料擴(kuò)散是有利的。

4.1.2 流線圖

圖6為煙氣入口處的氣體流線圖，曲線顏色顯示為氣流在爐體內(nèi)運(yùn)動(dòng)時(shí)間。從流線圖上可以看出，氣體在軸心兩側(cè)均勻分布，平均運(yùn)動(dòng)時(shí)間在5s內(nèi)。

從圖7可以看出，從物料口進(jìn)入的氣流由于風(fēng)速較小，被窯尾煙氣氣流吹散，基本在窯尾煙氣氣流兩側(cè)運(yùn)動(dòng)，經(jīng)過(guò)縮口，能量損失，然后貼近右側(cè)處運(yùn)動(dòng)。同時(shí)從放大圖可以看出，從物料口進(jìn)入的氣流在交界面處附近形成回旋運(yùn)動(dòng)。從物料口進(jìn)入的氣流在爐內(nèi)運(yùn)動(dòng)時(shí)間為13s內(nèi)。根據(jù)定量結(jié)果提取，得出全爐內(nèi)氣體的平均運(yùn)動(dòng)時(shí)間為8.2s。

4.1.3 速度圖（圖8）

從下端進(jìn)入的窯尾煙氣和右側(cè)進(jìn)入的物料氣流進(jìn)入爐內(nèi)匯合，由于窯尾煙氣速度較大，匯合后的氣流大部分射向物料進(jìn)口處，然后與爐壁面約成30°夾角以較大的流速向上運(yùn)動(dòng)，從Y=0面速度圖來(lái)看，由于只有一個(gè)入口的原因，速度場(chǎng)偏對(duì)稱分布。結(jié)合流線圖，可以看出物料氣流在窯尾煙氣較大的氣流速度帶動(dòng)下，向周圍區(qū)域迅速擴(kuò)散，有利于生料的充分混合和化學(xué)反應(yīng)。

從分解爐各段速度矢量圖可以看出，由于物料下料口速度較低，氣流在出交界面后靠近下料口端有較小漩渦。同時(shí)，在縮口處并沒(méi)有形成DD爐的噴騰效應(yīng)。分析其原因，可能是分解爐直筒段較長(zhǎng)，動(dòng)量變化梯度較小，使壓力與速度呈均勻分布，因此在縮口處未出現(xiàn)較明顯的噴騰效應(yīng)。

4.1.4 渦耗散圖（圖9）

從等值圖及各段等值圖可以看出，主要三段有漩渦：分解爐出口處、縮口處及交界面段。其中分解爐出口及縮口處的渦變化較小。在Y=-13m～Y=-10m處，漩渦強(qiáng)度最大。而此范圍正好是物料進(jìn)口位置所設(shè)之處，且開(kāi)口位置恰好離這幾個(gè)漩渦中心最近。預(yù)示若物料以一定的速度噴入爐內(nèi)，極有利于生料的擴(kuò)散并充分混合和反應(yīng)，進(jìn)一步表明該分解爐的生料進(jìn)口位置設(shè)計(jì)較合理。

4.2 氣固兩相流計(jì)算結(jié)果

按表2顆粒分布。

圖10為顆粒分布停留時(shí)間圖。截取其中10支流線來(lái)分析，物料從下料口進(jìn)入，在附近受氣流影響，做回旋運(yùn)動(dòng)，這是由于兩股氣流壓力差及速度場(chǎng)相互作用，產(chǎn)生漩渦能量，隨后引起的顆粒在交界面處回旋運(yùn)動(dòng)。

同時(shí)，物料在窯尾煙氣的作用下，非對(duì)稱地運(yùn)動(dòng)，結(jié)合圖9aY=-13m、Y=-12m、Y=-11m三個(gè)平面速度云圖可知，由于在軸向有較大渦流能量的存在，物料略微偏離生料口運(yùn)動(dòng)。其中，物料在中心處停留時(shí)間最長(zhǎng)，在18s左右。從這點(diǎn)也可以看出，窯尾煙氣的速度和物料下料口位置在設(shè)計(jì)中匹配合適。

5 DD分解爐優(yōu)化方案分析

由以上分析可知，該結(jié)構(gòu)縮口作用不大，沒(méi)有形成DD爐的噴騰效應(yīng)。分析原因，一方面是工況參數(shù)不適合，另一方面是結(jié)構(gòu)設(shè)計(jì)，因?yàn)榉纸鉅t長(zhǎng)度太長(zhǎng)，縮口角度及尺寸不匹配，沒(méi)有形成DD爐的壓力及速度變化梯度所致。

表2 顆粒質(zhì)量分布比例

表3 修改前后結(jié)果比對(duì)

為達(dá)到DD爐的噴騰效應(yīng)，有多種修改方案，本文從設(shè)計(jì)結(jié)構(gòu)尺寸上修改。將直徑縮小，從1410mm縮小到1350mm；直段部分高度縮小，從450mm更改成200mm。

計(jì)算結(jié)果見(jiàn)表3、圖11。停留時(shí)間增加，截面風(fēng)速減小，壓降增加。

從圖12可以看出，改進(jìn)后的結(jié)構(gòu)，能更好地促進(jìn)顆粒在中心位置的收集。因此，縮口處直段長(zhǎng)度縮短，能增加系統(tǒng)壓力，延長(zhǎng)顆粒在縮口處的運(yùn)動(dòng)時(shí)間。

從圖13可以看出，縮口處的湍流強(qiáng)度也呈現(xiàn)雙峰效應(yīng)，中心處湍流強(qiáng)度比較高，兩側(cè)比較低。這樣有助于顆粒在中心運(yùn)動(dòng)，不至于爐壁出現(xiàn)燒結(jié)現(xiàn)象。

從圖14a可以看出，此種修改后的結(jié)構(gòu)，其速度流場(chǎng)在縮口中心處值較高，兩側(cè)低，窯尾煙氣縮口處流場(chǎng)偏中心，因此，適當(dāng)提供窯尾煙氣的速度有助于中心反應(yīng)。從圖14b壓力場(chǎng)可以看出，分解爐接近縮口下端壓力較高，經(jīng)縮口，壓力擴(kuò)散，因此見(jiàn)圖13b，看出在出縮口處流場(chǎng)能形成雙峰現(xiàn)象。

由此，縮短縮口處直徑及直段距離，是可行的改進(jìn)方案之一。

6 結(jié)論

（1）所設(shè)計(jì)的結(jié)構(gòu)生料下料口位置設(shè)計(jì)合理，流場(chǎng)分布趨勢(shì)合理。

（2）原設(shè)計(jì)結(jié)構(gòu)定量分析壓降為386.6Pa。

（3）為形成該爐的雙噴騰效應(yīng)，進(jìn)行設(shè)計(jì)改進(jìn)，得出縮短縮口處直段直徑及高度是改進(jìn)方案之一，本次調(diào)整將直徑從1410mm縮小到1350mm，直段的高度從450mm縮小到200mm。

（4）由以上分析可知，采用模擬仿真技術(shù)，對(duì)水泥設(shè)備設(shè)計(jì)及優(yōu)化過(guò)程有指導(dǎo)性作用。

[1]黃來(lái).水泥預(yù)分解爐內(nèi)物理化學(xué)過(guò)程模擬和優(yōu)化設(shè)計(jì)研究.華中科技大學(xué),2006,4.

[2]肖國(guó)權(quán).分解爐內(nèi)伴有燃燒、化學(xué)反應(yīng)的兩相湍流流動(dòng)的數(shù)值模擬.武漢理工大學(xué).2005.3.

[3]王漢立，張振平.水泥熱工設(shè)備與測(cè)試技術(shù)[M].化學(xué)工業(yè)出版社，2010.2.

[4]周力行.湍流兩相流動(dòng)與燃燒的數(shù)值模擬[M].北京；清華大學(xué)出版社，1991.

[5]孫銳，李爭(zhēng)起，等.不同湍流模型對(duì)強(qiáng)旋流動(dòng)的數(shù)值模擬.動(dòng)力工程.第22卷第三期2002.6.

[6]Sommerfeld M，Kohnen G and Ruger M.Some open questions and inconsistencies of La?grangian particle dispersion models，Proc.Ninth Symp On Turbulent Shear Flows,Kyoto,Jan?pan，Paper 5.1，1993.

[7]Chen X Q and Pereira J C F.Efficient com?putation of Particle dispersion in turbulent flow with a stochastic-probabilistic method，Int.J.Heat Mass Transfer,1997,40:1727-1741.

潘麗萍1，周濤2

Simulation Analysis of Flow Field in DD Calciner

PAN Li-ping1,ZHOU Tao2(1.Wuhan Building Material Industry Design&Research Institute Co.,Ltd.,Wuhan 430071,China;
2.Pera Global Technology Beijing Co.,Ltd.,Wuhan 430071,China)

According to the type of DD calciner designed by my Insti?tute for some program which its size based on experiences，the owners request us to check the design structure reasonableness now.In order to reduce the investment cost and shorten the test cycle,this paper adopts the numerical simulation technique to establish the 3D model and theory model which replaced the traditional test method.Through the simulation of internal flow field,we get the results such as the rele?vant pressure field,temperature field,particle flow and so on.All that will be used to prove the reasonable of the designed structure and find optimization of design scheme.

Calciner;Three-dimensional;Optimization;Simulation

TQ172.622.29

A

1001-6171（2012）04-0038-05

通訊地址：1武漢建筑材料工業(yè)設(shè)計(jì)研究院有限公司，湖北 武漢 430071；

2安世亞太科技北京有限公司，湖北 武漢 430071；

2011-11-29；

趙 蓮