喻健良，高 遠(yuǎn)，閆興清，高 偉

(大連理工大學(xué)化工機(jī)械學(xué)院，遼寧 大連 116024)

?

高濃度氬氣稀釋氣體爆轟波臨界管徑和臨界間距關(guān)系

喻健良，高 遠(yuǎn)，閆興清，高 偉

(大連理工大學(xué)化工機(jī)械學(xué)院，遼寧 大連 116024)

建立圓管及環(huán)形管道系統(tǒng)研究臨近極限下爆轟波在管道內(nèi)傳播失效機(jī)理。選用C2H2+2.5O2+70%Ar氣體，采用光纖探針測(cè)量爆轟波在管道內(nèi)傳播速度，用煙跡法記錄管道內(nèi)爆轟波胞格結(jié)構(gòu)。結(jié)果表明：初始?jí)毫h(yuǎn)大于爆轟極限壓力時(shí)，爆轟波在管道內(nèi)以穩(wěn)定速度傳播；隨著初始?jí)毫Φ臏p小，爆轟波速度逐漸降低；當(dāng)初始?jí)毫σ欢〞r(shí)，爆轟波速度隨著管道尺寸的減小而逐漸減??；當(dāng)初始?jí)毫_(dá)到臨界壓力時(shí)，爆轟波在進(jìn)入到管道內(nèi)后其速度會(huì)逐漸衰減直至爆轟波完全失效。對(duì)于不同幾何尺寸的圓管與環(huán)管，通過(guò)引入無(wú)量綱參數(shù)d/λ及w/λ(d為圓管管徑，w為環(huán)管間距，λ為爆轟胞格尺寸)得出，爆轟波在管道內(nèi)傳播的臨界圓管直徑為環(huán)形間距的2倍，與理論模型結(jié)果相吻合，驗(yàn)證了穩(wěn)態(tài)氣體基于爆轟波波面曲率的失效機(jī)理。

爆炸力學(xué)；爆轟極限；臨界直徑；臨界間距；爆轟波波面曲率

爆轟波在管道內(nèi)傳播時(shí)需要滿(mǎn)足一定的初始條件(初始?jí)毫?、?dāng)量比等)及邊界條件(管道尺寸、幾何形狀)，若上述條件無(wú)法滿(mǎn)足，則爆轟波不能夠在管道內(nèi)自持穩(wěn)定傳播，并最終失效，此時(shí)達(dá)到爆轟極限[1]。爆轟極限的研究對(duì)深入理解爆轟波在管道內(nèi)的傳播機(jī)理具有重要意義。爆轟波波面是瞬態(tài)三維非穩(wěn)定結(jié)構(gòu)，且在傳播過(guò)程中會(huì)與管道壁面產(chǎn)生相互擾動(dòng)[2-4]，并導(dǎo)致爆轟波熱量、動(dòng)量及質(zhì)量的損失[5]。當(dāng)條件接近爆轟極限時(shí)，爆轟波的傳播變得非常復(fù)雜且不穩(wěn)定，在介于自持穩(wěn)定傳播與失效之間出現(xiàn)一系列的非穩(wěn)定傳播現(xiàn)象，如：?jiǎn)晤^螺旋爆轟、馳振爆轟、低速爆轟等[6-7]。近期研究表明，對(duì)于高濃度氬氣稀釋氣體，其化學(xué)反應(yīng)區(qū)具有層流結(jié)構(gòu)且胞格結(jié)構(gòu)非常規(guī)則，可用ZND模型描述其爆轟波結(jié)構(gòu)[8-10]。對(duì)于此類(lèi)氣體，爆轟波波面非穩(wěn)定性對(duì)于其傳播影響很小，因此，在研究此類(lèi)氣體的爆轟極限時(shí)，可以忽略爆轟波波面非穩(wěn)定性因素而僅從邊界條件對(duì)其傳播影響著手研究爆轟極限[11]。根據(jù)Fay-Dabora理論模型[12]，管道壁面的邊界層效應(yīng)導(dǎo)致爆轟波化學(xué)反應(yīng)區(qū)中的流場(chǎng)擴(kuò)散，使得爆轟波波陣面產(chǎn)生彎曲，導(dǎo)致熱量及動(dòng)量的損失，從而使得爆轟速度衰減。當(dāng)爆轟波波面曲率達(dá)到某個(gè)臨界值時(shí)，爆轟波波速特征值不再存在，此時(shí)達(dá)到爆轟極限，爆轟波失效。

