張?坤，Valiev Damir

微尺度細管內(nèi)熱損失對火焰動力學的影響

張?坤1, 2, 3，Valiev Damir1, 2, 3

(1. 清華大學能源與動力工程系，北京 100084；2. 清華大學燃燒能源中心，北京 100084；3. 清華大學熱科學與動力工程教育部重點實驗室，北京 100084)

針對微尺度細管內(nèi)火焰?zhèn)鞑栴}，通過數(shù)值模擬研究了點火能量和壁面熱損失對火焰?zhèn)鞑ツＪ降挠绊?數(shù)值模擬中觀察到熄火現(xiàn)象和低速穩(wěn)定傳播的火焰，但并未觀察到爆燃向爆轟轉(zhuǎn)變. 在以二維管道高度和初始已燃區(qū)長度構(gòu)成的坐標圖中，低速穩(wěn)定火焰?zhèn)鞑ツＪ匠尸F(xiàn)半島形，具有上、下兩個熄火極限，其中熄火下極限位于由初始點火區(qū)域面容比決定的雙曲線上方．對于低速穩(wěn)定火焰?zhèn)鞑ツＪ剑鹧婕舛说膫鞑ニ俣扰c層流火焰速度較為接近．火焰面前方的氣體幾乎靜止，火焰面后方有一個回流區(qū)，已燃氣體向閉口端運動．

爆燃；熄火；點火能量；熱損失

微尺度細管內(nèi)的火焰?zhèn)鞑ナ侨紵茖W中的一個基本問題．其在微型燃燒器、爆轟發(fā)動機以及火災(zāi)事故預(yù)防等領(lǐng)域有廣泛應(yīng)用．在物理尺度小于1mm的微型燃燒器和微型推進器中，成功點火和火焰穩(wěn)定是需要解決的關(guān)鍵問題[1]．在爆轟發(fā)動機中，一種起爆的方法是細管內(nèi)火焰加速及爆燃向爆轟轉(zhuǎn)變．這種方法可以縮短起爆的時間和距離[2]．成功實現(xiàn)火焰加速需要了解細管內(nèi)的火焰?zhèn)鞑ヒ?guī)律．阻火器廣泛應(yīng)用在石油天然氣輸送管道中．它的核心部件是一系列平行的細管，火焰穿過這些細管時會發(fā)生淬熄，進而阻止爆燃向爆轟轉(zhuǎn)變[3]．研究細管內(nèi)的火焰?zhèn)鞑タ梢暂o助阻火器的設(shè)計，在實現(xiàn)裝置輕量化的同時確保在火災(zāi)事故發(fā)生時保護油氣輸送管路．

在實際燃燒設(shè)備和阻燃裝置中通常存在壁面熱損失．壁面熱損失減小了已燃區(qū)氣體的熱膨脹，阻止火焰加速．而根據(jù)Schelkin火焰加速機理和Bychkov等[4]的火焰加速理論，已燃區(qū)氣體的熱膨脹對管內(nèi)火焰加速有重要作用．化學反應(yīng)釋熱引發(fā)已燃區(qū)氣體熱膨脹，壁面熱損失則會使其收縮．這兩種效應(yīng)的競爭對細管內(nèi)火焰的傳播模式的影響值得研究. Akkerman等[5]通過理論預(yù)測了穩(wěn)定傳播的低速燃燒波，但未在模擬中觀察到．Akkerman等[5]的理論將火焰面簡化成無限薄的間斷面，只考慮了氣體熱膨脹而忽略了詳細的火焰結(jié)構(gòu)．與之相比，反應(yīng)擴散模型考慮了詳細的火焰結(jié)構(gòu)但忽略了氣體熱膨脹[6]．如果采用該模型，理論上可以證明二維和三維行波解的存?在[7].數(shù)值模擬中也觀察到穩(wěn)定的低速燃燒波，例如Daou等[8]發(fā)現(xiàn)管道中央燃燒波持續(xù)傳播而壁面附近發(fā)生淬熄．Short等[9]拓展了Daou等[8]的模型，考慮了氣體密度變化，在模擬中觀察到穩(wěn)定傳播的低速燃燒波．Wang等[10]數(shù)值求解了包含化學反應(yīng)的可壓縮納維-斯托克斯方程，觀察到了熄火現(xiàn)象，壁面的熱損失阻止了爆燃向爆轟轉(zhuǎn)變．Han等[11]在數(shù)值模擬中觀察到火焰震蕩，但其并未轉(zhuǎn)變?yōu)楸Z．在管內(nèi)乙烯氧氣預(yù)混氣燃燒實驗中，Wu等[12]發(fā)現(xiàn)了爆燃向爆轟轉(zhuǎn)變，震蕩火焰，穩(wěn)定低速火焰和熄火現(xiàn)象．Du等[13]在毫米級乙烯氧氣預(yù)混火焰?zhèn)鞑嶒炛性黾赢斄勘龋来斡^察到穩(wěn)定低速火焰，忽停忽動火焰，震蕩火焰，震蕩火焰轉(zhuǎn)變?yōu)楸Z以及忽停忽動火焰轉(zhuǎn)變?yōu)楸Z．

