李玉蘭

摘要：根據微尺度燃燒室內氫氣和空氣的預混燃燒過程，利用STAR-CD計算軟件建立微燃燒的物理和數學模型，在試驗驗證的基礎上模擬壁面導熱系數、壁厚和外壁面?zhèn)鳠嵯禂档缺诿鎱祵錃夂涂諝忸A混合燃燒的影響。結果表明：靠近燃燒室入口處的火焰中心溫度隨導熱系數減小而增加，出口處混合氣體的溫度隨導熱系數減小而降低；隨著外壁面厚度增加，外壁面溫度降低，壁面溫度分布更均勻；隨著外壁面?zhèn)鳠嵯禂翟黾?，燃燒室內火焰溫度降低，出口處的氣體溫度和外壁面溫度相應降低。

關鍵詞：微燃燒室；壁面參數；數值模擬；燃燒特性；氫氣

中圖分類號：TK42 文獻標識碼：A 文章編號：1674-1161（2017）09-0044-04

微燃燒器是微熱光電系統(tǒng)的核心部件，其結構和性能對外壁面溫度分布的均勻性有重要影響，進而影響光電能量的轉換效率。近年來，國內外學者對微燃燒器的結構和性能進行了大量研究。受加工尺寸的限制，研究壁面參數對微燃燒的影響比較困難。為此，利用STAR-CD計算軟件建立微燃燒的物理和數學模型，在試驗驗證的基礎上模擬研究燃燒室壁面導熱系數、壁厚和外壁面?zhèn)鳠嵯禂档缺诿鎱祵ξ⒊叨热紵覂葰錃夂涂諝忸A混合燃燒的影響，旨在為微燃燒室結構設計和性能改善提供一定的參考。

1 數值模型和方法

1.1 物理數學模型

研究的模擬對象（如圖1所示），為內徑0.2 mm、壁厚0.05 mm、長1.6 mm的微型圓柱管。采用STAR-CD前處理軟件pro-star對微燃燒室進行網格劃分。燃燒室結構較為規(guī)則，因此采用結構化網格形式，在簡單計算區(qū)域內可以直接生成。利用網格特點可以構造各種算法。

1.2 數學模型

數學模型包括以下基本方程：質量守恒方程、動量守恒方程、化學組分守恒方程、流體能量守恒方程、壁面能量守恒方程。

1.3 流體流動模型

該燃燒室尺度較小，但研究的傳熱及流動問題仍屬于連續(xù)介質區(qū)，適用Navier-Stokes方程及Fourier熱傳導定律，可運用商業(yè)通用軟件STAR-CD進行數值模擬研究。在微圓管內通入氫氣/空氣預混氣體燃燒，以空氣物性參數作為整個混合氣流的物性參數，認為是層流粘性流，可用理想氣體定律計算流體密度。用物質性質的質量分數的加權平均數來計算流體混合物的比熱、粘度、熱傳導率。

1.4 燃燒模型

燃燒反應模型很多?？紤]到氫氣和氧氣的化學反應，選擇適于層流預混合燃燒的Chemical Kinetics模型，其化學反應率遵循Arrhenius定律。

1.5 傳熱模型

在微尺度燃燒室內燃燒過程中，由于表面積體積比增大，因此壁面熱損失不能忽略，必須綜合考慮導熱、對流和輻射散熱產生的熱損失。

1.5.1 導熱 壁面導熱損失采用Conduct and conjugate Heat Transfer模型計算。模擬計算時采用雙精度，這樣利于解決潛在的收斂問題。

1.5.2 對流換熱 在計算中需要確定對流換熱產生熱損失的傳熱系數，計算公式為：

1.5.3 輻射換熱 主要考慮燃燒室外壁面的輻射散熱對微尺度燃燒室內燃燒的影響。在STAR-CD模擬中，如果計算區(qū)域中沒有建立空氣網格，燃燒室外壁面輻射散熱可不設置輻射模型，但需要在固體外壁面邊界條件中給出輻射散熱的熱阻。輻射傳熱系數hra的計算公式為：

1.6 邊界條件和算法

在邊界條件中，入口條件選取速度邊界條件，需要給出預混合氣在入口的流速、溫度和氣體混合物質量百分比。燃燒室的結構屬于軸對稱，為全面分析微燃燒室的火焰?zhèn)鞑ヌ匦?，采用三維數值模擬。數值模擬采用1/2圓柱作為計算區(qū)域，計算區(qū)域的底部邊界為軸對稱邊界條件。出口壓力邊界條件固定為

1.013 25×105 Pa。燃燒室壁面采用無滑移、無物質流出表面壁面邊界條件。

燃燒室壁面的熱邊界條件與2個主要傳熱因素有關：壁面散失到外界環(huán)境的熱損失和墻壁內的熱傳導。壁面包括內壁面邊界和外壁面邊界。內壁面邊界熱條件選擇為傳導，不需要給出具體的導熱熱阻，STAR-CD會根據材料特性和壁面厚度自行計算。外壁面邊界條件只需給出綜合的傳熱熱阻就可以確定散熱量。

