一種燃燒器爐底布置的煤粉工業(yè)鍋爐的數(shù)值模擬研究

丘性通

福建省鍋爐壓力容器檢驗(yàn)研究院

一種燃燒器爐底布置的煤粉工業(yè)鍋爐的數(shù)值模擬研究

丘性通

福建省鍋爐壓力容器檢驗(yàn)研究院

煤粉工業(yè)鍋爐中燃燒器爐底布置方案具有燃燒效率高、燃燒器便于安裝與檢修、風(fēng)機(jī)耗電量少等優(yōu)點(diǎn)，具有較好的應(yīng)用前景。運(yùn)用商業(yè)CFD軟件STAR-CCM+對(duì)某燃燒器爐底布置的煤粉工業(yè)鍋爐進(jìn)行數(shù)值模擬研究，計(jì)算了爐內(nèi)的流場(chǎng)、溫度場(chǎng)、氧氣濃度分布與氮氧化物濃度分布，分析一次風(fēng)風(fēng)速、內(nèi)二次風(fēng)風(fēng)速與二次風(fēng)旋流強(qiáng)度對(duì)流場(chǎng)的影響，為后續(xù)燃燒器爐底布置的煤粉工業(yè)鍋爐的設(shè)計(jì)及優(yōu)化提供了指導(dǎo)。

煤粉工業(yè)鍋爐；燃燒器爐底布置；數(shù)值模擬

高效的煤粉工業(yè)鍋爐是傳統(tǒng)燃煤工業(yè)鍋爐的升級(jí)換代產(chǎn)品，其燃料燃盡率接近99%，鍋爐熱效率達(dá)到88%左右，具有高效節(jié)能、污染排放低、啟停便捷、負(fù)荷調(diào)節(jié)性好等優(yōu)點(diǎn)，受到不少用戶的好評(píng)[1]。尤其煤粉工業(yè)鍋爐的爐膛所具有的低溫運(yùn)行特點(diǎn)，如果能應(yīng)用燃燒器空氣分級(jí)和煙氣再循環(huán)等低氮燃燒技術(shù)，則可以在不使用復(fù)雜昂貴的非選擇性催化還原（SNCR）煙氣脫硝技術(shù)與選擇性催化還原（SCR）煙氣脫硝技術(shù)的前提下實(shí)現(xiàn)NOx的低排放，滿足日益提高的環(huán)保要求[2]。目前，我國已有大量煤粉工業(yè)鍋爐投入運(yùn)行，由此可見煤粉工業(yè)鍋爐具有良好的發(fā)展前景。

煤粉工業(yè)鍋爐容量在40 t/h以下時(shí)，從降低NOx排放的角度出發(fā)，在鍋爐爐頂或爐底布置單個(gè)燃燒器的方案較為合理，這樣可以保證燃燒器火焰能以爐膛中心為對(duì)稱，燃燒器各個(gè)風(fēng)道出來的風(fēng)與煤粉可以較好貫徹低溫低氧設(shè)計(jì)意圖，有效地組織還原性氣氛和降低火焰峰值溫度；同時(shí)充分利用整個(gè)爐膛高度，對(duì)燃盡有利。以長(zhǎng)期致力于煤粉工業(yè)鍋爐研發(fā)的德國為例，其多款工業(yè)鍋爐產(chǎn)品都采用了燃燒器爐膛頂部布置的方案。燃燒器爐膛頂部布置方案相比，燃燒器爐底布置的方案較少。從理論上分析，燃燒器爐底布置有下列優(yōu)點(diǎn)：煤粉燃燒時(shí)穿過整個(gè)爐膛，燃盡時(shí)間最長(zhǎng)，飛灰含碳量低；燃燒器離地面最近，便于安裝與檢修；一次風(fēng)管道短彎頭少，且位置靠下，可減少一次風(fēng)機(jī)耗電量；爐膛順向高溫拔風(fēng)增強(qiáng)，可減少引風(fēng)機(jī)耗電量等[3]。但爐底布置也可能帶來爐膛結(jié)渣等問題，因此目前投入使用的工業(yè)鍋爐不多，其具體布置方案與運(yùn)行參數(shù)設(shè)置尚有待深入研究。

