冷林濤,翁春生,白橋棟,林 玲

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

單組元粉末發(fā)動(dòng)機(jī)內(nèi)流場(chǎng)數(shù)值模擬研究

冷林濤,翁春生,白橋棟,林 玲

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

為了研究粉末發(fā)動(dòng)機(jī)的工作特性,建立了以火藥顆粒為燃料的單組元粉末發(fā)動(dòng)機(jī)理論模型,采用CE/SE方法對(duì)粉末發(fā)動(dòng)機(jī)內(nèi)流場(chǎng)進(jìn)行數(shù)值模擬。分析了初始顆粒粒徑、初始固相體積分?jǐn)?shù)以及堵蓋對(duì)粉末發(fā)動(dòng)機(jī)燃燒過程的影響。計(jì)算結(jié)果表明:顆粒初始粒徑的減小和固相初始體積分?jǐn)?shù)的提高均能提高發(fā)動(dòng)機(jī)燃燒室內(nèi)的壓力;固相體積分?jǐn)?shù)對(duì)發(fā)動(dòng)機(jī)內(nèi)溫度峰值影響不大;在點(diǎn)傳火階段,堵蓋能提高點(diǎn)火壓強(qiáng),堵蓋打開壓力越大,增壓越明顯。

粉末發(fā)動(dòng)機(jī);顆粒燃燒;單組元;內(nèi)流場(chǎng);堵蓋

粉末發(fā)動(dòng)機(jī)是一種以高能粉末為燃料的新概念發(fā)動(dòng)機(jī)[1],它使用的粉末燃料存儲(chǔ)在燃料箱中,工作時(shí)由專用的輸送系統(tǒng)送入燃燒室,具有推力可調(diào)、安全性好等優(yōu)點(diǎn),在航空航天領(lǐng)域有很好的應(yīng)用前景,特別是在太空探測(cè)方面,粉末發(fā)動(dòng)機(jī)一方面可通過自身攜帶的燃料/氧化劑在真空環(huán)境下工作,另一方面也可以只攜帶純高能粉末燃料,利用其他行星大氣中的氣體為氧化劑作動(dòng)力系統(tǒng)[2]。

早在20世紀(jì)60年代,美國(guó)就提出用金屬粉末作為發(fā)動(dòng)機(jī)燃料,并開展了多種金屬粉末在不同大氣環(huán)境下的燃燒特性試驗(yàn),國(guó)外對(duì)于粉末發(fā)動(dòng)機(jī)的研究大多以研究金屬顆粒的燃燒特性為主,Wang等[3]研究了不同粒徑的鎂粉在空氣-乙炔環(huán)境下的燃燒特性。目前對(duì)粉末發(fā)動(dòng)機(jī)的試驗(yàn)研究比較多[4-6],對(duì)于數(shù)值計(jì)算,則比較少。D’yachenko等[7]研究了粉末燃燒產(chǎn)物在噴管處的流動(dòng)情況;申慧君等[8-9]開展了粉末發(fā)動(dòng)機(jī)燃燒室數(shù)值模擬研究工作,利用吉布斯最小自由能法計(jì)算了給定壓強(qiáng)下粉末燃燒產(chǎn)物的組份,并用顆粒軌道模型分析了產(chǎn)物相態(tài)等因素對(duì)發(fā)動(dòng)機(jī)燃燒效率的影響;李悅[10]等分析了不同氧燃比對(duì)粉末發(fā)動(dòng)機(jī)燃燒性能的影響,并指出在一定范圍內(nèi)氧燃比較高時(shí),燃燒室溫度反而較低。

國(guó)內(nèi)外學(xué)者對(duì)粉末發(fā)動(dòng)機(jī)的研究大多以金屬粉末為燃料,水/氣體為氧化劑的雙組元發(fā)動(dòng)機(jī)為主,對(duì)于以火藥顆粒為燃料的單組元發(fā)動(dòng)機(jī)研究尚未見報(bào)道。火藥由于具有較高的能量密度和優(yōu)異的燃燒性能常被用作火炮的發(fā)射藥,為了研究粉末發(fā)動(dòng)機(jī)的工作特性,本文采用雙流體假設(shè)模型,對(duì)設(shè)計(jì)以小粒黑火藥為燃料的粉末發(fā)動(dòng)機(jī)進(jìn)行兩相流軸對(duì)稱數(shù)值模擬研究,采用高精度CE/SE算法研究不同初始顆粒粒徑、不同初始固相體積分?jǐn)?shù)和堵蓋的粉末發(fā)動(dòng)機(jī)燃燒過程,為粉末發(fā)動(dòng)機(jī)的設(shè)計(jì)及熱態(tài)試驗(yàn)提供參考。

