賈國(guó)海，李立君，高自成，李際平，陳喜龍

（中南林業(yè)科技大學(xué)機(jī)電工程學(xué)院，湖南長(zhǎng)沙410004）

　　摘要：以新型生物質(zhì)顆粒回轉(zhuǎn)燃燒器為研究對(duì)象，采用Pro/E軟件對(duì)生物質(zhì)顆粒回轉(zhuǎn)燃燒器內(nèi)流體進(jìn)行幾何建模，并用Gambit軟件對(duì)模型應(yīng)用非結(jié)構(gòu)化網(wǎng)格生成技術(shù)劃分網(wǎng)格，并進(jìn)行有限元前處理。在采用k-ε湍流模型的基礎(chǔ)上，應(yīng)用計(jì)算流體力學(xué)軟件Fluent模擬穩(wěn)定工況下回轉(zhuǎn)燃燒器內(nèi)流體的流動(dòng)特性，分析進(jìn)氣速度對(duì)回轉(zhuǎn)燃燒器內(nèi)的壓力分布、速度分布以及湍動(dòng)能分布的影響。結(jié)果表明：燃燒器內(nèi)部的氣體壓力、速度和湍動(dòng)能隨風(fēng)機(jī)引風(fēng)速度的增大逐漸增大，風(fēng)機(jī)出口到回轉(zhuǎn)燃燒室之間的氣體壓力、速度和湍動(dòng)能較大，在二次風(fēng)口處達(dá)到最大值，回轉(zhuǎn)燃燒室內(nèi)氣體的壓力、速度和湍動(dòng)能分布較為均勻。

　　生物質(zhì)顆粒燃料是一種典型的生物質(zhì)固體成型燃料，具有高效、潔凈、容易點(diǎn)火、CO₂近零排放等優(yōu)點(diǎn)，可替代煤炭等化石燃料應(yīng)用于炊事、供暖等民用領(lǐng)域和鍋爐燃燒、發(fā)電等工業(yè)領(lǐng)域，近幾年來(lái)在歐盟、北美、中國(guó)得到了迅速發(fā)展。文獻(xiàn)[1]分析了我國(guó)生物質(zhì)能源技術(shù)的現(xiàn)狀、發(fā)展趨勢(shì)和存在問(wèn)題，提出了發(fā)展對(duì)策和建議。文獻(xiàn)[2]對(duì)南方主要薪炭樹(shù)種的生物量、熱值進(jìn)行聚類分析，并結(jié)合生物生態(tài)學(xué)特性及種植推廣情況進(jìn)行了綜合評(píng)價(jià)。文獻(xiàn)[3]對(duì)13種生物質(zhì)原料的物理性能及燃燒特性進(jìn)行了研究，并分析了燃燒后的氣體排放成分、灰分組成等特性。文獻(xiàn)[4]對(duì)不同生物質(zhì)燃料的元素組成和燃燒特性進(jìn)行了對(duì)比分析。文獻(xiàn)[5]對(duì)木屑和城市固體廢物進(jìn)行燃燒特性的模擬和試驗(yàn)研究，分析了燃燒爐在不同風(fēng)量影響下的燃燒特性。文獻(xiàn)[6]為探討木質(zhì)顆粒燃料應(yīng)用的效益，對(duì)比分析了常規(guī)能源和木質(zhì)顆粒燃料的單位能價(jià)及典型鍋爐實(shí)際運(yùn)行成本。文獻(xiàn)[7]采用Fluent軟件對(duì)一種煤摻燒生物質(zhì)的旋流燃燒器進(jìn)行了冷態(tài)數(shù)值模擬，計(jì)算結(jié)果表明：安裝合適角度的旋流葉片可使得旋流燃燒器中的氣流擾動(dòng)更加強(qiáng)烈，從而可提高生物質(zhì)的燃燒效率。文獻(xiàn)[8]采用Fluent對(duì)燃燒器空氣流場(chǎng)情況進(jìn)行了數(shù)值模擬，分析了不同二次風(fēng)門開(kāi)度下旋流燃燒器空氣流場(chǎng)特性。文獻(xiàn)[9]采用渦流耗散（EDC）模型對(duì)一種預(yù)熱室旋流燃燒器的燃燒特性進(jìn)行了數(shù)值模擬，并對(duì)6種不同含氧量下的燃燒特性進(jìn)行了分析。