利用CFD技術(shù)對(duì)不同工況下鍋爐與燃燒器流場(chǎng)進(jìn)行模擬是目前缺乏實(shí)驗(yàn)數(shù)據(jù)情況下研究與優(yōu)化燃燒器的有效方法[4,5]，而且對(duì)于氮氧化物等較難在爐內(nèi)測(cè)量的成分而言[6-8]，數(shù)值模擬可以詳細(xì)地反映出其在爐內(nèi)的產(chǎn)生過程，因此該方法已普遍用于爐內(nèi)燃燒過程的研究[7-9]。近年來，英、美等國家均已對(duì)鍋爐爐膛內(nèi)三維兩相流動(dòng)和煤粉燃燒進(jìn)行模擬，探討將數(shù)值模擬用于爐膛的優(yōu)化的途徑[10、11]，國內(nèi)相關(guān)高校[12]也借助CFD軟件平臺(tái)模擬煤粉鍋爐的燃燒，模擬結(jié)果和實(shí)測(cè)值較為吻合，為鍋爐的優(yōu)化改造提供了參考。

鑒于上述研究現(xiàn)狀，本文擬運(yùn)用商業(yè)CFD軟件STAR-CCM+對(duì)某燃燒器爐底布置的煤粉工業(yè)鍋爐進(jìn)行數(shù)值模擬研究，計(jì)算爐內(nèi)的流場(chǎng)、溫度場(chǎng)、氧氣濃度分布與氮氧化物濃度分布，分析一次風(fēng)風(fēng)速、二次風(fēng)風(fēng)速與二次風(fēng)旋流強(qiáng)度對(duì)流場(chǎng)的影響，為后續(xù)燃燒器爐底布置的煤粉工業(yè)鍋爐的設(shè)計(jì)及優(yōu)化提供指導(dǎo)。

圖1 鍋爐結(jié)構(gòu)圖

1 模擬對(duì)象介紹

1.1 鍋爐介紹

數(shù)值模擬的對(duì)象為某DHS35-1.6-A III型煤粉工業(yè)鍋爐，該鍋爐是一款高效鍋爐，擁有高效節(jié)能，工作環(huán)境友好，鍋爐操作簡(jiǎn)單，潔凈排放，節(jié)約用地，性價(jià)比高等優(yōu)點(diǎn)，鍋爐結(jié)構(gòu)如圖1所示。其主要參數(shù)為：額定蒸發(fā)量35 t/h，額定蒸汽壓力1.6 Mpa，鍋爐出口蒸汽溫度204.2℃，設(shè)計(jì)熱效率88.31%。

該鍋爐為單鍋筒室燃爐，有一個(gè)旋流燃燒器，燃燒器為爐底布置。煤粉與空氣在爐膛內(nèi)燃燒，所生成的高溫?zé)煔庠跔t膛內(nèi)輻射與對(duì)流換熱后依次沖刷對(duì)流管排，省煤器與空氣預(yù)熱器后排至煙囪。其具體結(jié)構(gòu)與三維建模圖2所示。

圖2 鍋爐三維建模圖

1.2 燃燒器介紹

該燃燒器為35 t/h蒸發(fā)量的旋流煤粉燃燒器。燃燒器軸心線上為一次風(fēng)煤粉氣流通道，以直流射流的形式進(jìn)入爐膛；二次風(fēng)分為內(nèi)二次風(fēng)與外二次風(fēng)，皆為旋轉(zhuǎn)射流，通過環(huán)形通道口進(jìn)入爐膛。燃燒器結(jié)構(gòu)如圖3所示。

燃燒器設(shè)計(jì)燃料為III類煙煤煤粉，燃料收到基的碳、氫、氧、氮、硫、水分、灰分分別為55.82%、4.95%、8.77%、1.04%、0.51%、12.2%、16.71%。

燃料揮發(fā)分為46.04%，低位發(fā)熱量5 345 kcal/kg，燃料額定消耗量4 126.99 kg/h，理論空氣量5.999 m3/kg。

圖3 燃燒器三維建模圖

2 模擬方程與邊界條件

2.1 模擬方程