1 理論模型

粉末發(fā)動(dòng)機(jī)的燃料以噴注的方式進(jìn)入燃燒室,燃燒室內(nèi)的燃料以氣固兩相流的形式存在,顆粒速度不高且濃度較大。顆粒在粉末發(fā)動(dòng)機(jī)內(nèi)停留時(shí)間較短,使顆粒在較短的時(shí)間內(nèi)著火并迅速燃燒是提高發(fā)動(dòng)機(jī)燃燒效率的關(guān)鍵。初始時(shí)刻設(shè)立局部高溫區(qū)域進(jìn)行點(diǎn)火,顆粒達(dá)到著火溫度后在發(fā)動(dòng)機(jī)頭部首先燃燒,產(chǎn)生高溫燃?xì)?高溫燃?xì)馔ㄟ^對(duì)流輻射換熱將能量傳給未燃區(qū)域,最終使燃燒室內(nèi)顆粒全部被點(diǎn)燃,燃?xì)庠趪姽苤信蛎涀龉Α?/p>

本文計(jì)算模型如圖1,根據(jù)對(duì)稱性采用1/2模型計(jì)算,發(fā)動(dòng)機(jī)頭部設(shè)有局部高溫區(qū)域模擬點(diǎn)火過程,整個(gè)發(fā)動(dòng)機(jī)長(zhǎng)度L=0.2 m,其中噴管收斂段長(zhǎng)0.03 m,擴(kuò)張段長(zhǎng)0.02 m。燃燒室半徑R1=0.03 m,喉部半徑R2=0.015 m,噴管出口半徑R3=0.02 m。

1.1 基本假設(shè)

由于兩相間的作用十分復(fù)雜,為了研究方便,本模型研究的單一球形火藥顆粒均勻分布在燃燒室中,主要假設(shè)如下:

①顆粒由單一藥粒組成,其形狀和大小嚴(yán)格一致且均勻分布在燃燒室中,單顆粒火藥燃燒規(guī)律相同;

②將顆粒當(dāng)作擬流體處理,不考慮粘性,顆粒初始時(shí)刻溫度相同;

③顆粒燃燒產(chǎn)物組份不變,火藥氣體熱力學(xué)參數(shù)保持不變;

④顆粒不可壓縮,即固相物質(zhì)密度保持為常量;

⑤顆粒著火取表面溫度準(zhǔn)則,忽略壁面與外界換熱;

⑥顆粒噴入燃燒室后均勻分布;

⑦點(diǎn)火準(zhǔn)則采用Gough點(diǎn)火準(zhǔn)則,即顆粒表面溫度達(dá)到著火溫度時(shí)被點(diǎn)燃。

1.2 基本方程

根據(jù)以上假設(shè)可推導(dǎo)出粉末發(fā)動(dòng)機(jī)內(nèi)流場(chǎng)軸對(duì)稱兩相流方程組為[11]

(1)

(2)

(3)

(4)

(5)

(6)

(7)

式中:下標(biāo)g,s分別代表氣相與固相;x,r分別為軸向與徑向坐標(biāo);ρ為密度;u為軸向速度;v為徑向速度;p為氣相壓力;E為總能;φ為體積分?jǐn)?shù),滿足φg+φs=1;Fs為相間阻力;mc為燃?xì)馍伤俾?es為火藥化學(xué)潛能;Qd為相間傳熱,具體求解公式見文獻(xiàn)[11～12]。

1.3 計(jì)算方法

本文采用CE/SE方法進(jìn)行計(jì)算,CE/SE方法具有精度高、計(jì)算格式簡(jiǎn)單等優(yōu)點(diǎn),這種方法構(gòu)造思想是將時(shí)間和空間統(tǒng)一起來處理,在處理多維問題時(shí)不需要方向交替處理,計(jì)算實(shí)踐表明CE/SE方法具有很高的分辨率[13]。CE/SE的計(jì)算格式如下:

(8)

1.4 邊界條件與初始條件