　　文獻(xiàn)[10]對(duì)一種具有空氣分級(jí)裝置的小型生物質(zhì)燃燒系統(tǒng)進(jìn)行了CFD仿真分析，計(jì)算結(jié)果表明該裝置可以在較低的空氣比下實(shí)現(xiàn)生物質(zhì)的完全燃燒。文獻(xiàn)[11]對(duì)一種生物質(zhì)顆粒氣化爐的燃燒效率和污染物排放特性進(jìn)行了分析，得出在控制好空氣混合比的情況下可大大減少污染物排放并提高燃燒效率。文獻(xiàn)[12]對(duì)固定床鍋爐中的生物質(zhì)燃燒特性進(jìn)行了數(shù)值模擬和實(shí)驗(yàn)測(cè)試對(duì)比分析，得到了固定床中固相和鍋爐爐中氣相的主要變化曲線。文獻(xiàn)[13]采用計(jì)算流體動(dòng)力學(xué)（CFD）軟件對(duì)生物質(zhì)氣化爐進(jìn)行建模，并對(duì)不同生物質(zhì)顆粒的氣化特性進(jìn)行了仿真與試驗(yàn)研究。

　　以上分析可知，生物質(zhì)顆粒燃料廣泛應(yīng)用于生物質(zhì)鍋爐、熱風(fēng)爐、采暖爐中，但大部分燃燒器不具有抗結(jié)渣性的特點(diǎn)。針對(duì)這種情況，本研究團(tuán)隊(duì)設(shè)計(jì)一種可抗結(jié)渣性的生物質(zhì)回轉(zhuǎn)顆粒燃燒器。燃燒器內(nèi)流速分布特性對(duì)生物質(zhì)顆粒燃料的燃燒效率以及燃燒器的使用壽命均有重要的影響，筆者利用流體計(jì)算軟件對(duì)該新型生物質(zhì)顆粒回轉(zhuǎn)燃燒器內(nèi)流場(chǎng)特性進(jìn)行研究，能夠減少試驗(yàn)工作量，縮短設(shè)計(jì)周期。

　　1燃燒器結(jié)構(gòu)

　　為了提高生物質(zhì)顆粒回轉(zhuǎn)燃燒器熱效率，使燃燒更充分，減少煙塵量的生成，生物質(zhì)顆粒回轉(zhuǎn)燃燒器采用回轉(zhuǎn)燃燒室。此外，還采用螺旋進(jìn)料裝置控制進(jìn)料量，保證顆粒燃料充分燃燒，從而控制嚴(yán)重結(jié)渣。

　　生物質(zhì)顆粒回轉(zhuǎn)燃燒器結(jié)構(gòu)示意圖如圖1所示，主要由螺旋進(jìn)料機(jī)構(gòu)、回轉(zhuǎn)燃燒室、小燃燒室和聯(lián)接法蘭等部件組成。燃燒器在工作時(shí)，生物質(zhì)顆粒從進(jìn)料口進(jìn)入到進(jìn)料螺旋管中，在進(jìn)料螺旋機(jī)構(gòu)的作用下，快速、平穩(wěn)、均勻地推進(jìn)到回轉(zhuǎn)燃燒室。在回轉(zhuǎn)燃燒室內(nèi)，燃料被干燥，然后進(jìn)入燃燒區(qū)域，生物質(zhì)顆粒完成一次燃燒所需要的空氣通過(guò)鼓風(fēng)機(jī)將氣流從回轉(zhuǎn)燃燒室的下部和環(huán)繞在燃燒床內(nèi)的燃燒室上的小孔進(jìn)入，并提供適宜的燃料/空氣比，營(yíng)造理想的燃燒氛圍。生物質(zhì)顆粒燃料在環(huán)繞整個(gè)燃燒室外的二次風(fēng)口處吹入二次風(fēng)，這些二次空氣使得燃燒更加充分，保證所有的生物質(zhì)燃料在回轉(zhuǎn)燃燒室內(nèi)完全燃燒。

　　2理論模型和計(jì)算方法

　　鼓風(fēng)機(jī)首先將空氣噴射到燃燒區(qū)域，然后生物質(zhì)顆粒燃料在螺旋進(jìn)料機(jī)構(gòu)的輸送下進(jìn)入燃燒室，與空氣混合進(jìn)行燃燒。生物質(zhì)顆粒回轉(zhuǎn)燃燒器可以在較短的時(shí)間內(nèi)達(dá)到穩(wěn)態(tài)燃燒，故對(duì)燃燒器在穩(wěn)態(tài)狀態(tài)下的流動(dòng)特性進(jìn)行研究。

