李延偉， 胡善苗， 馬達(dá)夫

(1. 國(guó)家電力投資集團(tuán)河南電力有限公司沁陽發(fā)電分公司，河南焦作454150；2.東南大學(xué) 能源與環(huán)境學(xué)院，南京 210096；3. 上海發(fā)電設(shè)備成套設(shè)計(jì)研究院有限責(zé)任公司，上海 200240)

近年來，隨著國(guó)家各類節(jié)能環(huán)保政策的相繼出臺(tái)，煤電機(jī)組大力開展了脫硫脫硝工作。因此，對(duì)引風(fēng)機(jī)及其進(jìn)出口煙道進(jìn)行降阻尤為重要。同時(shí)，深入了解引風(fēng)機(jī)與其進(jìn)出口煙道流場(chǎng)間的相互影響，對(duì)提高引風(fēng)機(jī)性能、降低系統(tǒng)阻力、節(jié)約廠用電等具有重要的意義[1-2]。

目前，關(guān)于引風(fēng)機(jī)及其進(jìn)出口煙道流場(chǎng)與性能耦合的研究較少，并且多為單一方面的研究。已有的研究大都僅考慮引風(fēng)機(jī)本體，風(fēng)機(jī)邊界條件假設(shè)為風(fēng)機(jī)進(jìn)口速度或質(zhì)量流量均勻，未考慮實(shí)際運(yùn)行過程中引風(fēng)機(jī)進(jìn)出口煙道對(duì)風(fēng)機(jī)邊界條件的影響；大多研究?jī)H考慮進(jìn)出口煙道區(qū)域，忽略風(fēng)機(jī)內(nèi)部流場(chǎng)特征及其對(duì)整個(gè)煙道流場(chǎng)的影響[3-5]。對(duì)于引風(fēng)機(jī)進(jìn)出口煙道，其幾何特征不僅會(huì)影響鍋爐風(fēng)煙系統(tǒng)的阻力，而且會(huì)影響引風(fēng)機(jī)的運(yùn)行狀態(tài)。同時(shí)，引風(fēng)機(jī)進(jìn)出口煙道流場(chǎng)分布不均會(huì)影響引風(fēng)機(jī)的穩(wěn)定運(yùn)行，嚴(yán)重時(shí)會(huì)造成引風(fēng)機(jī)異常振動(dòng)，進(jìn)而對(duì)機(jī)組運(yùn)行的可靠性和安全性產(chǎn)生不良影響[6-7]。

將引風(fēng)機(jī)放入整個(gè)煙道系統(tǒng)中，對(duì)引風(fēng)機(jī)及其進(jìn)出口煙道的整體進(jìn)行全三維數(shù)值模擬計(jì)算，研究引風(fēng)機(jī)性能和煙道流場(chǎng)的變化，以使其更符合實(shí)際工程應(yīng)用情況。同時(shí)，以此為基礎(chǔ)可挖掘進(jìn)出口煙道外部幾何特征與煙道流場(chǎng)和引風(fēng)機(jī)性能之間的聯(lián)系，了解煙道流場(chǎng)紊亂的原因，具有實(shí)際指導(dǎo)意義與工程應(yīng)用價(jià)值。

筆者利用計(jì)算流體力學(xué)(CFD)技術(shù)對(duì)雙級(jí)動(dòng)葉可調(diào)軸流風(fēng)機(jī)及其進(jìn)出口煙道的流場(chǎng)與性能進(jìn)行數(shù)值模擬研究，為引風(fēng)機(jī)及其進(jìn)出口煙道流場(chǎng)優(yōu)化的數(shù)值模擬提供更加準(zhǔn)確可靠的理論依據(jù)，進(jìn)而為軸流引風(fēng)機(jī)的優(yōu)化設(shè)計(jì)提供改進(jìn)方向。

1 設(shè)備概況

某電廠600 MW機(jī)組配備2臺(tái)引風(fēng)機(jī)，引風(fēng)機(jī)采用雙級(jí)動(dòng)葉可調(diào)軸流風(fēng)機(jī)。引風(fēng)機(jī)設(shè)計(jì)參數(shù)見表1(TB工況為風(fēng)機(jī)最大工況，BMCR工況為鍋爐最大連續(xù)蒸發(fā)量工況)，電動(dòng)機(jī)額定功率為5 200 kW。

表1 引風(fēng)機(jī)設(shè)計(jì)參數(shù)

