張楊竣 逯紅梅

1.天津大學(xué)環(huán)境科學(xué)與工程學(xué)院 2.中國市政工程華北設(shè)計研究總院3.青島農(nóng)業(yè)大學(xué)建筑工程學(xué)院

天然氣部分預(yù)混燃燒一次空氣系數(shù)與燃燒器頭部溫度關(guān)系的模擬仿真與實驗研究

張楊竣1，2逯紅梅3

1.天津大學(xué)環(huán)境科學(xué)與工程學(xué)院 2.中國市政工程華北設(shè)計研究總院3.青島農(nóng)業(yè)大學(xué)建筑工程學(xué)院

針對部分預(yù)混燃燒過程一次空氣系數(shù)變化的研究，很少涉及燃燒器頭部溫度變化的影響，利用ICEM CFD對部分預(yù)混燃燒器進(jìn)行建模，通過模擬仿真計算不同天然氣組分在不同燃燒器溫度情況對引射一次空氣系數(shù)的影響，并結(jié)合實驗測試，對比模擬仿真結(jié)果和實驗測試結(jié)果，明確燃燒器頭部溫度與部分預(yù)混引射一次空氣系數(shù)之間的變化關(guān)系。模擬結(jié)果符合隨著天然氣沃泊指數(shù)的增加大氣式燃燒器引射一次空氣系數(shù)降低的變化規(guī)律，一次空氣系數(shù)實驗測試值與模擬值偏差在±0.02之內(nèi)。當(dāng)燃?xì)獍l(fā)生置換時，一次空氣系數(shù)的變化不僅與沃泊指數(shù)成反比，同時也受燃燒器頭部溫度的影響，應(yīng)在現(xiàn)有一次空氣系數(shù)計算公式基礎(chǔ)上對其進(jìn)行溫度修正。當(dāng)貧氣置換富氣時，建議對一次空氣系數(shù)變化的溫度修正系數(shù)k取值1.1；反之，當(dāng)富氣置換貧氣時，建議溫度修正系數(shù)k取值0.9；氣質(zhì)相近的兩兩置換，無需進(jìn)行溫度修正。

部分預(yù)混燃燒 一次空氣系數(shù) 頭部溫度 沃泊指數(shù) 溫度修正

部分預(yù)混燃燒器采用噴嘴引射，通過燃?xì)鈧?cè)壓力射流引射周圍空氣，達(dá)到預(yù)混效果。引射器作為大氣式燃燒器的關(guān)鍵部件之一，一直是燃燒器設(shè)計的研究重點(diǎn)，且由于缺乏系統(tǒng)的設(shè)計理論指導(dǎo)，通常依靠經(jīng)驗、半經(jīng)驗來設(shè)計［1］。近年來，隨著CFD技術(shù)的發(fā)展，國內(nèi)一些學(xué)者使用CFD模擬軟件，對引射器結(jié)構(gòu)內(nèi)部的流場進(jìn)行模擬分析，為初步設(shè)計燃燒器提供技術(shù)參數(shù)依據(jù)，從而方便燃燒器設(shè)計和結(jié)構(gòu)改進(jìn)［2－10］。馮良等［2］討論了利用Fluent軟件對大氣式燃燒器的引射器流動進(jìn)行數(shù)值模擬研究的可行性，利用三維穩(wěn)態(tài)流動數(shù)值模擬，獲得了大型的流體狀態(tài)參數(shù)，從而對燃燒器的工作狀態(tài)有一個全面的了解，數(shù)值模擬研究對引射器的設(shè)計有一定的指導(dǎo)意義。但由于模型建立結(jié)構(gòu)簡化，導(dǎo)致模擬結(jié)果在一定程度上存在偏差。方媛媛等［3］同樣利用Fluent軟件，對低壓大氣式燃燒器引射器內(nèi)部的流場、溫度場與濃度場進(jìn)行模擬分析。相比馮良等人的研究，方媛媛等人建立引射器模型更合理，網(wǎng)格劃分和邊界條件的設(shè)置較理想，所得模擬結(jié)果更可靠。且在此基礎(chǔ)上，方媛媛等［4］就噴嘴位置對引射器性能影響進(jìn)行了數(shù)值模擬，模擬結(jié)果顯示，噴嘴位置偏移對引射器一次空氣系數(shù)的影響極大。陳偉雄等［5］通過數(shù)值模擬研究，在滿足最大的引射比的情況下優(yōu)化了天然氣引射器的結(jié)構(gòu)參數(shù)，得出了引射器中混合段和擴(kuò)壓段的最佳張角、噴嘴最佳直徑等參數(shù)。此外，倪娟娟等［6－8］也對引射器的流場進(jìn)行了數(shù)值模擬研究。郭甲生等［9］利用Fluent軟件對上進(jìn)風(fēng)燃?xì)庠畹囊湫阅苓M(jìn)行了模擬研究，通過對比實驗結(jié)果，驗證數(shù)值模擬來分析燃?xì)庠顑?nèi)部氣流流動問題的可行性。鄧海燕等［10］通過對低壓引射大氣式燃燒器內(nèi)部流場進(jìn)行實驗和數(shù)值模擬，深入分析燃燒器各部位對引射能力、流場分布、燃燒工況的影響，為燃燒器的設(shè)計和優(yōu)化提供理論和實驗依據(jù)。

