一次風(fēng)速對(duì)預(yù)燃室式煤粉旋流燃燒器性能的影響

張 旭 芳

(1.煤科院節(jié)能技術(shù)有限公司，北京 100013； 2.煤炭資源高效開采與潔凈利用國家重點(diǎn)實(shí)驗(yàn)室，北京 100013；3.國家能源煤炭高效利用與節(jié)能減排技術(shù)裝備重點(diǎn)實(shí)驗(yàn)室，北京 100013)

0 引 言

煤粉工業(yè)鍋爐承擔(dān)著我國大部分地區(qū)工業(yè)用蒸汽及居民采暖的重任。據(jù)統(tǒng)計(jì)，我國目前使用的工業(yè)鍋爐數(shù)量約60余萬臺(tái)，其中80%以上均為燃煤鍋爐。煤粉工業(yè)鍋爐具有燃燒效率高、污染物排放水平低的優(yōu)勢(shì)，逐漸淘汰了落后的燃煤鏈條爐，成為工業(yè)鍋爐的“主力軍”[1-2]。煤粉燃燒器是煤粉工業(yè)鍋爐燃燒組織的核心設(shè)備，其結(jié)構(gòu)設(shè)計(jì)及運(yùn)行參數(shù)對(duì)于煤粉穩(wěn)定燃燒和高效燃盡具有重要的意義[3-4]。

由于煤質(zhì)及運(yùn)行負(fù)荷的變化，煤粉工業(yè)鍋爐易出現(xiàn)燃燒不穩(wěn)定的問題。為了強(qiáng)化煤粉燃燒火焰穩(wěn)定性，預(yù)燃室式煤粉燃燒器被廣泛研究并投入生產(chǎn)實(shí)踐中，目前已取得良好的穩(wěn)燃效果。由煤科院節(jié)能技術(shù)有限公司研發(fā)的雙錐逆噴旋流煤粉燃燒器具有濃相供粉、點(diǎn)火迅速、燃燒效率高的特點(diǎn)，已經(jīng)成功應(yīng)用于多種容量的煤粉工業(yè)鍋爐，獲得了較好的經(jīng)濟(jì)效益[5-6]；雙錐逆噴旋流煤粉燃燒器由一次風(fēng)管、點(diǎn)火器、回流帽、穩(wěn)燃腔以及加速腔構(gòu)成[7-8]。清華大學(xué)開發(fā)了1款煤粉濃縮預(yù)熱低NOx燃燒器(簡稱PRP)，通過卷吸爐內(nèi)高溫?zé)煔馐節(jié)饷悍蹥饬鬟M(jìn)入爐膛前被加熱至接近著火點(diǎn)的高溫狀態(tài)，同時(shí)保證燃燒前期處于還原性氣氛，實(shí)現(xiàn)了與傳統(tǒng)高溫空氣燃燒同等的氮氧化物排放水平。預(yù)熱室是PRP燃燒器的關(guān)鍵部件，由于偏心的高速濃粉氣流引射作用，預(yù)熱室上部空間壓力小于爐膛，造成爐內(nèi)熱煙氣回流至預(yù)熱室內(nèi)，使得流經(jīng)預(yù)熱室的濃相煤粉氣流與熱煙氣強(qiáng)烈混合后迅速升溫，實(shí)現(xiàn)煤粉的迅速升溫及穩(wěn)定著火。張海等[9]在12 MW試驗(yàn)鍋爐上研究了旋流式PRP燃燒器的穩(wěn)燃性能，燃燒器燃煤量為1 310 kg/h，一次風(fēng)溫90 ℃～175 ℃，研究結(jié)果表明，一次風(fēng)速為20 m/s時(shí)則預(yù)熱室內(nèi)升溫速率約為24 000 ℃/s，很好地實(shí)現(xiàn)了煤粉氣流的快速升溫，大幅改善低揮發(fā)煤著火穩(wěn)燃特性。朱霖等[10]提出了1種適用于煤粉工業(yè)鍋爐的預(yù)燃室式燃燒器，該燃燒器預(yù)燃室結(jié)構(gòu)與上述清華大學(xué)開發(fā)的燃燒器較為相似，不同點(diǎn)在于該燃燒器一次風(fēng)采用直流，且在一次風(fēng)管口設(shè)置有1個(gè)圓錐形鈍體用于穩(wěn)定火焰，冷態(tài)實(shí)驗(yàn)及熱態(tài)工業(yè)試驗(yàn)均驗(yàn)證了該燃燒器良好的火焰穩(wěn)定特性。王進(jìn)卿等[11]采用數(shù)值計(jì)算軟件對(duì)1種新開發(fā)的采用鈍體的預(yù)燃室燃燒器進(jìn)行數(shù)值計(jì)算，針對(duì)一次風(fēng)速、二次風(fēng)速以及二次風(fēng)旋流強(qiáng)度等對(duì)回流強(qiáng)度的影響進(jìn)行研究。煤粉工業(yè)鍋爐大多采用中儲(chǔ)式煤粉供料系統(tǒng)，一次風(fēng)中煤粉濃度可調(diào)，對(duì)于縮短著火時(shí)間、適應(yīng)不同煤種的穩(wěn)定燃燒具有重要的作用。煤粉濃度越高煤粉著火所需的點(diǎn)火熱低，但風(fēng)粉混合以及燃盡性能會(huì)變差，因此存在1個(gè)較佳的一次風(fēng)速，其對(duì)于煤粉工業(yè)鍋爐的安全、經(jīng)濟(jì)、穩(wěn)定運(yùn)行至關(guān)重要[12-14]。一次風(fēng)動(dòng)量過大或過小都不利于煤粉在預(yù)燃室內(nèi)穩(wěn)定著火及高效燃燒；另外，一次風(fēng)動(dòng)量過大還可能破壞預(yù)燃室內(nèi)形成的回流區(qū)，降低旋流火焰穩(wěn)定性。