邊界條件:壁面和入口處封閉端為固壁邊界條件。中心對(duì)稱軸為軸對(duì)稱邊界條件。噴管出口處根據(jù)堵蓋打開情況分別處理,堵蓋打開前為固壁邊界條件,堵蓋打開后若燃?xì)庖詠喡曀倭鞒龉芡鈩t管內(nèi)燃?xì)鈮毫Φ扔诔隹陟o壓,管內(nèi)流場(chǎng)受到管外壓力影響,壓力為出口靜壓;若燃?xì)庖猿曀倭鞒龉芡鈩t管內(nèi)流場(chǎng)不受管外流場(chǎng)影響。

初始條件:在發(fā)動(dòng)機(jī)頭部設(shè)置局部高溫區(qū)域以模擬點(diǎn)火過程,整個(gè)燃燒室內(nèi)顆粒均勻分布。

2 計(jì)算結(jié)果及其討論

2.1 發(fā)動(dòng)機(jī)內(nèi)流場(chǎng)分析

以顆粒半徑為300 μm,固相體積分?jǐn)?shù)為0.1,發(fā)動(dòng)機(jī)無堵蓋為例進(jìn)行發(fā)動(dòng)機(jī)內(nèi)流場(chǎng)分析。

圖2為在10.26 ms時(shí)刻發(fā)動(dòng)機(jī)內(nèi)壓力和溫度的瞬時(shí)分布云圖?？傮w上,這些參數(shù)結(jié)果分布符合發(fā)動(dòng)機(jī)內(nèi)流場(chǎng)的物理行為,壓力分布和溫度分布一致,燃料顆粒在發(fā)動(dòng)機(jī)頭部點(diǎn)火后產(chǎn)生高溫燃?xì)?在燃燒室內(nèi)積聚向右傳播,燃?xì)馔ㄟ^噴管膨脹做功,燃燒室內(nèi)的壓力梯度和溫度梯度不大,燃?xì)庠趪姽苤袎毫ο陆?溫度下降,噴管喉部壓力和溫度分別為2.3 MPa和1 182 K。圖3為中心軸線處不同時(shí)刻的壓力分布,由圖3可知,發(fā)動(dòng)機(jī)頭部壓力最大,在同一時(shí)刻燃燒室內(nèi)壓力基本不變,噴管處壓力下降明顯,11.05 ms時(shí)刻發(fā)動(dòng)機(jī)頭部壓力約為4.2 MPa,噴管出口壓力為0.3 MPa。在穩(wěn)定燃燒階段,燃燒室內(nèi)壓力平均每1 ms增加約1 MPa,可見粉末發(fā)動(dòng)機(jī)內(nèi)顆粒有較高的燃燒效率。

圖4為發(fā)動(dòng)機(jī)點(diǎn)火傳火過程中不同時(shí)刻發(fā)動(dòng)機(jī)內(nèi)溫度分布云圖。由圖可以清晰地看到發(fā)動(dòng)機(jī)內(nèi)的溫度變化情況,燃燒波從發(fā)動(dòng)機(jī)頭部傳向尾部,在3.05ms時(shí)刻頭部區(qū)域由初始溫度升高至1 163K(如圖4(a)所示)。顆粒燃燒過程中從燃燒波波前到波后依次分為未燃區(qū)、預(yù)熱區(qū)和已燃區(qū),高溫燃?xì)庖詫?duì)流和輻射的形式不斷將能量向軸向和徑向傳遞,在此狀態(tài)下已燃區(qū)的溫度是未燃區(qū)溫度的2～3倍(如圖4(b)所示)。由圖4(c)可知,在7.03ms時(shí)刻,燃燒波傳至發(fā)動(dòng)機(jī)出口,發(fā)動(dòng)機(jī)內(nèi)顆粒被全部點(diǎn)燃,在點(diǎn)傳火階段發(fā)動(dòng)機(jī)內(nèi)最高溫度達(dá)1 600K。

圖5為顆粒燃燒過程中不同時(shí)刻氣相速度分布流線圖。由圖可知,在點(diǎn)傳火初始階段,發(fā)動(dòng)機(jī)內(nèi)徑向流動(dòng)明顯,隨著燃燒波向出口處傳播,徑向流動(dòng)逐漸減弱,傳火過程發(fā)動(dòng)機(jī)內(nèi)速度場(chǎng)逐漸以軸向流動(dòng)為主。