國內(nèi)針對部分預(yù)混一次空氣系數(shù)的模擬，更多的是利用CFD軟件模擬仿真各結(jié)構(gòu)尺寸對引射一次空氣系數(shù)的影響，而很少研究涉及燃?xì)饨M分變化、燃燒過程變化等因素導(dǎo)致燃燒器頭部溫度變化對引射一次空氣系數(shù)的作用。本文利用ICEM CFD對部分預(yù)混燃燒器進(jìn)行建模，通過模擬仿真計算不同天然氣組分在不同燃燒器溫度情況下對引射一次空氣系數(shù)的影響，并結(jié)合實驗測試，對比模擬仿真結(jié)果和實驗測試結(jié)果，明確燃燒器頭部溫度與部分預(yù)混引射一次空氣系數(shù)之間的變化關(guān)系，為研究多氣源天然氣發(fā)生置換時燃燒器燃燒工況和火焰形態(tài)的變化做基礎(chǔ)，并為燃?xì)饣Q性［11］預(yù)測方法研究提供理論計算依據(jù)。

1 模型建立

本文設(shè)計加工了一個直立型的部分預(yù)混燃燒器用于實驗測試，并根據(jù)實際設(shè)計尺寸進(jìn)行模型建立仿真計算。根據(jù)目前國內(nèi)市場主流大氣式燃?xì)庠罹叩墓に囁?，噴嘴和火蓋選用銅材質(zhì)、分氣盤為鋁材質(zhì)、引射器為鑄鐵材質(zhì)。燃燒器各部分建模結(jié)構(gòu)尺寸及ICEM CFD非結(jié)構(gòu)化網(wǎng)格劃分情況，如圖1所示，考慮模型本身為軸對稱圓柱體結(jié)構(gòu)，故建立1／4軸對稱模型幾何結(jié)構(gòu)。

本文主要討論不同燃?xì)饨M分下燃燒器頭部溫度對引射一次空氣系數(shù)的變化影響。模擬進(jìn)行一定的簡化，只分析燃?xì)馀c空氣的物理混合過程，而不考慮預(yù)混燃?xì)馊紵幕瘜W(xué)反應(yīng)湍流不可壓縮流動問題。采用標(biāo)準(zhǔn)的k－ε計算模型，壓力和速度耦合采用SIMPLE算法，對流場采用二階迎風(fēng)格式，并采用非耦合穩(wěn)態(tài)隱式格式求解，定義收斂條件為殘差值小于10－5。流體定義為不可壓縮的牛頓流體，流動遵循動量守恒定律、能量守恒定律和質(zhì)量守恒定律。且在建立數(shù)學(xué)模型時，做如下假設(shè)：①流場已充分發(fā)展為穩(wěn)態(tài)湍流流動；②忽略體積力的影響；③流動在壁面上無滑移［2］。模擬用各氣體組分及物理特性參數(shù)，如表1所示。

