圓弧激波緩聚點(diǎn)火及后續(xù)燃燒傳播

侯自豪, 楊劍挺, 朱雨建,*, 楊基明

(1. 中國(guó)科學(xué)技術(shù)大學(xué) 近代力學(xué)系, 合肥 230027; 2. 昆明理工大學(xué) 民航與航空學(xué)院, 昆明 650504)

0 引 言

為了實(shí)現(xiàn)更經(jīng)濟(jì)、有效的高超聲速飛行，一些基于爆轟的推進(jìn)方式相繼被提出，如脈沖爆轟發(fā)動(dòng)機(jī)、斜爆轟發(fā)動(dòng)機(jī)、旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)等。在超聲速流動(dòng)環(huán)境中組織有效的燃燒，從而提供充足的動(dòng)力，是這些推進(jìn)系統(tǒng)共同的研究熱點(diǎn)和難點(diǎn)之一[1-2]。其中，可燃?xì)怏w經(jīng)氣動(dòng)壓縮處理(壓縮波、激波、壓縮構(gòu)型等)后發(fā)生點(diǎn)火以及燃燒的發(fā)展和傳播過(guò)程倍受關(guān)注。已有研究表明，燃料的點(diǎn)火過(guò)程和燃燒特性受流場(chǎng)的影響涉及到諸多因素[3]，而實(shí)際工程中的點(diǎn)火和燃燒現(xiàn)象更是復(fù)雜多變。因此，發(fā)展并采用合適的研究方法來(lái)解析這種復(fù)雜現(xiàn)象中的關(guān)鍵機(jī)理是相關(guān)研究的迫切需求。

運(yùn)動(dòng)激波誘導(dǎo)點(diǎn)火是研究點(diǎn)火和燃燒特性的一種重要方法。其中，平面激波誘導(dǎo)點(diǎn)火因?qū)嶒?yàn)上較容易實(shí)現(xiàn)且激波壓縮后的狀態(tài)相對(duì)明確而被廣泛應(yīng)用[4-9]。但采用平面入射激波誘導(dǎo)點(diǎn)火難以構(gòu)造多樣的點(diǎn)火環(huán)境(溫度梯度、組分濃度梯度、點(diǎn)火延遲時(shí)間梯度等)。此外，為產(chǎn)生足夠強(qiáng)的平面激波，常采用平面反射激波誘導(dǎo)點(diǎn)火，此時(shí)反射激波與入射激波波后邊界層的相互作用會(huì)對(duì)點(diǎn)火形成干擾[10-11]，破壞了該方法構(gòu)造點(diǎn)火環(huán)境的簡(jiǎn)潔性。與此相對(duì)，匯聚激波誘導(dǎo)點(diǎn)火則可構(gòu)造多樣的點(diǎn)火環(huán)境。Gelfand等詳細(xì)討論了平面激波進(jìn)入拋物形、圓弧形、楔形等不同形狀凹腔內(nèi)匯聚點(diǎn)火的現(xiàn)象[12-13]，王春、滕宏輝、何立明等也開(kāi)展了大量研究[14-15]。需要指出的是，上述激波匯聚過(guò)程往往伴隨復(fù)雜的激波反射、干涉等現(xiàn)象，這給有效評(píng)估點(diǎn)火環(huán)境以及分析點(diǎn)火發(fā)生與燃燒傳播的動(dòng)力學(xué)機(jī)理帶來(lái)不確定性。

