雙脈沖發(fā)動機點火過程三維數(shù)值模擬

嚴登超，陳 雄，李映坤，朱 亮，孫 姍

(1.南京理工大學(xué) 機械工程學(xué)院，南京 210094；2.江南機電設(shè)計研究所，貴陽 550009)

2016-10-10；

2016-12-20。

國家自然科學(xué)基金(51606098)；中央高校基本科研業(yè)務(wù)費專項資金資助(30915118805)。

嚴登超(1989—)，男，碩士，研究領(lǐng)域為雙脈沖發(fā)動機技術(shù)。E-mailydc203@163.com

雙脈沖發(fā)動機點火過程三維數(shù)值模擬

嚴登超1,2，陳 雄1，李映坤1，朱 亮1，孫 姍1

(1.南京理工大學(xué) 機械工程學(xué)院，南京 210094；2.江南機電設(shè)計研究所，貴陽 550009)

以噴射棒式雙脈沖發(fā)動機燃燒室、級間隔離裝置和噴管一體化為研究對象，采用數(shù)值仿真技術(shù)對Ⅱ脈沖點火過程三維流場特性進行分析研究。計算結(jié)果表明，點火初期燃氣壓力波峰超前于火焰峰到達級間隔離裝置，并以壓強沖擊波形式傳播，Ⅱ脈沖燃燒室相對高壓區(qū)位置不斷發(fā)生改變；級間孔打開過程對藥柱末端壓強影響較大，但對Ⅱ脈沖燃燒室壓強整體上升過程影響較小；級間孔打開后，燃氣經(jīng)級間孔加速后形成高度欠膨脹射流，并在Ⅰ脈沖燃燒室內(nèi)形成非對稱帶狀低壓區(qū)；級間孔分布的非對稱性，導(dǎo)致壓強及溫度在發(fā)動機燃燒室中呈現(xiàn)顯著的三維分布特性；高溫區(qū)出現(xiàn)在隔板附近，而在裝藥前端、裝藥末端及外圍級間孔軸線附近出現(xiàn)低溫區(qū)。

雙脈沖發(fā)動機；三維流場；數(shù)值模擬；點火過程

0 引言

雙脈沖發(fā)動機主要由兩個燃燒室、兩個點火具、級間隔離裝置(以下簡稱隔板)及噴管組件組成。當Ⅰ脈沖工作時，隔板承受Ⅰ脈沖燃燒室高溫高壓燃氣的作用，保證Ⅱ脈沖燃燒室不受影響。Ⅱ脈沖工作后，隔板可靠打開，Ⅱ脈沖工作產(chǎn)生的燃氣穿過隔板進入Ⅰ脈沖燃燒室，從噴管噴出從而產(chǎn)生推力。實現(xiàn)固體火箭發(fā)動機的多次關(guān)機、啟動，滿足現(xiàn)代戰(zhàn)爭對導(dǎo)彈高機動性、高可控性、大末端速度及遠射程等指標要求[1]。

雙脈沖發(fā)動機作為現(xiàn)有飛行系統(tǒng)的一種先進動力裝置，美國、日本、加拿大等國家從20世紀60年代開始，對雙脈沖固體火箭發(fā)動機技術(shù)進行了全面系統(tǒng)的研究[2]。其中，Naumann[3]、Stadler[4]、Schilling[5]、Harold[6]、Javed[7]等分別針對雙脈沖發(fā)動機的設(shè)計、制造、飛行實驗及內(nèi)流場特性進行了廣泛的研究。國內(nèi)針對雙脈沖發(fā)動機的研究雖然起步較晚，但仍取得了較大的進展，北京航空航天大學(xué)的王長輝[8]、西北工業(yè)大學(xué)的王偉、王春光[9-10]、中國航天科技集團第四研究院的劉雨等[11]對雙脈沖發(fā)動機陶瓷隔艙式、金屬膜片式及軟隔層式級間隔離裝置分別進行了有限元結(jié)構(gòu)強度分析、理論和試驗驗證；上海航天動力技術(shù)研究所的孫娜[12]、哈爾濱工程大學(xué)的朱衛(wèi)兵[13]、北京理工大學(xué)的楊春慶[14]、南京理工大學(xué)的李映坤等[15]針對雙脈沖發(fā)動機內(nèi)流場特征、級間隔離裝置對流場參數(shù)影響、Ⅱ脈沖點火延遲影響因素及燃氣對絕熱層燒蝕規(guī)律等進行了數(shù)值模擬研究。