圖6為中心軸線處距發(fā)動(dòng)機(jī)頭部不同位置溫度隨時(shí)間的變化,圖7為不同時(shí)刻發(fā)動(dòng)機(jī)中心軸線處的溫度變化。由圖6和圖7可以看出,發(fā)動(dòng)機(jī)在整個(gè)工作過程中頭部溫度最高,隨著燃燒時(shí)間的推進(jìn),更多的能量被釋放,在13.26ms時(shí)刻出現(xiàn)最高溫度,達(dá)2 246K。初始時(shí)刻在經(jīng)歷一個(gè)短暫的低溫區(qū)后溫度迅速升高,之后發(fā)動(dòng)機(jī)內(nèi)出現(xiàn)一段溫度平臺(tái)期,平臺(tái)期過后溫度繼續(xù)增加至最大。分析原因主要是由于在發(fā)動(dòng)機(jī)點(diǎn)火初期,室內(nèi)每一處顆粒均存在一個(gè)點(diǎn)火延遲期,在點(diǎn)火未完成階段由于該位置的顆粒沒有燃燒釋熱,所以初始階段表現(xiàn)為升溫緩慢,當(dāng)該位置顆粒達(dá)到著火溫度后,顆粒會(huì)在極短的時(shí)間內(nèi)迅速著火燃燒,能量大量釋放,溫度迅速升高。達(dá)到一定溫度后(約1 000K),在固相顆粒不斷燃燒,高溫燃?xì)庀虻蜏貐^(qū)傳遞能量并向出口排出的過程中,溫度會(huì)有一個(gè)平臺(tái)期,距發(fā)動(dòng)機(jī)頭部越遠(yuǎn),溫度平臺(tái)期越小,由圖6可知,溫度平臺(tái)在噴管處消失。從圖7中可以看出發(fā)動(dòng)機(jī)中心軸線處溫度隨時(shí)間的變化,在傳火過程中火焰峰溫度基本不變,燃燒波傳至發(fā)動(dòng)機(jī)出口時(shí)火焰峰消失,整個(gè)發(fā)動(dòng)機(jī)內(nèi)顆粒處于全面燃燒狀態(tài),溫度整體上升。

2.2 顆粒粒徑對(duì)發(fā)動(dòng)機(jī)的影響

以固相體積分?jǐn)?shù)為0.1且發(fā)動(dòng)機(jī)無堵蓋為例,分析不同顆粒粒徑下的發(fā)動(dòng)機(jī)內(nèi)流場(chǎng)。

圖8為中心軸線上距發(fā)動(dòng)機(jī)頭部0.1m處不同顆粒粒徑d下溫度隨時(shí)間的變化圖,圖9為中心軸線上距發(fā)動(dòng)機(jī)頭部0.1m處不同顆粒粒徑下壓力隨時(shí)間的變化圖。由圖8和圖9可知,顆粒粒徑越小,溫度平臺(tái)越低,顆粒著火時(shí)間越短,點(diǎn)火延遲期越短,發(fā)動(dòng)機(jī)內(nèi)該位置處的溫度峰值和壓力峰值越高,而到達(dá)溫度峰值和壓力峰值的時(shí)間越短,相同大小的顆粒燃燒到達(dá)的溫度峰值和壓力峰值的時(shí)間是一致的。300μm顆粒的點(diǎn)火延遲期比400μm顆粒的點(diǎn)火延遲期短約0.8ms,較小顆粒(300μm)的粉末發(fā)動(dòng)機(jī)在距發(fā)動(dòng)機(jī)0.1m處的溫度峰值和壓力峰值要比較大顆粒(400μm)粉末發(fā)動(dòng)機(jī)分別高約30%和60%,可見顆粒粒徑對(duì)粉末發(fā)動(dòng)機(jī)的影響顯著。分析原因:一方面由于顆粒粒徑越小,其比表面積越大,單位體積內(nèi)顆粒吸收對(duì)流和輻射的能量也就越多,所以表現(xiàn)為發(fā)動(dòng)機(jī)內(nèi)溫度平臺(tái)越低;另一方面由于顆粒粒徑越小,燃燒越充分,燃燒效率越高,所以表現(xiàn)為發(fā)動(dòng)機(jī)內(nèi)溫度峰值和壓力峰值越高,到達(dá)溫度峰值和壓力峰值的時(shí)間越快。由圖8可知,粒徑為300μm的顆粒燃燒時(shí),約在13.26ms時(shí)刻最大溫度和壓力分別為1 625K和5.41MPa。當(dāng)顆粒粒徑較大時(shí),變化規(guī)律相反。