模型自除塵器出口至脫硫塔進(jìn)口，包含引風(fēng)機(jī)和煙道系統(tǒng)，將其稱為全煙道系統(tǒng)。因?yàn)閱闻_(tái)機(jī)組兩側(cè)除塵器出口至脫硫塔進(jìn)口的煙道布置完全對(duì)稱，所以僅選取一側(cè)煙道進(jìn)行建模計(jì)算。

2 數(shù)值模擬

2.1 模型建立

采用三維掃描儀對(duì)引風(fēng)機(jī)葉片的外形等進(jìn)行數(shù)據(jù)采集，利用三維建模軟件建立引風(fēng)機(jī)三維整機(jī)全流道模型，煙道通流部分的幾何特征根據(jù)施工圖紙按照1∶1進(jìn)行建模[8-9]。引風(fēng)機(jī)模型和全煙道系統(tǒng)模型見圖1。在滿足工程要求的前提下，為便于建模與模擬，對(duì)全煙道系統(tǒng)進(jìn)行如下假設(shè)：(1)將煙氣視為不可壓縮牛頓流體，定常流動(dòng)；(2)假設(shè)除塵器出口煙氣速度分布均勻；(3)忽略煙道中支撐桿等對(duì)流場(chǎng)影響較小的內(nèi)部構(gòu)造的影響[10-11]。

圖1 引風(fēng)機(jī)模型和全煙道系統(tǒng)模型

采用六面體結(jié)構(gòu)化網(wǎng)格劃分引風(fēng)機(jī)進(jìn)出口煙道，根據(jù)引風(fēng)機(jī)的結(jié)構(gòu)特征，將其分為進(jìn)氣區(qū)、動(dòng)葉輪區(qū)、導(dǎo)葉輪區(qū)和擴(kuò)壓區(qū)等區(qū)域[12-13]?？紤]到葉片的結(jié)構(gòu)復(fù)雜性，首先對(duì)動(dòng)葉輪區(qū)、導(dǎo)葉輪區(qū)進(jìn)行單流道網(wǎng)格劃分，然后將得到的單流道模型網(wǎng)格繞中心軸旋轉(zhuǎn)復(fù)制生成全流道結(jié)構(gòu)網(wǎng)格。進(jìn)出口煙道與引風(fēng)機(jī)之間通過Interface面進(jìn)行信息的交換傳輸[14]。對(duì)于煙道壁面等速度梯度較大的區(qū)域采用網(wǎng)格加密處理，首層網(wǎng)格高度為5 mm，網(wǎng)格增長(zhǎng)率為1.2，模型最大網(wǎng)格尺寸為500 mm。將風(fēng)機(jī)全壓、效率和煙道壓力損失作為特征參數(shù)，進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證，最終確定網(wǎng)格數(shù)量為1 600萬(見圖2)。

圖2 全煙道系統(tǒng)網(wǎng)格圖

2.2 邊界條件

將除塵器出口、脫硫塔進(jìn)口分別作為模型進(jìn)口、出口，同時(shí)將邊界條件設(shè)置為均勻速度進(jìn)口、壓力出口；風(fēng)機(jī)內(nèi)部旋轉(zhuǎn)區(qū)和靜止區(qū)之間采用多參考系模型進(jìn)行耦合計(jì)算，不同區(qū)域之間通過Interface面進(jìn)行連接。選擇RNGk-ε湍流模型，采用標(biāo)準(zhǔn)壁面函數(shù)，應(yīng)用SIMPLE算法求解壓力-速度的耦合方程，方程源項(xiàng)中的壓為項(xiàng)采用PRESTO！差分格式，對(duì)流項(xiàng)、擴(kuò)散項(xiàng)及湍流黏性系數(shù)均采用一階迎風(fēng)格式離散[15-16]。當(dāng)殘差小于10-4，并且流場(chǎng)中監(jiān)測(cè)點(diǎn)參數(shù)波動(dòng)小于3%時(shí)，認(rèn)定計(jì)算收斂。

3 模擬驗(yàn)證

需要對(duì)數(shù)值模擬得到的計(jì)算結(jié)果進(jìn)行準(zhǔn)確性驗(yàn)證，根據(jù)電廠提供的設(shè)計(jì)參數(shù)，對(duì)引風(fēng)機(jī)動(dòng)葉葉頂安裝角為32.0°時(shí)的性能參數(shù)進(jìn)行模擬驗(yàn)證。圖3為引風(fēng)機(jī)運(yùn)行特性曲線。

圖3 引風(fēng)機(jī)運(yùn)行特性曲線