目前，針對雙脈沖發(fā)動機數(shù)值模擬研究主要集中在Ⅱ脈沖點火延遲和穩(wěn)態(tài)流場特性的研究。然而，對于噴射棒式雙脈沖發(fā)動機，隔板上開有流道截面積劇烈變化的孔(以下簡稱級間孔)，導(dǎo)致在Ⅱ脈沖點火升壓過程中隔板前后燃氣流動紊亂，并具有顯著的三維流動現(xiàn)象，對此現(xiàn)象仍沒有得到足夠深入研究。

針對上述問題，本文以噴射棒式雙脈沖發(fā)動機為研究對象，采用FLUENT商業(yè)軟件，進行發(fā)動機內(nèi)流場仿真計算。分析Ⅱ脈沖點火過程裝藥內(nèi)外表面升壓過程和級間孔分布對燃氣流動的影響。計算結(jié)果揭示了噴射棒式雙脈沖發(fā)動機Ⅱ脈沖點火過程的流場特性，對后續(xù)雙脈沖發(fā)動機點火過程的認識及設(shè)計具有一定的借鑒意義。

1 物理模型及計算方法

1.1 物理模型

物理模型是典型的噴射棒式雙脈沖發(fā)動機[16]，考慮到本文著重討論Ⅰ脈沖工作結(jié)束后，Ⅱ脈沖點火工作過程，在確保計算精確的條件下，對發(fā)動機做適當?shù)暮喕鐖D1所示。

Ⅰ脈沖燃燒室長度為56 mm，Ⅱ脈沖燃燒室長度為125 mm，Ⅰ、Ⅱ脈沖燃燒室內(nèi)徑均為85 mm，如圖1(b)所示；隔板上開有9個級間孔，級間孔均為兩臺階形式，小孔臺階長度為2 mm，大孔臺階長度為4 mm，如圖1(a)所示。將點火具簡化為等直徑圓管，長度L=7 mm，直徑D=40 mm，Ⅱ脈沖燃燒室推進劑采用兩端包覆、自由裝填的改性雙基推進劑。

1.2 計算模型及網(wǎng)格劃分

針對雙脈沖發(fā)動機及隔板結(jié)構(gòu)特點，為節(jié)省計算資源，選取模型的1/4作為計算域進行計算。采用Gambit專業(yè)網(wǎng)格劃分軟件，對計算域不同部分進行網(wǎng)格劃分，如圖2所示。其中，級間孔區(qū)域使用非結(jié)構(gòu)網(wǎng)格，其他區(qū)域使用結(jié)構(gòu)網(wǎng)格，通過Interface面將結(jié)構(gòu)網(wǎng)格與非結(jié)構(gòu)網(wǎng)格進行對接。為了獲得流場精確的計算結(jié)果，在點火具出口、級間孔區(qū)域、裝藥加質(zhì)面及噴管喉部等局部區(qū)域進行網(wǎng)格加密處理，網(wǎng)格總數(shù)為679 840。

為了更好地顯示流場計算結(jié)果，在計算域中作相關(guān)輔助面，如圖3(a)所示，Ⅱ脈沖燃燒室前封頭位于X軸坐標為0，各切面位于X軸坐標如表1所示。為得到點火升壓過程燃燒室內(nèi)壓強變化規(guī)律，在裝藥內(nèi)外壁面的前端、中端、末端、Ⅰ脈沖燃燒室及噴管收斂段等位置設(shè)置監(jiān)測點，如圖3(b)所示。各監(jiān)測點X軸坐標如表2所示。