2.3 固相體積分?jǐn)?shù)對(duì)發(fā)動(dòng)機(jī)的影響

以顆粒粒徑為300μm且發(fā)動(dòng)機(jī)無堵蓋為例,分析不同固相體積分?jǐn)?shù)下的發(fā)動(dòng)機(jī)內(nèi)流場(chǎng)。

圖10為中心軸線處發(fā)動(dòng)機(jī)頭部不同固相體積分?jǐn)?shù)下壓力隨時(shí)間的變化,圖11為不同固相體積分?jǐn)?shù)下發(fā)動(dòng)機(jī)內(nèi)溫度峰值Tmax與壓力峰值pmax的變化圖。由圖10可知,固相體積分?jǐn)?shù)越大,該位置處的壓力峰值越高,到達(dá)壓力峰值的時(shí)間越短,表明增大固相體積分?jǐn)?shù)能提高發(fā)動(dòng)機(jī)室內(nèi)壓力。在本模型的計(jì)算條件下,固相體積分?jǐn)?shù)每增加0.01,發(fā)動(dòng)機(jī)室內(nèi)壓力峰值平均增加30%,到達(dá)壓力峰值的時(shí)間提前0.85ms。圖11中線條是由實(shí)心點(diǎn)擬合所得,實(shí)心點(diǎn)為模型的計(jì)算值。由圖11可知,固相體積分?jǐn)?shù)對(duì)溫度峰值的影響較小,對(duì)壓力峰值影響較大,在固相體積分?jǐn)?shù)較大的情況下,發(fā)動(dòng)機(jī)內(nèi)所達(dá)到的壓力峰值和溫度峰值也越大。較大的固相體積分?jǐn)?shù)能顯著提高發(fā)動(dòng)機(jī)的室內(nèi)壓力,當(dāng)固相體積分?jǐn)?shù)增加時(shí),其燃燒釋放的熱量也將增大。楊朝晉等[14]通過研究粉末發(fā)動(dòng)機(jī)內(nèi)鎂顆粒的燃燒過程發(fā)現(xiàn),在顆粒濃度相對(duì)較小的情況下增加顆粒濃度有助于提高火焰溫度,而在顆粒濃度相對(duì)較大的情況下變化規(guī)律相反。由于本模型使用的是單組元燃料,不需要額外的氧化劑,而鎂粉燃燒過程需要額外添加氧化劑,顆粒濃度過大時(shí)會(huì)引發(fā)鎂顆粒對(duì)氣相中氧含量的爭(zhēng)奪,在氧化劑含量一定的情況下,過高的固相體積分?jǐn)?shù)反而會(huì)降低發(fā)動(dòng)機(jī)的燃燒性能。這也是單組元粉末發(fā)動(dòng)機(jī)的優(yōu)勢(shì)之一,因此在單組元粉末發(fā)動(dòng)機(jī)中應(yīng)盡可能地提高噴入燃燒室內(nèi)的燃料顆粒的含量。本模型在不同固相體積分?jǐn)?shù)下計(jì)算所得的黑火藥燃燒溫度約為2 100K,文獻(xiàn)[15]采用多光譜輻射測(cè)溫系統(tǒng)測(cè)得黑火藥燃燒溫度約1 950K,兩者比較接近。

2.4 堵蓋對(duì)發(fā)動(dòng)機(jī)的影響

在火箭發(fā)動(dòng)機(jī)設(shè)計(jì)中,若初始時(shí)刻室內(nèi)壓力過低,會(huì)對(duì)發(fā)動(dòng)機(jī)的點(diǎn)火傳火過程產(chǎn)生不利影響,因此常設(shè)計(jì)噴管堵蓋對(duì)發(fā)動(dòng)機(jī)進(jìn)行增壓,促使發(fā)動(dòng)機(jī)內(nèi)燃料顆粒的穩(wěn)定燃燒。文獻(xiàn)[16]通過調(diào)整堵蓋厚度來試驗(yàn)不同的堵蓋打開壓力,在不同的試驗(yàn)條件下測(cè)得堵蓋打開壓力為0.74～2.53MPa。本模型中考慮的堵蓋打開壓力在此范圍內(nèi),堵蓋設(shè)置在噴管出口。