以下詳細(xì)介紹以預(yù)燃室式煤粉燃燒器為基礎(chǔ)而新設(shè)計(jì)的1款預(yù)燃室式煤粉旋流燃燒器，并針對(duì)其相應(yīng)的數(shù)值計(jì)算進(jìn)行分析，通過改變一次風(fēng)速以分析預(yù)燃室內(nèi)流場(chǎng)、溫度及組分濃度場(chǎng)，探究一次風(fēng)動(dòng)量對(duì)預(yù)燃室式燃燒器內(nèi)流場(chǎng)以及燃燒特性的影響，對(duì)于揭示中心給粉預(yù)燃室式煤粉旋流燃燒器一次風(fēng)風(fēng)速對(duì)空氣動(dòng)力場(chǎng)和燃燒性能的影響具有一定的積極意義。

1 燃燒器結(jié)構(gòu)及網(wǎng)格劃分

預(yù)燃室式煤粉旋流燃燒器結(jié)構(gòu)如圖1所示。該燃燒器的工作原理簡述如下：煤粉空氣混合物由一次風(fēng)管1進(jìn)入；二次風(fēng)由風(fēng)道2進(jìn)入，經(jīng)過軸向葉片3進(jìn)入圓錐形預(yù)燃室前腔4；三次風(fēng)由三次風(fēng)入口5處進(jìn)入圓筒形預(yù)燃室6，燃燒后產(chǎn)生的熱煙氣由預(yù)燃室出口7噴入爐膛。

圖1 預(yù)燃室式煤粉旋流燃燒器結(jié)構(gòu)

預(yù)燃室式煤粉旋流燃燒器的設(shè)計(jì)思路：二次風(fēng)進(jìn)口處圓錐形進(jìn)口端蓋的存在，使得受限的旋流二次風(fēng)被引向圓筒形預(yù)燃室邊壁，從而在預(yù)燃室中央形成1個(gè)較大的回流區(qū)；同時(shí)，合理的一次風(fēng)粉噴射速度確保了煤粉在預(yù)燃室內(nèi)受熱著火階段的停留時(shí)間，同時(shí)，大量的高溫?zé)煔庥肯蛞淮物L(fēng)粉的根部，一次風(fēng)粉直接進(jìn)入回流區(qū)核心區(qū)域與高溫?zé)煔饣旌?，使得煤粉著火穩(wěn)定性大幅增強(qiáng)。