表1 模擬和實驗測試用氣體組分及特性參數(shù)（15℃，15℃，101.325 kPa）Table 1 Gas components and properties for simulation and experiment（15℃，15℃，101.325 kPa）

根據(jù)設(shè)計燃燒器的實驗測試結(jié)果，對模型各邊界條件進(jìn)行如下設(shè)置：①燃?xì)膺M(jìn)口設(shè)為壓力進(jìn)口，2 000 Pa，溫度為室溫298K，紊流強(qiáng)度設(shè)為4%，水力直徑設(shè)為2mm；②空氣進(jìn)口設(shè)為壓力進(jìn)口，0Pa，溫度為室溫298K，紊流強(qiáng)度設(shè)為4%，水力直徑設(shè)為16mm；③燃燒器出口設(shè)為壓力出口5Pa，溫度為室溫298K，紊流強(qiáng)度設(shè)為4%，水力直徑設(shè)為2.5mm；④各材質(zhì)表面壁面厚度2.5mm，模擬定義3種頭部溫度情況，各結(jié)構(gòu)常壁溫表面參數(shù)設(shè)置列于表2。

2 模擬仿真結(jié)果討論

定義模型頭部溫度為正常穩(wěn)定燃燒情況，模擬得到5種不同天然氣組分特性情況下的燃燒器引射參數(shù)，列于表3。表1中列出各氣源沃泊指數(shù)特性WPNG－R＞W(wǎng)LNG＞W(wǎng)CH4＞W(wǎng)PNG＞W(wǎng)PNG－P，表3模擬結(jié)果顯示，α＇LNG＜α＇PNG－R＜α＇CH4＜α＇PNG＜α＇PNG－P。模擬結(jié)果符合隨著天然氣沃泊指數(shù)的增加，大氣式燃燒器引射一次空氣系數(shù)降低，引射一次空氣系數(shù)與氣源沃泊指數(shù)呈反比關(guān)系，即：

表2 3種頭部溫度情況下各常壁溫表面參數(shù)設(shè)置Table 2 Temperature settings of constant temperature walls in three situations

表3 不同天然氣情況下燃燒器頭部溫度為正常穩(wěn)定燃燒時引射參數(shù)模擬結(jié)果Table 3 Inject simulation results of normal combustion temperature of gas burner head under different natural gases

在燃燒器頭部溫度為正常穩(wěn)定燃燒情況下的模擬基礎(chǔ)之上，本文對CH4、PNG和LNG 3種天然氣組分在頭部溫度分別為高溫和低溫的情況進(jìn)行一次空氣引射模擬，得到不同頭部溫度情況下的燃燒器引射參數(shù)，列于表4。

表4 不同頭部溫度下各燃?xì)饨M分燃燒器引射參數(shù)模擬結(jié)果Table 4 Inject parameters simulation results of different gas burner head temperatures

根據(jù)圖2模擬數(shù)據(jù)，在保證引射氣體特性不變的情況下，隨著頭部溫度的升高，3種天然氣下的燃燒器引射一次空氣系數(shù)均降低，且3條擬合線性斜率彼此接近。對于大氣式燃燒器，引射一次空氣系數(shù)的變化與頭部溫度呈反比關(guān)系，即：

3 模擬結(jié)果與實驗結(jié)果對比

實驗對CH4、PNG、PNG－R、PNG－P 4種氣源進(jìn)行一次空氣系數(shù)變化測試。實驗時，燃?xì)膺M(jìn)口壓力保持在2kPa不變，設(shè)計的部分預(yù)混燃燒器空氣閥開度保持不變。用注射器在燃燒器頭部火孔內(nèi)抽取混合氣樣本，在燃燒器噴嘴進(jìn)口前100mm管道處采集燃?xì)鈽颖荆缓髮⒓內(nèi)細(xì)鈽颖竞突旌蠚鈽颖痉謩e進(jìn)行氣相色譜分析。使用注射器抽氣采樣時，先沖洗注射器5～8次，取樣動作保持勻速緩慢進(jìn)行，保證采得的氣體樣本為實際氣體組分。通過氣相色譜分析法分析出純?nèi)細(xì)夂鸵淮慰諝猓細(xì)饣旌蠚獾慕M分，利用燃?xì)夂涂諝饣旌线^程中某組分i的平衡式（3），即可計算出一次空氣系數(shù)，見公式（4）。

