低塵生物質(zhì)顆粒燃燒機(jī)內(nèi)燃燒特性的數(shù)值模擬
摘要：介紹了一種用于中小型工業(yè)窯爐的新型低塵燃燒技術(shù)，利用計(jì)算機(jī)數(shù)值模擬考察了低塵旋流生物質(zhì)燃燒機(jī)的特性．在合理選擇氣相流動(dòng)、固相流動(dòng)、煤燃燒及NO的生成等模型的同時(shí)，針對(duì)旋流燃燒場(chǎng)中固體顆粒在壁面附近的碰撞及熔融特性，探討了在壁面處的運(yùn)動(dòng)模型，并以此為基礎(chǔ)考察了燃燒場(chǎng)的兩相流動(dòng)特性，模擬了生物質(zhì)燃燒機(jī)內(nèi)的燃燒過(guò)程及各物理量的分布．在與實(shí)驗(yàn)比較的基礎(chǔ)上，對(duì)生物質(zhì)燃燒機(jī)的結(jié)構(gòu)進(jìn)行了改進(jìn)．結(jié)果表明，在低化學(xué)計(jì)量比下，改進(jìn)后的生物質(zhì)燃燒機(jī)性能更好，顆粒在生物質(zhì)燃燒機(jī)內(nèi)充分燃盡，在保證液排渣效果的同時(shí)，NO的排放遠(yuǎn)低于常規(guī)液排渣旋風(fēng)器的NO排放量．
1  茼言
    在現(xiàn)有的燃煤技術(shù)中，燃煤直接加熱方式的火焰和煙氣中含有大量灰渣，污染嚴(yán)重，水煤漿的形態(tài)類(lèi)似油，但其燃燒特性基本類(lèi)似普通燃燒，火焰和煙氣的含灰量仍很高，目前在鋼鐵等大型行業(yè)中應(yīng)用較多的煤氣化技術(shù)，附屬設(shè)備較多，操作較復(fù)雜，成本較高，一般中小型工業(yè)企業(yè)難以承型1，2]．陳恩鑒等‘31在傳統(tǒng)的液排渣旋風(fēng)燃燒技術(shù)基礎(chǔ)上提出了一種“低塵燃燒技術(shù)”，使在一個(gè)圓柱形的燃燒室中進(jìn)行欠氧燃燒，煤中的灰渣在旋轉(zhuǎn)流場(chǎng)作用下被燃燒室壁面捕集，在高溫下呈液態(tài)渣除去，進(jìn)入爐膛的是潔凈的高溫還原性火焰，將這種燃燒技術(shù)應(yīng)用于燃油工業(yè)窯爐進(jìn)行以煤代油改造及對(duì)老式燃煤工業(yè)窯爐進(jìn)行技術(shù)改造，具有廣闊的前景，
    圖1為低塵生物質(zhì)燃燒機(jī)示意圖，一次風(fēng)r占10%)采用壓縮空氣濃相送粉方式，在端部的葉片中與二次風(fēng)混合后送入生物質(zhì)燃燒機(jī)內(nèi)；熱風(fēng)進(jìn)入生物質(zhì)燃燒機(jī)后，一路作為三次風(fēng)r占30%)進(jìn)入生物質(zhì)燃燒機(jī)出口處，另一路經(jīng)過(guò)生物質(zhì)燃燒機(jī)筒壁，冷卻筒壁，然后經(jīng)過(guò)環(huán)型葉片，作為二次風(fēng)進(jìn)入生物質(zhì)燃燒機(jī)內(nèi)，形成的液渣由燃燒室前部的出渣口排出，與傳統(tǒng)的液排渣生物質(zhì)燃燒機(jī)相比，這種低塵生物質(zhì)燃燒機(jī)在開(kāi)發(fā)構(gòu)思上有以下特點(diǎn)：(1)首次提出取消燃燒室側(cè)壁上的二坎風(fēng)進(jìn)口，采用端面預(yù)旋技術(shù)，在燃燒室內(nèi)部形成對(duì)稱(chēng)且無(wú)阻礙的旋轉(zhuǎn)流場(chǎng)；(2)把噴入外層旋轉(zhuǎn)氣流，延長(zhǎng)了的停留時(shí)間，消除了因“短路”而逸出燃燒器的機(jī)會(huì)；(3)采用分級(jí)燃燒技術(shù)，在生物質(zhì)燃燒機(jī)內(nèi)進(jìn)行局部欠氧燃燒，富含未燃成分的高溫燃?xì)庠谘a(bǔ)入三次風(fēng)后進(jìn)入窯爐進(jìn)行二次燃燒，一方面分級(jí)燃燒在生物質(zhì)燃燒機(jī)內(nèi)形成還原氣氛有利于降低渣熔點(diǎn)，便于排渣，另一方面可抑制氮氧化物的生成‘4]．蔣利橋等‘51對(duì)冷態(tài)流場(chǎng)進(jìn)行了測(cè)試，林伯川等‘61進(jìn)行了初步的熱態(tài)中間實(shí)驗(yàn)，證實(shí)了這種液排渣低塵生物質(zhì)燃燒機(jī)具有連續(xù)排渣、灰捕集率高、低氮氧化物排放等優(yōu)點(diǎn)，但這種燃燒方式的燃燒強(qiáng)度高于普通的生物質(zhì)燃燒機(jī)，尤其是為了液態(tài)排渣，壁面附近溫度較高，對(duì)材料的耐熱性和減少熱損失都提出了很高要求，為了實(shí)現(xiàn)在液排渣前提下優(yōu)化生物質(zhì)燃燒機(jī)的結(jié)構(gòu)和操作條件，把握爐內(nèi)的傳熱特性、的運(yùn)動(dòng)和燃燒行為是十分必要的，本研究對(duì)液排渣生物質(zhì)燃燒機(jī)的燃燒特點(diǎn)進(jìn)行了計(jì)算機(jī)模擬與解析．
