進(jìn)氣總壓對(duì)連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)爆轟影響的二維數(shù)值模擬

李寶星，翁春生

(南京理工大學(xué) 瞬態(tài)物理國家重點(diǎn)實(shí)驗(yàn)室，南京　210094)

?

進(jìn)氣總壓對(duì)連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)爆轟影響的二維數(shù)值模擬

李寶星，翁春生

(南京理工大學(xué) 瞬態(tài)物理國家重點(diǎn)實(shí)驗(yàn)室，南京210094)

為了研究?jī)上噙B續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)爆轟特性，基于二維守恒元和求解元的方法，對(duì)汽油/富氧空氣為燃料的發(fā)動(dòng)機(jī)爆轟過程進(jìn)行數(shù)值模擬，獲得了燃燒室內(nèi)部流場(chǎng)結(jié)構(gòu)和爆轟波傳播特性，對(duì)不同進(jìn)氣總壓條件下發(fā)動(dòng)機(jī)的爆轟性能進(jìn)行了計(jì)算。結(jié)果表明，入口處流場(chǎng)隨著爆轟波傳播呈周期性變化；爆轟波前的預(yù)混燃料層是形成穩(wěn)定連續(xù)旋轉(zhuǎn)爆轟的關(guān)鍵，當(dāng)進(jìn)氣總壓過小時(shí)，點(diǎn)火后入口處不能形成有效預(yù)混燃料層，致使爆轟熄滅；出入口處的壓力、溫度、密度均受到進(jìn)氣總壓的影響；發(fā)動(dòng)機(jī)的比沖與進(jìn)氣總壓成正比例關(guān)系。

氣液兩相；連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)；內(nèi)部流場(chǎng)；預(yù)混燃料層

0　引言

連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)(Continuous Rotating Detonation Engine，CRDE)是利用爆轟波在環(huán)形燃燒室內(nèi)連續(xù)旋轉(zhuǎn)傳播，爆轟產(chǎn)物從出口排出，而產(chǎn)生穩(wěn)定持續(xù)推力的新型發(fā)動(dòng)機(jī)。該發(fā)動(dòng)機(jī)是通過爆轟方式，其熱力循環(huán)效率比常壓燃燒發(fā)動(dòng)機(jī)要高20%～30%。與脈沖爆轟發(fā)動(dòng)機(jī)相比，連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)具有以下優(yōu)點(diǎn)：僅需進(jìn)行一次點(diǎn)火起爆，在燃燒室內(nèi)形成穩(wěn)定旋轉(zhuǎn)傳播的爆轟波，旋轉(zhuǎn)頻率可高達(dá)幾千赫茲；可連續(xù)性地填充高速燃料，獲得較大的穩(wěn)定推力。CRDE是極具吸引力的新型推進(jìn)系統(tǒng)。

實(shí)驗(yàn)研究方面，國外Bykovskii等[1-4]利用不同的燃料(如氫氣、煤油、丙酮等)，在不同尺寸及不同形狀的燃燒室進(jìn)行爆轟實(shí)驗(yàn)，發(fā)現(xiàn)旋轉(zhuǎn)爆轟的形成與持續(xù)性都受到預(yù)混燃料進(jìn)氣壓力、燃燒室形狀尺寸、周圍的環(huán)境等因素的影響。實(shí)驗(yàn)中，氣體燃料與氧氣或空氣都能實(shí)現(xiàn)連續(xù)爆轟，但使用液體燃料和空氣組合無法成功起爆，通過加入氧氣提高其活性后，才能實(shí)現(xiàn)連續(xù)起爆。國內(nèi)林偉、周進(jìn)等[5-6]利用H2和空氣進(jìn)行了不同模態(tài)下的連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)推力實(shí)驗(yàn)研究。