參考相關(guān)文獻(xiàn)[13-16]，在模擬過程中基于質(zhì)量守恒和能量守恒定律，氣相燃燒模擬的控制方程采用SIMPLE算法求解，燃燒反應(yīng)采用EBU模型控制，粘性模型為標(biāo)準(zhǔn)k-ε湍流模型，輻射傳熱采用DOM輻射模型，揮發(fā)份析出采用雙方程模型控制速率，焦炭燃燒采用動(dòng)力學(xué)/擴(kuò)散控制反應(yīng)速率模型。

2.2 邊界條件

入口邊界設(shè)置：模型中旋流燃燒器的一次風(fēng)為直流，因此將一次風(fēng)進(jìn)口面設(shè)置為速度進(jìn)口，直接設(shè)置軸向進(jìn)口速度。旋流燃燒器的二次風(fēng)道中裝有切向可動(dòng)葉片，為了獲得不同的二次風(fēng)旋流強(qiáng)度需要對(duì)二次風(fēng)旋流葉片角度進(jìn)行調(diào)整。在以往的旋流燃燒器模擬中，通常將二次風(fēng)進(jìn)口面設(shè)置為速度進(jìn)口,在進(jìn)口處同時(shí)設(shè)置軸向和切向兩個(gè)方向的速度，以實(shí)現(xiàn)旋流效果。這樣做雖然降低了模型構(gòu)建工作量與計(jì)算量，但是如果沒有實(shí)地測(cè)量旋流燃燒器的二次風(fēng)出口風(fēng)速，則模擬結(jié)果往往與實(shí)際情況相差較遠(yuǎn)。因此本文在研究中采用直接建立帶有不同角度的葉片的燃燒器模型的方法來進(jìn)行模擬，使得模擬更加接近實(shí)際。

出口邊界設(shè)置：出口邊界條件設(shè)置為壓力出口?？紤]到鍋爐爐膛一般處于負(fù)壓運(yùn)行，根據(jù)鍋爐煙風(fēng)阻力估算，出口表壓設(shè)置為-10 Pa。

壁面邊界設(shè)置：壁面設(shè)置為定壁溫條件，且無滑移，無質(zhì)量滲透。

3 計(jì)算結(jié)果與分析

通過分析模型，該結(jié)構(gòu)主要難點(diǎn)在燃燒器設(shè)計(jì)建模上，需要局部捕捉燃燒器的細(xì)致結(jié)構(gòu)。模擬中主要采用STAR-CCM+軟件提供的多面體網(wǎng)格來進(jìn)行劃分，并配合多層邊界層網(wǎng)格，網(wǎng)格質(zhì)量滿足計(jì)算要求。整體網(wǎng)格數(shù)量104萬，如圖4所示。

圖4 鍋爐與燃燒器的網(wǎng)格圖

針對(duì)煤粉燃燒問題的模擬，求解順序?yàn)椋孩倮鋺B(tài)流動(dòng)的模擬，此時(shí)，在流場(chǎng)中并不撒播拉格朗日相；②當(dāng)冷態(tài)流動(dòng)獲得基本正確的解后，此時(shí)流場(chǎng)中放入拉格朗日相，但可通過調(diào)低煤粉溫度等措施，使流場(chǎng)內(nèi)不發(fā)生燃燒現(xiàn)象，成為冷態(tài)的氣固兩相流動(dòng)模擬；③當(dāng)計(jì)算穩(wěn)定后，可調(diào)高煤粉溫度到設(shè)定值，發(fā)生預(yù)期的反應(yīng)，從而計(jì)算得到最終的流場(chǎng)、溫度場(chǎng)等信息。

3.1 冷態(tài)模擬

模擬首先針對(duì)冷態(tài)流動(dòng)進(jìn)行計(jì)算。根據(jù)設(shè)計(jì)工況，當(dāng)入口一次風(fēng)速為25.5 m/s，內(nèi)二次風(fēng)速為28 m/s，外二次風(fēng)速為25 m/s，氣流溫度為25℃時(shí)，其冷態(tài)流場(chǎng)如圖5所示，