近年來(lái)，作者實(shí)驗(yàn)室將一種采用激波動(dòng)力學(xué)方法反設(shè)計(jì)獲得的圓柱面匯聚激波技術(shù)[16]應(yīng)用于激波誘導(dǎo)點(diǎn)火研究[17]。這種經(jīng)設(shè)計(jì)的柱面匯聚激波在光滑連續(xù)的不斷增強(qiáng)過(guò)程中，在流場(chǎng)中誘導(dǎo)產(chǎn)生空間連續(xù)的準(zhǔn)一維的梯度場(chǎng)(無(wú)復(fù)雜的激波干擾現(xiàn)象)，并實(shí)現(xiàn)了梯度場(chǎng)中預(yù)混可燃?xì)怏w的點(diǎn)火與燃燒波傳播。該研究發(fā)現(xiàn)：自點(diǎn)火誘導(dǎo)的初始燃燒波其傳播速度遠(yuǎn)大于當(dāng)?shù)谻J(Chapman-Jouguet)爆轟波的速度；這種超高速的燃燒波事實(shí)上是一種化學(xué)反應(yīng)“相波”，即各質(zhì)點(diǎn)自發(fā)化學(xué)演變過(guò)程中某一階段的時(shí)空分布，或依據(jù)相關(guān)研究稱(chēng)之為“自發(fā)反應(yīng)波”(spontaneous reaction wave/spontaneous auto-ignition wave)。根據(jù)Zeldovich的理論[18]，這種自發(fā)反應(yīng)波在減速轉(zhuǎn)變?yōu)楸Z波的過(guò)程中一般不形成解耦的激波和燃燒波面結(jié)構(gòu)，這有別于一般的DDT方式和直接起爆方式。

本文作為上述研究的延續(xù)，采用實(shí)驗(yàn)與數(shù)值模擬相結(jié)合的方法，對(duì)這種匯聚激波誘導(dǎo)預(yù)混可燃?xì)怏w點(diǎn)火現(xiàn)象作進(jìn)一步的研究，著重關(guān)注點(diǎn)火建立的爆轟波在后續(xù)傳播中遭遇“相波”所導(dǎo)致的特殊現(xiàn)象。

1 研究方法

1.1 實(shí)驗(yàn)方法

圖1為本研究所采用的激波管實(shí)驗(yàn)平臺(tái)示意圖。驅(qū)動(dòng)段和被驅(qū)動(dòng)段為等截面方管，長(zhǎng)度分別為2.5 m、1.5 m，內(nèi)截面大小為40 mm×70 mm。圖2顯示了實(shí)驗(yàn)段內(nèi)流道構(gòu)型以及觀察窗口的位置，與實(shí)驗(yàn)段末端相連的是等截面方管的泄壓段，長(zhǎng)度為2.5 m，內(nèi)截面大小為40 mm×70 mm。激波管平臺(tái)運(yùn)行時(shí)，驅(qū)動(dòng)段充入高壓驅(qū)動(dòng)氣體，被驅(qū)動(dòng)段、實(shí)驗(yàn)段及泄壓段內(nèi)充入預(yù)混可燃?xì)怏w；首先采用電控破膜產(chǎn)生強(qiáng)度較弱的平面激波，不足以直接起爆；平面激波在等截面管道運(yùn)動(dòng)一段距離,達(dá)到較理想狀態(tài)后進(jìn)入實(shí)驗(yàn)段，并與依據(jù)激波動(dòng)力學(xué)反設(shè)計(jì)的內(nèi)流道型線(xiàn)(紅色實(shí)線(xiàn))相互作用，從而連續(xù)光滑地轉(zhuǎn)變?yōu)閳A弧激波；圓弧激波在楔形段內(nèi)維持相似波面形態(tài)進(jìn)一步緩慢匯聚增強(qiáng)，并誘導(dǎo)發(fā)生點(diǎn)火。采用曝光時(shí)間為1 μs、拍攝頻率超過(guò)250000 fps的高速紋影系統(tǒng)對(duì)觀察窗內(nèi)流場(chǎng)進(jìn)行觀測(cè)，捕捉流場(chǎng)的演變過(guò)程。

圖1 實(shí)驗(yàn)平臺(tái)示意圖

圖2 實(shí)驗(yàn)段示意圖

1.2 數(shù)值計(jì)算方法