理論研究方面，近年來，國內(nèi)外從不同角度對(duì)CRDE進(jìn)行了數(shù)值模擬，Zhdan等[7]對(duì)氫/氧預(yù)混燃料旋轉(zhuǎn)爆轟進(jìn)行二維數(shù)值模擬，揭示出爆轟波的傳播機(jī)理；Schwer等[8]忽略徑向作用，采用二維模型來模擬旋轉(zhuǎn)爆轟，研究入口壓力和出口背壓的壓比對(duì)爆轟流場(chǎng)和爆轟性能的影響；Davidenko等[9]分析了入口壓力(1～6 MPa)和燃料入口有效面積比對(duì)氫/氧連續(xù)旋轉(zhuǎn)爆轟特性和燃燒室比沖等影響。范寶春等[10-11]對(duì)連續(xù)旋轉(zhuǎn)爆轟進(jìn)行了二維和三維的數(shù)值模擬，研究了爆轟波的結(jié)構(gòu)及連續(xù)旋轉(zhuǎn)傳播的自持機(jī)理；邵業(yè)濤等[12-13]對(duì)連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)流場(chǎng)進(jìn)行了二維數(shù)值模擬，得到與試驗(yàn)匹配較好的流場(chǎng)結(jié)構(gòu)，在三維數(shù)值計(jì)算中，分析了不同進(jìn)氣總壓條件下爆轟推進(jìn)性能；劉世杰等[14]對(duì)連續(xù)旋轉(zhuǎn)爆轟波細(xì)致結(jié)構(gòu)和自持機(jī)理進(jìn)行了數(shù)值研究；馬虎等[15]以氫氣為燃料，研究了壓力條件對(duì)CRDE爆轟性能影響。

目前，國內(nèi)外有關(guān)液體燃料連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)的數(shù)值模擬研究較少，本文對(duì)汽油為燃料、富氧空氣為氧化劑的氣液兩相連旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)進(jìn)行二維數(shù)值模擬，采用守恒元和求解元方法(the Method of Conservation Element and Solution Element，CE/SE方法)進(jìn)行求解，分析了氣液兩相連續(xù)旋轉(zhuǎn)爆轟穩(wěn)定狀態(tài)時(shí)燃燒室內(nèi)流場(chǎng)結(jié)構(gòu)，研究填充總壓對(duì)氣液兩相CRDE爆轟過程以及出口參數(shù)性能的影響，為以后氣液兩相連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)研究提供理論指導(dǎo)。

1　理論模型

實(shí)際上連續(xù)旋轉(zhuǎn)爆轟過程為三維爆轟過程，為了簡(jiǎn)化問題，只考慮沿環(huán)形燃燒室的二維爆轟過程。圖1(a)是三維環(huán)形燃燒室物理模型，下端為入口，上端為出口；假設(shè)該圓環(huán)沒有厚度，將圓環(huán)沿著母線ab剪開，得到矩形計(jì)算域圖1(b)，左邊界ab和右邊界a′b′是通過周期邊界(兩邊物理參數(shù)相同)相連，下端為進(jìn)氣端面，上端為出口端面。

(a)物理域

(b)計(jì)算域

旋轉(zhuǎn)爆轟過程采用氣液兩相流的爆轟控制方程描述[16-17]：

(1)

燃料液滴剝離和蒸發(fā)對(duì)氣相質(zhì)量的貢獻(xiàn)率m21計(jì)算為[18]

(2)

其中，r為燃料液滴半徑。燃料液滴半徑變化率由氣動(dòng)剝離與蒸發(fā)(等式右邊第一部分為剝離，第二部分為蒸發(fā))2部分組成[19]：

假設(shè)汽油為辛烷單一組分，辛烷的一步化學(xué)反應(yīng)方程式為

C8H18+12.5O2→8CO2+9H2O

(4)

(5)

2　計(jì)算方法及初始條件和邊界條件

2.1計(jì)算方法