式中：ri，g為燃?xì)庵衖組分的體積分?jǐn)?shù)；ri，a為空氣中i組分的體積分?jǐn)?shù)；ri，m為混合氣中i組分的體積分?jǐn)?shù)；V0為15℃、101.235kPa下1m3的燃?xì)馊紵枰睦碚摽諝饬浚琺3／m3；α′為一次空氣系數(shù)。

實驗以火焰根部（即火孔出口處）溫度為對比參照溫度，分別測試3種頭部溫度下一次空氣系數(shù)變化情況，測試結(jié)果如表5所示。對于測試用的4種天然氣，實驗與模擬存在的差別在于兩者的頭部溫度不同。實驗過程中無法進(jìn)行溫度的控制，燃燒器所能達(dá)到的最高溫度與燃?xì)馓匦杂嘘P(guān)；而模擬過程中對燃燒器火孔表面初始條件溫度進(jìn)行簡化設(shè)定，統(tǒng)一為1 000K、750K和550K。但實驗測試結(jié)果顯示，燃燒器從點(diǎn)燃到穩(wěn)定燃燒，頭部溫度逐漸升高，引射一次空氣系數(shù)逐漸降低，表現(xiàn)出與模擬結(jié)果一樣的變化趨勢，即：α＇∝1／T。

表5 4種天然氣在不同頭部溫度下的一次空氣系數(shù)實驗測試結(jié)果Table 5 Experiment results of primary air ratio during 4 kinds of natural gases at different gas burner head temperatures

設(shè)計燃燒器點(diǎn)燃30min后，基本達(dá)到穩(wěn)定燃燒狀態(tài)，頭部溫度達(dá)到最高值。故將30min后的一次空氣系數(shù)測試結(jié)果認(rèn)定為穩(wěn)定燃燒情況下的引射一次空氣系數(shù)值。將該值與表3中的模擬值進(jìn)行對比，如圖3所示。由于模型建立過程中對實際引射器入口進(jìn)行了簡化，將風(fēng)門結(jié)構(gòu)簡化成整個圓形面的空氣入口，所以模擬結(jié)果相比于實際情況表現(xiàn)為引射能力較差、引射一次空氣系數(shù)較小。4種天然氣的一次空氣系數(shù)實驗測試值與模擬值偏差在±0.02之內(nèi)，兩者可比，模擬結(jié)果可信。

當(dāng)天然氣發(fā)生置換時，部分預(yù)混燃燒器的一次空氣系數(shù)變化不僅與氣源沃泊指數(shù)成反比關(guān)系，同時也受到燃燒器頭部溫度的影響（兩者成反比關(guān)系），即

而目前，燃?xì)獍l(fā)生置換時，采用的一次空氣系數(shù)變化計算公式為［1］：

因此，對天然氣置換后的一次空氣系數(shù)變化計算公式應(yīng)進(jìn)行溫度修正。

表5為設(shè)計的直立型部分預(yù)混燃燒器一次空氣系數(shù)隨天然氣沃泊指數(shù)和燃燒頭部溫度變化的實驗測試數(shù)據(jù)。由于燃燒器在PNG－P情況下出現(xiàn)離焰，即測得的頭部燃燒溫度并不能真實地反應(yīng)PNG－P燃燒穩(wěn)定狀態(tài)下的燃燒器頭部溫度，其測試結(jié)果不可信，因此，在進(jìn)行燃?xì)庵脫Q計算中不采用該組實驗數(shù)據(jù)。對剩余3種天然氣進(jìn)行兩兩置換后α′s／α′a、Wa／Ws、V0a／V0s、Hhs／Hha、Ta／Ts的計算，計算結(jié)果列于表6。