圖3為相同構(gòu)型的匯聚角度為8°的數(shù)值計(jì)算域。內(nèi)流道沿軸向分成3部分：平直段、弧形過(guò)渡段(500.0 mm

圖3 匯聚角度8°數(shù)值計(jì)算域

本文對(duì)以上構(gòu)型內(nèi)流動(dòng)開(kāi)展變截面準(zhǔn)一維無(wú)粘數(shù)值模擬。控制方程為帶化學(xué)反應(yīng)的歐拉方程。計(jì)算采用基于非均勻網(wǎng)格的MUSCL格式，其中數(shù)值通量的計(jì)算采用HLLC格式?；瘜W(xué)反應(yīng)采用Oran的氫-空氣基元反應(yīng)模型[19]。反應(yīng)源項(xiàng)采用分裂步法解耦處理，其中解耦后的源項(xiàng)常微分方程組采用VODE程序包[20]求解。上述數(shù)值模擬方法及程序的適用性和準(zhǔn)確性已在Yang[17]及Shi[21]等的研究中得到驗(yàn)證。

計(jì)算采用非均勻網(wǎng)格，在出口前400 mm的收縮段內(nèi)共布網(wǎng)格約1600個(gè)，其中最大網(wǎng)格尺度0.5 mm(左)，最小網(wǎng)格尺度0.1 mm(右)。當(dāng)自點(diǎn)火發(fā)生時(shí)，匯聚激波波面與著火點(diǎn)距離約29 mm，該距離內(nèi)網(wǎng)格數(shù)不少于200個(gè)。這一網(wǎng)格尺度可保證點(diǎn)火過(guò)程的正確模擬。

2 結(jié)果與討論

本文考察兩種不同配比的預(yù)混氣體，如表1所

表1 實(shí)驗(yàn)氣體組成和屬性Table 1 Compositions and properties of test gases

示。通過(guò)調(diào)配氣體組分使A、B兩種不同配比混合氣體的聲速保持一致，從而當(dāng)入射平面激波強(qiáng)度一致時(shí)，冷態(tài)流場(chǎng)和燃燒流場(chǎng)具有可比性。

2.1 實(shí)驗(yàn)結(jié)果

首先選取不可反應(yīng)混合氣體A研究點(diǎn)火環(huán)境的特性。圖4為冷態(tài)流動(dòng)紋影實(shí)驗(yàn)結(jié)果。

選定激波到達(dá)過(guò)渡段起始位置為t=0時(shí)刻，圓弧激波在楔形段匯聚過(guò)程中沒(méi)有出現(xiàn)激波反射、相干等現(xiàn)象，流場(chǎng)中無(wú)明顯間斷(t=0.296～0.392 ms)；激波到達(dá)出口后形成膨脹扇(t=0.392 ms)，而后經(jīng)歷足夠長(zhǎng)時(shí)間(約80 μs)，由于氣流壅塞形成反射激波往上游傳播。需要指出，圓弧激波匯聚增強(qiáng)的過(guò)程中，理論上必然掃過(guò)預(yù)混氣體的著火點(diǎn)，這一過(guò)程發(fā)生在反射激波形成前，因此點(diǎn)火環(huán)境基本不會(huì)受到反射激波與邊界層相互作用的干擾。

圖4 圓弧激波匯聚過(guò)程冷態(tài)流動(dòng)紋影實(shí)驗(yàn)結(jié)果

Fig.4Schlierenpicturesofthenonreactiveshockconvergingprocess

為了評(píng)估上述點(diǎn)火環(huán)境的有效性以及對(duì)自點(diǎn)火的基本物理現(xiàn)象有個(gè)清晰的認(rèn)識(shí)，考察了與冷態(tài)實(shí)驗(yàn)相同入射平面激波馬赫數(shù)和初始?jí)毫ο職怏wB的點(diǎn)火及燃燒傳播過(guò)程。圖5是氣體B燃燒流場(chǎng)的典型紋影實(shí)驗(yàn)結(jié)果，其中，t=0.392～0.416 ms為燃燒傳播初期的精細(xì)刻畫(huà)。