文獻中描述了細管內(nèi)的多種火焰?zhèn)鞑ツＪ?，但鮮見探究初始條件對管內(nèi)火焰?zhèn)鞑ツＪ接绊懙恼撐模疄榱颂骄奎c火能量和壁面熱損失對光滑細管內(nèi)火焰?zhèn)鞑ツＪ降挠绊?，筆者針對二維細管高度和初始點火高溫區(qū)的長度做了一系列參數(shù)研究．二維細管高度和初始點火高溫區(qū)的長度刻畫了初始時刻加入系統(tǒng)中的能量和壁面熱損失．

1?數(shù)值方法

筆者數(shù)值求解了如下包含化學反應(yīng)和輸運過程的可壓縮納維-斯托克斯方程：

本文模擬中使用的代碼被大量應(yīng)用于層流預(yù)混火焰和燃燒穩(wěn)定性等問題的研究[15]，也應(yīng)用于冷壁面熄火研究[16]．求解器采用基于結(jié)構(gòu)網(wǎng)格的有限體積法求解了二維可壓縮納維-斯托克斯方程，采用半離散的方法，利用二階龍格-庫塔法進行時間推進，空間方向上對流項的離散格式為四階精度，黏性項的離散格式為二階精度．

2?結(jié)果與討論

2.1?點火能量和壁面熱損失對火焰?zhèn)鞑ツＪ降挠绊?/h3>

圖2 二維管道高度和初始點火高溫區(qū)長度對火焰?zhèn)鞑ツＪ降挠绊?/p>

圖3?火焰面尖端的速度

表1?熄火下極限

Tab.1?Quenching limit for the lower branch

2.2?對低速穩(wěn)定火焰的分析

下面對火焰尖端的速度做出估計．火焰面形狀可以描述為：

圖5?低速穩(wěn)定燃燒波火焰面附近的速度場和流線

將公式(13)和公式(12)代入公式(11)，可得：

3?結(jié)?語

本文的研究工作受到國家自然科學基金(項目編號52176118)和青年千人項目的支持．Florian Schmidt副教授為數(shù)值模擬提供幫助，筆者與繆承希進行了討論，在此表示感謝．

［1］ Ju Yiguang，Maruta Kaoru. Microscale combustion：Technology development and fundamental research[J].，2011，37(6)：669-715.

［2］ Wang Zhiwu，Zhang Yang，Huang Jingjing，et al. Ignition method effect on detonation initiation characteristics in a pulse detonation engine[J].，2016，93：1-7.

［3］ Wang Luqing，Ma Honghao，Shen Zhaowu，et al. Flame quenching by crimped ribbon flame arrestor：A brief review[J].，2019，38 (1)：27-41.

［4］ Bychkov Vitaly，Petchenko Arkady，Akkerman V’yacheslav，et al. Theory and modeling of accelerating flames in tubes[J].，2005，72(4)：046307.