本文中研究高濃度氬氣稀釋氣體在2種不同幾何形狀(圓管、環(huán)管)管道內(nèi)的傳播現(xiàn)象及爆轟極限，對(duì)比分析臨界管徑與臨界間距之間的關(guān)系，進(jìn)一步驗(yàn)證和闡述高濃度氬氣稀釋氣體爆轟波的失效機(jī)理。

1 實(shí)驗(yàn)裝置

實(shí)驗(yàn)裝置見(jiàn)圖1。圓管爆轟管道由驅(qū)動(dòng)段和測(cè)試段組成(圖1(a))，驅(qū)動(dòng)段為長(zhǎng)1.3 m，內(nèi)徑為70 mm，壁厚4.5 mm的鋼管。測(cè)試段為透明的聚碳酸酯塑料管。環(huán)形爆轟管道由驅(qū)動(dòng)段、過(guò)渡段及測(cè)試段組成(圖1(b))，其中間距(w)為3.2及5.9 mm的環(huán)形通道通過(guò)在內(nèi)徑為50.8 mm的聚碳酸酯塑料管內(nèi)分別插入外徑為44.4與39.0 mm 光滑銅管構(gòu)成。圓管及環(huán)管幾何尺寸見(jiàn)表1～2，其中：L為管道長(zhǎng)度，t為管壁厚度，d為圓管內(nèi)徑。由于環(huán)管間距相對(duì)于環(huán)管直徑很小，故爆轟波在環(huán)管內(nèi)傳播時(shí)，其徑向彎曲曲率對(duì)其傳播的影響可以忽略，因此可以將其視為二維通道。采用環(huán)管而不采用矩形管道的優(yōu)點(diǎn)在于，環(huán)管能夠避免矩形管道的邊緣效應(yīng)[11,13]。采用高壓電火花起爆，并在驅(qū)動(dòng)段內(nèi)靠近起爆端添加一長(zhǎng)200 mm的Shchelkin螺旋用以增加對(duì)爆炸波的擾動(dòng)，從而更易在驅(qū)動(dòng)段內(nèi)形成穩(wěn)態(tài)的爆轟。采用分壓法配置C2H2+2.5O2+70%Ar氣體。研究表明C2H2+2.5O2+70%Ar混合氣體爆轟時(shí)具有規(guī)則胞格結(jié)構(gòu)，且化學(xué)反應(yīng)區(qū)具有層流結(jié)構(gòu)[1,4,9]，因此被稱(chēng)為穩(wěn)定氣體。

圖1 實(shí)驗(yàn)裝置Fig.1 Schematic of the experimental apparatus

實(shí)驗(yàn)前采用真空泵對(duì)管道抽真空，使管道內(nèi)的壓力低于10 Pa，然后通過(guò)閥門(mén)控制加入混合氣體，并通過(guò)電子壓力表控制進(jìn)氣量。在驅(qū)動(dòng)段上安裝PCB壓力傳感器(型號(hào)：113B24，響應(yīng)頻率：500 kHz，量程：0～6.9 MPa，輸出電壓：0～5 V)，用以測(cè)定爆轟波的壓力和到達(dá)各壓力傳感器的時(shí)間，從而計(jì)算在驅(qū)動(dòng)段內(nèi)是否形成并達(dá)到穩(wěn)態(tài)的爆轟波。在測(cè)試段管道上等間距布置光纖探針，通過(guò)光纖探針測(cè)量測(cè)試段內(nèi)爆轟波到達(dá)各測(cè)點(diǎn)的時(shí)間，并計(jì)算出爆轟波在相鄰兩測(cè)點(diǎn)間的實(shí)際傳播速度。通過(guò)CEA程序[14]計(jì)算混合氣體在該初始狀態(tài)下的理論爆轟CJ速度。采用煙跡法記錄爆轟波在管道內(nèi)胞格變化規(guī)律并測(cè)量對(duì)應(yīng)條件下胞格尺寸。

表1 圓管幾何尺寸

表2 環(huán)管幾何尺寸

2 實(shí)驗(yàn)結(jié)果及討論

2.1 爆轟極限臨界壓力的確定

圖2 圓管內(nèi)爆轟波傳播Fig.2 Propagation of detonation wave in the round tube

圖3 環(huán)管內(nèi)爆轟波傳播Fig.3 Propagation of detonation wave in the annular channel

2.2 初始?jí)毫肮艿缼缀纬叽鐚?duì)爆轟速度的影響