表1 各切面位于X軸坐標

表2 P0切面內(nèi)各監(jiān)測點X坐標

根據(jù)雙脈沖發(fā)動機工作條件，在物理模型中采取如下假設(shè)：

(1)點火具生成燃氣與推進劑燃燒產(chǎn)物性質(zhì)一致，均為純氣相流且服從理想氣體狀態(tài)方程，推進劑及燃氣參數(shù)如表3所示。

表3 推進劑及燃氣參數(shù)

(1)

(2)

(3)

其中

(3)考慮到噴射棒噴射時間極短，假設(shè)噴射棒瞬間從級間孔噴出，忽略噴射棒在Ⅰ脈沖燃燒室中運動對內(nèi)流場的影響。

1.3 控制方程

本文數(shù)值模型計算的流動控制方程為三維非定常可壓縮N-S方程[17]

(4)

式中φ為某一變量；Γφ為擴散系數(shù)；Sφ為源項。

源項如下：

連續(xù)方程：

(5)

動量方程：

(6)

能量方程：

(7)

式中ρ、ui、p、E和H分別為密度、速度分量、壓強、總能和總焓；Sm為推進劑燃燒生成氣體的質(zhì)量加質(zhì)源項；Sfi為推進劑燃燒生成氣體沿i方向的動量加質(zhì)源項；Sh為推進劑燃燒生成氣體的能量加質(zhì)源項；其他相關(guān)參數(shù)及物理意義見參考文獻[18]。

考慮點火過程發(fā)動機內(nèi)流場呈現(xiàn)湍流特性，數(shù)值計算中湍流模型選擇RNGk-ε兩方程模型。

1.4 邊界條件及初始條件

邊界條件如下：

(1)采用模擬自由容積法得到點火具出口質(zhì)量流率隨時間的變化關(guān)系[19](如圖4所示)，質(zhì)量流率隨時間的變化關(guān)系轉(zhuǎn)換成質(zhì)量通量的形式，通過用戶自定義函數(shù)寫入FLUENT中，作為求解器的質(zhì)量通量邊界條件。

(2)計算域各壁面均為無滑移絕熱壁面，參數(shù)法向梯度均為零，采用標準壁面函數(shù)處理；級間孔入口達到打開壓強前，滿足壁面條件，達到打開壓強后，壁面設(shè)置為內(nèi)部流場區(qū)域。

(3)推進劑裝藥加質(zhì)面邊界條件為當加質(zhì)網(wǎng)格單元未達到點火溫度時，裝藥燃面作為絕熱固壁面邊界處理；當加質(zhì)網(wǎng)格單元溫度達到點火溫度時，加質(zhì)層網(wǎng)格轉(zhuǎn)變?yōu)榧淤|(zhì)源項，向流場中加入動量、質(zhì)量及能量。

初始條件：初始流場壓強、溫度、密度及推進劑表面溫度與環(huán)境空氣相同，燃氣初始流速為零，即u=v=w=0，T0=300 K，p0=101 325 Pa。

2 計算結(jié)果與分析

2.1 級間孔打開前Ⅱ脈沖內(nèi)流場分析

2.1.1 Ⅱ脈沖燃燒室燃氣填充過程

圖5給出了不同時刻燃氣在Ⅱ脈沖燃燒室p0截面中填充過程壓力和溫度云圖分布情況。