本文采用的CE/SE方法[16，20]是求解強(qiáng)間斷問題的一種新的計(jì)算方法。其計(jì)算格式簡(jiǎn)單、精度高、捕捉爆轟波等強(qiáng)間斷能力強(qiáng)，在求解爆轟等強(qiáng)間斷面物理問題方面已有很多成功算例[17-18]。連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)內(nèi)旋轉(zhuǎn)爆轟波為強(qiáng)間斷，利用CE/SE方法計(jì)算有獨(dú)特的優(yōu)勢(shì)，其詳細(xì)計(jì)算格式見文獻(xiàn)[16]。

2.2初始條件和邊界條件

在數(shù)值模擬過程中，周向長(zhǎng)度L=0.3 m、軸向長(zhǎng)度H=0.1 m的矩形區(qū)域.

由于計(jì)算條件的限制，分別采用240×80，300×100和400×140網(wǎng)格數(shù)進(jìn)行計(jì)算，比較3種網(wǎng)格對(duì)爆轟波間斷面捕捉的情況，計(jì)算結(jié)果見圖2。

圖2　不同網(wǎng)格數(shù)條件下入口處壓力分布情況

結(jié)果表明，采用300×100網(wǎng)格能捕捉到爆轟波間斷面，進(jìn)一步增加網(wǎng)格數(shù)，對(duì)計(jì)算結(jié)果的影響不大，可認(rèn)為300×100網(wǎng)格滿足計(jì)算精度要求。因此，本文采用300×100網(wǎng)格。

初始條件：圖1(b)左下角的區(qū)域1為點(diǎn)火區(qū)域，大小為20×30個(gè)網(wǎng)格數(shù)，點(diǎn)火條件為高溫、高壓及高速氣相周向初速度(為了快速形成向單一方向穩(wěn)定傳播的連續(xù)旋轉(zhuǎn)爆轟波)；區(qū)域2為預(yù)混新鮮燃料層；3為富氧空氣(空氣中增加氧的含量)。

邊界條件：下端為入口邊界，與進(jìn)氣管道相連，總壓為p0，總溫為T0。假設(shè)氣相、液相是以相同的速度進(jìn)入燃燒室，設(shè)燃燒室內(nèi)邊界臨近處的計(jì)算壓力為p，進(jìn)氣邊界分3種情況：

(1)當(dāng)p≥p0，此時(shí)預(yù)混燃料不能進(jìn)入燃燒室內(nèi)，將入口邊界按照固壁邊界處理；

(2)當(dāng)pcr

(3)當(dāng)p

式中R為氣體常數(shù)。

計(jì)算域的上邊aa為出口邊界，使用無反射自由邊界條件；當(dāng)出口為超聲速時(shí)，出口邊界條件根據(jù)二階外推得到；當(dāng)出口為亞聲速時(shí)，出口壓力等于環(huán)境壓力(0.1 MPa)。左右邊界即為周期邊界。

3　計(jì)算結(jié)果與分析

3.1旋轉(zhuǎn)爆轟波及其流場(chǎng)分析

氣液兩相預(yù)混燃料填充總壓0.4 MPa，總溫288.15 K，其中富氧空氣中含氧質(zhì)量分?jǐn)?shù)30%，液相中液滴半徑統(tǒng)一取R=50 μm(實(shí)驗(yàn)中測(cè)量的汽油霧化后液滴半徑分布為40～60 μm，文獻(xiàn)[17，19，20]均采用液滴半徑為50 μm進(jìn)行數(shù)值計(jì)算)。

