郭培卿+葛冰

摘要： 針對合成氣非預(yù)混火焰結(jié)構(gòu)開展數(shù)值模擬和試驗驗證，分析天然氣改燒合成氣后燃燒特性的變化規(guī)律.結(jié)果表明，大渦模擬（Large Eddy Simulation，LES）在速度分布和溫度分布的預(yù)測中與試驗結(jié)果比較吻合，而在對燃燒過程產(chǎn)物（如OH自由基）預(yù)測中則與試驗結(jié)果有所差異.采用數(shù)值模擬與試驗測量相結(jié)合的方法，探討和分析合成氣燃燒特性的變化規(guī)律：與天然氣火焰相比，合成氣燃燒時高溫區(qū)域更大，火焰穩(wěn)定性較好；隨著當(dāng)量比提高，燃燒室熱負荷不斷增大，同時最高回流速度增大，火焰根部受到壓縮，逐漸呈現(xiàn)出推舉火焰特征.

關(guān)鍵詞： 合成氣； 燃燒； 火焰結(jié)構(gòu)； 大渦模擬； 平面激光誘導(dǎo)熒光

中圖分類號： TK16； TP391.9文獻標志碼： B

Abstract： The numerical simulation and test validation are used to study syngas nonpremixed flame structure， and the change rule of combustion characteristics of syngas which is instead of natural gas is analyzed. The results indicate that， the velocity and temperature distribution predicted by Large Eddy Simulation（LES） are in good consistency with the test results， but the prediction of substances generated during combustion process such as OH radical prediction is different from the test results. The method of numerical simulation combined with test measurement is used to discuss and analyze the change rule of syngas combustion. It shows that， comparing with natural gas flame， the high temperature zone is enlarged and the flame stability is better； with the increase of equivalence ratio， the heat load of combustion chamber increases， the top back flow velocity increases， and the flame shows pushup characteristic.

Key words： syngas； combustion； flame structure； large eddy simulation； planar laser induced fluorescence

0引言

整體煤氣化聯(lián)合循環(huán)（Integrated Gasification Combined Cycle， IGCC）是目前世界上公認的煤炭清潔化利用的重要發(fā)展方向.在IGCC系統(tǒng)中，燃氣輪機使用煤氣化產(chǎn)生的合成氣，其可燃成分為H2，CO和其他少量高階碳氫化合物，熱值普遍較低.為保持輸出功率一致，與天然氣相比，合成氣燃料流量將大幅提高，因此如何降低NOx排放是合成氣燃氣輪機燃燒室設(shè)計面臨的主要難點之一.作為低污染燃燒技術(shù)之一的合成氣燃燒正逐步成為相關(guān)燃燒領(lǐng)域的研究熱點.[14]隨著計算機計算能力的不斷提高，CFD在燃燒現(xiàn)象研究領(lǐng)域中得到越來越多的應(yīng)用，其中非預(yù)混燃燒物理過程受湍流摻混、分子擴散和化學(xué)反應(yīng)的綜合影響，受試驗數(shù)據(jù)庫容量、化學(xué)反應(yīng)動力學(xué)和燃燒模型等多種限制，通過CFD計算精確求解湍流燃燒過程中間產(chǎn)物的時空分布依然十分困難，但在預(yù)測湍流燃燒流場中的主要物理量作為試驗研究的補充和延伸方面，CFD計算發(fā)揮極其重要的作用.由于其成本較低并且具備全流場模擬和特殊工況模擬能力等優(yōu)勢，CFD是開展燃燒現(xiàn)象研究的重要手段之一.本文以合成氣非預(yù)混火焰為研究對象開展大渦模擬（Large Eddy Simulation， LES），并利用PIV和OHPLIF等測量數(shù)據(jù)進行試驗驗證，根據(jù)數(shù)值模擬與試驗結(jié)果，針對合成氣與天然氣燃燒特性開展研究.

1數(shù)值模擬方法