在本算例中，旋流燃燒器的內(nèi)二次風(fēng)道與外二次風(fēng)道里的氣流在出燃燒器之前是作螺旋運(yùn)動(dòng)的。當(dāng)它剛離開燃燒器時(shí)，沿著螺旋線的切線方向運(yùn)動(dòng)，產(chǎn)生輻射狀的旋轉(zhuǎn)射流，不僅具有軸向速度，而且有著較大的切向速度。此時(shí)氣流的擾動(dòng)非常強(qiáng)烈，見圖6所示。在后續(xù)運(yùn)動(dòng)中，射流不斷卷吸周圍氣體，其切向速度的旋轉(zhuǎn)半徑不斷增大，切向速度很快衰減，因而降低后擾動(dòng)的劇烈程度。同時(shí)軸向速度也由于卷吸周圍氣體而很快衰減，在離開出口一段距離后射流的軸向速度降為負(fù)值，即產(chǎn)生一個(gè)回流區(qū)。此回流區(qū)將有助于煤粉的燃燒，其速度矢量如圖7所示。

圖5 冷態(tài)爐膛流場(chǎng)速度分布

圖6 冷態(tài)時(shí)燃燒器出口處速度分布

圖7 冷態(tài)時(shí)回流區(qū)速度分布

圖8 一次風(fēng)速變化時(shí)的軸向風(fēng)速

通過改變一次風(fēng)速、內(nèi)二次風(fēng)速、內(nèi)二次風(fēng)旋流強(qiáng)度，來尋求流場(chǎng)，特別是回流區(qū)的分布規(guī)律，結(jié)果見圖8、圖9和圖10所示，分別給出了各個(gè)工況下距旋流燃燒器出口0.5 m橫截面上的軸向速度分布情況。圖中可以看出，軸向速度均呈開口向上的拋物線形狀分布。分析原因可知，盡管直流的一次風(fēng)向上噴出，但在燃燒器出口鈍體和旋流內(nèi)二次風(fēng)的共同作用下，使得在燃燒器出口處內(nèi)流場(chǎng)中心區(qū)域形成負(fù)壓，導(dǎo)致了回流區(qū)的產(chǎn)生，因而在該位置處軸向速度為負(fù)值。

通過對(duì)比觀察回流區(qū)的軸向速度分布，當(dāng)單獨(dú)改變一次風(fēng)速度時(shí)，隨著一次風(fēng)速增大，中心回流區(qū)軸向速度分別為-6.2 m/s、-5.6 m/s、-4.3m/s，速度絕對(duì)值減小明顯，導(dǎo)致回流區(qū)減弱。

當(dāng)單獨(dú)改變內(nèi)二次風(fēng)速度時(shí)，隨著內(nèi)二次風(fēng)速增大，中心回流區(qū)軸向速度分別為-5.0 m/s、-5.5 m/s、-5.9 m/s，速度分布情況較為相似，速度絕對(duì)值增加，增強(qiáng)回流效果。

圖9 內(nèi)二次風(fēng)速變化時(shí)的軸向風(fēng)速

圖10 內(nèi)二次風(fēng)速變化時(shí)的軸向風(fēng)速

當(dāng)單獨(dú)改變內(nèi)二次風(fēng)旋流強(qiáng)度時(shí)，隨著旋流強(qiáng)度增大，中心回流區(qū)軸向速度分別為-4.9 m/s、-5.6 m/s、-6.2 m/s，回流速度絕對(duì)值增加明顯，而且，回流區(qū)速度較大區(qū)域變寬，導(dǎo)致平均回流速度增大。這種方案對(duì)增強(qiáng)回流區(qū)效果最為明顯，由此說明增強(qiáng)內(nèi)二次風(fēng)旋流強(qiáng)度可使得中心回流區(qū)流速增大，增加回流量，加強(qiáng)了此處流場(chǎng)的卷吸高溫?zé)煔獾哪芰?，有助于煤粉的燃燒?/p>

3.2 熱態(tài)模擬

冷態(tài)模擬完成后，根據(jù)冷態(tài)分析的結(jié)果，并綜合考慮工況條件，選擇最優(yōu)的入口條件進(jìn)行熱態(tài)模擬。當(dāng)入口一次風(fēng)速為25.5 m/s，溫度為100℃；內(nèi)二次風(fēng)速為28 m/s，溫度為250℃；外二次風(fēng)速為25 m/s，溫度為250℃；煤粉質(zhì)量流量為0.415 kg/s時(shí)，其熱態(tài)的流場(chǎng)與溫度場(chǎng)分別如圖11和12所示，氧氣與氮氧化物分布如圖13與圖14所示。