圖3為計(jì)算域內(nèi)(x=0.2 mm，y=0.004 mm)點(diǎn)處的壓力和溫度隨時(shí)間的變化曲線。圖3中顯示，從點(diǎn)火之后，一共循環(huán)了7個(gè)周期，前兩個(gè)周期內(nèi)，該點(diǎn)爆轟壓力和溫度波動(dòng)較大，尚未形成穩(wěn)定爆轟，從第3個(gè)周期到第7個(gè)周期內(nèi)，爆轟壓力和溫度平穩(wěn)，表明達(dá)到穩(wěn)定爆轟狀態(tài)。從589.85～1 427.51 μs這段時(shí)間內(nèi)，該點(diǎn)處壓力和溫度隨著爆轟波傳播呈周期性變化，爆轟波波陣面處高壓間斷面與高溫間斷面相互耦合在一起，突顯出了爆轟波的基本特征。通過計(jì)算得到爆轟波壓力峰值、溫度峰值以及傳播速度的平均值分別為3.97 MPa、2 407.6 K、1 432.5 m/s，得到對(duì)應(yīng)的連續(xù)旋轉(zhuǎn)爆轟頻率為4 775 Hz。

圖3　(x=0.2 mm，y=0.004 mm)處壓力和溫度隨時(shí)間變化曲線

圖4為t=1 011 μs時(shí)刻燃燒室內(nèi)流場(chǎng)分布圖。其中，圖4(a)的1是連續(xù)旋轉(zhuǎn)爆轟波，2是爆轟波后的爆轟燃燒產(chǎn)物，3是連續(xù)旋轉(zhuǎn)爆轟的斜激波，4是斜激波和爆轟產(chǎn)物接觸間斷面，5是新鮮的預(yù)混燃料，6是爆轟燃燒產(chǎn)物與新鮮預(yù)混燃料的接觸面；與Bykovkii[4]利用煤油(液態(tài))為燃料，在環(huán)形燃燒室內(nèi)進(jìn)行爆轟實(shí)驗(yàn)中所揭示的旋轉(zhuǎn)爆轟流場(chǎng)定性一致。

(a)溫度

(b)壓力

(c)氣相周向速度

(d)氣相軸向速度

從圖4(a)中可看出，燃燒室內(nèi)最高溫度可高達(dá)2 550 K，出現(xiàn)的位置在爆轟波、斜激波和爆轟產(chǎn)物與新鮮預(yù)混燃料接觸面的交匯處；由于爆轟產(chǎn)物存在較高溫度，使得新鮮預(yù)混燃料層與爆轟產(chǎn)物的接觸面處發(fā)生了緩燃，沒有能充分利用于爆轟，對(duì)氣液兩相連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)爆轟能量有一定損失。圖4(b)為壓力流場(chǎng)分布云圖，在入口x=0.2 m處，爆轟波波陣面處壓力峰值高達(dá)3.97 MPa；而斜激波出口處的前端壓力值最低，只有0.2 MPa，由于爆轟產(chǎn)物通過斜激波向出口膨脹，距離斜激波越遠(yuǎn)，爆轟產(chǎn)物膨脹越充分。圖4(c)、(d)分別是氣流的周向和軸向速度分布云圖；最大周向氣流速度為650 m/s，相對(duì)以氫氣為燃料的連續(xù)旋轉(zhuǎn)爆轟的周向氣流速度(1 000 m/s左右)小很多，主要是由于燃料的活性、氣液兩相之間作用力等因素的影響；軸向速度最大值985 m/s，出現(xiàn)在斜激波出口處，表明斜激波可促進(jìn)爆轟產(chǎn)物的膨脹。

圖5為t=1 011 μs時(shí)刻液滴半徑在入口處x方向上分布曲線。從圖5可看出，此時(shí)連續(xù)旋轉(zhuǎn)爆轟波正好傳播到x=0.2 m處，爆轟波波前液滴半徑R=50 μm，為未反應(yīng)的預(yù)混燃料；爆轟波掃過之后，液滴在蒸發(fā)和剝離的作用下迅速變?yōu)槠驼羝麉⒓踊瘜W(xué)反應(yīng)，在x=0.136 m處，液滴半徑變?yōu)?。在x=0.102 m處，開始有預(yù)混燃料開始進(jìn)入燃燒室，與高溫的爆轟產(chǎn)物接觸時(shí)，在接觸面處，液滴燃料會(huì)產(chǎn)生蒸發(fā)，同時(shí)在接觸面處發(fā)生燃燒，使得接觸面處的參數(shù)發(fā)生波動(dòng)。從圖5可看出，在出口x=0.008 8 m和x=0.009 4 m位置，液滴半徑發(fā)生明顯的波動(dòng)。