LES將湍流瞬時運動通過濾波方法分解為大、小尺度渦2部分，其中大尺度渦在計算中直接求解，然后利用亞網(wǎng)格模型建立與大尺度渦的關(guān)系對其小尺度渦進行模擬.主要的亞網(wǎng)格模型包括Smagorinsky亞網(wǎng)格模型[5]和Germano動態(tài)亞網(wǎng)格模型[6]等.與直接數(shù)值模擬相比，LES計算量大大減少[7]，同時可以精確求解，獲得的大尺度渦湍流運動能夠真實反映湍流運動的主要瞬態(tài)特征[8]；另一方面，與傳統(tǒng)雷諾時均模擬方法相比計算量大得多仍是LES應(yīng)用于實際工程問題的最大障礙[9].

在模擬非預(yù)混燃燒時，用概率密度函數(shù)（βPDF分布）描述混合分數(shù)的脈動性質(zhì)，根據(jù)不同燃料生成對應(yīng)的PDF分布，應(yīng)用在相應(yīng)的燃燒過程計算中，燃料組分見表1.

5結(jié)論

利用LES針對合成氣非預(yù)混燃燒開展研究，并根據(jù)試驗驗證結(jié)果，對合成氣燃燒規(guī)律進行對比分析，得出以下結(jié)論：

（1）LES對燃燒室內(nèi)流場和溫度場分布的預(yù)測與試驗結(jié)果較為吻合，Germano亞網(wǎng)格模型對湍流脈動的預(yù)測更為精確；

（2）采用簡單化學(xué)反應(yīng)假設(shè)的PDF燃燒模型對燃燒中間產(chǎn)物如OH自由基分布的預(yù)測與試驗測量結(jié)果存在偏差；

（3）隨著當(dāng)量比的提高，合成氣火焰根部逐漸出現(xiàn)M型分布，形成典型的推舉火焰特征，高溫區(qū)向下游移動，火焰根部受到壓縮，旋流燃料與空氣間混合得到加強，最高回流速度增大.

參考文獻：

[1]STRAKEY P， SIDWELL T， ONTKO J. Investigation of the effects of hydrogen addition on lean extinction in a swirl stabilized combustor[J]. Proc Combustion Inst， 2007， 31（2）： 31733180.

[2]CASARSA L， MICHELI D， PEDIRODA V， et al. Investigations of pyrolysis syngas swirl flames in a combustor model[C]// Proc ASME Turbo Expo 2009， GT200959610. Orlando， 2009.

[3]張文興， 穆克進， 王岳， 等. 合成氣甲醇摻燒火焰研究[J]. 熱能動力工程， 2009， 24（2）： 236241.

[4]趙曉燕， 李祥晟， 豐鎮(zhèn)平. 燃氣輪機低熱值合成氣燃燒室內(nèi)三維湍流流動的數(shù)值模擬研究[J]. 動力工程， 2009， 29（4）： 330334.

[5]SMAGORINSKY J S. General circulation experiments with the primitive equations[J]. Mon Weather Rev， 1963， 91（3）： 99164.

[6]GERMANO M， PIOMELLI U， MOIN P， et al. A dynamic subgridscale eddy viscosity model[J]. Phys Fluids A， 1991， 3（7）： 17601765.

[7]LILLY D K. A proposed modification of the Germano subgridscale closure method[J]. Phys Fluids： A， 1992， 4（3）： 633635.

[8]趙堅行. 燃燒的數(shù)值模擬[M]. 北京： 科學(xué)出版社， 2002.

[9]周力行， 胡礫元， 王方. 湍流燃燒大渦模擬的最近研究進展[J]. 工程熱物理學(xué)報， 2006， 27（2）： 331334.

[10]寧晃， 高歌. 燃燒室氣動力學(xué)[M]. 2版. 北京： 科學(xué)出版社， 1987.

[11]WISSINK J G. DNS of separating， low Reynolds number flow in a turbine cascade with incoming wakes[J]. Int J Heat & Fluid Flow， 2003， 24（4）： 626635.

（3）隨著當(dāng)量比的提高，合成氣火焰根部逐漸出現(xiàn)M型分布，形成典型的推舉火焰特征，高溫區(qū)向下游移動，火焰根部受到壓縮，旋流燃料與空氣間混合得到加強，最高回流速度增大.

參考文獻：

[1]STRAKEY P， SIDWELL T， ONTKO J. Investigation of the effects of hydrogen addition on lean extinction in a swirl stabilized combustor[J]. Proc Combustion Inst， 2007， 31（2）： 31733180.