圓弧激波在匯聚過(guò)程中維持相似波面穩(wěn)定前行(t=0.356～0.368 ms)；匯聚增強(qiáng)到一定程度后，波后流場(chǎng)中發(fā)生自點(diǎn)火(t=0.380～0.392 ms)，反應(yīng)區(qū)向上游和下游分別擴(kuò)張。紅框內(nèi)紋影結(jié)果顯示初始燃燒波傳播階段并未形成間斷特征明顯的激波，而后壓縮波累積，激波逐漸凸顯(t=0.408 ms)，在激波后深黑色區(qū)域逐漸增加表明化學(xué)反應(yīng)強(qiáng)度增強(qiáng)(t=0.408～0.416 ms)。往上游擴(kuò)張的燃燒波面先形成激波-火焰復(fù)合結(jié)構(gòu)(t=0.416 ms)，然后激波和火焰形成相對(duì)穩(wěn)定的強(qiáng)耦合，演變?yōu)樽猿謧鞑サ谋Z波(t=0.420 ms)?？梢钥吹剑y影圖中波面的深黑色區(qū)域逐漸追趕上激波的淺色區(qū)域，表明反應(yīng)導(dǎo)致的流場(chǎng)密度降低趨勢(shì)陡峭，此時(shí)反應(yīng)較劇烈。同時(shí)，出口附近壅塞堆積產(chǎn)生反射激波，跟隨左行爆轟傳播。隨后，自持傳播的爆轟波在傳播一定距離后逐漸弱化為激波(t=0.480～0.500 ms)，紋影圖中顯示為燃燒波面逐漸由黑色變淡，直至僅剩下跨激波密度躍升所致的淺色區(qū)域，表明此時(shí)緊隨激波的化學(xué)反應(yīng)正逐漸減弱，而在激波的上游也逐漸凸顯出新的燃燒波面。

圖5 圓弧匯聚激波誘導(dǎo)自點(diǎn)火紋影實(shí)驗(yàn)結(jié)果

Fig.5Schlierenpicturesofself-ignitioninducedbyimplodingarc-shapedshockwave

上述點(diǎn)火及燃燒傳播發(fā)生在反射激波前，不同于反射激波匯聚點(diǎn)火，上述自點(diǎn)火的發(fā)生不受反射激波擾動(dòng)，體現(xiàn)了圓弧匯聚激波構(gòu)造點(diǎn)火環(huán)境的簡(jiǎn)潔性。

2.2 圓弧匯聚激波誘導(dǎo)的點(diǎn)火環(huán)境

為了對(duì)冷態(tài)實(shí)驗(yàn)現(xiàn)象進(jìn)行討論分析，開(kāi)展了與冷態(tài)實(shí)驗(yàn)相同工況下的準(zhǔn)一維數(shù)值模擬。從圖6中密度梯度的x-t圖與實(shí)驗(yàn)中各激波位置的對(duì)比結(jié)果來(lái)看，無(wú)論是入射激波還是反射激波，實(shí)驗(yàn)測(cè)得的軌跡雖稍滯后于數(shù)值結(jié)果(實(shí)驗(yàn)中粘性、壁面?zhèn)鳠嵋约坝^察窗玻璃與內(nèi)楔的間隙等均傾向削弱反射激波強(qiáng)度)，但在流動(dòng)特征上兩者定性一致。

圖6中的等值線(xiàn)是相應(yīng)的準(zhǔn)一維流動(dòng)x-t圖的溫度分布。結(jié)果表明，由于激波強(qiáng)度逐漸增強(qiáng)，波后氣體溫度連續(xù)升高，且緊隨激波面位置的等值線(xiàn)越來(lái)越密集(ΔT=50 K)。匯聚激波作用于靜止氣體使其獲得沿激波運(yùn)動(dòng)方向的誘導(dǎo)速度，對(duì)于某一特定流體質(zhì)點(diǎn)(如帶箭頭流動(dòng)跡線(xiàn)所示，兩質(zhì)點(diǎn)分別位于過(guò)渡段的首尾)，由于波后的非定常運(yùn)動(dòng)導(dǎo)致其不斷往高溫區(qū)運(yùn)動(dòng)。溫度的連續(xù)升高，導(dǎo)致各流體質(zhì)點(diǎn)的反應(yīng)誘導(dǎo)時(shí)間也在不斷縮小。