為了探究一次風(fēng)速對(duì)該燃燒器預(yù)燃室內(nèi)氣流動(dòng)力場(chǎng)特性及燃燒性能的影響，建立了1∶1燃燒器結(jié)構(gòu)模型，并采用ICEM軟件對(duì)其進(jìn)行網(wǎng)格劃分。預(yù)燃室式煤粉旋流燃燒器的物理結(jié)構(gòu)模型及對(duì)應(yīng)的計(jì)算域網(wǎng)格如圖2所示。

圖2 燃燒器的物理結(jié)構(gòu)模型及計(jì)算域網(wǎng)格劃分示意

2 煤粉燃燒計(jì)算模型及網(wǎng)格無關(guān)性驗(yàn)證

煤粉燃燒過程伴隨著強(qiáng)烈的質(zhì)量傳遞、能量傳遞、動(dòng)量傳遞及復(fù)雜化學(xué)反應(yīng)。數(shù)值計(jì)算作為1種重要的研究煤粉燃燒的方法，已得到廣泛的應(yīng)用。煤粉燃燒數(shù)值計(jì)算需要確定合適的計(jì)算模型，主要包括煤粉燃燒數(shù)學(xué)模型和煤粉燃燒化學(xué)反應(yīng)模型。

2.1 煤粉燃燒數(shù)學(xué)模型

流體流動(dòng)過程在歐拉坐標(biāo)系下進(jìn)行求解。連續(xù)相方程為質(zhì)量、動(dòng)量、能量的連續(xù)性方程和時(shí)間均值的納維斯托克斯方程。同時(shí)，在連續(xù)相方程中添加組分輸運(yùn)方程以求解反應(yīng)流。

湍流數(shù)值模擬可以分為直接數(shù)值模擬和非直接數(shù)值模擬，其中直接數(shù)值模擬對(duì)計(jì)算機(jī)要求很高，且運(yùn)算量大、成本高，很少被采用。非直接數(shù)值模擬方法被廣泛使用，其中以大渦模擬方法和雷諾平均模擬方法為典型代表[15]。此次研究采用Realizablek-ε雙方程湍流模型。

基于FLUENT模擬軟件，當(dāng)光學(xué)深度>1，可選用P1和Rossland模型。P1模型被證明適合用于顆粒燃料燃燒時(shí)遠(yuǎn)離火焰周邊的研究，計(jì)算效率較高。因此筆者選擇P1輻射模型，煤粉顆粒使用隨機(jī)顆粒軌道模型進(jìn)行模擬。

2.2 煤粉燃燒化學(xué)反應(yīng)模型

煤粉燃燒均相反應(yīng)在高溫、高速的湍流化學(xué)反應(yīng)條件下，此次選擇組分輸運(yùn)模型(Species Transport Model)研究燃燒器及爐膛內(nèi)的化學(xué)反應(yīng)過程。組分輸運(yùn)模型主要包括渦耗散模型(Eddy-dissipation model)、有限速率模型(Finite Rate Kinetic Model)、有限速率/渦耗散模型(Finite Rate/Eddy Dissipation Model)3種化學(xué)反應(yīng)模型。其中，渦耗散模型假設(shè)氣相化學(xué)反應(yīng)時(shí)間尺度遠(yuǎn)小于湍流作用下氣相物質(zhì)之間的混合時(shí)間尺度，認(rèn)為氣相化學(xué)反應(yīng)動(dòng)力學(xué)控制步驟為擴(kuò)散；有限速率模型考慮溫度對(duì)化學(xué)反應(yīng)的影響，但忽略湍流對(duì)化學(xué)反應(yīng)的影響；有限速率/渦耗散模型分別基于層流有限速率模型和渦耗散模型發(fā)展而來。該模型按實(shí)際情況，對(duì)比分別由有限速率和渦耗散計(jì)算而得的2個(gè)反應(yīng)速率，選用小者作為反應(yīng)速率?；诖?，此次研究采用組分輸運(yùn)模型計(jì)算化學(xué)反應(yīng)過程，揮發(fā)份析出時(shí)為雙反應(yīng)競(jìng)爭(zhēng)模型，揮發(fā)份燃燒時(shí)的氣相反應(yīng)采用有限速率/渦耗散模型(Finite-Rate/Eddy-Dissipation Model)。