2  教學(xué)模型
    對(duì)有反應(yīng)的湍流氣一固兩相流動(dòng)及的燃燒過(guò)程采用拉格朗日一歐拉混合模型，在歐拉坐標(biāo)系下求解氣相的質(zhì)量守恒、動(dòng)量守恒和能量守恒方程，顆粒運(yùn)動(dòng)采用拉格朗日法計(jì)算，忽略氣體密度脈動(dòng)、阻力脈動(dòng)及重力項(xiàng)，對(duì)氣一固相流動(dòng)采用k- 8/RNG模型和拉格朗日法顆粒隨機(jī)軌道模型；煤的燃燒考慮熱解、焦炭燃燒、氣相燃燒三部分，分別采用雙平行競(jìng)爭(zhēng)反應(yīng)模型、擴(kuò)散一動(dòng)力燃燒模型、EBU-Arrhenius模型；熱力型NO的生成采用Zeldovich模型，燃料型和快速型NO的生成采用De Soete模型；輻射模型采用離散傳播法(DT)．孫學(xué)信‘71對(duì)上述模型的機(jī)理進(jìn)行了詳細(xì)介紹，燃燒場(chǎng)視為穩(wěn)態(tài)軸對(duì)稱(chēng)，在二維柱坐標(biāo)系下，各物理量的控制方程可表示為以下通用形式：式中妒代表速度、溫度、濃度等不同的物理量，/表示廣義擴(kuò)散系數(shù)，S(p是氣相湍流引入的源項(xiàng)，S(pp是考慮顆粒影響的源項(xiàng)，其具體含義及k- 8/RNG湍流模型參數(shù)見(jiàn)文獻(xiàn)[8]．
3  數(shù)值解法及計(jì)算條件
3.1數(shù)值解法
    由于低塵生物質(zhì)燃燒機(jī)的出口條件較復(fù)雜，為便于給出出口邊界條件，在生物質(zhì)燃燒機(jī)后部加上了爐膛的簡(jiǎn)化模型，二者同場(chǎng)求解，圖2為整場(chǎng)的計(jì)算網(wǎng)格，其中燃燒器部分設(shè)置r方向網(wǎng)格數(shù)為56，z方向網(wǎng)格數(shù)為60，在近壁面物理量變化急刷的地方采用加密網(wǎng)格，由于爐膛不是重點(diǎn)考察對(duì)象，所以計(jì)算網(wǎng)格較稀，計(jì)算平臺(tái)用通用CFD軟件Star-CD3.10，對(duì)氣相控制方程組用有限體積數(shù)值解法穩(wěn)態(tài)求解，應(yīng)用MARS二階差分格式，在計(jì)算域的各控制微元體內(nèi)積分各控制方程，得到變量在節(jié)點(diǎn)處的差分方程，差分方程的求解采用p-v修正的SIMPLE算法，用代數(shù)多重網(wǎng)格法(Algebraic Multi-gridMethod)加速收斂‘9]，考慮顆粒和氣流間的相互作用，兩相耦合循環(huán)迭代，
3.2進(jìn)出口邊界條件
    在實(shí)際運(yùn)行中，的進(jìn)料點(diǎn)在環(huán)形葉柵附近r圖1)，且一次風(fēng)量較小(占總風(fēng)量10%)，因此計(jì)算中將一次風(fēng)和二次風(fēng)合并在一起考慮，本研究重點(diǎn)考察生物質(zhì)燃燒機(jī)內(nèi)的燃燒情況，所以暫時(shí)未考慮窯爐爐膛內(nèi)引入的三次風(fēng)，二次風(fēng)進(jìn)風(fēng)通道為口徑20 mm的環(huán)形葉柵通道r圖2)，葉片偏轉(zhuǎn)角約85.5。，z=0的其他區(qū)域?yàn)楸诿妫諝夂偷钠渌M(jìn)口條件如表1所示，爐膛出口取定壓邊界條件，絕對(duì)壓力設(shè)為0.1 MPa．