圖6 實(shí)驗(yàn)和數(shù)值紋影中入射激波和反射激波軌跡及溫度等值分布

Fig.6Experimentalandnumericaltrajectoriesofincidentshockandreflectedshockandx-tdiagramofthetemperaturecontourdistribution

圖7為t=0～0.60 ms時(shí)段沿程壓力、密度及溫度分布。入射激波往下游運(yùn)動(dòng)時(shí)，在激波波面處，壓力、密度、溫度均陡然上升，且隨著激波強(qiáng)度的增強(qiáng)，各激波波面處的參數(shù)值連續(xù)上升。同時(shí)可以看到，入射激波所誘導(dǎo)的溫度峰值位于激波波面處；而由于波后運(yùn)動(dòng)的非定常性，壓力、密度峰值均不在激波面位置，而是在波后一定距離處。在大約t=0.38 ms時(shí)，入射激波到達(dá)高度為4 mm的出口，氣流壅塞導(dǎo)致壓力、密度、溫度峰值連續(xù)攀升，在反射激波形成后，各峰值往上游前傳且逐步下降。

進(jìn)一步分析沿程溫度分布，認(rèn)為該點(diǎn)火環(huán)境由3部分構(gòu)成。入射激波在平直段運(yùn)動(dòng)時(shí)(如t=0 ms)，其波后溫度處處均勻；當(dāng)激波進(jìn)入弧線(xiàn)過(guò)渡段，激波強(qiáng)度連續(xù)增強(qiáng)，其波后溫度場(chǎng)存在梯度(如t=0.10 ms)；匯聚激波在過(guò)渡段末完成向圓弧激波的轉(zhuǎn)變，而后在楔形管道內(nèi)進(jìn)一步匯聚增強(qiáng)，此時(shí)所誘導(dǎo)的溫度分布呈現(xiàn)中間平緩、頭部尖銳的趨勢(shì)(如t=0.35 ms)。溫度沿程分布表明：在匯聚管道內(nèi)，圓弧激波誘導(dǎo)的溫度梯度沿軸向具有復(fù)雜關(guān)系，且由于流體質(zhì)點(diǎn)的運(yùn)動(dòng)使得溫度梯度存在非定常性。這區(qū)別于已有理論和數(shù)值模擬研究中人為給定的定常線(xiàn)性溫度梯度環(huán)境[22-23]。

圖7 空間壓力、密度、溫度分布的演變過(guò)程

Fig.7Evolutionprocessesofspatialpressure、densityandtemperaturedistribution

2.3 典型現(xiàn)象的分析

借助數(shù)值模擬可以對(duì)上述點(diǎn)火環(huán)境中自點(diǎn)火及點(diǎn)火之后往上游傳播燃燒波的發(fā)展過(guò)程進(jìn)行詳細(xì)刻畫(huà)。圖8表示燃燒流場(chǎng)中激波、燃燒波傳播軌跡實(shí)驗(yàn)和數(shù)值紋影之間的比較。圓點(diǎn)為實(shí)驗(yàn)中入射激波的位置，方形深色點(diǎn)為爆轟波的位置，方形淺色點(diǎn)為弱化激波的位置，三角形點(diǎn)為反射激波的位置；對(duì)于數(shù)值紋影中燃燒波的發(fā)展演變過(guò)程，淺色區(qū)域反映了激波壓縮所致的流場(chǎng)密度上升趨勢(shì)，深色區(qū)域反映了燃燒所致的流場(chǎng)密度下降趨勢(shì)。對(duì)比結(jié)果表明，盡管實(shí)驗(yàn)與數(shù)值紋影結(jié)果中激波和燃燒波的傳播有一定差異，但在波系結(jié)構(gòu)和流動(dòng)特征上兩者定性一致。