　　2.1三維網(wǎng)格模型的建立

　　在進(jìn)行數(shù)值仿真計(jì)算前，采用Pro/E軟件建立生物質(zhì)顆粒回轉(zhuǎn)燃燒器的三維模型，如圖2所示。

　　為了使仿真模型的計(jì)算收斂速度更快，在不影響計(jì)算精度的情況下對(duì)燃燒器表面倒角和細(xì)小結(jié)構(gòu)進(jìn)行適當(dāng)簡(jiǎn)化。采用Gambit軟件對(duì)燃燒器內(nèi)流體模型進(jìn)行幾何處理并使用多面體網(wǎng)格功能對(duì)燃燒器內(nèi)流體模型進(jìn)行體網(wǎng)格的劃分，燃燒器流體的三維體網(wǎng)格模型如圖3所示，其節(jié)點(diǎn)數(shù)為1104429個(gè)，網(wǎng)格數(shù)為338087個(gè)。

　　2.2三維流動(dòng)的控制方程

　　在燃燒器的流動(dòng)仿真計(jì)算中，將其視為穩(wěn)態(tài)的湍流運(yùn)動(dòng)，其流動(dòng)和傳熱過(guò)程都遵從質(zhì)量守恒、能量守恒定律和動(dòng)量守恒[10]。

　　2.3湍流模型

　　反映湍流脈動(dòng)量對(duì)流場(chǎng)影響的湍流動(dòng)能方程和湍流應(yīng)力方程可通過(guò)k-ε方程[12]得到：

　　2.4邊界條件及求解設(shè)置

　　因?yàn)槿紵鞯囊L(fēng)口為矩形，所以其等效水力直徑[8]為

　　根據(jù)式（7）-（9）得到不同引風(fēng)速度下的引風(fēng)口和燃料入口的等效水力直徑、雷諾數(shù)和湍流強(qiáng)度，如表1所示。

　　以Fluent軟件為基礎(chǔ)，采用k-ε湍流模型模擬穩(wěn)定工況下該型生物質(zhì)回轉(zhuǎn)燃燒器內(nèi)流體的流動(dòng)特性，入口邊界采用速度入口邊界條件，燃燒器出口處邊界條件選擇出口壓力邊界，依照實(shí)際工作要求，設(shè)置為大氣壓力。

　　為了模擬絕熱燃燒過(guò)程，壁面設(shè)為絕熱壁面，沒(méi)有熱量和質(zhì)量流量通過(guò)，壁面為無(wú)滑移邊界條件，近壁區(qū)采用標(biāo)準(zhǔn)壁面函數(shù)。離散格式中壓力插補(bǔ)格式采用PRESTO格式，對(duì)流項(xiàng)采用QUICK格式。設(shè)置收斂為平均殘差小于10-5，并且當(dāng)進(jìn)風(fēng)口、出風(fēng)口幾個(gè)預(yù)設(shè)監(jiān)測(cè)點(diǎn)速度、壓力都趨于穩(wěn)定時(shí)，可認(rèn)為收斂并停止求解。

　　3仿真結(jié)果分析

　　3.1燃燒器內(nèi)壓力分布

　　生物質(zhì)回轉(zhuǎn)燃燒器在不同引風(fēng)速度下的氣相壓力分布如圖4所示。仿真結(jié)果表明：在生物質(zhì)旋流燃燒器整體壓力最大區(qū)域主要分布在風(fēng)機(jī)引風(fēng)處到回轉(zhuǎn)燃燒室二次風(fēng)口之間，其他區(qū)域壓力值分布較為均勻；氣體進(jìn)入生物質(zhì)回轉(zhuǎn)燃燒室上的二次風(fēng)口后，壓力值逐漸減小，在回轉(zhuǎn)燃燒室內(nèi)部壓力基本保持一致。

　　從圖4可以看出：燃燒器內(nèi)部的整體壓力隨風(fēng)機(jī)引風(fēng)速度的增大逐漸增大。當(dāng)風(fēng)機(jī)引風(fēng)風(fēng)速為25m·s^-1時(shí)，二次風(fēng)口的最大壓力為541.0Pa；當(dāng)風(fēng)機(jī)引風(fēng)風(fēng)速為35m·s^-1時(shí)，二次風(fēng)口的最大壓力為1059.5Pa；當(dāng)風(fēng)機(jī)引風(fēng)風(fēng)速為45m·s^-1時(shí)，二次風(fēng)口的壓力增大為1750.7Pa。