3.3壁面邊界條件
    熱態(tài)下這種生物質(zhì)燃燒機(jī)“壁面”處的邊界祭件較復(fù)雜，在燃燒過(guò)程中，在靠近生物質(zhì)燃燒機(jī)側(cè)壁附近投入，壁面附近的燃燒強(qiáng)度高，在壁面處形成高粘度、高粗糙度的渣層，對(duì)于氣相，粗糙度的提高增加了流動(dòng)的阻力，對(duì)于顆粒運(yùn)動(dòng)而言，顆粒運(yùn)動(dòng)到壁面時(shí)，會(huì)由于熔渣層的粘附作用而被捕集難以再回到氣相，沒(méi)有燃盡的將附在渣層上燃燒，現(xiàn)階段的計(jì)算還沒(méi)有考慮渣層的運(yùn)動(dòng)特性，對(duì)這種復(fù)雜的氣固壁面流動(dòng)和燃燒現(xiàn)象，本計(jì)算針對(duì)實(shí)際熔渣的表面特性假設(shè)計(jì)算中的“壁面”條件，采用改變近壁處的壁面函數(shù)和顆粒碰撞條件的方法來(lái)進(jìn)行合理的模擬．
    (1)氣相流動(dòng)采用無(wú)滑移的壁面邊界條件，近壁面區(qū)域采用以下完全粗糙條件下的壁面函數(shù)‘101：
    Launder等[10]將特征厚度D定義為：當(dāng)與壁面的距離降到y(tǒng)o+D時(shí)，氣流的速度減為0（如圖3所示）．考慮到與顆粒直接作用的表面是粗糙的熔渣表面，參數(shù)的取值如表2所示，增加壁面粗糙度的目的在于模擬氣流由于熔渣層的粘滯造成的減速，顆粒進(jìn)入這個(gè)“滯止區(qū)”以后速度將大幅度降低，
    (2)對(duì)生物質(zhì)燃燒機(jī)部分，氣相壁面?zhèn)鳠岵捎玫谌?lèi)邊界條件，外壁面的冷卻風(fēng)溫580 K，冷卻介質(zhì)、金屬管壁、生物質(zhì)燃燒機(jī)襯料層、渣層的總熱阻根據(jù)實(shí)驗(yàn)選取0.021112．K/W；對(duì)爐膛部分，壁面?zhèn)鳠岵捎玫谝活?lèi)邊界條件，根據(jù)設(shè)計(jì)溫度，取1500K．
    (3)固相顆粒和壁面間碰撞采用完全非彈性碰撞條件，即顆粒與壁面發(fā)生碰撞后速度變?yōu)?，但可被氣流再次攜帶，再次碰壁，如此循環(huán)，直至顆粒燃盡，以此近似模擬在壁面附近的運(yùn)動(dòng)情況．