圖4為不同管道內(nèi)，初始?jí)毫εc爆轟速度關(guān)系。隨著初始?jí)毫Φ慕档停Z波速度逐漸減小，這是由于C2H2+2.5O2+70%Ar穩(wěn)定氣體化學(xué)反應(yīng)區(qū)具有典型的ZND結(jié)構(gòu)。當(dāng)混合氣體初始?jí)毫档蜁r(shí)，爆轟波誘導(dǎo)區(qū)長(zhǎng)度增加，邊界層擴(kuò)散導(dǎo)致的能量損失逐漸增大，從而導(dǎo)致爆轟波的化學(xué)反應(yīng)變緩，并使得爆轟速度降低[12]。在臨界壓力處，當(dāng)爆轟波失效時(shí)，最小爆轟速度與管道直徑、環(huán)形間距無(wú)關(guān)，其值約為0.8vCJ。同時(shí)由圖4得到，當(dāng)初始?jí)毫σ欢〞r(shí)，爆轟速度隨管道尺寸減小而降低。這是由于在小內(nèi)徑管道中，邊界層所占比例增加，邊界層效應(yīng)增大，從而導(dǎo)致爆轟速度的降低。

圖4 爆轟速度隨著初始?jí)毫Φ淖兓疐ig.4 Detonation velocity versus initial pressure

對(duì)于圓形管，當(dāng)管徑(d)為1.5、3.2、12.7、31.7、50.8 mm時(shí)，對(duì)應(yīng)的爆轟極限極限壓力(pc)分別為20、9.8、3.2、1.5、1.1 kPa；對(duì)于環(huán)形管，當(dāng)間距(w)為3.2和5.9 mm時(shí)，對(duì)應(yīng)的爆轟極限極限壓力分別為7.0和3.5 kPa。通過(guò)比較不同幾何尺寸管道所對(duì)應(yīng)的臨界壓力可以看出，爆轟極限臨界壓力隨著管徑或間距的減小而顯著增大，說(shuō)明隨著管徑或間距的減小，管道內(nèi)壁邊界層所占的比例不斷增加，在小尺寸管道內(nèi)管道邊界層效應(yīng)對(duì)爆轟波在管道內(nèi)的傳播起到主導(dǎo)作用。

2.3 臨界圓管內(nèi)徑與環(huán)管間距的關(guān)系

當(dāng)混合氣初始?jí)毫σ欢〞r(shí)，圓管臨界管徑尺寸為環(huán)管臨界間距的2倍[12]，這是因?yàn)楸Z波在環(huán)管內(nèi)傳播時(shí)，由于環(huán)管可視為二維通道，爆轟波波面彎曲只發(fā)生在一個(gè)平面上，而對(duì)于圓管，爆轟波波面彎曲則發(fā)生在兩個(gè)相互垂直的平面上。本實(shí)驗(yàn)中共采用了5種不同內(nèi)徑的圓管及2種間距的環(huán)管，并得到其對(duì)應(yīng)的臨界壓力。由于臨界壓力值不同，為了消除初始?jí)毫Σ煌瑢?duì)結(jié)果的影響，引入無(wú)量綱參數(shù)dc/λ及wc/λ，通過(guò)無(wú)量綱化的方法消除初始?jí)毫?duì)結(jié)果的影響，并得到一個(gè)普適規(guī)律。無(wú)量綱參數(shù)dc/λ及wc/λ中，dc、wc和λ分別為臨界管道直徑、臨界環(huán)形間距及胞格尺寸，其中dc和wc可反映爆轟波在管道內(nèi)傳播的邊界條件，而胞格尺寸作為爆轟動(dòng)態(tài)參數(shù)[1]，可反映爆轟氣體的敏感度。對(duì)于給定的氣體，胞格尺寸大小與初始?jí)毫ο嚓P(guān)，因此胞格尺寸可間接反映爆轟氣體初始?jí)毫Α?/p>

圖5為典型的爆轟波胞格結(jié)構(gòu)隨初始?jí)毫ψ兓?guī)律(圓管，d=50.8 mm)。由圖可知，C2H2+2.5O2+70%Ar具有規(guī)則的胞格結(jié)構(gòu)，且隨著初始?jí)毫Φ闹饾u降低，爆轟波胞格尺寸逐漸增加。通過(guò)實(shí)驗(yàn)測(cè)定了不同管道及初始?jí)毫ο卤Z波胞格尺寸(如圖6所示)，并采用了最小二乘法擬合出胞格尺寸與初始?jí)毫χg的關(guān)系式：

(1)