采用只用于穩(wěn)態(tài)計算的SIMPLE算法求解上述控制方程。在計算過程中，收斂條件規(guī)定各個控制方程迭代計算的余量至少小于10-3。速度的差分格式采用重構格式（MARS，Monotone Advection and Reconstruction Scheme），這樣不僅可提高梯度劇烈變化時的計算精度，還能減少對網格的依賴性。

2 模型驗證

所模擬的研究對象尺寸較小，已經超出了目前國內能加工的最小尺寸。為驗證計算模型的準確性，與黃俊等開展的試驗進行對比，選用具有耐高溫性能的SiC陶瓷薄壁柱型微燃燒室。通過試驗發(fā)現，氫氧混合氣在薄壁柱型微燃燒室內易實現穩(wěn)定燃燒。圖2是不同氫氧混合比下的壁面溫度分布，圖3是模擬燃燒室外壁面的溫度分布。燃燒室外壁面溫度變化規(guī)律基本吻合，驗證了模擬結果的正確性。

3 結果及分析

3.1 微燃燒室壁面材料對微尺度燃燒的影響

燃燒室內徑0.2 mm，壁面厚度0.05 mm，外壁面散熱系數100 W/（m2·K），以化學當量比下氫氣和空氣預混合方式通入燃燒室，入口流速固定為1.5 m/s，入口溫度為300 K，選擇3種燃燒室壁面材料進行試驗，用以研究壁面熱傳導對微尺度燃燒室中燃燒特性的影響。第一種是高導熱率材料（如金屬和硅），導熱率為200 W/（m2·K）。第二種是中等導熱率材料（如陶瓷），導熱率為15 W/（m2·K）。第三種是低導熱率材料（如絕熱材料），導熱率在1 W/（m2·K）級別。不同壁面材料燃燒室內火焰溫度沿無量綱軸向距離分布如圖4所示。endprint

由圖4可知：隨著壁面導熱系數減小，火焰中心最大溫度增加，火焰中心（高溫區(qū)）稍微向入口下游移動。可見，靠近燃燒室入口處的火焰中心溫度隨導熱系數減小而增加。但在燃燒室出口處，混合氣體溫度隨導熱系數減小而降低。這是因為導熱系數減小，通過壁面的導熱量降低，壁面向環(huán)境的散熱量也減少。

另外，熱點和較高壁面溫度是壁面材料熔化或變脆的主要原因。壁面材料導熱率對燃燒室壁面的溫度分布有很大影響，對影響燃燒室特性的材料應力也有顯著影響。Bi準則數（Bi=hδ/λ）可以用來考察壁面溫度的分布特性。在研究中，綜合傳熱系數為100 W/（m2·K），壁面厚度為0.05 mm，Bi準則數隨燃燒室壁面導熱率的變化而變化。

不同壁面材料燃燒室外壁面溫度沿無量綱軸向距離的分布見圖5。

由圖5可見，低導熱率（Bi=5×10-3）的壁面溫度分布有較大梯度，對壁面材料的特性有較高要求，要求其能夠承受高溫度和溫度梯度（能夠導致高的材料應力）。然而，高導熱率（Bi=2.5×10-5）壁面的溫度分布較均勻，降低了壁面熱梯度，有利于降低燃燒室壁面的熱壓。

3.2 微燃燒室壁面厚度對微尺度燃燒的影響

燃燒室內徑0.2 mm，外壁面散熱系數100 W/（m2·K），以化學當量比下氫氣和空氣預混合方式通入燃燒室，入口流速固定為1.5 m/s，入口溫度300 K，分別選擇0.05，0.1，0.2 mm壁厚來研究壁面厚度對微尺度燃燒室燃燒特性的影響。

壁面厚度對燃燒室外壁面溫度分布有重要影響，詳見圖6。

由圖6可見，隨外壁面厚度增加，外壁面溫度降低，壁面溫度分布更均勻。薄壁面厚度和氣體溫度分布類似，在靠近火焰中心有一個150 K的峰值溫度，入口處溫度達960 K。壁面越厚，壁面溫度分布越均勻，這是因為它們有較大的面積產生軸向熱傳導。而均勻的溫度使它們有較低的壁面峰值溫度（大約為750～900 K）和較低的入口壁面溫度。較低的入口壁面溫度減弱了火焰穩(wěn)定性。為降低壁面內導熱產生的熱損失，燃燒室的壁厚越小越好?？紤]目前機電裝置加工技術的不足，控制壁面厚度有利于實現微尺度燃燒的火焰穩(wěn)定性。

3.3 外壁面?zhèn)鳠嵯禂祵ξ⒊叨热紵挠绊?/p>

