王 斌

（國(guó)電豐城發(fā)電有限公司，江西豐城 331100）

0 引言

鍋爐冷、熱態(tài)性能試驗(yàn)是火力發(fā)電廠的重要試驗(yàn)之一，但其工作量較大、工況調(diào)整困難和測(cè)量方法復(fù)雜等。并且要獲得爐內(nèi)各物理量的細(xì)節(jié)信息方面有很大難度，而計(jì)算流體力學(xué)CFD軟件為鍋爐研究提供了一個(gè)虛擬平臺(tái)。本文將針對(duì)該爐進(jìn)行爐內(nèi)冷態(tài)空氣動(dòng)力場(chǎng)的數(shù)值模擬，得到在該鍋爐結(jié)構(gòu)和燃燒器布置方式下爐內(nèi)空氣動(dòng)力特性和結(jié)焦原因。

1 鍋爐本體結(jié)構(gòu)介紹

本文所模擬的對(duì)象是某電廠2號(hào)鍋爐系哈爾濱鍋爐廠有限責(zé)任公司自主開發(fā)研制的HG-1964/25.4-YM17型墻式切圓燃燒超臨界鍋爐。該鍋爐為一次中間再熱、超臨界壓力變壓運(yùn)行帶內(nèi)置式再循環(huán)泵啟動(dòng)系統(tǒng)的直流鍋爐、單爐膛、平衡通風(fēng)、固態(tài)排焦、全鋼架、全懸吊結(jié)構(gòu)、露天布置的π型鍋爐。采用中速磨直吹制粉系統(tǒng)，每爐配6臺(tái)HP1003碗式中速磨煤機(jī)，燃用設(shè)計(jì)煤種時(shí)5臺(tái)運(yùn)行，1臺(tái)備用。其結(jié)構(gòu)示意圖見圖1。

2 冷態(tài)空氣動(dòng)力場(chǎng)的數(shù)值模擬

2.1 數(shù)值方法

圖1 爐膛的結(jié)構(gòu)尺寸和燃燒器切圓示意圖（單位：mm）

由于四墻切圓鍋爐爐膛結(jié)構(gòu)復(fù)雜，為了使數(shù)值計(jì)算可行、準(zhǔn)確，需對(duì)實(shí)際爐膛結(jié)構(gòu)進(jìn)行合理簡(jiǎn)化處理。本文采用了非均勻六面體網(wǎng)格，由于該鍋爐所劃分的網(wǎng)格線與流動(dòng)方向垂直，因此無(wú)需考慮模擬中產(chǎn)生的偽擴(kuò)散現(xiàn)象，共生成網(wǎng)格數(shù)1553536個(gè)。在通用方程的離散操作中,對(duì)流項(xiàng)采用二階迎風(fēng)格式，擴(kuò)散項(xiàng)采用二階中心差分格式。對(duì)。對(duì)壁面附近粘性層中的模擬計(jì)算,本文采用壁面函數(shù)法予以處理，壓力速度耦合采用SIMPLE算法處理

出口壓力設(shè)為-90Pa。

在固體壁面上，氣流按無(wú)滑移條件且無(wú)湍流運(yùn)動(dòng)處理，即

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

表2 流場(chǎng)微分方程組通用形式

2.3 邊界條件及模擬工況

模擬計(jì)算中燃燒器噴口處各物理量取均勻入口條件，入口條件設(shè)為速度入口，即一次風(fēng)風(fēng)速V1=25m/s，二次風(fēng)風(fēng)速V2=56m/s，Vsofa=60m/s，爐膛出口設(shè)為壓力出口。具體分3個(gè)模擬工況：

工況1：設(shè)計(jì)工況，即V1=25m/s，V2=56m/s，Vsofa=60m/s。

工況2：1、3號(hào)角速度同時(shí)減小，設(shè)為V1=20m/s，V2=48m/s，Vsofa=60m/s。2、4號(hào)角保持不變。

工況3：1號(hào)角速度減小，設(shè)為V1=20m/s，V2=48m/s，Vsofa=60m/s。2、3、4號(hào)角保持不變。

3 模擬結(jié)果及分析

圖2中，從工況1即設(shè)計(jì)工況可以看出，氣流切向運(yùn)動(dòng)形成明顯的切圓燃燒工況，切圓中心在爐膛中央，未發(fā)生氣流偏斜現(xiàn)象。前、后、左、右最大貼壁速度均小于6.5m/s,不會(huì)造成嚴(yán)重的氣流刷墻現(xiàn)象和結(jié)焦問(wèn)題。