3.4煤種特性
    現(xiàn)場(chǎng)實(shí)驗(yàn)中采用的大同混煤屬于煙煤，參照普通煙煤在高溫及低溫下的熱解動(dòng)力學(xué)參數(shù)，煤中揮發(fā)分的析出采用雙平行一級(jí)反應(yīng)模型‘7]．在焦炭燃燒中，假定炭表面以生咸CO為主；在氣相燃燒中，考慮了揮發(fā)分的兩步燃燒過(guò)程，揮發(fā)分的元素組成可由元素分析結(jié)果求出，本計(jì)算中將其描述為C。Hf，其中e=0.6758，f=0.6483.反應(yīng)的動(dòng)力學(xué)參數(shù)見(jiàn)文獻(xiàn)[11]．煤樣的元素和工業(yè)分析如表3所示，煤樣的粒度分布如表4所示，煤顆粒密度1200 kg/m3．本計(jì)算根據(jù)實(shí)際的粒度分布，設(shè)置了0.025，0.05，0.075，0.125，0.175 mm等5個(gè)粒徑，每個(gè)粒徑根據(jù)投入位置的不同平均分為8組，
4  結(jié)果與討論
4.1計(jì)算值與實(shí)驗(yàn)值的比較
    為了檢驗(yàn)?zāi)Ｐ驮O(shè)置的準(zhǔn)確性，圖4給出了a=0.77，聊。。。i=0.051 kg/s條件下實(shí)測(cè)溫度和計(jì)算溫度的比較（測(cè)圖5生物質(zhì)燃燒機(jī)尾都NO濃度隨化學(xué)計(jì)量比的變化Fig.5 NO concentration at the rear of combustor點(diǎn)位于生物質(zhì)燃燒機(jī)內(nèi)壁面約25 mm，即r=175 mm處，沿軸向設(shè)置5個(gè)測(cè)量熱電偶）．圖5給出了不同化學(xué)計(jì)量比下生物質(zhì)燃燒機(jī)尾部NO標(biāo)態(tài)濃度的變化
    從比較結(jié)果可以看出，計(jì)算值和實(shí)驗(yàn)值吻合較好，最大相對(duì)誤差在10%以?xún)?nèi)，說(shuō)明關(guān)于顆粒近壁條件的設(shè)置是比較合理的，數(shù)值模擬結(jié)果也能夠反映生物質(zhì)燃燒機(jī)內(nèi)部的燃燒特性，以此為基礎(chǔ)，對(duì)生物質(zhì)燃燒機(jī)進(jìn)行設(shè)計(jì)優(yōu)化是可行的．