以顆粒粒徑為300μm且固相體積分?jǐn)?shù)為0.1為例,分析堵蓋對(duì)發(fā)動(dòng)機(jī)內(nèi)流場(chǎng)的影響。

2.4.1 有/無堵蓋發(fā)動(dòng)機(jī)流場(chǎng)分析

本節(jié)在計(jì)算中設(shè)置堵蓋打開壓力為1MPa。

圖12為發(fā)動(dòng)機(jī)噴管喉部壓力隨時(shí)間的變化圖,圖13為發(fā)動(dòng)機(jī)頭部壓力隨時(shí)間的變化圖。由圖12可知,堵蓋在7.25ms時(shí)打開,堵蓋打開后發(fā)動(dòng)機(jī)尾部出現(xiàn)短暫壓力波動(dòng)。由圖13可以看出,堵蓋對(duì)發(fā)動(dòng)機(jī)頭部壓力影響較小。對(duì)比圖12和圖13可知,設(shè)計(jì)有堵蓋的粉末發(fā)動(dòng)機(jī)能顯著提高室內(nèi)壓力,加速顆粒燃燒,本模型計(jì)算下有堵蓋的發(fā)動(dòng)機(jī)比無堵蓋的發(fā)動(dòng)機(jī)燃燒室內(nèi)壓力增加約20%,出現(xiàn)壓力峰值的時(shí)間提前約11%。

圖14為點(diǎn)火過程中發(fā)動(dòng)機(jī)有/無堵蓋中心軸線處壓力分布圖。由圖14可知,有堵蓋發(fā)動(dòng)機(jī)比無堵蓋發(fā)動(dòng)機(jī)在點(diǎn)火階段的室內(nèi)壓力要高,在5.12ms時(shí)刻高約23%,增設(shè)堵蓋會(huì)在發(fā)動(dòng)機(jī)點(diǎn)火階段形成一個(gè)較高的壓力峰值,有利于點(diǎn)火的進(jìn)行。由圖14可以看出,在5.12ms時(shí)距發(fā)動(dòng)機(jī)頭部0.1m處存在點(diǎn)火壓力峰值為0.185MPa。

圖15為有/無堵蓋發(fā)動(dòng)機(jī)內(nèi)氣相速度分布云圖。由圖15可以看到發(fā)動(dòng)機(jī)內(nèi)速度的變化,軸向流動(dòng)明顯,徑向流動(dòng)微弱。有堵蓋的發(fā)動(dòng)機(jī)燃?xì)馑俣扔捎诒诿娴姆瓷渥饔迷趪姽芴幮纬伤俣然亓鲄^(qū),在x=0.185處形成速度為0的交界面,整個(gè)燃燒室中部速度較大。無堵蓋發(fā)動(dòng)機(jī)在噴管處有最大速度約62m/s,流場(chǎng)相對(duì)穩(wěn)定。

2.4.2 堵蓋打開壓力對(duì)發(fā)動(dòng)機(jī)的影響

圖16為3種堵蓋打開壓力pd下中心軸線上距發(fā)動(dòng)機(jī)頭部0.1m處的壓力分布。

由圖16可知,堵蓋壓力越大,發(fā)動(dòng)機(jī)內(nèi)所能達(dá)到的壓力越大,當(dāng)堵蓋壓力為0.8MPa時(shí),該位置處的最大壓力為6.15MPa,本模型計(jì)算條件下堵蓋壓力每增加0.2MPa,該位置處的最大壓力增加約0.4MPa,可見粉末發(fā)動(dòng)機(jī)增設(shè)堵蓋有利于提高發(fā)動(dòng)機(jī)的室內(nèi)壓力,增大顆粒的燃燒效率,提高顆粒在燃燒室內(nèi)的駐留時(shí)間。但是,堵蓋壓力過高可能引起初始時(shí)刻壓力升高異常,嚴(yán)重時(shí)還會(huì)導(dǎo)致發(fā)動(dòng)機(jī)爆炸,所以噴管堵蓋的設(shè)計(jì)應(yīng)結(jié)合粉末發(fā)動(dòng)機(jī)的具體使用范圍和條件綜合考慮。

3 結(jié)論