2.3 網(wǎng)格無關(guān)性驗(yàn)證

為了排除網(wǎng)格數(shù)量對(duì)計(jì)算結(jié)果的影響，須開展網(wǎng)格無關(guān)性驗(yàn)證，將網(wǎng)格數(shù)量劃分為10萬、30萬以及60萬3種網(wǎng)格。開展冷態(tài)數(shù)值計(jì)算檢驗(yàn)的方法，得到不同網(wǎng)格數(shù)量下燃燒器中心軸線上20個(gè)監(jiān)測(cè)點(diǎn)上的速度分布結(jié)果，如圖3所示。結(jié)果表明，網(wǎng)格數(shù)量大于30萬后，燃燒器中心軸線上速度分布趨勢(shì)保持不變。因此，同時(shí)兼顧計(jì)算精度和計(jì)算效率，選用數(shù)量為30萬的網(wǎng)格開展數(shù)值計(jì)算。

記憶中，在我4歲那年，由中國國際電視總公司出品了一部41集古裝神話劇——改編自明代小說家吳承恩同名文學(xué)古典名著《西游記》。1986年春節(jié)一經(jīng)播出便轟動(dòng)全國，可謂老少皆宜，獲得了極高評(píng)價(jià)。至今仍是寒暑假期間被重播最多的電視劇之一，百看不厭，成為一部公認(rèn)的、無法超越的經(jīng)典。該劇講述的是孫悟空、豬八戒、沙僧輔保大唐高僧玄奘(唐僧)去西天取經(jīng)的故事，師徒四人一路搶灘涉險(xiǎn)，降妖伏怪，歷經(jīng)九九八十一難，取回真經(jīng)，終修成正果的故事。

圖3 不同網(wǎng)格數(shù)計(jì)算結(jié)果對(duì)比

3 數(shù)值模擬工況與煤質(zhì)參數(shù)

筆者數(shù)值模擬所采用的煤種與模型準(zhǔn)確性驗(yàn)證試驗(yàn)中所使用的相同，所有模擬工況均在煤粉供料量為 700 kg/h的條件下。為了探究一次風(fēng)動(dòng)量的影響，保持二次風(fēng)風(fēng)量、三次風(fēng)風(fēng)量以及二次風(fēng)旋流數(shù)不變，開展3個(gè)工況下的數(shù)值計(jì)算，工況主要參數(shù)見表1。

表1 數(shù)值模擬工況主要參數(shù)

筆者研究所使用的煤種為神木煙煤，其工業(yè)及元素分析等煤質(zhì)參數(shù)見表2。

表2 神木煙煤的煤質(zhì)分析結(jié)果

4 結(jié)果與討論

4.1 速度場(chǎng)分析

預(yù)燃室中心截面軸向速度(u)分布云圖如圖4所示，可知該燃燒器的流場(chǎng)特點(diǎn)如下：

圖4 當(dāng)一次風(fēng)速v=12.43 m/s時(shí)中心截面軸向速度分布云圖

(1)一次風(fēng)粉以直流形式噴入預(yù)燃室內(nèi)，在預(yù)燃室中心形成具有一定剛性的直流氣流；二次風(fēng)為旋轉(zhuǎn)氣流，旋流數(shù)為2.16，該值大于0.6，因此二次風(fēng)屬強(qiáng)旋流，在預(yù)燃室內(nèi)形成明顯的回流區(qū)。

(2)當(dāng)一次風(fēng)速v=12.43 m/s時(shí)，預(yù)燃室內(nèi)形成了1個(gè)較大的且為整體的中心回流區(qū)，當(dāng)一次風(fēng)速v=18.64 m/s和24.86 m/s時(shí)，一次風(fēng)穿透原中心回流區(qū)，形成2個(gè)對(duì)稱的環(huán)形回流區(qū)，且一次風(fēng)速越大，將2個(gè)環(huán)形回流區(qū)“擠壓”的越小；