4.2燃燒特性分析
4.2.1流場(chǎng)分布
    圖6為a=l.0時(shí)生物質(zhì)燃燒機(jī)內(nèi)的縱截面速度矢量分布，圖7是生物質(zhì)燃燒機(jī)的截面流線圖，不同化學(xué)計(jì)量比下的流場(chǎng)形態(tài)相似，大小有所差異，與冷態(tài)測(cè)試結(jié)果‘51不同，在燃燒條件下，中心氣流受熱膨脹，中心回流區(qū)消失，環(huán)室回流區(qū)發(fā)生阻斷，主氣流區(qū)擴(kuò)張，在生物質(zhì)燃燒機(jī)的中后部和中心氣流區(qū)合并，環(huán)室回流區(qū)大大減小，因此在實(shí)際燃燒過(guò)程中，“環(huán)室回流區(qū)”將不足以使顆粒產(chǎn)生往復(fù)運(yùn)動(dòng)，絕大部分顆粒受強(qiáng)旋流作用集中在壁面附近，壁面附近局部氧濃度很低，顆粒很快運(yùn)動(dòng)到“無(wú)氧區(qū)”f圖8)，炭粒的燃燒效果不好，在化學(xué)計(jì)量比較低∞=0.7，0.8)時(shí)，燃燒率僅為70.1%和80.7%．本燃燒系統(tǒng)若采取分級(jí)燃燒、而生物質(zhì)燃燒機(jī)內(nèi)采用低化學(xué)計(jì)量比的燃燒方式，如何最大程度地提高炭粒的燃燒率是一個(gè)很重要的問(wèn)題．
4.2.2濕度場(chǎng)分布
    圖9為a=0.7～1.1時(shí)生物質(zhì)燃燒機(jī)內(nèi)溫度場(chǎng)的分布，隨著燃燒的進(jìn)行，溫度沿軸向先升高后略有降低，比較不同的化學(xué)計(jì)量比下的溫度場(chǎng)分布，低a(0.7，0.8)比高a(l.0，1.1)時(shí)的高溫區(qū)更靠前，在生物質(zhì)燃燒機(jī)端面觀察孔處的肉眼觀測(cè)結(jié)果也證明了這一點(diǎn)，從總體來(lái)看，燃燒區(qū)過(guò)分靠后，出口截面溫度都很高(1866～2032 K)，對(duì)NO的控制很不利，在實(shí)際設(shè)計(jì)過(guò)程中，希望主燃燒區(qū)能提前，以利于降低出口截面的溫度，這樣當(dāng)爐膛內(nèi)三次風(fēng)補(bǔ)入，NO的總生成量不會(huì)再次大幅度升甜12]．因此，燃燒器結(jié)構(gòu)需要進(jìn)一步改進(jìn)，
4.2.3低塵生物質(zhì)燃燒機(jī)灰渣的捕集
    投入生物質(zhì)燃燒機(jī)后，強(qiáng)烈的旋流使幾乎所有的顆粒均在壁面附近運(yùn)動(dòng)，旋流的離心分離作用使剩余顆粒難以隨氣流運(yùn)動(dòng)而被壁面熔渣層捕集，雖然實(shí)驗(yàn)測(cè)到的捕渣率在85%以上，也證實(shí)了生物質(zhì)燃燒機(jī)能實(shí)現(xiàn)液排渣，產(chǎn)生接近重油燃燒的清潔燃?xì)?，但在化學(xué)計(jì)量比較低r0.8）時(shí)，渣中的含碳量達(dá)10%以上，捕渣率的提高不應(yīng)以燃燒效率的降低為代價(jià)．