點火具工作后，從點火具噴出的高溫燃氣在下游形成高度欠膨脹射流，填充Ⅱ脈沖燃燒室自由容積，如圖5(a)所示。對比圖5(b)中的壓力、溫度云圖可知，燃氣壓力波峰明顯超前于火焰峰，在火焰峰前面，燃燒室空氣受壓力波峰擾動而壓強升高，溫度卻未升高，此時點火具出口燃氣壓力與溫度相對于流場其它區(qū)域最高。在0.15 ms時，裝藥內(nèi)孔燃面前端首次被引燃，引燃區(qū)開始向流場加質(zhì)，此時流場高壓區(qū)位于裝藥前端，而高溫區(qū)仍位于點火具出口與裝藥前端之間，如圖5(c)所示；隨著燃面加質(zhì)的進行，壓力沖擊波沿內(nèi)外通道向下游傳播速度加快，燃氣在隔板前集聚，使得隔板前區(qū)域壓力快速升高，此時流場中高壓區(qū)位于級間隔板處，而高溫區(qū)由之前的點火具出口區(qū)域轉(zhuǎn)移到裝藥加質(zhì)填充區(qū)域，如圖5(d)所示；隨著時間的推移，壓力沖擊波在隔板處產(chǎn)生反射波，并沿裝藥外通道繼續(xù)傳播，外通道空氣受前后來流壓縮，壓力快速上升，隨之Ⅱ脈沖燃燒室高壓區(qū)位于裝藥外通道的中間段，而其它區(qū)域壓強分布趨向平穩(wěn)。而高溫區(qū)同時出現(xiàn)于裝藥燃氣加質(zhì)的填充區(qū)，即內(nèi)外通道，且溫度沿軸向逐漸降低，如圖5(e)所示。

圖6給出了不同時刻切面P0上內(nèi)壁面監(jiān)測線(P1與P3相連線段)和外壁面(P4與P6相連線段)上的溫度分布情況。計算結(jié)果顯示，在點火初期，隨著點火具燃氣的填充，裝藥表面沿軸向溫度依次上升。0.1 ms時，裝藥內(nèi)壁面溫度開始上升；0.5 ms后，裝藥內(nèi)壁面完全點燃；0.3 ms時，裝藥外壁面溫度開始上升；0.6 ms后，裝藥外壁面也相繼完全點燃。

2.1.2 點火初期Ⅱ脈沖燃燒室升壓過程

圖7給出了裝藥內(nèi)外壁面監(jiān)測點(P1～P6)壓強時間曲線。從圖7可知，燃氣進入燃燒室后沿內(nèi)外通道依次傳播，且外通道壓力沖擊波傳播速度明顯滯后于內(nèi)通道。

從圖7(a)可知，0.58 ms前，內(nèi)壁面各監(jiān)測點出現(xiàn)首次壓強驟升過程，升壓速率沿軸向依次增加，壓強沿軸向分布依次遞減；0.58～0.82 ms之間，裝藥內(nèi)壁面壓強沿軸向分布為依次遞增，這主要是燃氣在隔板前的集聚作用，導(dǎo)致裝藥末端壓力上升速率大于前端；0.82 ms到級間孔打開，隨著裝藥燃燒加質(zhì)產(chǎn)生的燃氣大量填充，內(nèi)通道壓力快速上升，內(nèi)壁面壓強沿軸向呈遞減分布，各監(jiān)測點壓強差減小，且增壓速率保持一致。相比于內(nèi)壁面， 0.6 ms前，外壁面監(jiān)測點壓強沿軸向呈遞減分布；0.6 ms到級間孔打開，外壁面壓強沿軸向分布先增加、后減小，但前后壓強差逐步減小，且增壓速率保持一致，這與圖5的分析結(jié)果一致，如圖7(b)所示。

圖8給出了各級間孔達到打開壓強前的壓強時間曲線。由于各級間孔分布不同，致使各級間孔入口壓強出現(xiàn)一定的時間延遲，且壓強沿中心級間孔向外依次遞減。在0.53 ms前后，各孔出現(xiàn)不同程度的衰減，隨后穩(wěn)定上升。這主要是由于隔板對壓力沖擊波的反射作用所引起的，孔1受影響最大，孔4受影響最小，由于孔2與孔3分布的對稱性，孔2與孔3的壓力時間曲線重合度較好，與預(yù)期一致。