式中：λ單位為mm，p0單位為kPa。通過(guò)(1)式，可計(jì)算得到各管道在臨界壓力處對(duì)應(yīng)的胞格尺寸，并計(jì)算出相應(yīng)的dc/λ、wc/λ，結(jié)果見(jiàn)表3～4。由表3～4可以得出，dc/λ和wc/λ的平均值分別為0.185和0.095 4。對(duì)于不同內(nèi)徑的圓管及環(huán)管，在極限壓力處，dc/λ與wc/λ的平均比值為1.95，這與理論結(jié)果(dc/wc=2)相吻合，證明對(duì)于穩(wěn)態(tài)氣體，可采用理論分析預(yù)測(cè)其爆轟極限。

圖5 爆轟胞格結(jié)構(gòu)隨初始?jí)毫Φ淖兓?d=50.8 mm)Fig.5 Stuctures of the detonation cellulars vary with the initial pressure (d=50.8 mm)

圖6 爆轟胞格尺寸與初始?jí)毫Φ年P(guān)系Fig.6 Detonation cellular size versus initial pressure

對(duì)于高濃度氬氣稀釋的穩(wěn)態(tài)氣體，其胞格結(jié)構(gòu)非常規(guī)則，因此可采用經(jīng)典的ZND模型進(jìn)行分析，故除了采用胞格尺寸所建立無(wú)量綱參數(shù)外，亦可采用ZND誘導(dǎo)區(qū)長(zhǎng)度(ΔZND)建立無(wú)量綱參數(shù)，以此分析臨界管徑與臨界間距間關(guān)系。ZND誘導(dǎo)區(qū)長(zhǎng)度可通過(guò)Chemkin程序[15]和San Diego化學(xué)反應(yīng)動(dòng)力學(xué)原理[16]計(jì)算得到，結(jié)果見(jiàn)表5和6。B.Bazathazajan等[17]證明了San Diego化學(xué)反應(yīng)機(jī)理對(duì)于計(jì)算C2H2+2.5O2+70%Ar氣體爆轟波ZND誘導(dǎo)區(qū)長(zhǎng)度的有效性。由表5和6可知，dc/ΔZND與wc/ΔZND的平均值分別為24.07和12.25。dc/wc的平均值為1.96。當(dāng)采用ZND誘導(dǎo)區(qū)長(zhǎng)度建立無(wú)量綱參數(shù)時(shí)，結(jié)果同樣與理論結(jié)果吻合。進(jìn)一步證明了采用爆轟波波面曲率失效機(jī)理魚(yú)刺爆轟極限的正確性。

表3 圓管實(shí)驗(yàn)結(jié)果

表4 環(huán)管實(shí)驗(yàn)結(jié)果

表5 圓管ZND化學(xué)反應(yīng)區(qū)長(zhǎng)度

表6 環(huán)管ZND化學(xué)反應(yīng)區(qū)長(zhǎng)度

3 結(jié) 論

(1) 當(dāng)初始?jí)毫h(yuǎn)大于爆轟極限臨界壓力時(shí)，爆轟波在管道內(nèi)以一穩(wěn)定速度傳播。隨著預(yù)混氣初始?jí)毫Φ臏p小，爆轟波在管道內(nèi)傳播的速度逐漸降低。當(dāng)初始?jí)毫_(dá)到臨界壓力時(shí)，爆轟波在進(jìn)入到管道內(nèi)并傳播一段距離后，其速度將逐漸衰減直至爆轟波衰變成爆燃波并完全失效。通過(guò)比較不同幾何尺寸管道內(nèi)爆轟波傳播速度可知，爆轟波速度隨著管道幾何尺寸的減小而減小，表明爆轟波在管道內(nèi)傳播時(shí)受到邊界條件及初始條件的影響。

(2) 對(duì)于不同幾何尺寸管道，通過(guò)無(wú)量綱參數(shù)分析研究得出，爆轟波在管道內(nèi)傳播時(shí)臨界圓管直徑為臨界環(huán)管間距的2倍，這與理論預(yù)測(cè)相符合。表明高濃度氬氣稀釋氣體爆轟波在管道內(nèi)傳播時(shí)，管道壁面的邊界層效應(yīng)所導(dǎo)致的爆轟波化學(xué)反應(yīng)區(qū)內(nèi)的動(dòng)量、熱量損失是導(dǎo)致此種高濃度氬氣稀釋氣體爆轟波在管道內(nèi)傳播失效的主要原因。通過(guò)對(duì)爆轟波在圓形及環(huán)形管道內(nèi)的失效分析,進(jìn)一步驗(yàn)證了此種高濃度氬氣稀釋氣體基于爆轟波波面曲率的失效機(jī)理。

[1] Lee J H S. The detonation phenomenon[M]. Cambrige, UK: Cambridge University Press, 2008.