根據(jù)前文分析，燃?xì)庵脫Q后一次空氣系數(shù)不僅與燃?xì)馕植粗笖?shù)成反比，同時與燃燒器頭部溫度也成反比，因此，在原燃?xì)庵脫Q后一次空氣系數(shù)變化計算公式（6）的基礎(chǔ)上，加入溫度影響因素，即一次空氣系數(shù)變化計算公式改寫成公式（7）。將各兩兩置換氣按公式（7）分別計算公式兩側(cè)數(shù)值，并進(jìn)行比值統(tǒng)計分析，列于表7。進(jìn)過3次獨(dú)立實驗測試，3種天然氣之間兩兩置換后，測得數(shù)據(jù)按公式（7）計算得到的結(jié)果線性平均值均為1，即燃?xì)庵脫Q后一次空氣系數(shù)變化應(yīng)按公式（7）進(jìn)行計算。

表6 4種天然氣兩兩置換后各參數(shù)比值Table 6 Parameters ratio when natural gases interchanged with each other

表7 公式（7）兩側(cè)實驗數(shù)據(jù)計算值比值Table 7 Specific value of two sides of equation（7）calculated by experiment results

由于進(jìn)行燃燒器頭部溫度測試，將增加實際工作量，實際操作不易行。在實際應(yīng)用中，評價燃?xì)庵脫Q后一次空氣系數(shù)變化時，一般直接根據(jù)氣質(zhì)特性進(jìn)行計算。因此，為了便于實際應(yīng)用，在公式（7）的基礎(chǔ)之上，將燃?xì)庵脫Q后一次空氣系數(shù)變化公式改寫成公式（8），其中k為溫度修正系數(shù)，即將燃燒器頭部溫度對一次空氣系數(shù)的影響簡化為溫度修正。根據(jù)3種天然氣兩兩置換實驗測試數(shù)據(jù)，分成較貧燃?xì)庵脫Q較富燃?xì)狻⑤^富燃?xì)庵脫Q較貧燃?xì)夂蜌赓|(zhì)相近燃?xì)庀嗷ブ脫Q3種情況，按公式（8）計算3種情況下各自的溫度修正系數(shù)k，所得結(jié)果列于表8。當(dāng)貧氣置換富氣時，燃燒器頭部溫度將從原先的高溫變成低溫，而一次空氣系數(shù)與頭部溫度成反比，根據(jù)實驗結(jié)果計算，建議此類情況對一次空氣系數(shù)變化的溫度修正系數(shù)k取值1.1；反之，當(dāng)富氣置換貧氣時，建議溫度修正系數(shù)k取值0.9；氣質(zhì)相近的兩兩置換，無需進(jìn)行溫度修正。

表8 按公式（8）進(jìn)行溫度修正實驗測得的溫度修正系數(shù)kTable 8 Experiment results of temperature correction coefficient value used in equation（8）

4 結(jié)論

通過模擬仿真計算不同天然氣組分在不同燃燒器溫度情況下對引射一次空氣系數(shù)的影響，并進(jìn)行實驗測試研究，得出以下結(jié)論：

（1）模擬結(jié)果符合隨著天然氣沃泊指數(shù)的增加大氣式燃燒器引射一次空氣系數(shù)降低的變化規(guī)律，一次空氣系數(shù)實驗測試值與模擬值偏差在±0.02之內(nèi)，兩者可比，證明可通過CFD對部分預(yù)混燃燒在不同燃燒器頭部溫度下不同氣質(zhì)特性的引射情況進(jìn)行模擬仿真，結(jié)果可信，類似的實驗測試可省略，節(jié)約實驗費(fèi)用。

（2）在保證引射氣體特性不變的情況下，隨著頭部溫度的升高，燃燒器引射一次空氣系數(shù)均降低，且模擬得到的3種天然氣下一次空氣系數(shù)變化擬合線性斜率彼此接近；實驗測試結(jié)果顯示，燃燒器從點(diǎn)燃到穩(wěn)定燃燒，頭部溫度逐漸升高，引射一次空氣系數(shù)逐漸降低，表現(xiàn)出與模擬結(jié)果一樣的變化趨勢。

（3）當(dāng)天然氣發(fā)生置換時，部分預(yù)混燃燒器的一次空氣系數(shù)變化不僅與氣源沃泊指數(shù)成反比關(guān)系，同時也受到燃燒器頭部溫度的影響（兩者成反比關(guān)系），故一次空氣系數(shù)變化計算采用公式（7）或公式（8），其中，貧氣置換富氣k建議取值1.1、富氣置換貧氣k建議取值0.9、氣質(zhì)相近的兩兩置換k＝1。