進(jìn)一步考察數(shù)值紋影中燃燒波的發(fā)展演變過(guò)程。圖8表明，自點(diǎn)火發(fā)生于入射激波后某一點(diǎn)，初始傳播階段并未形成激波和燃燒波面結(jié)構(gòu)；隨著反應(yīng)區(qū)的擴(kuò)張速度降低，淺色區(qū)域的激波逐漸形成，同時(shí)波后的化學(xué)反應(yīng)加強(qiáng)，導(dǎo)致深色區(qū)域也逐漸增加，兩者互為促進(jìn)，共同構(gòu)成激波-火焰復(fù)合結(jié)構(gòu)(如圖8中放大圖所示)；化學(xué)反應(yīng)的逐步增強(qiáng)推動(dòng)深色區(qū)域的反應(yīng)面與淺色區(qū)域的激波強(qiáng)耦合在一起，發(fā)展為向上游自持傳播的爆轟波；與實(shí)驗(yàn)現(xiàn)象一致，爆轟波傳播一定距離后開(kāi)始慢慢解耦，逐漸弱化為激波，同時(shí)波面上游凸顯出速度較快的燃燒波面。數(shù)值紋影預(yù)測(cè)到實(shí)驗(yàn)觀察窗口以外的現(xiàn)象，可以看出，新的燃燒波面處化學(xué)反應(yīng)強(qiáng)度較弱，隨著氣動(dòng)過(guò)程和化學(xué)反應(yīng)的耦合增強(qiáng)，又再次形成激波-火焰復(fù)合結(jié)構(gòu)，此后，反應(yīng)面和激波重新強(qiáng)耦合演變?yōu)樾碌谋Z波。

圖8 實(shí)驗(yàn)和數(shù)值紋影中激波和燃燒波軌跡

Fig.8Experimentalandnumericaltrajectoriesofshockwavesandcombustionwaves

從數(shù)值結(jié)果進(jìn)一步分析這種獨(dú)特現(xiàn)象的具體流場(chǎng)特性，圖9為流場(chǎng)壓力、密度、溫度及化學(xué)反應(yīng)中間產(chǎn)物(自由基OH)的沿程分布。

自由基OH的分布一定程度上表征了反應(yīng)區(qū)，自點(diǎn)火(約t=0.366 ms)后自由基OH含量迅速增加，并向上下游蔓延。初期燃燒時(shí)，高溫反應(yīng)區(qū)內(nèi)呈現(xiàn)出壓力均勻上升的特征；而從沿程密度分布來(lái)看(t=0.366～0.370 ms)，燃燒流場(chǎng)(B氣體)與冷態(tài)流場(chǎng)(A氣體)的空間密度分布基本沒(méi)有差別，說(shuō)明初期燃燒近似等容燃燒過(guò)程，氣體動(dòng)力學(xué)過(guò)程被有效地“關(guān)閉”。直到t=0.374 ms，流場(chǎng)中密度才出現(xiàn)明顯的強(qiáng)間斷，此時(shí)，激波-火焰復(fù)合結(jié)構(gòu)逐漸形成，這也解釋了圖5實(shí)驗(yàn)紋影中初期燃燒沒(méi)有形成激波的原因。隨后，左行壓縮波迅速累積為激波并與燃燒耦合傳播(t=0.434 ms)。從各參數(shù)峰值形狀來(lái)看，隨后從間斷波面處凸出新的燃燒波面，波面平緩并越過(guò)激波波面，后面的間斷燃燒波面逐漸弱化為在燃燒產(chǎn)物中傳播的激波，且激波強(qiáng)度進(jìn)一步衰減。此后，燃燒波面處的各參數(shù)峰值逐漸變得尖銳，預(yù)示激波和燃燒耦合的爆轟波再次形成(t=0.554 ms)。

圖9 空間壓力、密度、溫度、自由基OH分布演變過(guò)程