[2]CASARSA L， MICHELI D， PEDIRODA V， et al. Investigations of pyrolysis syngas swirl flames in a combustor model[C]// Proc ASME Turbo Expo 2009， GT200959610. Orlando， 2009.

[3]張文興， 穆克進， 王岳， 等. 合成氣甲醇摻燒火焰研究[J]. 熱能動力工程， 2009， 24（2）： 236241.

[4]趙曉燕， 李祥晟， 豐鎮(zhèn)平. 燃氣輪機低熱值合成氣燃燒室內(nèi)三維湍流流動的數(shù)值模擬研究[J]. 動力工程， 2009， 29（4）： 330334.

[5]SMAGORINSKY J S. General circulation experiments with the primitive equations[J]. Mon Weather Rev， 1963， 91（3）： 99164.

[6]GERMANO M， PIOMELLI U， MOIN P， et al. A dynamic subgridscale eddy viscosity model[J]. Phys Fluids A， 1991， 3（7）： 17601765.

[7]LILLY D K. A proposed modification of the Germano subgridscale closure method[J]. Phys Fluids： A， 1992， 4（3）： 633635.

[8]趙堅行. 燃燒的數(shù)值模擬[M]. 北京： 科學(xué)出版社， 2002.

[9]周力行， 胡礫元， 王方. 湍流燃燒大渦模擬的最近研究進展[J]. 工程熱物理學(xué)報， 2006， 27（2）： 331334.

[10]寧晃， 高歌. 燃燒室氣動力學(xué)[M]. 2版. 北京： 科學(xué)出版社， 1987.

[11]WISSINK J G. DNS of separating， low Reynolds number flow in a turbine cascade with incoming wakes[J]. Int J Heat & Fluid Flow， 2003， 24（4）： 626635.

（3）隨著當(dāng)量比的提高，合成氣火焰根部逐漸出現(xiàn)M型分布，形成典型的推舉火焰特征，高溫區(qū)向下游移動，火焰根部受到壓縮，旋流燃料與空氣間混合得到加強，最高回流速度增大.

參考文獻：

[1]STRAKEY P， SIDWELL T， ONTKO J. Investigation of the effects of hydrogen addition on lean extinction in a swirl stabilized combustor[J]. Proc Combustion Inst， 2007， 31（2）： 31733180.

[2]CASARSA L， MICHELI D， PEDIRODA V， et al. Investigations of pyrolysis syngas swirl flames in a combustor model[C]// Proc ASME Turbo Expo 2009， GT200959610. Orlando， 2009.

[3]張文興， 穆克進， 王岳， 等. 合成氣甲醇摻燒火焰研究[J]. 熱能動力工程， 2009， 24（2）： 236241.

[4]趙曉燕， 李祥晟， 豐鎮(zhèn)平. 燃氣輪機低熱值合成氣燃燒室內(nèi)三維湍流流動的數(shù)值模擬研究[J]. 動力工程， 2009， 29（4）： 330334.

[5]SMAGORINSKY J S. General circulation experiments with the primitive equations[J]. Mon Weather Rev， 1963， 91（3）： 99164.

[6]GERMANO M， PIOMELLI U， MOIN P， et al. A dynamic subgridscale eddy viscosity model[J]. Phys Fluids A， 1991， 3（7）： 17601765.

[7]LILLY D K. A proposed modification of the Germano subgridscale closure method[J]. Phys Fluids： A， 1992， 4（3）： 633635.

[8]趙堅行. 燃燒的數(shù)值模擬[M]. 北京： 科學(xué)出版社， 2002.

[9]周力行， 胡礫元， 王方. 湍流燃燒大渦模擬的最近研究進展[J]. 工程熱物理學(xué)報， 2006， 27（2）： 331334.

[10]寧晃， 高歌. 燃燒室氣動力學(xué)[M]. 2版. 北京： 科學(xué)出版社， 1987.

[11]WISSINK J G. DNS of separating， low Reynolds number flow in a turbine cascade with incoming wakes[J]. Int J Heat & Fluid Flow， 2003， 24（4）： 626635.