［1］同濟(jì)大學(xué).燃?xì)馊紵c應(yīng)用［M］.上海：中國建筑工業(yè)出版社，2011.

［2］馮良，劉鯤，韓國園.大氣式燃?xì)馊紵饕淦鞯腃FD研究［J］.上海煤氣，2003（2）：13－16.

［3］方媛媛，郭全，傅忠誠.低壓大氣式燃燒器中引射器的數(shù)值模擬［J］.北京建筑工程學(xué)院學(xué)報，2006，22（4）：55－58.

［4］方媛媛，郭全.噴嘴位置對引射器性能影響的數(shù)值模擬［J］.煤氣與熱力，2007，27（7）：42－44.

［5］CHEN W X，CHONG D T，YAN J J，et al.Numerical optimization on the geometrical factors of natural gas ejectors［J］.International Journal of Thermal Sciences，2011，50（4）：1554－1561.

［6］倪娟娟，張海軍.大氣式燃燒器內(nèi)引射器的數(shù)值模擬與實驗研究［J］.上海煤氣，2012（2）：22－25.

［7］游超林，陳迪龍，蔡國漢，等.多射流引射器的數(shù)值模擬［J］.煤氣與熱力，2010，30（6）：B12－B15.

［8］林興偉.大氣式燃燒器適應(yīng)性實驗及模擬研究［D］.哈爾濱：哈爾濱工業(yè)大學(xué)，2011：40－57.

［9］郭甲生，秦朝葵，孫秀麗.上進(jìn)風(fēng)燃?xì)庠钜湫阅艿臄?shù)值模擬［J］.工業(yè)加熱，2008，37（2）：12－15.

［10］鄧海燕，仇明貴.低壓引射大氣式燃燒器內(nèi)部流場分析［C］／／中國土木工程學(xué)會燃?xì)夥謺?昆明，2012：113－119.

［11］劉鵬君，秦朝葵，楊賢潮.工業(yè)燃?xì)馊紵骰Q性實驗測試方法研究［J］.石油與天然氣化工，2013，42（5）：487－491.

Simulation and experiment research on the relationship between primary air ratio and burner head temperature of partially premixed combustion burning natural gas

Zhang Yangjun1，2，Lu Hongmei3
（1.School of Environmental Science and Engineering，Tianjin University，Tianjin300072，China）
（2.North China Municipal Engineering Design &Research Institute，Tianjin300074，China）
（3.School of Architectural Engineering，Qingdao Agricultural University，Qingdao 266109，China）

There are rarely studies on the influence of gas burner head temperature caused by gas components change or combustion process variation with primary air ratio.In this paper，it simulated the relationship between primary air ratio and gas burner head temperature of partially premixed combustion by using ICEM CFD，and compared with the experiment results，it confirmed the temperature correction coefficient to calculate the specify value of primary air ratio when gas was substituted.The simulation results conform to the characteristic of primary air ratio decreasing with Wobbe Index increasing，the deviation of simulation value and experiment value is within±0.02.When gas was interchanged，the primary air ratio of partially premixed combustion is not only inversed with Wobbe Index，but also affected by gas burner head temperature.From the experiment results，when rich gas is substituted by poor gas，the suggestion temperature correction coefficient k is 1.1；and when poor gas is substituted by rich gas，the suggestion temperature correction coefficient k is 0.9；but when the properties of interchange gases are similar，it is not necessary to make temperature correction.

partially premixed combustion，primary air ratio，burner head temperature，Wobbe index，temperature correction

TE624.9

A

10.3969／j.issn.1007－3426.2015.05.010

張楊竣（1986－），男，浙江浦江，2014年7月畢業(yè)于同濟(jì)大學(xué)供熱、供燃?xì)狻⑼L(fēng)及空調(diào)工程專業(yè)，博士學(xué)歷；天津大學(xué)環(huán)境科學(xué)與工程在站企業(yè)博士后；主要從事燃?xì)馊紵c應(yīng)用、燃?xì)饣Q性研究工作。E－mail：zyjtongji＠163.com

2015－05－18；編輯：康 莉