圖5　t=1 011 μs時(shí)入口處液滴半徑在x方向上的變化曲線

3.2不同進(jìn)氣壓力對(duì)連續(xù)旋轉(zhuǎn)爆轟參數(shù)的影響

在其他條件保持不變的情況下，研究不同進(jìn)氣總壓力0.25、0.3、0.4、0.5、0.6 MPa條件對(duì)氣液兩相連續(xù)旋轉(zhuǎn)爆轟特性的影響。爆轟達(dá)到穩(wěn)定狀態(tài)時(shí)，取爆轟波傳播到周向上x=0.2 m處各參數(shù)，研究爆轟變化規(guī)律。

圖6為不同進(jìn)氣總壓條件下，爆轟波的壓力峰值和溫度峰值變化情況。爆轟波壓力峰值隨著進(jìn)氣總壓的增加而增大，增大趨勢(shì)較明顯；但爆轟波溫度峰值隨著進(jìn)氣總壓的增加而有所減少，從2 300 K降低到2 100 K，變化趨勢(shì)較小，保持相對(duì)平緩。

表1為不同進(jìn)氣總壓時(shí)，燃燒室入口處的相關(guān)參數(shù)。從表1可看出，入口處的平均壓力和密度隨進(jìn)氣總壓的增加而呈線性增加；但入口處的平均溫度則隨進(jìn)氣總壓的增加而降低。入口處的軸向速度是以燃料的進(jìn)氣速度占主導(dǎo)部分，只有小部分是由于爆轟波波后爆轟產(chǎn)物側(cè)向膨脹引起與進(jìn)氣速度相反的氣流，通過計(jì)算得到入口處的平均軸向速度隨著進(jìn)氣總壓增加而略有增加。入口處的周向速度由3部分組成的：(1)為正在爆轟的周向速度；(2)為爆轟產(chǎn)物的側(cè)向膨脹引起的周向速度；(3)為新鮮預(yù)混燃料層內(nèi)的周向速度。通過計(jì)算得到入口處的周向平均速度方向朝x的反方向，大小隨進(jìn)氣壓力的增加而略有減小，說明爆轟產(chǎn)物側(cè)向膨脹引起的周向速度占主導(dǎo)地位。進(jìn)氣總壓對(duì)爆轟波速度幾乎沒有影響。

圖6　不同進(jìn)氣總壓條件下爆轟波的壓力峰值和溫度峰值變化情況

進(jìn)氣壓力/MPa入口平均壓力/MPa入口平均溫度/K入口平均密度/(kg/m3)入口平均軸向速度/(m/s)入口平均周向速度/(m/s)爆轟波波速/(m/s)0.25——————0.300.5912102.1987.2-62.91452.20.400.7211102.9889.3-41.81431.70.500.8410103.8092.0-29.41421.80.600.989644.5695.4-28.21419.1

當(dāng)進(jìn)氣總壓為0.25 MPa時(shí)，不同時(shí)刻燃燒室溫度分布云圖如圖7所示，從圖7(a)、(b)可看出，點(diǎn)火之后，t=294.9 μs到t=323.8 μs爆轟波前的新鮮預(yù)混燃料層區(qū)域逐漸變小。這是由于新鮮預(yù)混燃料層右端的爆轟產(chǎn)物側(cè)向膨脹不夠充分，在入口處的壓力仍大于進(jìn)氣總壓，導(dǎo)致新鮮預(yù)混燃料不能順利進(jìn)入燃燒室，沒有新鮮預(yù)混燃料支持爆轟波繼續(xù)傳播；圖7(c)爆轟燃燒已經(jīng)開始衰減直到消失，使得整個(gè)燃燒室入口處的壓力都小于進(jìn)氣總壓；圖7(d)所示，整個(gè)入口處都有低溫預(yù)混燃料進(jìn)入燃燒室，導(dǎo)致整個(gè)爆轟過程熄滅。