［5］ Akkerman V’yacheslav，Law Chung K，Bychkov Vitaly，et al. Analysis of flame acceleration induced by wall friction in open tubes[J].，2010，22(5)：053606.

［6］ Matalon M，Matkowsky B J. Flames as gas dynamic discontinuities[J].，1982，124(1)：239-259.

［7］ 步真會. 點火型和單穩(wěn)型反應(yīng)擴散方程的多維傳播現(xiàn)象[D]. 蘭州：蘭州大學，2018.

Bu Zhenghui. Multidimensional Propagation Phenomena of Reaction-Diffusion Equations with Combustion and Monostable Nonlinearities[D]. Lanzhou：Lanzhou University，2018(in Chinese).

［8］ Daou J，Matalon M. Influence of conductive heat-losses on the propagation of premixed flames in channels[J].，2002，128(4)：321-339.

［9］ Short Mark，Kessler David. Asymptotic and numerical study of variable-density premixed flame propagation in a narrow channel[J].，2009，638：305-337.

［10］ Wang Cheng，Zhao Yongyao，Han Wenhu. Effect of heat-loss boundary on flame acceleration and deflagration-to-detonation transition in narrow channels[J].，2017，189(7-9)：1605-1623.

［11］ Han Wenhu，Huang Jin，Gu Gongtian，et al. Surface heat loss and chemical kinetic response in deflagration-to-detonation transition in microchannels[J].，2020，5(5)：053201.

［12］ Wu Minghsun，Wang Chanyu. Reaction propagation modes in millimeter-scale tubes for ethylene/oxygen mixtures[J].，2011，33(2)：2287-2293.

［13］ Du Ning，Kong Wenjun，Han Wenhu，et al. Flame propagation in millimeter-scale tubes for lean ethylene-oxygen mixtures[J].，2019，58(4)：1-11.

［14］ Travnikov O Yu，Liberman M A，Bychkov V V. Stability of a planar flame front in a compressible flow[J].，1997，9(12)：3935-3937.

［15］ Valiev Damir，Bychkov Vitaly，Akkerman V’ yacheslav，et al. Different stages of flame acceleration from slow burning to Chapman-Jouguet deflagration[J].，2009，80(3)：036317.

［16］ Dion Claude，Valiev Damir，Akkerman V’yacheslav，et al. Dynamics of flame extinction in narrow channels with cold walls：Heat loss versus acceleration[J].，2021，33(3)：033610.

Influence of Heat Loss on Flame Dynamics in Microscale Channels

Zhang Kun1, 2, 3，Valiev Damir1, 2, 3

(1. Department of Energy and Power Engineering，Tsinghua University，Beijing 100084，China；2. Center for Combustion Energy，Tsinghua University，Beijing 100084，China；3. Key Laboratory for Thermal Science and Power Engineering of Ministry of Education，Tsinghua University，Beijing 100084，China)

The numerical study is focused on the influence of ignition energy and heat loss on the flame propagation regimes in microscale narrow channels with non-slip walls. Within the range of channel heights adopted in simulations，extinction and quasi-steady deflagration are observed while deflagration-to-detonation transition is not. Flame quenching and propagation modes are mapped in a diagram against channel height and initial burned zone length. The peninsula-shaped diagram exhibits two quenching limits. The criterion for the lower quenching limit is formulated based on surface-to-volume ratio. For the quasi-steady deflagration，the flame tip velocity is around laminar flame speed. The flow field upstream of the flame is nearly motionless，and the gas in the burned zone is moving backwards towards the closed end.

deflagration；extinction；ignition energy；heat loss

TK401

A

1006-8740(2022)04-0397-06

10.11715/rskxjs.R202206006

2021-04-07.

國家自然科學基金資助項目(52176118).

張?坤（1996—??），男，博士研究生，zhang-k18@mails.tsinghua.edu.cn.

Valiev Damir，男，博士，副教授，dvaliev@tsinghua.edu.cn.

(責任編輯：梁?霞)