工況2中切圓中心偏向較為明顯，即偏向射流剛性弱的1號(hào)角。從工況2中可以看出前墻和右墻的貼壁速度明顯大于后墻和左墻的貼壁速度，貼壁速度最大為18m/s，易造成氣流刷墻的現(xiàn)象，引起爐內(nèi)產(chǎn)生局部超溫和結(jié)焦的現(xiàn)象。

工況3中由于1、3號(hào)角氣流動(dòng)量比2、4號(hào)角動(dòng)量小，射流易發(fā)生偏斜，導(dǎo)致前、后、左、右墻都有明顯的貼壁現(xiàn)象，貼壁速度大概為12m/s。由于是均勻減小對(duì)角一、二次風(fēng)風(fēng)速，所以爐內(nèi)切圓整體未發(fā)生明顯的偏離，但切圓直徑有所增大，這是導(dǎo)致貼壁速度增大的主要原因。

圖2 同層一、二次風(fēng)噴嘴截面上的速度矢量

4 模擬結(jié)果與試驗(yàn)結(jié)果對(duì)比分析

為較好地驗(yàn)證數(shù)值模擬的正確性及可信度,將數(shù)值模擬結(jié)果與冷態(tài)試驗(yàn)數(shù)據(jù)進(jìn)行定量對(duì)比分析,以驗(yàn)證誤差是否在允許范圍內(nèi)。以設(shè)計(jì)工況為模擬和試驗(yàn)工況，將其結(jié)果進(jìn)行對(duì)比及分析，得出的結(jié)論如圖3所示。

由圖3可知，爐內(nèi)氣流速度沿逆時(shí)針?lè)较蛐D(zhuǎn)，模擬切圓比實(shí)際切圓稍大，且兩切圓橢圓度都較好，切圓中心基本與爐膛中心吻合，符合設(shè)計(jì)工況的要求。模擬值與實(shí)測(cè)值的曲線分布趨勢(shì)較相互較為吻合，但存在一定程度的偏差，大約是3～5m/s。誤差的原因主要是由于三維模型的簡(jiǎn)化和儀器的系統(tǒng)誤差引起的。

圖3 爐內(nèi)測(cè)量風(fēng)速得到的實(shí)際切圓和模擬切圓

5 結(jié)論

1）從工況1可知，模擬結(jié)果與試驗(yàn)結(jié)果基本一致，說(shuō)明該數(shù)值模型和方法能較好地反映爐內(nèi)真實(shí)的流場(chǎng)特性，更好的保證結(jié)論的可靠性和準(zhǔn)確性，將更好的指導(dǎo)實(shí)際生產(chǎn)工作。

2）從數(shù)值模擬結(jié)果來(lái)看，射流的剛性與氣流的動(dòng)量有很大的關(guān)系，動(dòng)量越大則越不易偏斜，動(dòng)量小則較易偏斜。氣流偏斜較嚴(yán)重時(shí)，則容易引起氣流刷墻，從而導(dǎo)致水冷壁的高溫腐蝕和結(jié)焦等問(wèn)題。

3）從三個(gè)工況來(lái)看，以不平衡供風(fēng)工況與正常工況對(duì)比，得出不均等配風(fēng)工況對(duì)爐內(nèi)空氣動(dòng)力場(chǎng)的影響，可以通過(guò)結(jié)合數(shù)值模擬，改變配風(fēng)來(lái)優(yōu)化爐內(nèi)流場(chǎng)，保證鍋爐的安全經(jīng)濟(jì)運(yùn)行。

[1]陶文銓.數(shù)值傳熱學(xué)[M].(第2版).西安:西安交通大學(xué)出版社,2001.

[2]劉勇,魏風(fēng),唐必光.四墻切圓鍋爐冷態(tài)空氣動(dòng)力場(chǎng)流動(dòng)特性的實(shí)驗(yàn)研究[J].武漢大學(xué)學(xué)報(bào)(工學(xué)版),2002,35(6):52～55.

[3]王鵬，翔楊威，陸慧林.200MW四墻切向煤粉鍋爐空氣動(dòng)力場(chǎng)冷態(tài)?；蛿?shù)值模擬[J].電站系統(tǒng)工程學(xué)報(bào)，2002，18（4）：22～25.

[4]宋國(guó)良,周俊虎,翁衛(wèi)國(guó).75t/h四墻切圓鍋爐冷態(tài)空氣動(dòng)力場(chǎng)的試驗(yàn)研究[J].電站系統(tǒng)工程學(xué)報(bào)，2005，21（1）：62～65.

[5]潘衛(wèi)國(guó),林泉,孫堅(jiān)榮,等.一次風(fēng)噴口速度分布對(duì)四墻切圓燃燒鍋爐爐內(nèi)空氣動(dòng)力場(chǎng)影響的數(shù)值研究[J].動(dòng)力工程,2004,24(4):495-500.