由圖3可得：數(shù)值模擬得到的曲線與廠家提供的曲線的變化趨勢(shì)基本一致，并且兩者較為吻合。但是，在模擬中，建立模型與選擇數(shù)值計(jì)算方法時(shí)均進(jìn)行了一定程度的簡(jiǎn)化，因此數(shù)值模擬計(jì)算結(jié)果和實(shí)驗(yàn)結(jié)果之間存在一定偏差。在BMCR工況下，數(shù)值模擬計(jì)算得到的風(fēng)機(jī)全壓為7 620 Pa，效率為85%，而實(shí)驗(yàn)中風(fēng)機(jī)全壓為7 350 Pa，效率為88.8%，風(fēng)機(jī)全壓和效率的模擬值與實(shí)驗(yàn)值的平均相對(duì)誤差分別為3.6%與4.3%，誤差較小。因此，認(rèn)為模擬結(jié)果能較準(zhǔn)確地反映引風(fēng)機(jī)的實(shí)際運(yùn)行性能，可采用該數(shù)值計(jì)算方法對(duì)引風(fēng)機(jī)運(yùn)行特性及內(nèi)部流場(chǎng)特征進(jìn)行分析研究。

4 結(jié)果分析

4.1 引風(fēng)機(jī)運(yùn)行特性與內(nèi)部流場(chǎng)特征分析

圖4為數(shù)值模擬得到的引風(fēng)機(jī)運(yùn)行特性曲線。由圖4可得：全煙道系統(tǒng)與引風(fēng)機(jī)本體的風(fēng)機(jī)運(yùn)行特性曲線的趨勢(shì)一致。在一定的動(dòng)葉葉頂安裝角下，全壓隨體積流量的增大而減小，效率隨體積流量的增大先增大后減?。粚?duì)比全煙道系統(tǒng)與引風(fēng)機(jī)本體的風(fēng)機(jī)運(yùn)行特性曲線，全煙道系統(tǒng)在小流量工況下風(fēng)機(jī)的全壓與效率均有所降低，而全煙道系統(tǒng)在大流量工況下風(fēng)機(jī)的全壓與效率幾乎重合。這說明引風(fēng)機(jī)進(jìn)出口煙道對(duì)于小流量工況下的風(fēng)機(jī)性能具有一定的影響。

圖4 引風(fēng)機(jī)運(yùn)行特性曲線對(duì)比

為研究引風(fēng)機(jī)本體及全煙道系統(tǒng)的流場(chǎng)特征，選取引風(fēng)機(jī)部分截面進(jìn)行分析[17-18]，各特征截面位置見圖5，表2為各特征截面的軸向速度不均勻系數(shù)。

圖5 引風(fēng)機(jī)本體及全煙道系統(tǒng)的特征截面位置

表2 各特征截面軸向速度不均勻系數(shù)

圖6為BMCR工況下集流器中間截面的軸向速度分布云圖。由圖6可得：對(duì)于引風(fēng)機(jī)本體與全煙道系統(tǒng)，集流器中間截面的軸向速度分布呈上下相反趨勢(shì)。對(duì)于引風(fēng)機(jī)本體，氣流經(jīng)進(jìn)氣箱90°折角彎頭后，受到一定的慣性作用，在下側(cè)有所積聚，因此軸向速度較高。對(duì)于全煙道系統(tǒng)，集流器中間截面軸向速度在上、下兩側(cè)分布不均勻，上側(cè)近機(jī)殼的大部分區(qū)域軸向速度較高，這主要是因?yàn)闅饬魇艿揭L(fēng)機(jī)進(jìn)口煙道紊亂流場(chǎng)分布與葉輪旋轉(zhuǎn)作用的雙重影響，流場(chǎng)產(chǎn)生相互干涉。軸向速度不均勻系數(shù)從引風(fēng)機(jī)本體的4.53%增加至全煙道系統(tǒng)的5.92%，軸向速度分布均勻性降低。

圖6 集流器中間截面的軸向速度分布云圖