(a)t=294.9 μs

(b)t=323.8 μs

(c)t=411.0 μs

(d)t=503.2 μs

圖8為不同進(jìn)氣總壓爆轟波傳播到x=0.2 m處，液滴半徑在入口處x方向上分布情況。從圖8可看出，液滴被爆轟波掃過之后發(fā)生蒸發(fā)和剝離，液滴半徑迅速減??；隨進(jìn)氣總壓增大，入口處開始填充燃料的位置距離爆轟波波陣面的位置越近，同時(shí)周向方向可填充新鮮燃料區(qū)域越長(zhǎng)，使得燃料填充比(在入口處周向方向上可填充燃料的長(zhǎng)度與周向總長(zhǎng)度的比)增加。當(dāng)進(jìn)氣總壓為0.25 MPa時(shí)，液滴半徑仍為50 μm，表明液滴在燃燒室內(nèi)并沒發(fā)生蒸發(fā)和剝離；由圖7可知，進(jìn)氣總壓為0.25 MPa，在t=294.9 μs之后，爆轟波前的燃料填充區(qū)域變得越來越小，無法支持爆轟波繼續(xù)傳播，表明進(jìn)氣總壓過小，無法成功形成連續(xù)旋轉(zhuǎn)爆轟波。

圖8　在x方向上入口處的液滴半徑的分布情況

3.3不同進(jìn)氣總壓對(duì)出口流場(chǎng)參數(shù)及其性能影響

衡量發(fā)動(dòng)機(jī)性能的重要指標(biāo)為單位質(zhì)量的推進(jìn)劑比沖，計(jì)算公式如下：

(4)

式中ρ為密度；v為軸向速度；pb為出口環(huán)境壓力；g為重力加速度，取9.8 m/s2；wf為燃料的質(zhì)量分?jǐn)?shù)。

燃料比沖與出口燃?xì)鈮毫?、環(huán)境壓力、出口燃?xì)饷芏纫约俺隹谔幍妮S向燃?xì)馑俣扔嘘P(guān)，從而研究不同進(jìn)氣總壓對(duì)出口處相關(guān)參數(shù)影響。圖9為在不同進(jìn)氣壓總壓條件下，出口處壓力、密度以及軸向速度的變化曲線。由圖9(a)可看出，出口處的壓力參數(shù)整體隨著進(jìn)氣總壓的增加而增加，壓力存在較大的波動(dòng)，這些波動(dòng)主要是斜激波傳到出口處所引起的；隨進(jìn)氣壓力的增加，出口壓力峰值在周向方向出現(xiàn)的位置會(huì)滯后，且出口壓力均高于環(huán)境壓力，由出口邊界條件可得到出口處均為超聲速氣流。圖9(b)為出口密度變化曲線，密度的變化規(guī)律與壓力變化規(guī)律類似，由于爆轟產(chǎn)物在斜激波和側(cè)向膨脹波的作用下，不斷向出口處膨脹，達(dá)到出口處時(shí)密度偏??；當(dāng)進(jìn)氣總壓力為0.6 MPa時(shí)，出口密度峰值為1.5 kg/m3，而進(jìn)氣總壓力為0.3 MPa時(shí)，此時(shí)出口密度峰值只有0.75 kg/m3，出口處的密度分布情況受進(jìn)氣總壓影響較大。圖9(c)為不同進(jìn)氣總壓出口處軸向速度在周向方向的分布曲線，出口軸向速度受進(jìn)氣總壓影響不大，軸向速度峰值范圍950～1 050 m/s，在斜激波出口位置的前后，都出現(xiàn)較明顯的波動(dòng)。計(jì)算得到基于燃料的比沖與進(jìn)氣總壓力的關(guān)系如圖10所示。