(3)三次風(fēng)以直流風(fēng)的形式緊貼預(yù)燃室壁面進(jìn)入，在圓柱段預(yù)燃室壁面附近形成了1個(gè)高速的冷卻空氣層，可有效避免預(yù)燃室壁面高溫、沾污以及腐蝕的發(fā)生。

中心軸線上的軸向速度分布對(duì)比曲線如圖5所示。

圖5 中心軸線上的軸向速度分布

分析圖5可知，一次風(fēng)速v=12.43 m/s時(shí)，中心軸線上速度先基本保持不變，隨后逐漸下降，直至軸向距離為0.33 m時(shí)，軸向速度衰減為0，此時(shí)預(yù)燃室中心均為由旋流產(chǎn)生的“逆流”氣流，即進(jìn)入回流區(qū)。而當(dāng)一次風(fēng)速v=18.64 m/s和24.86 m/s時(shí)，在預(yù)燃室中心軸線上未出現(xiàn)軸向速度小于0的情況，即中心回流區(qū)被一次風(fēng)穿透而消逝。

預(yù)燃室中心截面的軸向速度矢量圖如圖6所示。對(duì)比預(yù)燃室內(nèi)不同一次風(fēng)速下的流場(chǎng)分布可知，一次風(fēng)速較小而未穿透中心回流區(qū)時(shí)，預(yù)燃室內(nèi)回流區(qū)的面積顯著大于其他2種情況，即隨著一次風(fēng)速增大，環(huán)形回流區(qū)向壁面附近有所移動(dòng)；由于一次風(fēng)的風(fēng)速高、風(fēng)量大，導(dǎo)致一次風(fēng)粉的動(dòng)量顯著增大，對(duì)環(huán)形回流區(qū)造成明顯的“擠壓”現(xiàn)象。

4.2 溫度場(chǎng)分析

預(yù)燃室中心截面上溫度分布云圖如圖7所示。由圖7可知，一次風(fēng)粉以常溫(300 K)噴入燃燒器中，其在預(yù)燃室內(nèi)經(jīng)歷了預(yù)熱、燃燒的過程，隨一次風(fēng)速增大至v=24.86 m/s時(shí)，預(yù)燃室內(nèi)中心高溫區(qū)徹底被一次風(fēng)穿透，形成了中心低、外側(cè)高的分布，且高溫區(qū)有向前移動(dòng)的趨勢(shì)。

圖7 中心截面溫度分布云圖

中心軸線上溫度分布如圖8所示。

圖8 中心軸線上的溫度分布

分析圖8中心軸線上溫度分布曲線可知：v=12.43 m/s、中心軸線上x≤0.23 m時(shí)，溫度基本保持不變，隨后溫度迅速上升直至x=0.42 m時(shí)達(dá)到最高并基本保持不變。一次風(fēng)速v=18.64 m/s和24.86 m/s時(shí)，隨著風(fēng)速增大，中心軸線上溫升更慢，且溫升速率也較低，最終預(yù)燃室燃燒器出口中心點(diǎn)上一次風(fēng)速越大溫度越低。主要由于隨著中心回流區(qū)的破壞，靠近中心軸線上的高溫回流煙氣占比降低，且一次風(fēng)粉動(dòng)量越強(qiáng)，氣流剛性越大，極大地限制了一次風(fēng)粉與周圍煙氣的混合與對(duì)流傳熱，因此溫升速率低，溫度相對(duì)較低。

考察預(yù)燃室燃燒器出口平均溫度以及焦炭質(zhì)量流率，可以判斷煤粉在預(yù)燃室內(nèi)著火以及燃盡情況，數(shù)值計(jì)算結(jié)果見表3。由表3可知，3種不同的一次風(fēng)速下，預(yù)燃室出口截面平均溫度均大于1 200 K，可實(shí)現(xiàn)煤粉的穩(wěn)定著火。煤粉進(jìn)入預(yù)燃室初始質(zhì)量流率是700 kg/h(即0.194 kg/s)，根據(jù)預(yù)燃室出口截面焦炭質(zhì)量流率推算焦炭轉(zhuǎn)化率，3種情況下焦炭的轉(zhuǎn)率均高于99.99%，煤粉在該預(yù)燃室式燃燒器內(nèi)可以實(shí)現(xiàn)高效燃燒。