本文建立了粉末發(fā)動(dòng)機(jī)內(nèi)流場(chǎng)理論模型,研究表明采用CE/SE方法數(shù)值模擬火藥顆粒在粉末發(fā)動(dòng)機(jī)內(nèi)的燃燒過程是合理的,為粉末發(fā)動(dòng)機(jī)的研究提供了一種有力的數(shù)值工具。得出的主要結(jié)論如下:

①顆粒粒徑和固相體積分?jǐn)?shù)對(duì)粉末發(fā)動(dòng)機(jī)內(nèi)流場(chǎng)影響顯著。粒徑越小,固相體積分?jǐn)?shù)越大,發(fā)動(dòng)機(jī)內(nèi)的溫度峰值與壓力峰值越大,固相體積分?jǐn)?shù)對(duì)發(fā)動(dòng)機(jī)內(nèi)溫度峰值影響較小。數(shù)值計(jì)算表明,減小顆粒粒徑與增大固相體積分?jǐn)?shù)均能提高粉末發(fā)動(dòng)機(jī)的溫度與壓力。

②堵蓋對(duì)粉末發(fā)動(dòng)機(jī)的壓力影響明顯,在發(fā)動(dòng)機(jī)出口處增設(shè)堵蓋可有效提高燃燒室內(nèi)的壓力,提高點(diǎn)火性能,堵蓋打開前,增設(shè)堵蓋1MPa能使發(fā)動(dòng)機(jī)燃燒室內(nèi)壓力增加約20%。堵蓋打開壓力越大,點(diǎn)傳火階段增壓越明顯。有堵蓋時(shí)發(fā)動(dòng)機(jī)出口處存在一定的壓力震蕩,但對(duì)發(fā)動(dòng)機(jī)整體的影響很小,可以忽略。

[1] SHAFIROVICH E,VARMA A.Metal-CO2propulsion for mars missions:current status and opportunities[J].Journal of Propulsion and Power,2008,24(3):385-394.

[2] SZABO J,MILLER T,HERR J,et al.Metal bipropellant rockets for martian ascent vehicles:AIAA 2011-5834[R].USA:AIAA,2011.

[3] WANG S,CORCORAN A L,DREIZIN E L.Combustion of magnesium powders in products of an air/acetylene flame[J].Combustion and Flame,2015,162(4):1 316-1 325.

[4] MILLER T F,HERR J D.Green rocket propulsion by reaction of Al and Mg powders and water:AIAA 2004-4037[R].USA:AIAA,2004.

[5] XIA Z X,SHEN H J,HU J X,et al.Experimental investigation of powdered metals fuel ramjet:AIAA 2008-5131[R].USA:AIAA,2008.

[6] 鄧哲,胡春波,盧子元,等.Metal/N2O粉末火箭發(fā)動(dòng)機(jī)實(shí)驗(yàn)研究[J].固體火箭技術(shù),2015,38(2):220-224. DEGN Zhe,HU Chun-bo,LU Zi-yuan,et al.Experimental research on metal/N2O powder rocket engine[J].Journal of Solid Rocket Technology,2015,38(2):220-224.(in Chinese)

[7] D’YACHENKO N N,D’YACHENKO L I,GUROVA V S,et al.Numerical investigation of the flow of combustion products containing high-dispersive aluminum powder in the solid-fuel rocket engine nozzles[J].Vestn Tomsk Gos Univ Mat Mekh,2016,2(40):63-70.

[8] 申慧君,夏智勛,胡建新,等.粉末燃料沖壓發(fā)動(dòng)機(jī)理論性能分析[J].推進(jìn)技術(shù),2007,28(2):181-185. SHEN Hui-jun,XIA Zhi-xun,HU Jian-xin,et al.Theoretical performance analysis of the powdered fuel ramjet[J].Journal of Propulsion Technology,2007,28(2):181-185.(in Chinese)

[9] 申慧君,夏智勛,胡建新,等.粉末燃料沖壓發(fā)動(dòng)機(jī)燃燒室兩相流數(shù)值模擬[J].固體火箭技術(shù),2007,30(6):474-477. SHEN Hui-jun,XIA Zhi-xun,HU Jian-xin,et al.Numerical simulation on two-phase flow in the powder fuel ramjet combustor[J].Journal of Solid Rocket Technology,2007,30(6):474-477.(in Chinese)