如圖9所示，結(jié)合圖7可知壓強沖擊波沿內(nèi)壁面?zhèn)鞑ニ俾蚀笥谕獗诿妫瑑?nèi)壁面升壓速率大于外壁面，燃氣遇級間隔板發(fā)生反射沿內(nèi)外通道向上游傳播。隨著外通道燃氣的填充，裝藥內(nèi)外壁面壓差隨之上升。0.4 ms后，藥柱(P2～P5)、(P3～P6)位置內(nèi)外壓差急劇下降。由于燃氣在孔1附近集聚效應(yīng)大于其他級間孔，使得在裝藥末端(P3～P6)內(nèi)外壓差經(jīng)短暫降低后快速上升，而裝藥中段(P2～P5)因外通道壓力上升，內(nèi)外壓差持續(xù)下降。0.8 ms后，內(nèi)外通道壓差逐漸趨于穩(wěn)定。此時，裝藥內(nèi)外壁面平均壓差在0.05 MPa左右。

2.2 級間孔打開后壓力參數(shù)分析

假設(shè)各級間孔打開壓強為各孔入口壓強的平均值，設(shè)定該壓強為1.5 MPa。圖10給出了各級間孔入口打開前后壓強時間曲線。由圖10可知，在1.06 ms時，級間孔打開。由圖10放大圖可知，由于各級間孔瞬間打開，導(dǎo)致高溫高壓燃氣突然釋放，級間孔入口壓強均發(fā)生不同程度的衰減，孔1～孔4壓強降低幅度分別為26.5%、28.2%、28.1%、25.7%。級間孔打開后，隨著燃氣不斷填充，壓強出現(xiàn)振蕩上升，在7.16 ms時，級間孔入口壓強達到最大值8.16 MPa，隨后壓強開始緩慢下降。

圖11為級間孔打開前后Ⅱ脈沖燃燒室監(jiān)測點壓強時間曲線。從放大圖可知，級間孔打開過程對燃燒室不同位置壓強影響程度不同，對靠近級間孔處影響較大，對內(nèi)通道影響大于外通道，對裝藥末端影響大于前端。由圖11可知，外通道前部監(jiān)測點(P4～P5)幾乎不受影響，且級間孔打開過程對整個壓強建立過程的影響有限，壓強經(jīng)短暫振蕩后穩(wěn)定上升，在7 ms時達到峰值。

圖12為Ⅰ脈沖燃燒室監(jiān)測點壓強時間曲線。從局部放大圖可知，燃氣在級間孔釋放初期壓強振蕩劇烈，且中部監(jiān)測點(P7)比收斂段監(jiān)測點(P8)壓強振蕩明顯。這是因為高壓燃氣壓縮Ⅰ脈沖低壓氣體形成激波沿軸向朝尾部噴管傳播，激波依次掃過監(jiān)測點P7、P8，使得各監(jiān)測點處壓強先后劇增。監(jiān)測點(P8)受燃氣在收斂段的集聚作用，并沒有出現(xiàn)監(jiān)測點(P7)的壓力衰減現(xiàn)象，而是持續(xù)上升。Ⅰ脈沖燃燒室壓強在8.3 ms前后達到最大值6.4 MPa，各監(jiān)測點壓強時間曲線變化趨勢基本保持一致。

2.3 級間孔打開后流場特性分析

在發(fā)動機內(nèi)流場中，級間孔分布的非對稱性導(dǎo)致Ⅰ脈沖燃燒室流場分布具有明顯的三維特性；同時，Ⅱ脈沖燃燒室流場也形成非對稱分布。

圖13為1.2 ms時輔助面上壓強和溫度分布云圖。從圖13(a)可看出，級間孔打開后，級間孔入口壓強為1.5 MPa，遠大于級間孔出口大氣壓強，致使燃氣經(jīng)級間孔形成高度欠膨脹超聲速射流[20]，燃氣經(jīng)級間孔后持續(xù)膨脹，在各級間孔軸線下游，壓強急速下降出現(xiàn)帶狀低壓區(qū)。