　　3.2燃燒器內(nèi)流速分布

　　燃燒器在不同風(fēng)機(jī)引風(fēng)風(fēng)速下空氣流場(chǎng)速度的數(shù)值模擬結(jié)果如圖5所示，當(dāng)風(fēng)機(jī)引風(fēng)風(fēng)速為25m·s^-1時(shí)，二次風(fēng)口的最大速度為25.20m·s^-1；當(dāng)風(fēng)機(jī)引風(fēng)風(fēng)速為35m·s^-1時(shí)，二次風(fēng)口的最大速度為35.32m·s^-1；當(dāng)風(fēng)機(jī)引風(fēng)風(fēng)速為45m·s^-1時(shí)，二次風(fēng)口的最大速度為45.40m·s^-1。從總體上來(lái)說(shuō)，燃燒器內(nèi)部的氣體流速隨風(fēng)機(jī)引風(fēng)風(fēng)速的增大逐漸增大，但燃燒室內(nèi)的空氣速度流場(chǎng)分布情況略有差異。以風(fēng)機(jī)引風(fēng)風(fēng)速為45m·s^-1為例，對(duì)燃燒器內(nèi)部的風(fēng)速進(jìn)行分析。靠近風(fēng)機(jī)引風(fēng)附近的二次風(fēng)口區(qū)域的風(fēng)速達(dá)到最大值44.30m·s^-1，空氣流速?gòu)娘L(fēng)機(jī)入口處沿著回轉(zhuǎn)燃燒室不斷下降，至燃燒器的下方，達(dá)到最小值2.33m·s^-1；經(jīng)過(guò)回轉(zhuǎn)燃燒室的二次風(fēng)口時(shí)，風(fēng)速又快速上升，平均速度為21.00m·s^-1，回轉(zhuǎn)燃燒室內(nèi)部的流速比較穩(wěn)定，平均流速為9.93m·s^-1，攪動(dòng)筋板對(duì)氣流的流動(dòng)起到均勻的導(dǎo)向作用。

　　3.3燃燒器內(nèi)湍動(dòng)能分布

　　燃燒器的湍動(dòng)能分布云圖如圖6所示。湍動(dòng)能是與流體速度相關(guān)的參數(shù)，它是脈動(dòng)動(dòng)能而非平均流動(dòng)能。

　　從圖6可以看出：湍動(dòng)能變化較大的區(qū)域主要集中在風(fēng)機(jī)引風(fēng)到對(duì)應(yīng)的二次風(fēng)口之間，經(jīng)過(guò)二次風(fēng)口處的湍動(dòng)能達(dá)到最大值，在其他的區(qū)域，湍動(dòng)能值較小，基本保持穩(wěn)定不變；隨著風(fēng)機(jī)引風(fēng)速度的增大，湍動(dòng)能也相應(yīng)地增大：當(dāng)風(fēng)機(jī)引風(fēng)風(fēng)速為25m·s^-1時(shí)，二次風(fēng)口的最大湍動(dòng)能為77.58m²·s^-2；當(dāng)風(fēng)機(jī)引風(fēng)風(fēng)速為35m·s^-1時(shí)，二次風(fēng)口的最大湍動(dòng)能為152.04m²·s^-2；當(dāng)風(fēng)機(jī)引風(fēng)風(fēng)速為45m·s^-1時(shí)，二次風(fēng)口的湍動(dòng)能增大為251.32m²·s^-2。湍流具有擴(kuò)散能力，能夠耗散氣體的保勢(shì)動(dòng)能，形成時(shí)而膨脹時(shí)而壓縮的氣體介質(zhì)，也會(huì)引起氣體的回流和渦流，阻礙氣體的流動(dòng)，增加各部分的靜壓力，產(chǎn)生壓力損失。比較圖4，6可以看出，湍動(dòng)能大的地方壓力變化也較大。

　　4結(jié)論