(a)壓力

(b)密度

(c)軸向速度

圖10　發(fā)動(dòng)機(jī)燃料比沖與進(jìn)氣總壓關(guān)系曲線

在其他條件不變的情況下，進(jìn)氣壓力越高，發(fā)動(dòng)機(jī)基于燃料的比沖越大，與進(jìn)氣壓力幾乎成線性關(guān)系。

4　結(jié)論

利用守恒元和求解元方法，計(jì)算CRDE氣液兩相爆轟流場(chǎng)。從計(jì)算結(jié)果可看出，守恒元和求解元方法能夠有效地捕捉到爆轟波等強(qiáng)間斷面，并獲得了穩(wěn)定爆轟狀態(tài)時(shí)燃燒室的流場(chǎng)結(jié)構(gòu)，計(jì)算結(jié)果與文獻(xiàn)中結(jié)果定性的一致。

(1)分析了液滴半徑在入口處周向方向的變化規(guī)律，得出爆轟波前填充的燃料層是形成穩(wěn)定連續(xù)旋轉(zhuǎn)爆轟的關(guān)鍵。當(dāng)進(jìn)氣壓力過小，在入口處不能形成有效的燃料填充比，致使爆轟熄滅。

(2)當(dāng)出口環(huán)境壓力不變時(shí)，爆轟波壓力、溫度受進(jìn)氣總壓力的影響，出口處壓力和密度隨著進(jìn)氣總壓力增加而增大，進(jìn)氣總壓力對(duì)出口處的周向速度以及軸向速度影響很小。

(3)在進(jìn)氣總壓力0.3～0.6 MPa范圍內(nèi)，CRDE達(dá)到穩(wěn)定爆轟狀態(tài)時(shí)，出口處氣流速度均為超音速氣流；當(dāng)出口環(huán)境壓力不變時(shí)，CRDE的比沖與進(jìn)氣總壓力成正比例關(guān)系。

[1]Bykovskii F A,Zhdan S A and Adernikov E F.Continuous spin detonation in annular combustion[J].Combustion,Explosion and Shock Waves,2005,41(4):449-459.

[2]Bykovskii F A,Zhdan S A and Adernikov E F.Continuous spin detonation[J].Journal of Propulsion and Power,2006,22(6):1204-1216.

[3]Bykovskii F A,Zhdan S A and Adernikov E F.Continuous spin detonation of fuel-air mixtures[J].Combustion,Explosion,and Shock Waves,2006,42(4):463-471.

[4]Bykoviskii F A,Zhdan S A and Adernikov E F.Continuous spin detonation of hydrogen-oxygen mixture.1.Annular cylindrical combustor[J].Combustion,Explosion,and Shock Waves,2008,44(2):150-162.

[5]林偉,周進(jìn),林志勇,等.H2/Air連續(xù)旋轉(zhuǎn)爆震發(fā)動(dòng)機(jī)推力測(cè)試(Ⅰ):單波模態(tài)下的推力[J].推進(jìn)技術(shù),2015,36(4):497-503.

[6]林偉,周進(jìn),林志勇,等.H2/Air連續(xù)旋轉(zhuǎn)爆震發(fā)動(dòng)機(jī)推力測(cè)試(Ⅱ):雙波模態(tài)下的推力[J].推進(jìn)技術(shù),2015,36(5):641-649.

[7]Zhdan S A,Bykovskii F A,Vedernikov E F.Mathematical modeling of a rotating detonation wave in a hydrogen-oxygen mixture[J].Combustion,Explosion and Shock Waves,2007,43(4):449-459.[8]Schwer D,Kailasanath K.Numerical investigation of the physics of rotating detonation engines[J].Proceedings of the Combustion Institute.2011,33(2):2195-2202.

[9]Davidenko D M,Gokalp I and Kudryavtsev A N.Numerical study of the continuous detonation wave rocket engine[R].AIAA 2008-2680.