微尺度燃燒室有較高的表面積體積比，壁面熱損失對燃燒穩(wěn)定性有重要影響。取燃燒室內徑0.2 mm、壁厚0.05 mm、長1.6 mm的微硅管，以化學當量比下氫氣與空氣預混合方式通入燃燒室，入口溫度為300 K，入口流速為1.5 m/s，模擬研究外壁面?zhèn)鳠嵯禂捣謩e為50，100，200，250 W/（m2·K）時的情況。硅的輻射發(fā)射率取0.85，環(huán)境溫度取300 K，管壁面材料選擇導熱系數較大的硅，這樣可使燃燒室外壁面溫度分布比較均勻，有利于壁面材料保護。導熱產生的熱損失較大，有利于研究不同傳熱系數條件下燃燒室內火焰的熄滅情況。不同外壁面?zhèn)鳠嵯禂禇l件下對燃燒的影響見表1。

由表1可知，隨燃燒室外壁面?zhèn)鳠嵯禂翟黾?，燃燒室內的火焰最高溫度降低，火焰中心區(qū)域減小，并向燃燒室入口處移動，出口處氣體和燃燒室外壁面溫度也降低。外壁面?zhèn)鳠嵯禂颠^大，導致管內最高溫度低于點火溫度，使火焰不能持續(xù)。另外，隨燃燒室外壁面?zhèn)鳠嵯禂翟黾?，化學反應區(qū)域增加，出口處氫氣的質量分數增加。這是因為外壁面?zhèn)鳠嵯禂翟黾?，使熱能通過壁面?zhèn)鬟f給外部環(huán)境的熱損失增加。

外壁面?zhèn)鳠崾贡诿鏈囟鹊陀谧罡弑诿娌牧系目稍试S溫度。當外壁面?zhèn)鳠嵯禂禐?0 W/（m2·K）時，燃燒室內的火焰溫度最高，燃燒效率最高，但其外壁面溫度在1 100 K以上，比硅材料可允許最高溫度900 K高很多，所以外壁面?zhèn)鳠嵯禂挡⒉皇窃叫≡胶?。綜合考慮表1中各項數據，選用100 W/（m2·K）的傳熱系數。

4 結論

模擬試驗結果表明：靠近燃燒室入口處的火焰中心溫度隨導熱系數減小而增加，但燃燒室出口處混合氣體的溫度隨導熱系數減小而降低；隨燃燒室外壁面厚度增加，外壁面溫度降低，且壁面溫度分布更均勻；隨燃燒室外壁面?zhèn)鳠嵯禂翟黾?，燃燒室內的火焰溫度降低，燃燒室出口氣體溫度、外壁面溫度也相應降低。

參考文獻

[1] YUNFEI YAN， HAIBO WANG， WENLI PAN .Numerical study of effect of wall parameters on catalytic combustion characteristics of CH4/air in a heat recirculation micro combustor[J].Energy Conversion and Management，2016（118）：474–484.

[2] A.DI STAZIO，C.CHAUVEAU，G.DAYMA. Oscillating flames in micro-combustion[J].Combustion and Flame，2016（167）：392–394.

[3] WAN JL， YANG W， FAN AW，et al. A numerical investigation on combustion characteristics of H2/air mixture in a micro-combustor with wall cavities[J]. Hydrogen Energy，2014（39）：8 138-8 145.

[4] 潘劍鋒，陳琳琳，盧青波，等.微燃燒室內截面突變對氫氧預混燃燒的影響[J].江蘇大學學報（自然科學版），2017，38（2）：150-154.

[5] 冉景煜，秦昌雷，吳晟.壁面參數對甲烷微尺度催化燃燒的影響[J].重慶大學學報，2012（2）：1-9.

[6] 何正，朱林.微燃燒器外墻壁熱損失對微尺度燃燒特性的影響分析[J].機械工程，2015（9）：74-77.endprint

Abstract： According to the premixed combustion process between hydrogen and air in micro-scale combustion chamber internal， the computing software STAR - CD was used to establish micro combustion of the physical and mathematical model. On the basis of experimental verification the wall parameters such as wall thermal conductivity， wall thickness and outside wall heat transfer coefficient was simulated to the effect on hydrogen and air premixed combustion. The results showed that： The close to the flame center at the entrance to the combustion chamber temperature increased with the decrease of the coefficient of thermal conductivity， but the mixed gas temperature at the combustion chamber exit reduced with the decrease of the coefficient of thermal conductivity； With the increase of burning outdoor wall thickness， outer wall temperature was reduced， and the wall temperature distribution more uniform； With the increase of burning outside wall heat transfer coefficient， flame temperature inside the combustion chamber was reduced， and gas temperature at the combustion chamber exit and outside wall temperature was correspondingly reduced.

Key words： micro combustor； wall parameter； numerical simulation； combustion characteristic； hydrogenendprint