　　1）燃燒器內(nèi)壓力主要分布在2個(gè)區(qū)域：回轉(zhuǎn)燃燒室的內(nèi)部和外部。燃燒室內(nèi)部的壓力分布較為均勻，燃燒室外部的壓力大致呈現(xiàn)遞減趨勢(shì)，在風(fēng)機(jī)引風(fēng)到二次風(fēng)口處的壓力差最大，主要因?yàn)槎物L(fēng)口阻力造成回轉(zhuǎn)燃燒室的阻力損失。

　　2）燃燒器內(nèi)部的氣體流速隨風(fēng)機(jī)引風(fēng)風(fēng)速的增大逐漸增大，風(fēng)機(jī)引風(fēng)和二次風(fēng)口附近的流速比較大，風(fēng)機(jī)引風(fēng)和二次風(fēng)口附近的渦流現(xiàn)象比較明顯，回轉(zhuǎn)燃燒室內(nèi)氣體的流速分布基本均勻。

　　3）燃燒器湍動(dòng)能變化較大的區(qū)域主要集中在風(fēng)機(jī)引風(fēng)到對(duì)應(yīng)的二次風(fēng)口之間，在其他的區(qū)域，湍動(dòng)能較小，基本保持穩(wěn)定不變。隨著風(fēng)機(jī)引風(fēng)速度的增大，湍動(dòng)能也相應(yīng)增大，燃燒器的壁面和回轉(zhuǎn)燃燒室上的小孔對(duì)流體的限制和擾動(dòng)是其內(nèi)部產(chǎn)生湍流的主要原因。

　　參考文獻(xiàn)

　　[1]吳創(chuàng)之，周肇秋，陰秀麗，等.我國(guó)生物質(zhì)能源發(fā)展現(xiàn)狀與思考[J].農(nóng)業(yè)機(jī)械學(xué)報(bào)，2009，40(1)：91-98.

　　[2]李際平，陳喜龍，朱寧華，等.南方薪炭林樹(shù)種選擇及其木質(zhì)顆粒燃燒特性研究[J].中南林業(yè)科技大學(xué)學(xué)報(bào)，2013，33(12)：126-129.

　　[3]GARCÍA R，PIZARRO C，LAVÍN A G，et al.Biomass sources for thermal conversion techno-economical over-view[J].Fuel，2017，195：182-189.

　　[4]DEMIRABAS A.Combustion characteristics of different biomass fuel[J].Progress in Energy and Combustion Science，2004，30：219-230.

　　[5]YANG Y B，SHARIFI V N，SWITHENBANK J.Effect of air flow rate and fuel moisture on the burning beha-viors of biomass and simulated municipal solid wastes in packed beds[J].Fuel，2004，83：1553-1562.

　　[6]陳喜龍，李際平，王義強(qiáng)，等.木質(zhì)顆粒燃料鍋爐替代燃油燃?xì)忮仩t效益分析[J].農(nóng)業(yè)工程學(xué)報(bào)，2011，27(增刊2)：131-134.

　　[7]劉偉，陳琪.煤摻燒生物質(zhì)旋流燃燒器流場(chǎng)的數(shù)值模擬[J].工業(yè)加熱，2011，40(5)：24-27.

　　[8]楊培凱，莊智勇，陳偉博.旋流燃燒器空氣流場(chǎng)特性CFD 數(shù)值模擬[J].機(jī)電工程技術(shù)，2015，44(1)：23-28.

　　[9]CUI K，LIU B，WU Y X，et al.Numerical simulation of oxy-coal combustion for a swirl burner with EDC mo-del[J].Chinese Journal of Chemical Engineering，2014，22(2)：193-201.

　　[10]BUCHMAYR M，GRUBER J，HARGASSNER M，et al.A computationally inexpensive CFD approach for small-scale biomass burners equipped with enhanced air staging[J].Energy Conversion and Management，2016，115：32-42.

　　[11]CHEN Y C，SHEN G F，SU S，et al.Efficiencies and pollutant emissions from forced-draft biomass-pellet semi-gasifier stoves：comparison of international and Chinese water boiling test protocols[J].Energy for Sus-tainable Development，2016，32：22-30.

　　[12]GÓMEZ M A，PORTEIRO J，PATIÑO D，et al.Eule-rian CFD modelling for biomass combustion transient simulation of an underfeed pellet boiler[J].Energy Conversion and Management，2015，101：666-680.

　　[13]JEONG H J，HWANG I S，PARK S S，et al.Investiga-tion on co-gasification of coal and biomass in Shell gasi-fier by using a validated gasification model[J].Fuel，2017，196：371-377.