根據(jù)圖11與圖12可知，燃燒器出口存在一個(gè)1 800 K左右的高溫區(qū)，該高溫區(qū)域在燃燒器預(yù)燃室內(nèi)及預(yù)燃室出口后段呈叉狀形態(tài)，其原因是旋流燃燒器采用了強(qiáng)烈的旋流來強(qiáng)化燃燒，在旋流的作用下，燃燒器預(yù)燃室中心會(huì)產(chǎn)生一個(gè)回流區(qū)將煤粉燃燒形成的高溫?zé)煔庀蚧鼐砦己筮M(jìn)入的煤粉顆粒，因此高溫區(qū)中間呈向鈍體方向凹陷的形狀。在回流區(qū)的邊界部位熱、質(zhì)量交換最為激烈，且該部分氧量充足,所以煤粉燃燒放熱充分，產(chǎn)生接近2 000 K的高溫。從圖11與圖12的模擬可見，火焰作為煤粉穩(wěn)定燃燒的熱源能深入整個(gè)爐膛，為燃燒提供了一個(gè)較好的爐內(nèi)溫度環(huán)境，并最終達(dá)到較好的燃盡水平。

圖13中的氧氣濃度分布情況可以反映出空間中各部分的煤粉燃燒程度。由圖13可見從燃燒器預(yù)燃室起延伸到爐膛內(nèi)很長(zhǎng)一段距離內(nèi)有一個(gè)低氧濃度區(qū)域，尤其在燃燒器出口截面處，氧氣濃度分布呈中間接近零兩邊略高的形態(tài)。這是由于圖12中燃燒器出口回流區(qū)的存在，使得爐內(nèi)煤粉燃燒產(chǎn)生的高溫低氧煙氣被引回到預(yù)燃室出口，從而令出口截面中心處氧氣含量降低，形成缺氧燃燒。這種缺氧燃燒可以有效抑制煤粉燃燒過程中NOx的生成，這一點(diǎn)在圖14的氮氧化物分布中可以被觀察到。

由上面的分析可見，合理選擇工況條件可以使得燃燒器爐底布置的煤粉工業(yè)鍋爐在高效燃燒的前提下實(shí)現(xiàn)NOx低排放，實(shí)現(xiàn)節(jié)能低污染的良好效果。

圖11 熱態(tài)爐膛溫度場(chǎng)分布

圖12 熱態(tài)爐膛流場(chǎng)分布

4 結(jié)論

（1）工業(yè)鍋爐中燃燒器爐底布置的方案具有煤粉燃盡時(shí)間最長(zhǎng)、飛灰含碳量低、燃燒器便于安裝與檢修、風(fēng)機(jī)耗電量少等優(yōu)點(diǎn)，但是其也可能存在爐膛結(jié)渣等問題，具體布置方案與運(yùn)行參數(shù)設(shè)置尚有待深入研究。

（2）利用CFD工具可以模擬燃燒器爐底布置的煤粉工業(yè)鍋爐的冷態(tài)與熱態(tài)工況，得到的模擬結(jié)果較為可靠，可將其用于燃燒器爐底布置的煤粉工業(yè)鍋爐的設(shè)計(jì)與優(yōu)化中。

（3）根據(jù)CFD模擬可知，旋流燃燒器出口產(chǎn)生回流區(qū)。回流區(qū)加強(qiáng)了流場(chǎng)卷吸高溫?zé)煔獾哪芰?，有助于煤粉的燃燒。一次風(fēng)速度的減小，內(nèi)二次風(fēng)速度的增加以及內(nèi)二次風(fēng)旋流強(qiáng)度的增加都可以增強(qiáng)回流效果，其中內(nèi)二次風(fēng)旋流強(qiáng)度對(duì)回流效果的影響最為顯著。在煤粉燃燒時(shí)，回流區(qū)使得燃燒器出口截面中心處氧氣含量降低，有效抑制煤粉燃燒過程中NOx的生成，實(shí)現(xiàn)了NOx的低排放。

圖13 熱態(tài)爐膛內(nèi)氧氣分布

圖14 熱態(tài)爐膛內(nèi)氮氧化物濃度分布