[10] 李悅,胡春波,孫?？?等.粉末火箭發(fā)動(dòng)機(jī)燃燒室燃燒流動(dòng)特性研究[J].固體火箭技術(shù),2014,37(6):792-796. LI Yue,HU Chun-bo,SUN Hai-jun,et al.Study on flow and combustion characteristics in the thrust chamber of powder rocket[J].Journal of Solid Rocket Technology,2014,37(6):792-796.(in Chinese)

[11] 翁春生,王浩.計(jì)算內(nèi)彈道學(xué)[M].北京:國(guó)防工業(yè)出版社,2006. WENG Chun-sheng,WANG Hao.Computational interior ballistics[M].Beijing:National Defense Industry Press,2006.(in Chinese)

[12] 曾卓雄,胡春波,姜培正,等.密相、湍流、可壓、氣固兩相流理論及在管道中的應(yīng)用[J].空氣動(dòng)力學(xué)學(xué)報(bào),2001,19(1):107-118. ZENG Zhuo-xiong,HU Chun-bo,JIANG Pei-zheng,et al.Analysis of dense,turbulent,compressible,gas-particle two-phase flow and its application to pipe flow[J].Acta Aerodynamica Sinica,2001,19(1):107-118.(in Chinese)

[13] WANG G,WANG J T,LIU K X.New numerical algorithms in super CE/SE and their applications in explosion mechanics[J].Science China Physics,Mechanics & Astronomy,2010,53(2):237-243.

[14] 楊晉朝,夏智勛,胡建新.鎂顆粒群非穩(wěn)態(tài)著火過程數(shù)值模擬[J].物理學(xué)報(bào),2012,61(16):306-315. YANG Jin-chao,XIA Zhi-xun,HU Jian-xin.Numerical studies of unsteady ignition of pulverized magnesium particle cloud[J].Acta Physica Sinica,2012,61(16):306-315.(in Chinese)

[15] 李占英,席蘭霞,陳軍,等.多光譜輻射測(cè)溫技術(shù)測(cè)量火工煙火藥劑燃燒溫度[J].光譜學(xué)與光譜分析,2010,30(8):2 062-2 064. LI Zhan-ying,XI Lan-xia,CHEN Jun,et al.Study on flame temperature measurement of pyrotechnics using multi-spectral thermometer[J].Spectroscopy and Spectral Analysis,2010,30(8):2 062-2 064.(in Chinese)

[16] 屠小昌,王占利,王小娟.硬質(zhì)泡沫噴管堵蓋試驗(yàn)及研究[J].固體火箭技術(shù),2004,27(1):81-83. TU Xiao-chang,WANG Zhan-li,WANG Xiao-juan.Test and research of nozzle closure made of hard foamed plastics[J].Journal of Solid Rocket Technology,2004,27(1):81-83.(in Chinese)

Numerical Simulation on Internal Flow-field of Monopropellant Powder Engine

LENG Lin-tao,WENG Chun-sheng,BAI Qiao-dong,LIN Ling

(National Key Laboratory of Transient Physics,Nanjing University of Science and Technology,Nanjing 210094,China)

To study the operation performance of monopropellant powder engine,a theoretical model of the monopropellant powder engine using grains as fuel was established,and the internal flow-field of powder engine was numerically simulated by the CE/SE method.The effects of the initial particles size,the initial solid volume fraction and the cover on the combustion process of powder engine were analyzed.The results show that the decrease of the particles size and the increase of the solid volume fraction can improve the chamber pressure significantly;the solid volume-fraction has little effect on the peak temperature of the powder engine.The opening pressure of cover has great influence on ignition process;the greater the opening pressure of cover,the higher the combustion chamber pressure.

powder engine;particle combustion;monopropellant;internal flow field;cover

2016-12-14

國(guó)家自然科學(xué)基金項(xiàng)目(11372141,11472138);國(guó)防預(yù)研基金項(xiàng)目(9140C300205140C30137);裝備預(yù)研基金項(xiàng)目(9140C300202120C30)

冷林濤(1992- ),男,碩士研究生,研究領(lǐng)域?yàn)榉勰┌l(fā)動(dòng)機(jī)。E-mail:lenglintaollt@163.com。

翁春生(1964- ),男,教授,博士生導(dǎo)師,研究領(lǐng)域?yàn)橥七M(jìn)技術(shù)。E-mail:wengcs@126.com。

V435.12

A

1004-499X(2017)02-0058-07