第35卷第6期中國電機工程學報Vol.35No.6Mar.20,20152015年3月20日Proceedings of the CSeeC2015 Chin. Soc. for Elec Eng. 1429DOI:10.133410258-8013, pose201506.018文章編號:0258-8013(2015)06-1429-07中圖分類號:TK472合成氣燃氣輪機燃燒室CFD模擬的模型選擇及優(yōu)化王翰林,雷福林,邵衛(wèi)衛(wèi),熊燕,張哲巔,肖云漢(中國科學院先進能源動力重點實驗室(工程熱物理研究所),北京市海淀區(qū)100190)Screening and Modification of CFD Models for Syngas Turbine CombustorWANG Hanlin. LEI Fulin SHAO Weiwei. XIONG Yan. ZHANG Zhedian. XIAO Yunhan(Key Laboratory of Advanced Energy and power (Institute of Engineering Thermophysics, Chinese Academy of Sciences)Haidian District, Beijing 100190, China)ABSTRACT: A standard swirling diffusion flame and a scaled關(guān)鍵詞:合成氣燃氣輪杋燃燒室;計算流體力學(CFD);湍syngas turbine were numerically studied using the method of流模型;火焰結(jié)構(gòu);NO,排放RANSto analyze the effects of 4 turbulence models onsimulations. Compared with experimental profiles of velocity, 0 Iftemperature and mixture fraction, effects of models and model整體煤氣化聯(lián)合循環(huán)( integrated gasificationparameters on simulation results were discussed. The mostsuitable realizable K-8 model was chosen, and the constant of combined cycle,lGCC)利用燃氣-蒸汽聯(lián)合循環(huán)實現(xiàn)C2 in the a equation was modified. Coupled with the proper了煤的清潔高效利用,成為未來電力的重要發(fā)展方combustion models and chemical mechanisms, an appropriate向之一,因而發(fā)展合成氣燃氣輪機有重要意義。燃method for CFD simulation of syngas combustor was燒室的壁面溫度對于燃燒室的設(shè)計和使用壽命具d. A scaled syngas turbine combustor experiment wasut. Using the former method, the surface temperatu有重要意義。此外,隨著污染物排放標準的日益嚴O emissions of the scaled syngas turbine combustor格,燃氣輪機燃燒室NO污染物排放控制受到挑were calculated. The computational results show good戰(zhàn)。雷諾平均( Reynolds average Navier-Stocksagreements with experiments. The method was verifiedRANS)計算流體力學( computational fluid dynamics,KEY WORDS: syngas turbine; computational fluid dynamics CFD模擬為預測燃氣輪機燃燒室燃燒及污染物排(CFD); turbulence models; flame structure: NO, emissions放特性提供了一種基礎(chǔ)手段,從而為合理設(shè)計燃燒摘要:采用雷諾平均的計算流體力學方法,通過模擬與燃燒室提供了一條有效途徑。室模態(tài)相似的標準旋流擴散燃燒實驗以及合成氣模型旋流RANS方法的燃燒模擬結(jié)果和計算速度受到所燃燒室實驗,分析了4種湍流模型對旋流燃燒及合成氣燃燒室的數(shù)值模擬結(jié)果的影響。對比標準實驗的速度、溫度以及采用的湍流燃燒模型、化學反應(yīng)模型、湍流模型等混合分數(shù)等實驗數(shù)據(jù),討論不同模型以及參數(shù)取值的模擬效因素影響。合成氣燃燒室內(nèi)的燃燒通常屬于旋流擴果,選取最合適燃燒室內(nèi)復雜流型預測的 realizable kε湍流散燃燒模態(tài)口,這一模態(tài)中常采用的幾種湍流燃燒模型,并對位于湍流動能耗散率(a)方程源項中的模型常數(shù)模型特點不同:渦耗散( eddy dissipation,ED)模型C2做岀修正。耦合相應(yīng)較為準確的燃燒模型以及反應(yīng)機理采用“混合即燃燒”假設(shè),計算快速,然而 Nemitallah形成針對合成氣燃燒室性能預測的較為理想的模擬方法開等采用ED模型模擬非預混火焰的研究表明,此展合成氣模型旋流燃燒室燃燒實驗,利用形成的數(shù)值模擬方法,對燃燒室壁溫及NO排放進行了預測,結(jié)果與實驗較模型會高估火焰尾部的燃燒溫度;概率密度方程吻合,驗證模擬方法可靠。( Probability Density Function,PDF)計算速率相對較快, Ziani等中國煤化工燒火焰的數(shù)基金項目:國家863高技術(shù)基金項目(2008AAO5A302)The National High Technology Research andt of china值研究證明CNMHG預測火焰結(jié)863 Program(2008AA05A302)構(gòu);渦耗散概念( eddy dissipation concept,EDC)模1430中國電機工程學報第35卷型相對精細但計算速度較慢, Fukumoto等的研究1模擬對象結(jié)果表明EDC模型可以準確的計算出燃燒火焰結(jié)1悉尼大學標準旋流燃燒器構(gòu)。根據(jù)燃燒模型的不同特點,耦合相應(yīng)的化學反根據(jù)合成氣燃燒室內(nèi)包括旋流射流等復雜流應(yīng)模型可實現(xiàn)燃燒室工作過程以及污染物排放模型的流動特點,選取模態(tài)與之相似并測量出火焰結(jié)擬預測。燃燒及污染物NO的生成過程化學反應(yīng)復雜,更容易受到化學反應(yīng)動力學的影響,因此化構(gòu)的悉尼大學標準旋流擴散燃燒實驗進行分析,其結(jié)構(gòu)及邊界條件在文獻[18-19中有詳細說明。燃燒學反應(yīng)模型需要足夠精細,才能使模擬結(jié)果更可信。器包含一個直徑60mm環(huán)形腔體,可供空氣圍繞模擬湍流燃燒常用的湍流模型中, standard K-個直徑50m的鈍體進行旋流;燃料噴口位于鈍體模型對旋流和射流都有較理想的效果因此應(yīng)用較圓心,直徑3.6mm;整個燃燒器置于較大的風筒中,多6;RNGκ-模型對旋流的描述更為準確,因而風筒中的空氣流速保持20m/s在旋流燃燒的模擬中有所應(yīng)用; realizableκ-ε模型針對湍流黏度以及ε方程做出修正,對多種流型分別選取文獻[18-19]中冷態(tài)工況N29S054和都有理想效果,在燃燒室設(shè)計以及預測上應(yīng)用廣熱態(tài)工況SM1進行研究。其中冷態(tài)試驗,燃料孔泛0;雷諾應(yīng)力模型( Reynolds stress model,RSM出口速度(U)66ms,旋流切速度(U,)19.lms,旋流數(shù)0.54。熱態(tài)實驗采用甲烷燃料,U為327m/s,是理論上最為精密的模型,對旋流模擬效果較為理U為19.1m/s,旋流數(shù)0.5。該實驗在本研究中只用想Curci等對合成氣湍流燃燒的數(shù)值研究表明,火焰結(jié)構(gòu)的模擬結(jié)果對采用的湍流模型敏感于分析湍流模型對計算結(jié)果的影響,因此燃料并不影響結(jié)論程度高:由于模擬中所采用的湍流模型對于速度耗1.2合成氣模型燃燒室散的描述不同,因此不同的湍流模型適用于不同的流動形態(tài)。此外,為了獲得更理想的預測效果,不本實驗室合成氣模型旋流燃燒室,主要由旋流噴嘴、旋流器和火焰筒組成,置于同軸機匣內(nèi)?？胀瑢W者針對不同的流型對湍流模型做岀了修正,氣由外側(cè)機匣進入燃燒室。采用K型熱電偶測量燃如:Pope等將k-g模型中的模型常數(shù)與流場中的燒室的外壁面溫度,燃燒室外形結(jié)構(gòu)及溫度測點分應(yīng)變率關(guān)聯(lián),取得對射流模擬效果的改進; Khaleghi布如圖1所示。采用PEC-304FT紅外氣體分析儀測等對RNGκ-g模型中E方程源項進行了修正,改量煙氣中O2、NO、NO2、CO、H2O和未燃碳氫(以善了模型對于旋流的模擬效果; Hossain等1l的研CH4、C3H3表征)等組分濃度究表明將 standardκ-g模型中的C1由144修正為1.6可以更好的模擬鈍體火焰結(jié)構(gòu)等等。合成氣燃測點2測點3測點6測點9測點11氣輪機燃燒室具有復雜的幾何結(jié)構(gòu),燃燒室的頭部通常采用旋流結(jié)構(gòu)而壁面上包含主燃孔及摻混孔等結(jié)構(gòu)。因此合成氣旋流燃燒室內(nèi)的流型結(jié)構(gòu)屬于射流、旋流以及交叉射流等幾種流型的復雜結(jié)合,所以有必要針對旋流合成氣燃燒室數(shù)值模擬中測點1測點4測點5測點7測點8湍流模型的影響進行深入分析。圖1燃燒室壁面熱電偶測點分布示意圖本研究在選定較為準確的詳細化學反應(yīng)機理和Fig. 1 Distribution of thermocouples on combustor燃燒模型基礎(chǔ)上,通過模擬與燃燒室模態(tài)相似的悉surfaces尼大學標準旋流擴散燃燒實驗189,對比各種實驗實驗采用熱值12.2MJkg合成氣作為燃料,其數(shù)據(jù),討論前文提到的4種不同湍流模型及部分模組分體積比H2:CO:N2為0.37:0.48:0.15。實驗工型參數(shù)對于模擬結(jié)果的影響;選擇最適合預測合成況設(shè)置如表1,空氣以及燃料進氣溫度均為293K,氣旋流燃燒室燃燒的湍流模型,并對模型參數(shù)做出總過量空氣系數(shù)保持3.83。優(yōu)化,形成較為準確的模擬方法;開展合成氣模型旋流燃燒室實驗研究,對比壁溫以及NO-排放數(shù)據(jù)2數(shù)值模擬中國煤化工對方法進行驗證和應(yīng)用。結(jié)論對合成氣燃氣輪機燃2.1湍流模型CNMHG燒室流場及燃燒性能預估具有一定的參考意義。RANS模擬方法中的動量控制方程可表示為第6期王翰林等:合成氣燃氣輪機燃燒室CFD模擬的模型選擇及優(yōu)化1431表1旋流合成氣模型燃燒室實驗工況反應(yīng)動力學模型,總反應(yīng)模型包含32種組分和176Tab. 1 Operating conditions for the scaled syngas步化學反應(yīng)。combust計算方法進氣壓力MPa燃料流量/(kgh)氣流量/(kg/h)對于標準燃燒實驗,將整個燃燒器燃燒區(qū)域作為計算域,共4×105結(jié)構(gòu)化網(wǎng)格單元。對于模型燃338.737.0燒室,根據(jù)其幾何對稱性,選取燃燒室的四分之0.3546.2564.6作為計算域,并對火焰筒壁面上的冷卻孔進行簡化ap a處理,把每排冷卻孔改為等面積的窄縫,使得計算域具有周期對稱性。整個計算域1.27×10°結(jié)構(gòu)化網(wǎng)格單元。壁面設(shè)置為輻射邊界條件。[(+-5)+(-pu2)模擬計算求解采用 ANSYS FLUENT13.0商用方程中的雷諾應(yīng)力項-p需要本構(gòu)建模,其軟件,其中輻射換熱采用了離散坐標(DO)射模型壓力速度解耦采用SⅠMPLE算法,壓力迭代采用中一種常用的假設(shè)是1877年 Boussinesq提出的渦PRESTO!算法,其他物理量的空間離散格式采用旋黏性系數(shù)的概念,將有效流剪切應(yīng)力pm用二階迎風格式。近壁面采用加強壁面處理 (enhance湍流黏度H以及平均速度的梯度表示出來。為了使 wall treatmen),提高壁面換熱的計算精確度得整個控制方程組封閉,還需要添加其他變量控制物性利用相應(yīng)反應(yīng)模型中的輸運和熱物性數(shù)方程,使得未知量的數(shù)量與方程數(shù)量相等。由此假據(jù)包,采用可壓縮理想氣體方法計算混合氣體密設(shè)出來變量湍流動能k及其耗散率e,并將它們與度;混合定律計算定壓比熱(p);理想氣體混合定湍流黏度利用假設(shè)關(guān)聯(lián),形成 standard k-e模型。律計算熱導系數(shù)()和運動黏度();動力學理論計在此基礎(chǔ)上,后來的研究利用不同方法對方程以算質(zhì)量擴散率(D)和熱質(zhì)量擴散系數(shù)(D1),以此反映及G方程進行了發(fā)展,添加各向異性的假設(shè),形成壓力對物性參數(shù)的影響。了RNGκ-ε模型、 realizableκ-←g模型以及RSM模型。本研究對比RANS法中針對旋流燃燒室常用的3結(jié)果及討論standardκ-ε模型、旋流主導的RNGk-E模型3.1標準旋流燃燒realizableκ-<ε模型以及RSM模型的模擬效果,選3.1.1·不同模型的模擬結(jié)果擇最合適的湍流模型,并進一步做出修正在采用 FLUENT中湍流模型參數(shù)默認值情況2.2燃燒及化學反應(yīng)模型下,使用不同湍流模型模擬分別模擬標準旋流冷態(tài)標準實驗中燃料為甲烷,其燃燒反應(yīng)過程所包實驗N29S054的流場以及熱態(tài)燃燒實驗SM軸線含的反應(yīng)機理步數(shù)較多0,若采用有限反應(yīng)速率的上的軸向速度及溫度分布,對比實驗結(jié)果示于圖2燃燒模型(如EDC模型)耦合詳細反應(yīng)模型計算時圖2中,X和r分別為距燃料出口的軸向高度和徑間過長,因此采用較為快速且準確的PDF燃燒模向距離,D和R分別代表鈍體的直徑和半徑。型3,13耦合GR2.11骨干化學反應(yīng)機制。N29S054的計算結(jié)果表明,4種模型在頭部考慮到合成氣詳細化學反應(yīng)模型組分數(shù)較少,(XD=0.2)的預測效果差別不大,然而對于旋流以能滿足CFD模擬計算軟件的要求,并根據(jù)燃燒室及流速耗散速率的描述不同,使得流場后部(XNO排放過程的特點,模型合成氣燃燒室選用合成06)模擬結(jié)果差異較大。旋流主導的RNGk=g模型氣燃燒及NO生成的詳細化學反應(yīng)機理,耦合EDC和RSM模型對速度衰減預測偏差比較大; standard燃燒模型。已有的硏究表明,EDC模型耦合詳細κ-ε模型在流場中心的流速預測上較接近實驗,但的化學反應(yīng)機理相較于其他模型能夠較為準確的速度最低值及低速區(qū)位置不夠準確;綜合流場中心模擬岀合成氣湍流火焰的火焰結(jié)構(gòu)。反應(yīng)機理采用以及低速區(qū)的計算結(jié)果,采用 realizableκ-ε模型更在多種條件下優(yōu)化后的H/CO燃燒詳細化學反應(yīng)適合預測旋流中國煤化工。計算時間模型2;NO2生成及NO2燃燒相互關(guān)系子模型采用方面,旋流主CNMHGM模型時間Faravelli等人23提出并在多種實驗下驗證的詳細較長, realizable k-8模型和 standardκ-模型計算成1432中國電機工程學報第35卷standard K-8XD=0.2關(guān)系,可以得出結(jié)論,采用 realizable k-ε模型更適realizalbe k-8RNG K-8合模擬旋流燃燒的火焰結(jié)構(gòu),且計算量較小。針對旋流燃燒實驗,其他學者也得到了相似的結(jié)果◇24-25。然而,此模型仍然低估了速度的耗散速率3.1.2 realizableκ-ε模型修正E方程源項中的常數(shù)取值會影響速度耗散的計standard K-g算結(jié)果,對其進行一定的修正優(yōu)化2627,可以提高realizalbe k-8模型對不同流型下速度耗散描述的準確性。選取較為理想的 realizableκ-ε模型做進一步的分析,其ε方程如下為aaat0.4徑向位置rRpCSe-C2pE/(K+vve)+C: CE Gne/K+S(a)N29s054不同軸向高度半徑上的軸向速度分布式中:C1=max[0.43,n/(7+5);n=SK/e;S為應(yīng)變realizalbe K-8率,S=√2SS;Gn為浮力引起的湍流動能:;SRNGRSM為其他源項;σ2為湍流普朗特數(shù):ν為動力黏度,實驗m2/s;C2、Ca1、Ca3為常數(shù)。模型常數(shù)Ca1、Ca3在考慮浮力影響時通常根據(jù)實際情況由經(jīng)驗關(guān)系式計算得出,因此可以針對參數(shù)C2進行優(yōu)化校正,其默認值為1.9。根據(jù)控制方程,增大參數(shù)C2會使得方程計算所得出的ε減小,相應(yīng)的κ隨之增大,湍流動能耗散standard的時間尺度s/κ縮小。這使得計算出的湍流動能更0RNG K-8易維持,湍流擾動增大,速度耗散相應(yīng)加快。因此,RSM增大C2理論上可以優(yōu)化模型對旋流的模擬效果。150將C2分別設(shè)置為1.9、20和2.1,分別計算N29S054軸向位置mm(b)SM火焰中心線速度及溫度分布的流場以及SMI軸線上的軸向速度、混合分數(shù)以圖2不同模型對標準旋流實驗的及溫度分布,其計算結(jié)果及實驗結(jié)果示于圖3模擬結(jié)果與實驗數(shù)據(jù)對比根據(jù)N29S054的計算結(jié)果,修正后的模型對于Fig 2 Comparison between computational results b頭部的影響并不明顯;流場后部的的計算結(jié)果說different models and standard experimental datas明,增大C2使模型預測的速度耗散加快,同一軸本較低向位置上的中心流速降低,低速區(qū)向外移動,徑向SM!的結(jié)果表明,不同模型計算出的低速及回上的速度最低值增大。C2=1.9時,中心速率偏差較流區(qū)分布差別很大,因此溫度分布計算結(jié)果差異較大;C2=2.1時,低速區(qū)的偏差較大;綜合比較,大。 standardκ-ε模型模擬軸向上速度耗散過快,對C2=2.0較為合適。速度的模擬相對其他模型偏差較大。由于在火焰頭SMⅠ的計算結(jié)果也表明,增大C2可以使得整部旋流較強,對于旋流有精細描述的旋流主導的個熱態(tài)流玚的速度衰減加快,回流區(qū)更接近頭部RNGκ-ε模型和RSM模型相對比較準確;而在火回流區(qū)縮小;混合分數(shù)的衰減冋樣加快,說明摻混焰尾部,旋流減弱,射流特性更明顯,因此這2種更強;溫度峰值位置更靠近頭部。軸線上的軸向速模型對火焰尾部的溫度分布、回流區(qū)位置以及大小度分布、回流I凵中國煤化工布以及火焰的預測上相差較大。綜合速度分布、回流區(qū)位置及峰值位置與CNMH,在選擇大小以及火焰溫度分布的峰值位置與實驗值之間的 realizableκ-ε模型模擬與燃燒室內(nèi)模態(tài)接近的火焰第6期王翰林等:合成氣燃氣輪機燃燒室CFD模擬的模型選擇及優(yōu)化1433XD=0.2室燃燒實驗及數(shù)值模擬硏究。對比4種湍流模型耦合前文合成氣燃燒及反應(yīng)模型對模型旋流燃燒室▲實驗壁面溫度的模擬效果。保留默認參數(shù)設(shè)置條件下,模擬進口壓力0.10MPa工況,計算與實驗結(jié)果示于圖4 standardκ-ε模型在頭部的模擬效果相對其他模型偏差較大,壁面溫度最高點更靠前。采用旋流主導的RNGx-E模型和RSM模型在頭部預測效果C2=2.0C2=2,1上比較理想,而尾部的壁面溫度預測上相較于實驗值高出很多,計算出的高溫面積較大,溫度最高點更靠近頭部。 realizableκ-ε模型計算的壁面溫度最大偏差相對較小,壁面溫度值以及壁溫隨軸向位置的變化趨勢都更接近于實驗結(jié)果。因此與前文一0.4徑向位置rR致,相較于其他模型,對射流以及旋流都有較為準(a)N29s054不同軸向高度半徑上的軸向速度分布確的描述的 realizableκ-ε模型對于壁面溫度的預測比較理想,更適于旋流合成氣燃燒室內(nèi)的燃燒過C2程的模擬實驗匾0.4△-C2=2.0測點編號圖4不同湍流模型對燃燒室壁溫的模擬結(jié)果與實驗數(shù)據(jù)對比Fig 4 Comparison of surface temperature between2.0simulations by different models and model combustorexperiments選取燃燒室進口壓力0.10MPa以及0.35MPaC2=2.1這2個工況,采用C2修正為2.0的 realizable k-模型模擬旋流合成氣模型燃燒室壁面溫度,計算結(jié)軸向位置/m果與實驗測量值對比如圖5。模擬結(jié)果表明,將(b)SMl火焰中心線速度、混合分數(shù)及溫度分布realizableκ-E模型C2修正為2.0后,計算所得的壁圖3C2不同取值的 realizableκg模型對標準旋流實驗的面溫度相較于前文的結(jié)果明顯更接近實驗值。只有擬結(jié)果與實驗數(shù)據(jù)對比在0.35MPa工況條件下測點9處偏差相對略大,這Fig.3 Comparison between simulations by realizable K-E主要是由于在計算中,采用空氣逆流進入機匣的簡model with different C2 and standard experiments化,而實驗中則是在機匣的壁面開孔,開孔位置臨結(jié)構(gòu)時,將C2修正為2.0,可以得到較理想的模擬近測點9位置,壓力增大之后,空氣流量增大,因效果,而避免使用計算量較大的精確模擬方法此使得這個區(qū)域的壁面換熱增強,9號測點的溫度3.2合成氣模型旋流燃燒室偏小,偏差增」中國煤化工用C2修正3.2.1壁面溫度為2.0的reaiCNMHG燒模型及相為了驗證前文結(jié)論,開展合成氣模型旋流燃燒應(yīng)的詳細反應(yīng)模擬合成氣模型旋流燃燒室的1434中國電機工程學報第35卷3■1atm實驗要的影響。湍流模型對旋流及速度耗散描述的不口1atm模擬1.0}435am實驗同,使得流場及火焰結(jié)構(gòu)計算結(jié)果存在差異。其中△3.5atm模擬realizableκ-ε模型的計算結(jié)果與標準實驗結(jié)果最接近,最適用于旋流擴散燃燒模擬。島0.42)ε方程中的常數(shù)取值會影響流場速度耗散作用,本研究針對 realizableκ-ε模型常數(shù)C2修正。C2的取值增大可以增大速度耗散率,將C2由1.9測點編修正為2.0,計算標準實驗所得的速度及溫度分布圖5修正后的 realizableκ-ε模型對燃燒室壁溫的模擬結(jié)于實驗最為接近,改進了模型的模擬效果。果與實驗數(shù)據(jù)對比3)合成氣模型旋流燃燒室燃燒實驗驗證了Fig 5 Comparison of surface temperature betweenrealizableκ-ε模型最適用于燃燒室模擬。修正后的simulations by modified realizable K-g model and scaledrealizableκ-ε模型可以改進模擬效果,并能較為理combustor experiments想的預測出NO排放。壁面溫度及火焰結(jié)構(gòu)具有理想的效果322NO2排放4)采用 realizableκ-E湍流模型,將C2修正為利用已經(jīng)形成的合成氣燃燒室計算模型耦合2.0,耦合EDC燃燒模型及文中相應(yīng)的化學反應(yīng)模前文NO2反應(yīng)模型,對不同壓力條件下NO2排放進型,適用于合成氣燃氣輪機燃燒室的數(shù)值模擬,可行了預測,結(jié)果折合為干基15%O2濃度條件下的摩以對旋流合成氣燃氣輪機燃燒室燃燒排放特性預爾分數(shù),對比實驗如圖6。結(jié)果顯示,在實驗工況測及低污染物燃燒室設(shè)計和改進提供指導。條件下,模型計算燃燒室出口NO排放摩爾分數(shù)與參考文獻實驗吻合較為理想。根據(jù)NO排放的特點,除NO[]張永生,穆克進,張哲巔,等.同向和反向合成氣旋流化學反應(yīng)模型外,NO排放的計算準確性很大程度擴散燃燒研究.中國電機工程學報門.2009,29(14):上取決于燃燒流場及溫度分布計算的準確性。由于NO2反應(yīng)模型在不同研究中均有驗證,因此進一步Zhang Yongsheng, Mu Kejin, Zhang Zhedian, et al. Studyon the co-swirl and counter-swirl combustion of the證明模型對于燃燒室內(nèi)的流動及燃燒過程模擬syngas diluted diffuse flame[J]. Proceedings of the Csee可靠。[2] Nemitallah M A, Habib MA. Experimetal and numerical▲模擬Investigationsatmospheric十△實驗oxy-combustion flame in a gas turbine monbustor[J. Applied ene[3] Ziani L, Chaker A, Chetehouna K, Malek A. Numericalsimulations of non-premixed turbulent combustion ofCHa- H mixtures using the PDF approach. InternationalJournal of Hydrogen Energy, 2013, 38(1):8597-86030.070.140.210.28[4] Fukumoto操作壓力MPfor reducing chemical mechanisms using look-up table ofchemical equilibrium calculations Application to圖6燃燒室NO排放的模擬結(jié)果與實驗數(shù)據(jù)對比CO/H/air turbulent non-premixed flame[J]. ComputersFig 6 Comparison of NO, emissions between simulationsFluids,2012,66(1):98-106[5]Rebola A, Coelho P J, Costa M. Assessment of the結(jié)論performance of several turbulence and combustion modelsin the numerical simulation of a flameless combustor[J]本研究采用4種不同湍流模型,耦合相應(yīng)的燃Combustion Science and Technology, 2013, 185(1)燒及化學反應(yīng)模型,對與燃燒室模態(tài)相似的標準旋流燃燒實驗以及合成氣旋流模型燃燒室的進行了mmacharacterist中國煤化工me[J]. Journal數(shù)值分析,并以實驗作為驗證,結(jié)果表明CNMHG3,135(4)1)湍流模型對于燃燒過程的模擬結(jié)果具有重422041-0422048第6期王翰林等:合成氣燃氣輪機燃燒室CFD模擬的模型選擇及優(yōu)化1435[7 Brizuela E A, Roudsari M z Comparison of RANS/CMCturbulent non-premixed flames of methane: flow field andmodeling of Flame D with conventional and withcompositional structure[J]. Proceedings of the Combustionpresumed mapping function statistics[J]. CombustionInstitute,2002,29(1):1913-1919Theory and Modeling, 2011, 15(5): 671-690[20]姚強,李水清,王宇,燃燒學導論(第二版)M.北京[8 Khaleghi M, Hosseini S E, Wahid M A. Eission and華大學出版社,2009:126-134combustion characteristics of hydrogen in vortexYao Qiang, Li Shuiqing, Wang Yu. An introduction toflame]. Science Engineering, 2014, 66(2): 47-51[9〕劉富強,穆勇,房愛兵,等,雙級旋流環(huán)形燃燒室流場edition)[IM]. Beijing: Tsinghua University Press, 2009的數(shù)值分析與實驗研究叮.工程電機工程學報,2012,126-134(in Chinese)32(8):28-36[21]李曉宇,游小清.在可行域內(nèi)優(yōu)化H2/CO燃燒化學反Liu Fuqiang, Mu Yong, Fang Aibing, et al. Numerical應(yīng)模型(C]y中國工程熱物理學會燃燒學學術(shù)會議論文analysis and experiment research on flow field of an集:北京,2012.combustor with a dual-stage swirler[JLi Xiaoyu, You Xiaoqing. Optimize kinetic model ofProceeding of the CSsEe, 2012, 32(8):28-36(H,/CO combustion on feasible range[C]/Chinese NationalConference of Engineer Thermophysics, Beijing, 2012(in[10 Arghode v K, Gupta A K. Effect of flow field forChinese)colorless distributed combustion(CDC)for gas turbine[22] Faravelli T, Frassoldati A, Ranzi e. Kinetic modeling ofcombustion[J]. Applied Energy, 2010,87(1): 1631-1640the interactions between NO and hydrocarbons in the[ll] Khalil A E, Gupta AK. Swirling flowfield for colorlessoxidation of hydrocarbons at low temperatures[J]distributed combusion[J]. Applied Energy, 2014, 113(1)Combustion and Flame, 2003, 132(1): 188-207223 Faravelli T, Frassoldati A, Ranzi e. Kinetic modeling of208-218[1】]田曉晶,崔玉峰,房愛兵,等.預混段結(jié)構(gòu)與氫燃料旋the interactions between NO and hydrocabons at high流預混燃燒誘導渦破碎回火極限影響的數(shù)值研究凹.中temperature[J]. Combustion and Flame, 2003, 135(1)國電機工程學報,2014,34(8):1276-128497-112[24 Vondal J, Hajek J. Swirling flow prediction in modelTian Xiaojing, Cui Yufeng, Fang Aibing, et al. Numericalcombustor with axial guide vane swirler[J]. Chemicinestigation on the effects of mixing zone structure onEngineering Transactions, 2012, 29(1):1069-1074combustion induced nortex breakdown flashback limits of25]張智博,鄭泆濤,趙煜.湍流燃燒模型對某燃氣輪杋燃a swirl-premixed hydrogen flame[J]. Proceeding of the燒室性能計算的影響叮.燃氣輪機技術(shù),2013,26(4CSSEE,2014,34(8):1276-1284( in Chinese)[13 Cuoci A, Frassoldati A, Buzzi Ferraris G, FaravelliTZhang Zhibo, Zheng Hongtao, Zhao Yu. Effects ofRanzi e. The ignition, combustion andstructure ofturbulence combustion models on performance of the gascarbon monoxide/hydrogen mixturesturbine combustor]. Gas Turbine Technology, 2013dynamics and kinetic aspects of syngas combustion[J]26(4): 17-22(in Chinese)ernational Journal of Hydrogen Energy, 2007, 32(15)26Lilleberg B, Christ D, Ertevag I S, Rian K e, Kneer R3486-3500Numerical simulation with an extinction database for use[14]Pope S B. An explanation of the round jet/plane jetith the eddy dissipation concept for turbulentAnomaly[J]. AIAA Journal, 1978, 16(3): 279-81combustion[J]. Flow Turbulence Combust, 2013, 91(1)[15] Khaleghi M, Hosseini S E, Mohsin M, Wahid M A319-346Computational and experimental investigations of [27] Georgiadis N J, Yoder D A, Engblom WA. Evaluation ofturbulent asymmetric vortex flames J. Internationlmodified two-equation turbulence models for jet flowCommunications in Heat and Mass Transfer, 2011, 38(1)predictions[c/The 44th AIAA Aerospace Sciences353-362Meeting and Exhibit, Reno, Nevada, 2006[16 Hossain M, Jones J C, Malalasekera W. Modeling of abluff-body non-premixed flame using a coupled收稿日期:2014-09-15。radiation/flamelet combustion model[J]. Flow, Turbulence作者簡介and Combustion, 2001, 67(1): 217-234王翰林(1989),男,博士研究生,主要[17 Lefebvre A H, Ballal D R. Gas turbine combustion M]事燃氣輪機燃燒室方面的研究New York: CRC Press, 2010: 1-33ranghanlin@iet.cn[18] Al-Abdeli Y, Masri A. Recirculation and flowfield雷福林(1975),男,研究員,博士生導regimes of unconfined non-reacting swirling flows[王翰林Experimental Thermal and Fluid Science, 2003, 27(1)HE中國煤化工室、媒氣化、- nlaiet cnCNMHG[19]Kalt P, Al-Abdeli Y, Masri A, Barlow R. Swirling(責任編輯王慶霞)