圖13(b)給出了輔助面上溫度云圖分布，由于隔板對燃氣的阻礙作用，導(dǎo)致在隔板附近出現(xiàn)高溫區(qū)，后封頭及噴管收斂段直接承受高溫氣流沖擊，而在裝藥前端、裝藥末端及外圈級間孔軸線附近出現(xiàn)低溫區(qū)。在各級間孔橫截面壓強分布云圖上，前臺階截面壓強沿徑向向外降低，后臺階截面壓強分布在軸心附近與壁面之間存在低壓環(huán)，高壓區(qū)均位于各級間孔軸線處，這與Ⅰ脈沖燃燒室中各級間孔軸線下游壓強分布形成強烈反差，如圖13(c)所示。

3 結(jié)論

(1)點火初期Ⅱ脈沖燃燒室燃氣填充過程，燃氣火焰峰明顯滯后于壓力波峰。隨著燃燒室壓強不斷上升，Ⅱ脈沖燃燒室相對高壓區(qū)位置不斷發(fā)生轉(zhuǎn)變。從點火具出口到裝藥內(nèi)壁面前端，再到級間孔入口位置，高壓區(qū)最后穩(wěn)定在外通道中間位置，而在其他區(qū)域，壓強趨于均勻。

(2)點火升壓過程裝藥內(nèi)外壁面壓強差變化較大，首次升壓過程裝藥內(nèi)外壁面壓差逐漸上升；隨著外通道二次升壓，裝藥內(nèi)外壁面壓差急劇下降；臨近級間孔打開，裝藥內(nèi)外壁面沿軸向平均壓差趨于穩(wěn)定。

(3)級間孔打開對隔板附近區(qū)域壓強影響較大。由于級間孔流道面積小，對Ⅱ脈沖燃燒室整體升壓過程影響較小。級間孔打開初期，在級間孔下游區(qū)域，形成非對稱高速膨脹區(qū)。

(4)由于級間孔分布的非對稱性，導(dǎo)致壓強及溫度在燃燒室中分布呈現(xiàn)明顯的三維特性；在各級間孔中心軸線下游，均出現(xiàn)帶狀低壓區(qū)，隔板附近出現(xiàn)高溫區(qū)；而在裝藥前端、裝藥末端及外圈級間孔軸線附近，出現(xiàn)低溫區(qū)。

[1] 葉定友,王敬超. 固體發(fā)動機多次啟動技術(shù)及應(yīng)用[J]. 推進技術(shù),1989,10(4):71-74.

[2] 龔士杰,戴耀松. 多脈沖固體火箭發(fā)動機述評[J]. 飛航導(dǎo)彈,1993,23(1):32-39.

[3] Naumann K W,Stadler L. Double-pulse solid rocket motor technology-applications and technical solutions[R].AIAA 2010-6754.

[4] Stadler L J,Hoffmann S,Niedermaier H,et al. Testing and verification of the LFK-NG dual pulse motor[R]. AIAA 2006-4765.

[5] Schilling S,Trouillot P,Weigand A. On the development and testing of a 120 mm caliber double pulse motor[R]. AIAA 2004-3387.

[6] Harold D,Barry J. Demonstration of solid propellant pulse motor technologies[R]. AIAA 96-3157.

[7] Javed A,Manna P,Chakraborty D. Numerical simulation of a dual pulse solid rocket motor flow field[J]. Defence Science Journal,2012,62(6):369-374.

[8] 王長輝,劉亞冰,林震,等. 脈沖發(fā)動機陶瓷隔板的設(shè)計和試驗研究[J]. 固體火箭技術(shù),2010,33(3):327-331.