4.3低塵生物質(zhì)燃燒機(jī)的改進(jìn)
    燃燒效率、液排渣效果及污染物的排放情況是本燃燒器燃燒性能的主要考察因素，但3個(gè)因素間又相互影響，燃燒效率的提高在很大程度上依賴(lài)溫度水平的增加，這對(duì)液排渣有利，但對(duì)NO的排放不利，要控制NO的排放，就要適當(dāng)降低燃燒室出口的溫度水平，營(yíng)造室內(nèi)還原氣氛，這顯然會(huì)影響燃燒效率及液排渣的效果，因而如何合理地平衡三者間的關(guān)系，是低塵生物質(zhì)燃燒機(jī)設(shè)計(jì)的關(guān)鍵問(wèn)題，
    表5給出了不同化學(xué)計(jì)量比下生物質(zhì)燃燒機(jī)內(nèi)一些特性參數(shù)的計(jì)算值比較，在化學(xué)計(jì)量比較低時(shí)，壁面附近更易處于缺氧狀態(tài)，不利于炭粒燃盡，化學(xué)計(jì)量比較高∞=1.0，1.1，1.2)時(shí)，炭粒的燃燒率較高，但出口溫度的提高或出口氧濃度的提高都木利于NO的控制，
    針對(duì)以上一些矛盾和缺點(diǎn)，本研究將生物質(zhì)燃燒機(jī)的尺寸加大，直徑由400 mm改為500 mm，長(zhǎng)度由600 mm改為1000 mm，出渣口的位置延至z=800 mm處，表6給出了a=0.8，聊。。。i=0.026 kg/s時(shí)改進(jìn)前后的部分計(jì)算值比較，從表可以看出，由于環(huán)形進(jìn)風(fēng)通道面積及生物質(zhì)燃燒機(jī)長(zhǎng)度的增加，在同樣質(zhì)量流量下，進(jìn)風(fēng)速度減小，顆粒運(yùn)動(dòng)減慢，停留時(shí)間變長(zhǎng)，顆粒的燃燒率得到提高；主燃燒區(qū)相對(duì)提前，出口截面溫度下降，NO的排放進(jìn)一步降低，在保證順利排渣的前提下，出渣口溫度有所下降，同時(shí)由于尺寸增大，燃煤量可適當(dāng)增加，目前最新一代生物質(zhì)燃燒機(jī)的燃煤負(fù)荷已達(dá)到0.056 kg/s．初步實(shí)驗(yàn)結(jié)果證實(shí)[3]，在a=0.8，塒。。。1=0.056 kg/s時(shí)，炭粒的燃燒率達(dá)99%以土，捕渣率達(dá)90%，渣中的含碳量小于0.1%，出口溫度降到1800 K，NO的排放僅為400 mg/Nm3，遠(yuǎn)低于常規(guī)液排渣爐約850～1150 mg/Nm3的排放量‘13]．進(jìn)一步的研究和開(kāi)發(fā)仍在進(jìn)行中．
5  結(jié)論
    將k-E/RNG湍流流動(dòng)模型、顆粒軌道模型、氣相EBU-Arrhenius模型、焦炭的擴(kuò)散動(dòng)力燃燒模型和各種NO的生成模型相結(jié)合，通過(guò)對(duì)壁面條件的處理，探討了液排渣生物質(zhì)燃燒機(jī)模擬中在壁面附近行為的模型化方法，模擬了低塵生物質(zhì)燃燒機(jī)內(nèi)的燃燒過(guò)程和顆粒運(yùn)動(dòng)特性，為進(jìn)一步建立和完善的熔渣及附壁模型奠定了基礎(chǔ)，結(jié)果表明：
    (1)計(jì)算模型和計(jì)算方法較好地模擬了欠氧狀態(tài)下生物質(zhì)燃燒機(jī)內(nèi)的溫度分布、濃度分布和氮氧化物的生成，采用改變近壁處的壁面函數(shù)和顆粒碰撞條件模擬壁面處的行為是可行的．
    (2)由于燃燒實(shí)驗(yàn)中壁面條件的改變及氣流受熱膨脹加速，冷態(tài)下強(qiáng)旋流造成的多重回流流場(chǎng)結(jié)構(gòu)在熱態(tài)下減弱甚至部分消失．
    (3)燃燒大部分集中在生物質(zhì)燃燒機(jī)壁面附近，大部分灰分可在燃燒過(guò)程中以熔渣的形式除去．

    (4)針對(duì)燃燒效率、液排渣效果及污染物排放三者間的平衡關(guān)系，將生物質(zhì)燃燒機(jī)加大加長(zhǎng)，一次燃燒仍采用欠氧燃燒條件，解決了低化掌計(jì)量比下顆粒燃燒率不高及出口截面溫度過(guò)高的問(wèn)題，同時(shí)保證了液排渣效果和NO的低排放，

生物質(zhì)氣化站，http://www.598jx.com