表3 預(yù)燃室燃燒器的數(shù)值計(jì)算結(jié)果

4.3 組分場(chǎng)分析

預(yù)燃室內(nèi)組分分布可直觀地反映燃燒特性，通過對(duì)關(guān)鍵組分的分析對(duì)整體燃燒情況進(jìn)行預(yù)測(cè)。預(yù)燃室中心截面揮發(fā)分質(zhì)量分?jǐn)?shù)分布如圖9所示。

圖9 中心截面揮發(fā)分濃度分布圖(體積分?jǐn)?shù)/%)

由圖9可知，隨一次風(fēng)速的增加，揮發(fā)分在預(yù)燃室內(nèi)分布區(qū)域出現(xiàn)逐漸后移，且高濃度區(qū)域面積逐漸增大的趨勢(shì)。一次風(fēng)速增大，一次風(fēng)與煤粉的預(yù)熱速度變緩，煤粉受熱析出的揮發(fā)分推遲是導(dǎo)致后移以及后區(qū)濃度較高的主要原因。因此，一次風(fēng)速增大不利于煤粉與周圍高溫?zé)煔獾幕旌希档兔悍蹨厣俣?、延遲揮發(fā)分的析出可能會(huì)造成煤粉整體燃盡效果變差。

由預(yù)燃室內(nèi)氧氣濃度分布可直觀地判斷預(yù)燃室內(nèi)貧燃以及富燃區(qū)域，此舉對(duì)于控制污染物NO的排放具有重要的意義。在貧燃區(qū)域，過量空氣系數(shù)高于1，揮發(fā)分氮以及焦炭氮容易通過氧化反應(yīng)生成大量的NO，導(dǎo)致較高的NO排放。而合理地控制預(yù)燃室內(nèi)過量空氣系數(shù)，利于創(chuàng)造適宜的還原性氣氛，揮發(fā)分中NO的前驅(qū)物HCN和NH3等更傾向于發(fā)生還原反應(yīng)并生成N2，進(jìn)而降低NO的生成。同時(shí)，在強(qiáng)還原性氣氛下，焦炭表面還存在較強(qiáng)的活性位點(diǎn)，通過非均相還原反應(yīng)將已生成的NO也可還原成N2，對(duì)于降低NO排放均具有重要的意義。預(yù)燃室中心軸線上O2體積分?jǐn)?shù)分布曲線如圖10所示。

圖10 中心軸線上氧氣體積分?jǐn)?shù)分布

由圖10可知，隨一次風(fēng)速的增加，中心軸線上氧氣開始顯著下降的位置顯著地向后平移，且氧氣濃度衰減的速率也逐漸放緩，一次風(fēng)速v=12.43 m/s時(shí)在預(yù)燃室內(nèi)形成了最大的低氧氛圍，其較強(qiáng)的還原性氣氛對(duì)于控制NO的生成具有重要意義。

4 結(jié) 論

(1)一次風(fēng)速過高則會(huì)穿透預(yù)燃室內(nèi)的中心回流區(qū)，從而形成2個(gè)對(duì)稱的環(huán)形回流區(qū)，且一次風(fēng)速越大，其將2個(gè)環(huán)形回流區(qū)“擠壓”得越小。

(2)隨一次風(fēng)速的增大，中心軸線上溫升速率越小，一次風(fēng)粉氣流與旋流二次風(fēng)的混合過程被削弱，限制了其與周圍高溫?zé)煔獾膫鳠?。一次風(fēng)速過大，預(yù)燃室內(nèi)中心高溫區(qū)被穿透，形成了中心低、外側(cè)高的分布，且高溫區(qū)有向前移動(dòng)的趨勢(shì)。

(3)一次風(fēng)速增大，煤粉溫升速率顯著降低，揮發(fā)分的析出存在明顯的延遲，可能會(huì)造成煤粉整體燃盡效果變差；隨一次風(fēng)速的增加，中心軸線上氧氣開始顯著下降的位置明顯地向后平移，且氧氣濃度衰減的速率也逐漸放緩，預(yù)燃室內(nèi)較強(qiáng)的還原性氣氛區(qū)域減小，不利于控制NO的生成。