[10]張旭東,范寶春,潘振華,等.旋轉(zhuǎn)爆轟自持機(jī)理的數(shù)值研究[J].彈道學(xué)報(bào),2011,23(1):1-4.

[11]歸明月,范寶春,張旭東,等.旋轉(zhuǎn)爆轟的三維數(shù)值模擬[J].推進(jìn)技術(shù),2010,31(1):82-86.

[12]邵業(yè)濤,王健平.連續(xù)爆轟發(fā)動(dòng)機(jī)的二維數(shù)值模擬研究[J].航空動(dòng)力學(xué)報(bào),2009,24(5):980-986.

[13]邵業(yè)濤,王健平,唐新猛,等.連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)流場(chǎng)三維數(shù)值模擬[J].航空動(dòng)力學(xué)報(bào),2010,25(8):1717-1722.

[14]劉世杰,譚慧,林志勇,等.連續(xù)旋轉(zhuǎn)爆震波細(xì)致結(jié)構(gòu)及自持機(jī)理[J].推進(jìn)技術(shù),2011,32(3):431-436.

[15]馬虎,封峰,武曉松,等.壓力條件對(duì)旋轉(zhuǎn)爆震發(fā)動(dòng)機(jī)的影響[J].彈道學(xué)報(bào),2012,24(4):94-98.

[16]翁春生,王浩.計(jì)算內(nèi)彈道學(xué)[M].北京:國防工業(yè)出版社,2006.

[17]馬丹花,翁春生.爆震管內(nèi)擾流片對(duì)爆震波影響的數(shù)值分析[J].推進(jìn)技術(shù),2011,32(3):425-430.

[18]洪滔,秦承森.氣體-燃料液滴兩相系統(tǒng)爆轟的數(shù)值模擬[J].爆炸與沖擊,1999,19(4):335-342.

[19]彭振,翁春生.脈沖爆轟發(fā)動(dòng)機(jī)中等離子體點(diǎn)火的數(shù)值計(jì)算[J].工程力學(xué),2012,29(5):242-250.

[20]Chang S C.A new approach for constructing highly stable high order CE/SE schemes[R].AIAA 2010-543.

(編輯：崔賢彬)

Two-dimensional numerical simulation of the inlet stagnation pressure influence on the continuous rotating detonation engine

LI Bao-xing，WENG Chun-sheng

(National Key Lab of Transient Physics，Nanjing University of Science and Technology，Nanjing210094，China)

To study the detonation characteristics of two-phase continuous rotating detonation engine(CRDE),based on the method of two-dimensional Conservation Element and Solution Element,the detonation process of gasoline and oxygen-enriched air engine was simulated,the inner flow field of combustion chamber and the propagation characteristic of the detonation wave were obtained, and the detonation performance was calculated at different inlet stagnation pressure. The results show that the inlet flow field changes periodically with the propagation of detonation wave,and the premixed fuel layer of the front of detonation was the key to the formation of stable continuous rotating detonation,the effective premixed fuel layer could not be formed at entrance after ignition without enough inlet stagnation pressure,resulting in the extinguish of detonation.The pressure,temperature,density at the inlet and exit were related to the inlet stagnation pressure,and the engine specific impulse was proportional to the inlet stagnation pressure.

gas-liquid two-phase；continuous rotating detonation engine；inner flow field；premixed fuel layer

2015-08-31；

2016-02-01。

國家自然科學(xué)基金(11472138)；中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)基金(30920130112007)；國防預(yù)研基金(9140c300202120c30)。

李寶星(1990—)，男，博士生，研究方向?yàn)楸Z推進(jìn)技術(shù)。E-mail：bestlibaoxing@163.com

翁春生(1964—)，男，教授/博導(dǎo)，研究方向?yàn)橥七M(jìn)技術(shù)。E-mail：wengcs@126.com

V435

A

1006-2793(2016)05-0612-07

10.7673/j.issn.1006-2793.2016.05.003