圖7為BMCR工況下引風(fēng)機(jī)出口截面的軸向速度分布云圖。由圖7可得：與引風(fēng)機(jī)本體相比，全煙道系統(tǒng)的引風(fēng)機(jī)出口截面的軸向速度分布更不均勻，軸向速度不均勻系數(shù)從引風(fēng)機(jī)本體的56.76%增加至全煙道系統(tǒng)的76.28%。全煙道系統(tǒng)沿周向近機(jī)殼區(qū)域存在大面積的低速渦流區(qū)，氣流充滿度不高，主要集中在旋轉(zhuǎn)軸附近；同時(shí)，氣流先經(jīng)過動(dòng)葉輪旋轉(zhuǎn)加速，雖然后續(xù)經(jīng)過導(dǎo)葉輪的整流作用，但是仍存在較大速度梯度，使得擴(kuò)壓器區(qū)域的氣流沿軸心有所偏轉(zhuǎn)，在A側(cè)存在部分高速區(qū)。綜上可知，引風(fēng)機(jī)內(nèi)部流場(chǎng)特征受引風(fēng)機(jī)進(jìn)出口煙道流場(chǎng)分布的影響，在對(duì)引風(fēng)機(jī)內(nèi)部流場(chǎng)特征進(jìn)行分析時(shí)，應(yīng)充分考慮引風(fēng)機(jī)及其進(jìn)出口煙道流場(chǎng)的干涉作用。

圖7 引風(fēng)機(jī)出口截面的軸向速度分布云圖

4.2 進(jìn)出口煙道流場(chǎng)分析

通過模擬BMCR工況下的無引風(fēng)機(jī)進(jìn)出口煙道和全煙道系統(tǒng)的流場(chǎng)，分析引風(fēng)機(jī)內(nèi)部流動(dòng)特征與進(jìn)出口煙道流場(chǎng)間的干涉。圖8、圖9分別為無引風(fēng)機(jī)進(jìn)口、出口煙道流線，圖10為全煙道系統(tǒng)流線。

圖8 無引風(fēng)機(jī)進(jìn)口煙道流線圖

圖9 無引風(fēng)機(jī)出口煙道流線

圖10 全煙道系統(tǒng)流線

對(duì)比圖8與圖10可得：在不考慮引風(fēng)機(jī)的情況下，引風(fēng)機(jī)進(jìn)口煙道除塵器出口兩側(cè)氣流對(duì)流匯沖，在水平煙道形成大范圍的旋渦區(qū)，流場(chǎng)較為紊亂，導(dǎo)致進(jìn)出豎直煙道的氣流相互纏繞，靠近后墻側(cè)的氣流充滿度較低(見圖8中圓圈部分)；對(duì)于全煙道系統(tǒng)，引風(fēng)機(jī)進(jìn)口煙道經(jīng)對(duì)流匯沖后的流場(chǎng)依然較為紊亂，但是水平煙道的旋渦區(qū)明顯減小，豎直煙道的氣流充滿度得到提高。因此，在無引風(fēng)機(jī)進(jìn)口煙道流場(chǎng)模擬中，不考慮后續(xù)引風(fēng)機(jī)的模擬結(jié)果比實(shí)際流場(chǎng)分布情況更差，煙道與引風(fēng)機(jī)串聯(lián)后，進(jìn)口煙道的流場(chǎng)得到一定程度的改善。

對(duì)比圖9與圖10可得：在不考慮引風(fēng)機(jī)的情況下，假設(shè)引風(fēng)機(jī)出口氣流分布均勻，引風(fēng)機(jī)出口煙道氣流充滿度較高，流線較為平滑流暢；而對(duì)于全煙道系統(tǒng)，氣流經(jīng)引風(fēng)機(jī)動(dòng)葉輪的旋轉(zhuǎn)加速作用到擴(kuò)壓器出口，流場(chǎng)的不均勻性得到進(jìn)一步發(fā)展，經(jīng)90°折角彎頭后氣流主要集中在A側(cè)煙道，B側(cè)煙道存在部分空區(qū)，出口煙道氣流充滿度低，形成大面積低速旋渦區(qū)(見圖10中圓圈部分)。因此，在實(shí)際運(yùn)行中，引風(fēng)機(jī)內(nèi)部流場(chǎng)對(duì)出口煙道具有較大的影響，在對(duì)引風(fēng)機(jī)出口煙道流場(chǎng)進(jìn)行優(yōu)化改造時(shí)，應(yīng)充分考慮引風(fēng)機(jī)出口速度的分布情況。

為研究進(jìn)出口煙道的流場(chǎng)特性，選取進(jìn)出口煙道內(nèi)部分截面進(jìn)行分析[19-20]，各特征截面位置見圖11。

圖11 特征截面位置