[9] 王偉,李江,王春光,等. 隔艙式雙脈沖發(fā)動機金屬膜片設(shè)計與實驗研究[J]. 推進技術(shù),2013,34(8):1115-1120.

[10] 王春光,田維平,任全彬,等. 脈沖發(fā)動機中隔層工作過程的數(shù)值分析及試驗[J]. 推進技術(shù),2012,33(5):790-795.

[11] 劉雨,利鳳祥,李越森,等. 多脈沖固體火箭發(fā)動機陶瓷艙蓋結(jié)構(gòu)分析[J]. 固體火箭技術(shù),2008,31(2):179-183.

[12] 孫娜,婁永春,孫長宏,等. 某雙脈沖發(fā)動機燃燒室兩相流場數(shù)值分析[J]. 固體火箭技術(shù),2012,35(3):335-338.

[13] 朱衛(wèi)兵,張永飛,陳宏,等. 雙脈沖發(fā)動機內(nèi)流場研究[J]. 彈箭與制導(dǎo)學(xué)報,2012,32(1):114-118.

[14] 楊春慶,魏志軍,張雷,等. 隔層式多脈沖發(fā)動機點火延遲數(shù)值仿真研究[J]. 推進技術(shù),2014,35(4):514-522.

[15] 李映坤,韓珺禮,陳雄,等. 級間通道構(gòu)型對雙脈沖發(fā)動機燃燒室局部受熱的影響[J]. 推進技術(shù),2014,35(11):1503-1510.

[16] Nishii S,Fukuda K,Kubota N. Combustion Tests of Two-Stage Pulse Rocket Motors[R]. AIAA 89-2426.

[17] 徐學(xué)文,王連生,牟俊林,等. 固體火箭發(fā)動機藥柱裂紋腔內(nèi)三維流場瞬態(tài)特性分析[J]. 固體火箭技術(shù),2007,30(1):39-43.

[18] 武曉松,陳軍,王棟. 固體火箭發(fā)動機原理[M]. 兵器工業(yè)出版社,2010.

[19] 張為華. 固體火箭發(fā)動機點火器流量的預(yù)示[J]. 宇航學(xué)報,1996,17(4):39-43.

[20] 李軍,曹從詠,徐強. 固體火箭燃氣射流近場形成與發(fā)展的數(shù)值模擬[J]. 推進技術(shù),2003,24(5):410-413.

Three-dimensional(3D)numericalsimulationfortheignitionprocessofadual-pulsemotor

YAN Deng-chao1,2, CHEN Xiong1, LI Ying-kun1, ZHU Liang1, SUN Shan1

(1.School of Mechanical Engineering,Nanjing University of Science and Technology, Nanjing 210094, China; 2. Jiangnan Design Institute of Machinery & Electricity，Guiyang 550009, China)

Taking the chamber, pulse separation device (PSD) and the nozzle of dual pulse motor as research objectives, the three dimensional (3D) flow field characteristics of fuel gas were studied and analyzed during the ignition process of the second pulse chamber by means of numerical simulations. The results showed that at the beginning of ignition, the ignition gas pressure wave arrived at the PSD prior to the flame and spread in the form of pressure shock. The position of the high pressure zone was constantly changed with time. The opening of interstage hole had a greater effect on pressure at the end of the propellant, but less effect on the whole process of pressure increasing in the second chamber. A highly under-expanded jet and a non-symmetric band high-speed expansion zone were formed in the first chamber after the gas-flow through the interstage hole. The pressure and temperature contours were found to be significant 3D flow field paramters in the chamber because of the asymmetry distribution of the PSD. The higher temperature zone was located in the front cavity of clapboard, but the lower temperature zone was located in the rear end of the propellant charge, as well as the axial neighborhood region of the interstage hole.

dual pulse motor；3D flow field；numerical simulation；ignition process

V435

A

1006-2793(2017)05-0545-07

10.7673/j.issn.1006-2793.2017.05.002

(編輯：崔賢彬)