氮氣稀釋富氫合成氣高壓燃燒特性數值模擬

第43卷第11期熱力發(fā)電Vol 43 No2014年11月THERMAL POWER GENERATIONNov.2014氮氣稀釋富氫合成氣高莊燃嬈痔性數值模擬付忠廣,盧可,周揚,朱一鳴,劉雪琦(華北電力大學(xué)國家火力發(fā)電工程技術(shù)研究中心,北京102206)摘要]為研究高壓下氮氣稀釋富氫合成氣的微混合燃燒規律,利用常壓燃燒試驗校正模型對不同壓力、功率和當量比下的溫度場(chǎng)、火焰形態(tài)以及污染物排放進(jìn)行了計算分析。結果表明,隨著(zhù)壓力的升高,燃燒器火焰溫度和岀口溫度先升髙后略微下降,火焰寬度逐漸變小,火焰高度有增大趨勢,NO,排放指數升高,但在壓力>14MPa時(shí),趨勢變緩。[關(guān)鍵詞]lGCC;氮氣稀釋;富氫合成氣;高壓燃燒;微預混;NO排放;數值模擬中圖分類(lèi)號]TK16[文獻標識碼]A[文章編號]1002-3364(2014)11-0019-05[DOⅠ編號]10.3969/j.issn.1002-3364.2014.11.019Numerical study on combustion characteristicsof nitrogen diluted hydrogen-rich syngas at high pressuresFU Zhongguang, LU Ke, ZHOU Yang, ZHU Yiming, LIU XueqiNational Thermal Power Engineering & Technology Research Center, North China Electric Power University, Beijing 102206,ChinaAbstract Aiming at investigating the micro-mixing combustion characteristics of nitrogen diluted hydrogen-rich syngas at high pressures, the combustion model corrected at atmospheric pressure was adopted to analyze the temperature field flame shape and pollution emissions under conditions with different pressurespowers and equivalent ratios. The results show that, with an increase in pressure, the flame temperatureand outlet temperature of the burner rose first and then dropped slightly; the flame width decreased gradually while its height grew; the No emission indexes increased and tended to be smooth when the pressureincreased to higher than 14 MPKey words: IGCC, nitrogen dilution, hydrogen-rich syngas, high-pressure combustion, micro-mixing, NO,emission, numerical simulation整體煤氣化聯(lián)合循環(huán)(lGCC)是未來(lái)能源清潔硏究了在1~60個(gè)大氣壓下壓力和重力對甲烷-空高效利用的重要發(fā)展方向。燃氣輪機作為lGCC系氣同向層流擴散火焰結構的影響。 McCrain l i統中最重要的一環(huán),成為目前研究的重點(diǎn)。由于材等利用激光誘導熾光法(LI測量甲烷/乙烯燃燒料的限制,依靠提髙燃氣初溫提升燃氣輪機及其聯(lián)的層流擴散火焰中壓力提高對煙灰產(chǎn)量的影響,并合循環(huán)效率的方法具有很大的局限性。提高壓力成觀(guān)測到了壓力對火焰形態(tài)的影響。 Kevin a thom為另一條途徑。而釆用氮氣稀釋的微混合燃燒也是son等利用光譜煙灰排放(SSE)和視線(xiàn)衰減降低富氫燃料燃燒排放的一般手段。 Marc j(L.OSA)測量法對0.5~4MPa內非預混的甲烷/空Charest等(通過(guò)數值模擬求解未修正全耦合方程,氣中國煤化工度進(jìn)行了研究。煙灰收稿日期:2014-06-30THCNMHG基金項目:中央高?；究蒲袠I(yè)務(wù)費專(zhuān)項資金資助(13XS09,2014ZZD0)作者簡(jiǎn)介:付忠廣(1963—),男,漢族,北京人,博士,教授,主要從事潔凈燃料發(fā)電技術(shù)的教學(xué)與研究。Fuze熱力發(fā)電2014年溫度測量結果顯示整體溫度隨壓力提高而増大,但為煙氣岀口,呈圓柱形狀。在圓柱側面設置了水平火焰高度越大時(shí),溫度增長(cháng)速率會(huì )變慢。 Hyun I.取樣孔,用來(lái)收集尾部煙氣以及測量煙溫Joo等研究了在寬壓力范圍中壓力對同流甲烷空1.2計算模型氣層流擴散火焰溫度場(chǎng)結構和煙灰形成特性1)首先采用 Gambit軟件建立燃燒室模型并進(jìn)為了進(jìn)一步提高現有IGCC系統效率,采用高行網(wǎng)格劃分(圖2)。其中三維矩形燃燒室的尺寸為壓力、提高機組參數也是一個(gè)非常有意義的發(fā)展方200mm×200mm×300mm,底部由25組旋流噴向。通過(guò)提高進(jìn)入燃燒室合成氣及空氣的壓力,從嘴組成,按照5×5矩陣分布,每組噴嘴包含8個(gè)空而提高進(jìn)入透平燃氣的壓力,可以有效地增大透平氣噴口,均傾斜20°。網(wǎng)格劃分采用結構化六面體的輸出功率。然而,采用高壓力則意味著(zhù)高壓燃燒,網(wǎng)格,在燃燒器附近采用加密處理,網(wǎng)格總數量約為高壓力下合成氣的燃燒將不同于常壓燃燒,相同工106萬(wàn)況下燃燒更為劇烈,燃燒溫度更高且更難以控制。由于富氫合成氣燃燒本身存在NO污染問(wèn)題,因而在高壓下污染物產(chǎn)生可能會(huì )更多。對此,本文通過(guò)數值模擬研究了高壓力下氮氣稀釋富氫合成氣的燃燒特性及燃燒規律。數值計算模型1.1燃燒器結構以圖1所示的25孔微混合燃料噴射旋流燃燒器為例,對其常壓下的燃燒進(jìn)行試驗研究。該燃燒器共設置了25組射流噴嘴,成5×5矩陣型排列,且每組噴嘴中的空氣噴口均傾斜20°,以產(chǎn)生旋流效圖2燃燒室模型及網(wǎng)格果。每個(gè)噴嘴均由1個(gè)中心燃料小噴口和周?chē)h(huán)繞Fig 2 Mesh generation of the combustor model的8個(gè)空氣噴口組成。燃料噴口直徑為2mm,空2)將25個(gè)中心噴口設置為燃料速度進(jìn)口,8×氣噴口直徑為1.22mm,燃料噴口與空氣噴口的距25個(gè)旋流噴口設置為空氣速度進(jìn)口,矩形燃燒室的離為6.6mm最上端面為出口。此外,采用壁面函數法對流場(chǎng)邊界壁面進(jìn)行處理,流場(chǎng)中存在的非線(xiàn)性問(wèn)題采用欠取樣口燃料噴口松弛方法。最后對計算的收斂標準進(jìn)行設定,取能量方程的殘差值小于10-6,其余各項殘差值小于空氣噴口200mm單個(gè)噴嘴10石英窗F3)由于流場(chǎng)具有很強的旋流特性,因而根據文獻[7-8]采用標準kε湍流模型進(jìn)行湍流流動(dòng)模擬計算穹氣你4)采用渦耗散模型③對燃燒進(jìn)行數值模擬,則反應r中物質(zhì)i的產(chǎn)生速率R為@你Ri.r= min vi. Ap mingY噴嘴分布Ap,M、ABp2Y(1)圖1噴嘴微混合燃料噴射燃燒器結構式中仟何,種產(chǎn)物的化學(xué)組分;YR為某種產(chǎn)Fig. 1 Structure of the micro-mixing fuel物中國煤化分子量;v,和v"分injecting combustor and injector別CNMH反應物和生成物化學(xué)燃燒室為立方體結構,在中間部分的3個(gè)面設恰當比系數;A、B為經(jīng)驗常數,分別取4.0和0.5。有耐高溫石英玻璃觀(guān)察窗,以便觀(guān)察和測量。上部5)采用P-1輻射模型進(jìn)行合成氣燃燒模擬計http://www.rlfd.comcnhttp:/rlfd.periodicals.net.cn第11期付忠廣等氮氣稀釋富氫合成氣高壓燃燒特性數值模擬算。采用灰氣體加權平均模型( weighted-sum ofgray gases model, WSGGM)計算氣體的吸收系數,1300該模型根據壓力以及對輻射有吸收能力的氣體組分的濃度來(lái)計算氣體的吸收系數10口丑100015 kW15kw20kW10020kW1.3模型驗證25kW25kw利用現有的合成器常壓燃燒試驗臺對常壓燃燒0.50.60.70.80.91.00.50.60.70.80.91.0當量比當量比進(jìn)行試驗,監測不同當量比下的出口溫度,并與根據a)壓力5MPab)壓力14MP上述模型計算出的出口溫度進(jìn)行對比,結果見(jiàn)圖3圖4燃燒器出口溫度與當量比的關(guān)系Fig.+ The influence of equivalent ratio on outlet temperature2.2溫度與壓力的關(guān)系在一定的當量比下,出口溫度和火焰溫度隨壓力的變化規律如圖5所示。從圖5可見(jiàn):在0.1--15kW計算值l0MPa范圍內,出口溫度和火焰溫度隨著(zhù)壓力的增20kW計算值加而迅速升高;在10~14MPa范圍內,壓力的提高-25kW計算值使復合反應增強,輻射損失增大,出口溫度及火焰溫◆25kW實(shí)測度的升高速率減慢,溫度幾乎保持不變甚至略有下當量比降;當壓力大于14MPa時(shí),燃燒室內輻射損失快速圖3常壓下出口溫度與當量比關(guān)系增大,火焰溫度顯著(zhù)下降。這是因為根據澤爾多維Fig 3 Relationship between the outlet temperature奇一弗蘭克一卡門(mén)尼茨基提出的預混火焰傳播速度and equivalent ratio at atmospheric pressure與壓力的關(guān)系取決于化學(xué)反應級數,對于大多碳氫從圖3可見(jiàn),在15、20、25kW3個(gè)功率下,出燃料,反應級數小于2,其火焰速度隨壓力增加而降溫度隨當量比的變化規律一致,即隨當量比的提高,低。輻射的影響會(huì )隨著(zhù)速度梯度的降低越來(lái)越明出口溫度升高。從數值上看,本文的數值模型在常顯,因此在考慮熱輻射損失的條件下,隨著(zhù)壓力的升壓下的計算結果比試驗結果略高,這是因為數值計高,火焰溫度先升高然后逐漸平緩下降算未計算環(huán)境的熱損失?？梢?jiàn),數值計算與試驗1200結果基本吻合,因此認為本文建立的數值模型是合1000理的。15kW15kw2結果分析-20kw+25kW+25kW利用數值模擬 Fluent軟件對當量比分別為1、0.83、0.71、0.63、0.5,功率分別為15、20、25kW,壓壓力/MPa壓力MPaa)當量比0.63力分別為0.1、1、5、10、14、20MPa工況下的燃燒器b)當量比0.71溫度場(chǎng)、火焰形態(tài)以及NO排放進(jìn)行模擬計算,分析在高壓下氮氣稀釋富氫合成氣的燃燒特性165016002.1溫度與當量比的關(guān)系圖4為壓力分別在5、14MPa下燃燒器出口溫15kW綏150015kW20kW度與當量比的關(guān)系。從圖4可見(jiàn),在不同壓力下,出-25 kw25kW1350口溫度均隨著(zhù)當量比的增加而升高,這是因為隨著(zhù)當量比的増加空氣流量減少,旋流射流速度降低,使中國煤化工壓力/MPa得氣體在燃燒器內的停留時(shí)間增加,反應時(shí)間變長(cháng)CNMHGd)當量比0.71ml火焰溫度的影響同時(shí),空氣流量減少也會(huì )使熱量隨空氣擴散的速度Fig. 5 The influence of pressure on outlet temperature變慢,因此溫度升高and flame temperaturehttp://www.rlfd.comcnhttp:/rlfd.periodicals.net.cn熱力發(fā)電2014年2.3火焰半徑較大貢獻,因此NO生成速率和排放指數仍然呈上圖6為壓力對燃燒器火焰寬度的影響。從圖6升態(tài)勢,但趨勢逐漸變緩?？梢?jiàn),壓力越大,火焰的寬度越小。這是因為高壓會(huì )1400使火焰的弗勞德數變小,浮力作用增大,浮力對火焰當量比0.71,功率25kW當量比0.83,功率20kW的影響作用超過(guò)重力,因此在浮力的作用下火焰向中心線(xiàn)靠攏并向高處延伸,火焰半徑逐漸變小。800+15kw204001620kW日1815kW起蘭25kW日16-20KW5kW圖8NO4排放指數與壓力關(guān)系壓力/MPaFig8 The relationship between NO emission index and pressurea)當量比0.71B)當量比0.83當量比0.71,功率25kW圖6壓力對火焰寬度的影響當量比0.71,功率25kW,熱力型NO25kW快速型NOFig. 6 The influence of pressure on flame radiusE604-+當里比3率kW熱力0火焰高度圖7為延燃燒器軸線(xiàn)的溫度分布。從圖7可見(jiàn),火焰最高溫度點(diǎn)的高度隨著(zhù)壓力的增大而變高。這說(shuō)明在一定當量比和功率下,火焰高度隨著(zhù)壓力的提高而增大。這是因為壓力越大,弗勞德數越小浮力對火焰的作用力越強,火焰向上延伸距離增大,壓力/MPa火焰高度則越高。本文計算時(shí)未考慮煙灰,因此火圖9NO2生成速率與壓力關(guān)系焰高度未做精確性統計。Fig 9 The relationship between NO, formation rate and pressure3結論1600·“sa1)提高壓力,火焰溫度和出口溫度先升高后略微下降。其中壓力低于10MPa時(shí),溫度升高速率1000較快;壓力在10~14MPa時(shí),溫度升高速率降低IMPa75MPa溫度基本不變甚至略有下降;壓力高于14MPa時(shí)b 10 MPa414 MPa溫度有所下降。在高壓下,提高當量比,出口溫度升●20Mpa火焰徑向高度/m2)提高壓力,火焰寬度逐漸變小,火焰高度有增圖7不同壓力下火焰徑向溫度分布大趨勢,這是由弗勞德數隨壓力提高而變小所致Fig. 7 The flame radial temperature distribution3)隨著(zhù)壓力的升高,NO的排放指數均升高。at different pressures在壓力高于14MPa時(shí),熱力型NO2的生成速率降低,而快速型NO,依然升高較快,因此壓力高于2.5NO,生成特性14MPa時(shí)NO生成速率和排放指數仍然呈上升態(tài)圖8、圖9分別為NO2排放指數及NO生成速率與壓力關(guān)系。從圖8、圖9可見(jiàn),隨著(zhù)壓力的升中國煤化工高,NO的排放指數均升高;在壓力>14MPa時(shí),CNMHGL爹又獻]熱力型NO)的生成速率降低,快速型NO依然升[1]趙曉燕,李祥晟,豐鎮平,等,燃氣輪機低熱值合成氣燃高較快,而快速型NO,對總的NO生成速率有著(zhù)燒室內三維湍流流動(dòng)的數值模擬研究[].動(dòng)力工程,http://www.rlfd.comcnhttp:/rlfd.periodicals.net.cn第11期付忠廣等氮氣稀釋富氫合成氣高壓燃燒特性數值模擬2009,29(4):330-334NO, emission of a multi-stage self- preheatingrnerZHAO Xiaoyan, LI Xiangsheng, Feng Zhenping. N[J. Power Engineering, 2012. 32(10): 765-769lerical simulation research on three dimensional tur-[7]王福珍,王能,黃帆,等.微型燃氣輪機富氧燃燒室數值bulence flow in a gas turbine combustor burning low研究[J].熱力發(fā)電,2014,43(4):76-heat value syngasLJ. Power Engineering. 2009. 29(4)WANG Fuzhen, WANG Neng, HUANG Fan, Numeri[2] Charest M J, Groth C, Gulder O L. Effects of gravitymicro gas turbines [J].Thermal Power Generationand pressure on laminar coflow methane-air diffusion2014,43(4):76-81flames at pressures from 1 to 60 atmospheres[J]. Com- [8 Shih T H, Liou WW, Shabbir A, et al. A new k-e eddybustion and Flame, 2011, 158(5):860-875viscosity model for high reynolds number turbulent[3 McCrain L L, Roberts W L Measurements of the sootows[J]. Computers & Fluids, 1995,24(3): 227-238volume field in laminar diffusion flames at elevated[9]張永生.同向和反向合成氣旋流擴散燃燒研究[J.中pressures[J]. Combustion and Flame, 2005.140(1):60-國電機工程學(xué)報,2009,29(14):63-68ZHANG Yongsheng. Study on the co-swirl and count-[4 Thomson K A Gilder O L, Weckman E J, et al. Sooter-swirl combustion of the syngas diluted diffuse flarconcentration and temperature measurements in co-an[J. Proceedings of the CSEE, 2009. 29(14):63-68.nular, nonpremixed CH,/air laminar flait pres-[10]張引弟,周懷春,謝明亮,等.湍流擴散火焰的大渦模擬sures up to 4 MPa[J]. Combustion and Flame, 2005及驗證[].動(dòng)力工程,2009,29(5):451-455140(3):222-232ZHANG Yindi, ZHoU Huaichun, XIE[5 Joo H I. Gulder O L.ootLarge eddy simulation and verification for turbulentfield structure in co-flow laminar methane-air diffusiondiffusion flame [JI. Power Engineering, 2009.29(5flames at pressures from 10 to 60atmLJ. Proceedings451-455of the Combustion Institute. 2009.32(1):769-775[1]朱紅鈞,林元華. Fluent12流體分析及工程仿真LM]6]陳冬林,成珊,贠英,等.預熱室結構對多段式自預熱燃北京:清華大學(xué)出版社,20燒器內燃燒及NO,排放特性的影響[J].動(dòng)力工程學(xué)ZHU Hongjun, LIN Yuanhua Fluent 12 fluid Analysis報,2012,32(10):765-769,785and engineering simulation[ M]. Beijing Tsinghua UniCHEN Donglin, CHENG Shan, YUN Ying. Effects ofversity press, 2011(in Chinese)the structure of preheating chamber on combustion and00。0000000000000。00000000000000000(上接第11頁(yè)6]張秋生,梁華,胡曉花,基于區域并網(wǎng)的AGC和一次調7]康靜秋,楊振勇,李衛華,主蒸汽溫度控制系統前饋的頻控制回路的改進(jìn)[].熱力發(fā)電,2010,39(7):69-71化[].熱力發(fā)電,2013,42(1):79-8ZHANG Qiusheng, LIANG Hua, HU Xiaohua. RetroKANG Jingqiu, YANG Zhenyong, LI Weihua, et alfit of control system for AGC and primary frequencyOptimization on feed forward controller in main steamconversion based on regional synchronizing and closingmperature control system[J]. Thermal Power Generinto power grid[J]. Thermal Power Generation, 2010ation,2013,42(1):79-8339(7):69-71H中國煤化工CNMHGhttp://www.rlfd.comcnhttp:/rlfd.periodicals.net.cn