一次風速對預燃室式煤粉旋流燃燒器性能的影響

張 旭 芳

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

0 引 言

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

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

以下詳細介紹以預燃室式煤粉燃燒器為基礎而新設計的1款預燃室式煤粉旋流燃燒器，并針對其相應的數值計算進行分析，通過改變一次風速以分析預燃室內流場、溫度及組分濃度場，探究一次風動量對預燃室式燃燒器內流場以及燃燒特性的影響，對于揭示中心給粉預燃室式煤粉旋流燃燒器一次風風速對空氣動力場和燃燒性能的影響具有一定的積極意義。

1 燃燒器結構及網格劃分

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

圖1 預燃室式煤粉旋流燃燒器結構

預燃室式煤粉旋流燃燒器的設計思路：二次風進口處圓錐形進口端蓋的存在，使得受限的旋流二次風被引向圓筒形預燃室邊壁，從而在預燃室中央形成1個較大的回流區；同時，合理的一次風粉噴射速度確保了煤粉在預燃室內受熱著火階段的停留時間，同時，大量的高溫煙氣涌向一次風粉的根部，一次風粉直接進入回流區核心區域與高溫煙氣混合，使得煤粉著火穩定性大幅增強。

為了探究一次風速對該燃燒器預燃室內氣流動力場特性及燃燒性能的影響，建立了1∶1燃燒器結構模型，并采用ICEM軟件對其進行網格劃分。預燃室式煤粉旋流燃燒器的物理結構模型及對應的計算域網格如圖2所示。

圖2 燃燒器的物理結構模型及計算域網格劃分示意

2 煤粉燃燒計算模型及網格無關性驗證

煤粉燃燒過程伴隨著強烈的質量傳遞、能量傳遞、動量傳遞及復雜化學反應。數值計算作為1種重要的研究煤粉燃燒的方法，已得到廣泛的應用。煤粉燃燒數值計算需要確定合適的計算模型，主要包括煤粉燃燒數學模型和煤粉燃燒化學反應模型。

2.1 煤粉燃燒數學模型

流體流動過程在歐拉坐標系下進行求解。連續相方程為質量、動量、能量的連續性方程和時間均值的納維斯托克斯方程。同時，在連續相方程中添加組分輸運方程以求解反應流。

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

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

2.2 煤粉燃燒化學反應模型

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

2.3 網格無關性驗證

為了排除網格數量對計算結果的影響，須開展網格無關性驗證，將網格數量劃分為10萬、30萬以及60萬3種網格。開展冷態數值計算檢驗的方法，得到不同網格數量下燃燒器中心軸線上20個監測點上的速度分布結果，如圖3所示。結果表明，網格數量大于30萬后，燃燒器中心軸線上速度分布趨勢保持不變。因此，同時兼顧計算精度和計算效率，選用數量為30萬的網格開展數值計算。

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

圖3 不同網格數計算結果對比

3 數值模擬工況與煤質參數

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

表1 數值模擬工況主要參數

筆者研究所使用的煤種為神木煙煤，其工業及元素分析等煤質參數見表2。

表2 神木煙煤的煤質分析結果

4 結果與討論

4.1 速度場分析

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

圖4 當一次風速v=12.43 m/s時中心截面軸向速度分布云圖

(1)一次風粉以直流形式噴入預燃室內，在預燃室中心形成具有一定剛性的直流氣流；二次風為旋轉氣流，旋流數為2.16，該值大于0.6，因此二次風屬強旋流，在預燃室內形成明顯的回流區。

(2)當一次風速v=12.43 m/s時，預燃室內形成了1個較大的且為整體的中心回流區，當一次風速v=18.64 m/s和24.86 m/s時，一次風穿透原中心回流區，形成2個對稱的環形回流區，且一次風速越大，將2個環形回流區“擠壓”的越?。?/p>

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

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

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

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

預燃室中心截面的軸向速度矢量圖如圖6所示。對比預燃室內不同一次風速下的流場分布可知，一次風速較小而未穿透中心回流區時，預燃室內回流區的面積顯著大于其他2種情況，即隨著一次風速增大，環形回流區向壁面附近有所移動；由于一次風的風速高、風量大，導致一次風粉的動量顯著增大，對環形回流區造成明顯的“擠壓”現象。

4.2 溫度場分析

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

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

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

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

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

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

表3 預燃室燃燒器的數值計算結果

4.3 組分場分析

預燃室內組分分布可直觀地反映燃燒特性，通過對關鍵組分的分析對整體燃燒情況進行預測。預燃室中心截面揮發分質量分數分布如圖9所示。

圖9 中心截面揮發分濃度分布圖(體積分數/%)

由圖9可知，隨一次風速的增加，揮發分在預燃室內分布區域出現逐漸后移，且高濃度區域面積逐漸增大的趨勢。一次風速增大，一次風與煤粉的預熱速度變緩，煤粉受熱析出的揮發分推遲是導致后移以及后區濃度較高的主要原因。因此，一次風速增大不利于煤粉與周圍高溫煙氣的混合，降低煤粉溫升速度、延遲揮發分的析出可能會造成煤粉整體燃盡效果變差。

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

圖10 中心軸線上氧氣體積分數分布

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

4 結 論

(1)一次風速過高則會穿透預燃室內的中心回流區，從而形成2個對稱的環形回流區，且一次風速越大，其將2個環形回流區“擠壓”得越小。

(2)隨一次風速的增大，中心軸線上溫升速率越小，一次風粉氣流與旋流二次風的混合過程被削弱，限制了其與周圍高溫煙氣的傳熱。一次風速過大，預燃室內中心高溫區被穿透，形成了中心低、外側高的分布，且高溫區有向前移動的趨勢。

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