Fig.9Evolutionprocessofspatialpressure、densityandtemperaturefreeradicalOHdistribution

為了分析這一現(xiàn)象產(chǎn)生的機(jī)理，基于Zeldovich的自發(fā)傳播機(jī)制[18]， 從考慮純粹的化學(xué)反應(yīng)過(guò)程的角度，對(duì)圓弧匯聚激波所誘導(dǎo)的非定常點(diǎn)火環(huán)境中的“自發(fā)反應(yīng)波”[22-23]進(jìn)行數(shù)值計(jì)算。Zeldovich指出當(dāng)預(yù)混可燃?xì)怏w沿空間位置x的點(diǎn)火敏感性——點(diǎn)火時(shí)間ti——連續(xù)變化時(shí)，繼最小ti點(diǎn)發(fā)生點(diǎn)火后，速度為usp的 “自發(fā)反應(yīng)波”出現(xiàn)，并在點(diǎn)火時(shí)間ti較長(zhǎng)的相鄰位置通過(guò)因果無(wú)關(guān)的連續(xù)自燃而“傳播”，其中：

(1)

由圖10中質(zhì)點(diǎn)軌跡可知，在圓弧匯聚激波所誘導(dǎo)的特定點(diǎn)火環(huán)境中，預(yù)混氣體的點(diǎn)火時(shí)間ti包括兩部分：激波到達(dá)時(shí)間tshock和激波預(yù)處理狀態(tài)下的反應(yīng)誘導(dǎo)時(shí)間tind，其中反應(yīng)誘導(dǎo)時(shí)間是通過(guò)追蹤各質(zhì)點(diǎn)軌跡計(jì)算絕熱燃燒直到放熱反應(yīng)速率達(dá)到最大值來(lái)獲得。然后可以確定計(jì)算所得點(diǎn)火時(shí)間ti的軌跡，也即上述“自發(fā)反應(yīng)波”的路徑。更具體的，如果定義化學(xué)反應(yīng)過(guò)程的進(jìn)度為λ(λ=0：反應(yīng)開(kāi)始，λ=1：反應(yīng)完全)，則這種“自發(fā)反應(yīng)波”實(shí)際上是流場(chǎng)中各點(diǎn)的λ取某一相同進(jìn)度值時(shí)所處時(shí)空點(diǎn)的集合，其本質(zhì)上是一種化學(xué)反應(yīng)“相波”。區(qū)別于一般的動(dòng)力學(xué)波(壓縮波、稀疏波、激波)，所謂化學(xué)反應(yīng)“相波”的發(fā)展并非由前后介質(zhì)分子間質(zhì)量、動(dòng)量和能量的傳遞所驅(qū)動(dòng)，而是僅僅取決于當(dāng)?shù)亟橘|(zhì)自發(fā)的反應(yīng)進(jìn)程。從這個(gè)意義上，ZND爆轟波可視作前導(dǎo)激波波后的燃燒反應(yīng)相波。圖10顯示了以上描述的“相波”(選取λ=0.3、0.5)與實(shí)際燃燒波傳播軌跡的比較。對(duì)于該算例富氫的情況，有：

(2)

圖11展示了圖10中“相波”速度usp與數(shù)值紋影中實(shí)際燃燒波速度u(均相對(duì)流體質(zhì)點(diǎn))對(duì)比結(jié)果。結(jié)合圖10、11可知，對(duì)于往上游傳播的燃燒波，在初始著火點(diǎn)(x1,t1)，“相波”速度遠(yuǎn)超CJ速度，實(shí)驗(yàn)和數(shù)值結(jié)果均顯示初期傳播未形成激波，燃燒波后流動(dòng)呈現(xiàn)一定的弱爆轟波特征。隨著反應(yīng)區(qū)迅速擴(kuò)張，“相波”速度急劇衰減，反應(yīng)區(qū)內(nèi)壓縮波逐漸形成，燃燒波轉(zhuǎn)變?yōu)榧げê突鹧娴膹?fù)合結(jié)構(gòu)，初期燃燒近似在等容環(huán)境中進(jìn)行(t1～t2)，基本未受到氣動(dòng)過(guò)程的影響，因此燃燒波軌跡與“相波”軌跡基本重合。當(dāng)“相波”速度衰減至CJ速度以下時(shí)，左行壓縮波累積形