圖12為進(jìn)氣箱進(jìn)口截面的軸向速度分布云圖。由圖12可得：無引風(fēng)機(jī)進(jìn)口煙道數(shù)值模擬時(shí)，進(jìn)氣箱進(jìn)口截面的軸向速度分布不均勻，軸向速度不均勻系數(shù)為29.68%，在后墻側(cè)軸向速度相對(duì)較小，這主要是因?yàn)樨Q直煙道后墻側(cè)氣流的充滿度不高；對(duì)于全煙道系統(tǒng)，進(jìn)氣箱進(jìn)口截面A、B側(cè)的軸向速度呈對(duì)稱分布，從后墻側(cè)到前墻側(cè)呈梯度減小。氣流經(jīng)進(jìn)氣箱進(jìn)入水平集流器時(shí)，由于慣性作用，在后墻側(cè)中間形成局部高速區(qū)；最高軸向速度達(dá)到50 m/s，軸向速度不均勻系數(shù)為19.21%，速度分布均勻性有所提高。全煙道系統(tǒng)的引風(fēng)機(jī)內(nèi)部流場(chǎng)特征對(duì)引風(fēng)機(jī)進(jìn)口煙道的流場(chǎng)具有一定的改善作用，在引風(fēng)機(jī)進(jìn)口煙道流場(chǎng)進(jìn)行優(yōu)化改造時(shí)，忽略引風(fēng)機(jī)的影響會(huì)使計(jì)算結(jié)果與實(shí)際情況存在一定偏差。

圖12 進(jìn)氣箱進(jìn)口截面的軸向速度分布云圖

圖13為出口煙道截面軸向速度分布云圖。由圖13可得：出口煙道截面軸向速度分布極不均勻，在煙道壁面附近存在大面積的空區(qū)，氣流充滿度較低。對(duì)于無引風(fēng)機(jī)煙道，氣流在煙道中部分布得較為均勻，但是在煙道四角存在空區(qū)。對(duì)于全煙道系統(tǒng)煙道，氣流主要集中在煙道中間近A側(cè)區(qū)域，B側(cè)空區(qū)范圍擴(kuò)大，并且在靠近煙道下側(cè)存在局部高速區(qū)，一方面是受引風(fēng)機(jī)出口速度分布的影響，另一方面是因?yàn)槌隹跓煹乐袣饬鲝乃綗煹懒鹘?jīng)90°折角彎頭進(jìn)入豎直煙道爬升時(shí)，慣性作用造成主流煙氣與上壁面分離，氣流主要分布在煙道下側(cè)。引風(fēng)機(jī)出口軸向速度分布不均勻，會(huì)使不均勻氣流向煙道下游繼續(xù)發(fā)展，進(jìn)而擾亂引風(fēng)機(jī)出口煙道的流場(chǎng)。

圖13 出口煙道截面的軸向速度分布云圖

因此，在對(duì)引風(fēng)機(jī)出口煙道進(jìn)行優(yōu)化改造時(shí)，假設(shè)模型進(jìn)口速度分布均勻是不合理的，需要充分考慮引風(fēng)機(jī)內(nèi)部流場(chǎng)對(duì)出口煙道流場(chǎng)的影響。

計(jì)算無引風(fēng)機(jī)和全煙道系統(tǒng)的進(jìn)出口煙道的壓力損失，得到的結(jié)果見表3。全煙道系統(tǒng)的進(jìn)口煙道的壓力損失沒有明顯變化，而出口煙道壓力損失大幅增加，由219.75 Pa增加至421.21 Pa，說明引風(fēng)機(jī)內(nèi)部流場(chǎng)對(duì)進(jìn)口煙道的影響較小，但對(duì)出口煙道會(huì)產(chǎn)生顯著的影響。

表3 進(jìn)出口煙道壓力損失

5 結(jié)語

(1) 引風(fēng)機(jī)的內(nèi)部流場(chǎng)與其進(jìn)出口煙道的流場(chǎng)相互影響。小流量工況下，相比于引風(fēng)機(jī)本體，全煙道系統(tǒng)下全壓與效率均有所降低；引風(fēng)機(jī)及其進(jìn)出口煙道的流場(chǎng)干涉作用會(huì)在一定程度上降低引風(fēng)機(jī)內(nèi)部流場(chǎng)的均勻性。

(2) BMCR工況下，引風(fēng)機(jī)內(nèi)部流場(chǎng)特征有利于改善進(jìn)口煙道流場(chǎng)的均勻性，但是會(huì)加劇出口煙道的能量損耗，擾亂出口煙道的流場(chǎng)。因此，在對(duì)引風(fēng)機(jī)進(jìn)出口煙道流場(chǎng)進(jìn)行分析時(shí)，應(yīng)充分考慮引風(fēng)機(jī)內(nèi)部流場(chǎng)的影響。