[1] 俞建洪, 王善武. 煤粉工業(yè)鍋爐及其發(fā)展前景[J]. 工業(yè)鍋爐, 2014(1). 1-5.

[2] 秦裕琨, 李爭(zhēng)起, 吳少華. 旋流煤粉燃燒技術(shù)的發(fā)展[J]. 熱能動(dòng)力工程, 1997(12). 241-245.

[3] 張譽(yù)鏵. 一種燃燒器爐底布置的新型高效煤粉工業(yè)鍋爐[J]. 工業(yè)鍋爐, 2013年(1). 15-18.

[4] Purimetla A, Cui J. CFD studies on burner seconda-ry airflow[J]. Applied Mathematical Modelling, 2009, 33(2). 1126-1140.

[5] 朱霖,池作和,孫公鋼,等.一種適用于煤粉工業(yè)鍋爐燃燒器的冷態(tài)流場(chǎng)特性研究[J].中國計(jì)量學(xué)院學(xué)報(bào). 2011, 22(1). 39-43.

[6] Ji C C, Cohen R D．An investigation of the combustion of pulverizedcoal-air mixture in different combustor geometries[J]．Combustion and Flame. 1992, 90(3-4). 307-343．

[7] Anagnostopoulos J S,Sargianos N P．The prediction of pulverized Greek lignite combustion in axisymmetric furnaces[J]．Combustion and Flame. 1993,92(3). 209-221.

[8] Faltsi-Saravelou O，Wild P. Detailed modeling of a swirling coal flame[J]．Combust Sci. and Tech. 1997, 123(1-6). 1-22．

[9] 王磊. 煤粉燃燒器出口旋轉(zhuǎn)氣固兩相流場(chǎng)的數(shù)值模擬[J].機(jī)械工程學(xué)報(bào), 2001,37(9). 9-11.

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

[11] Srdjan B, Miroslav S, Simeon O. Three-dimensional modeling of utility boiler pulverized coal tangentially fired furnace[J]. International Journal of Heat and Mass Transfer. 2006, 49 (3). 3371-3378.

[12] 劉霞, 葛新鋒. FLUENT 軟件及其在我國的應(yīng)用[J]. 能源研究與利用. 2003 (2). 36-38.

[13] Jones J M, Pourkashanian M, Williams A. Modeling the combustion of pulverized biomass in an industrial combustion test furnace [J]. Fuel. 2007, 10 (5). 16-19.

[14] Hu Z J, Lu J D, Huang L Numerical simulation study on gas-solid twophase flow in proclaimer [J]. Communicationsin Nonlinear Science and Numerical Simulation. 2006, 11 (4). 440-451.

[15] Zhou L X, Li L, Li R X. Simulation of 3-D gas-particle f lows and coal combustion in a tangentially fired furnace using a two-fluid-trajectory model[J]. Powder Technology. 2002, 125 (3). 226-233.

[16] Xie J，He F，Yuan R. Relationship between coal powder and its combustibility [J]. Journal of Wuhan University of Technology，Materials Science Edition. 2002. 17 (1), 72-75.

Numerical Simulation Study on One Kind of Industrial Pulverized Coal Boiler with Burner Bottom Arranged

Qiu Xingtong
Fujian Province Boiler and Pressured Vessel Inspection and Research Institute

Industrial pulverized coal boiler with burner bottom arranged has several advantages, such as high burning efficiency, easy installation and maintenance for burner, low power consumption for fan and better application prospect. The author uses commercial CFD software STAR-CCM+ to carry out numerical simulation study on industrial pulverized coal boiler with burner bottom arranged. Flow field, temperature field, oxygen concentration distribution and nitrogen oxide concentration distribution inside boiler are calculated. Impact on flow field of primary wind speed, inner secondary wind speed and secondary wind cyclone intensity are analyzed, which provides reference to design and optimization of industrial pulverized coal boiler with burner bottom arranged in the future.

Industrial Pulverized Coal Boiler, Burner Bottom Arranged, Numerical Simulation

10.13770/j.cnki.issn2095-705x.2016.12.006

丘性通：（1969-）,男，福建省鍋爐壓力容器檢驗(yàn)研究院節(jié)能中心副主任，高級(jí)工程師。