成激波，氣動(dòng)過(guò)程與化學(xué)反應(yīng)過(guò)程強(qiáng)耦合導(dǎo)致能量釋放和燃燒釋熱，激波-火焰復(fù)合結(jié)構(gòu)加速演變?yōu)樽猿謧鞑サ腃J爆轟波(t=t2)，使得實(shí)際燃燒波的速度超過(guò)“相波”的速度，因此“相波”滯后于實(shí)際燃燒波軌跡，此后，燃燒波速度維持在CJ速度附近(t2～t3)。

圖10 “相波”和數(shù)值紋影中燃燒波傳播軌跡的對(duì)比

Fig.10Comparisonoftrajectoriesof“spontaneousreactionwave”anddetonationwave

圖11 流動(dòng)系統(tǒng)中“相波”與實(shí)際燃燒波速度對(duì)比

Fig.11Comparisonofspeedsof“spontaneousreactionwave”anddetonationwave

從圖10可以看出，該工況下，由于上游特定的點(diǎn)火環(huán)境，梯度場(chǎng)中的“相波”的軌跡超過(guò)CJ爆轟波，成為新的燃燒波前。由于“相波”迅速擴(kuò)張并消耗爆轟波前反應(yīng)物，爆轟波失去化學(xué)反應(yīng)的加持逐漸弱化為在反應(yīng)產(chǎn)物中傳播的激波(t3～t4)，此時(shí)，燃燒波以“相波”的速度傳播，數(shù)值紋影中燃燒波軌跡與“相波”軌跡再次重合。而后，“相波”的速度減小至CJ速度以下，再次滯后于CJ爆轟波(t4～)，深黑色區(qū)域反應(yīng)面和淺色區(qū)域激波又重新耦合,象征爆轟波再次形成，演變?yōu)樾碌娜紵ㄇ啊?/p>

3 結(jié) 論

本文采用激波管實(shí)驗(yàn)和準(zhǔn)一維數(shù)值模擬的方法，獲得主要結(jié)論如下：

(1) 圓弧匯聚激波所誘導(dǎo)的非定常梯度環(huán)境由入射激波在不同形態(tài)下所誘導(dǎo)的3個(gè)梯度區(qū)共同構(gòu)成，該環(huán)境為自點(diǎn)火、“自發(fā)反應(yīng)波”的傳播以及爆轟波的形成研究提供了獨(dú)特的條件。

(2) 在圓弧匯聚激波所誘導(dǎo)的非定常梯度環(huán)境下，“自發(fā)反應(yīng)波”在爆轟起始過(guò)程中起著關(guān)鍵作用。從自點(diǎn)火出現(xiàn)至激波-火焰復(fù)合結(jié)構(gòu)(準(zhǔn)爆轟復(fù)合體)形成，一個(gè)近似等容爆炸過(guò)程出現(xiàn)，自圓弧激波后某處形成往上游傳播的、無(wú)激波耦合的反應(yīng)面，初期傳播速度遠(yuǎn)超CJ爆轟波速，呈現(xiàn)一定的弱爆轟波特征，其本質(zhì)是“自發(fā)反應(yīng)波”。

(3) 這種“自發(fā)反應(yīng)波”與爆轟傳播過(guò)程也存在著緊密聯(lián)系。在本文條件下，自點(diǎn)火后實(shí)現(xiàn)了“相波”趕超爆轟波的過(guò)程，分析認(rèn)為梯度環(huán)境中的“相波”速度usp與CJ爆轟波速的競(jìng)爭(zhēng)關(guān)系主導(dǎo)了這一過(guò)程的出現(xiàn)。