[2] Dupre G, Peraldi O, Lee J H S, et al. Propagation of detonation waves in an acoustic absorbing-walled tube[J]. Progress in Astronautics and Aeronautics, 1988,114:248-263.

[3] Teodoczyk A, Lee J H S. Detonation attenuation by foams and wire meshes lining the walls[J]. Shock Waves, 1995,4(4):225-236.

[4] Raduledscu M I, Lee J H S. The failure mechanism of gaseous detonations: Experiments in porous wall tubes[J].Combustion and Flame, 2002,131:29-46.

[5] Zeldovich Y B. On the theory of the propagation of detonation in gaseous systems[R]. Soviet Union: Soviet Physics-JETP, 1940.

[6] Lee J J, Dupre G, Knystautas R, et al. Doppler interferometer study of unstable detonations[J]. Shock Waves, 1995,5(3):175-181.

[7] Camargo A, Ng H D, Chao J, et al. Propagation of near-limit gaseous detonations in small diameter tubes[J]. Shock Waves, 2010,20(6):499-508.

[8] Radulescu M I. The propagation and failure mechanism of gaseous detonations: Experiments in porous-walled tubes[D].Montreal, Canada: McGill University, 2003.

[9] Radulescu M I, Ng H D, Lee J H S, et al. The effect of argon dilution on the stability of acetylene-oxygen detonations[J]. Proceedings of the Combustion Institute, 2002,29(2):2825-2831.

[10] Ng H D, Radulescu M I, Higgins A J, et al. Numerical investigation of the instability for one-dimensional Chapman-Jouguet detonations with chain-branching kinetics[J]. Combustion Theory and Modeling, 2005,9:385-401.

[11] Chao J, Ng H D, Lee J H S. Detonation limits in thin annular channels[J]. Proceedings of the Combustion Institute, 2009,32(2):2349-2354.

[12] Fay J A. Two-dimensional gaseous detonations: Velocity deficit[J]. Physics of Fluids, 1959,2(3):283-289.

[13] Meredith J, Ng H D, Lee J H S. Detonation diffraction from an annular channel[J]. Shock Waves, 2010,20(6):449-455.

[14] Mcbride B J, Gordon S. Computer program for calculation of complex chemical equilibrium compositions and applications[R]. NASA, 1996.

[15] Kee R J, Rupley F M, Millerja. Chemkin-Ⅱ: A fortran chemical kinetics package for the analysis of gas-phase chemical kinetics[R]. Sandia National Laboratories, 1989.

[16] The San Diego Mechanism. Chemical-kinetic mechanisms for combustion applications[EB/OL]. http://web.eng.ucsd.edu/mae/groups/combustion/mechanism.html

[17] Varatharajan B, Williams F A. Chemical-kinetic descriptions of high-temperature ignition and detonation of acetylene-oxygen diluents systems[J]. Combustion and Flame, 2001,124(4):624-645.

(責(zé)任編輯 王小飛)

Correlation between the critical tube diameter and annular interval for detonation wave in high-concentration argon diluted mixtures

Yu Jian-liang, Gao Yuan, Yan Xing-qing, Gao Wei

(SchoolofChemicalMachinery,DalianUniversityofTechnology,Dalian116024,Liaoning,China)

Detonation tube including driver section and test section was built to investigate the failure mechanism of detonation wave near the limits. The mixture of C2H2+2.5O2+70%Ar was investigated experimentally. Fiber optics was used to measure detonation velocity. Smocked foils were used to record the detonation cellular structure. The results show that, with the initial pressure far lager than the critical pressure, detonation wave propagates at a constant value in the tubes. Detonation velocity decreases with the decreasing initial pressure. With a given initial pressure, the detonation velocity decreased as the tube diameter (or channel interval) decreased. Under the critical pressure, the detonation velocity propagated a short distance in the tubes and then decreased gradually until complete failure. For different geometries tubes and channels, by introducing dimensionless parameterd/λandw/λ(dthe dameter of the round tube,wthe interval of the annular channel andλthe size of detonation cellular), the results show that the critical thickness is half of the critical diameter. Good agreement is found between the experimental measurements in both geometries which supports this conclusion and theoretical mode. The failure mechanisms based on the detonation front curvature for stable detonation in mixtures that are highly argon diluted are well defined.

mechanics of explosion; detonation limits; critical diameter; critical thickness; eetonation front curvature

10.11883/1001-1455(2015)04-0603-06

2013-12-21；

2014-05-04

國(guó)家自然科學(xué)基金項(xiàng)目(50974027)

喻健良(1963- )，男，博士，教授，yujianliang@dl.cn。

O382 國(guó)標(biāo)學(xué)科代碼： 13035

A