肖 華

江南機(jī)電設(shè)計(jì)研究所 貴州 貴陽(yáng) 550000

引言

當(dāng)今為了提高導(dǎo)彈生存能力,強(qiáng)調(diào)防區(qū)外發(fā)射,強(qiáng)調(diào)彈道可變性；為了對(duì)付高機(jī)動(dòng)性目標(biāo),要求導(dǎo)彈具有足夠大的末速度；為了在最短時(shí)間攻擊目標(biāo),要求導(dǎo)彈的平均速度最大。從這些角度出發(fā),雙脈沖固體火箭發(fā)動(dòng)機(jī)成為足以代替助推—主發(fā)動(dòng)機(jī)的另一種導(dǎo)彈發(fā)動(dòng)機(jī),雙脈沖固體火箭發(fā)動(dòng)機(jī)在提高導(dǎo)彈的有效射程,靈活控制導(dǎo)彈速度以及諸多方面可以充分發(fā)揮其優(yōu)越性。高速防空反導(dǎo)導(dǎo)彈由于高機(jī)動(dòng)要求,發(fā)動(dòng)機(jī)在大過(guò)載條件下內(nèi)部流場(chǎng)發(fā)生改變,會(huì)造成一脈沖燃燒室絕熱層的異常燒蝕而失效,從而造成發(fā)動(dòng)機(jī)在大過(guò)載條件下結(jié)構(gòu)完整性失效。發(fā)動(dòng)機(jī)在導(dǎo)彈飛行狀態(tài)與地面試車(chē)狀態(tài)差別較大,造成發(fā)動(dòng)機(jī)飛行安全性降低,因此對(duì)發(fā)動(dòng)機(jī)處于大過(guò)載條件下的內(nèi)部流場(chǎng)特性分析十分必要,為一脈沖燃燒室絕熱層設(shè)計(jì)提供依據(jù)。

本文采用商用CFD軟件,分析了發(fā)動(dòng)機(jī)二脈沖工作時(shí),不同過(guò)載條件下發(fā)動(dòng)機(jī)內(nèi)流場(chǎng)特性。

2 計(jì)算方法

2.1 基本假設(shè) 由于雙脈沖固體火箭發(fā)動(dòng)機(jī)燃燒室內(nèi)的流動(dòng)是一個(gè)相當(dāng)復(fù)雜的物理和化學(xué)問(wèn)題,為簡(jiǎn)化計(jì)算,本文模型基于如下假設(shè):

(1)在燃燒室內(nèi)流動(dòng)過(guò)程中,均不再發(fā)生化學(xué)反應(yīng),也忽略熱輻射的作用；

(2)忽略燃燒產(chǎn)物中的固相顆粒,將其視為純氣相的完全氣體混合物；

(3)不考慮侵蝕燃燒效應(yīng)的影響和內(nèi)壁面燒蝕作用。

2.2 計(jì)算模型 典型雙脈沖固體火箭發(fā)動(dòng)機(jī)結(jié)構(gòu)示意圖如圖1所示。

圖1 雙脈沖固體火箭發(fā)動(dòng)機(jī)示意圖

2.3 數(shù)值計(jì)算方法 采用CFD商業(yè)軟件Fluent,通過(guò)求解Reynolds平均Navier-Stokes方程數(shù)值模擬了發(fā)動(dòng)機(jī)內(nèi)流場(chǎng),控制方程的離散采用有限容積法,空間離散采用二階迎風(fēng)格式,采用標(biāo)準(zhǔn)k-ε湍流模型。采用結(jié)構(gòu)網(wǎng)格進(jìn)行網(wǎng)格劃分,對(duì)隔艙和噴管壁面處的網(wǎng)格進(jìn)行加密。

2.4 數(shù)值仿真結(jié)果

3.1 典型雙脈沖發(fā)動(dòng)機(jī)內(nèi)部流場(chǎng)特性 在無(wú)過(guò)載條件下對(duì)發(fā)動(dòng)機(jī)進(jìn)行了流場(chǎng)模擬,發(fā)動(dòng)機(jī)內(nèi)部流線圖如圖2所示。雙脈沖發(fā)動(dòng)機(jī)燃燒室內(nèi)流場(chǎng)特征:二脈沖工作時(shí),產(chǎn)生的燃?xì)饬鹘?jīng)過(guò)級(jí)間通道流入一脈沖燃燒室,此時(shí)一脈沖燃燒室絕熱層壁面完全暴露在高溫燃?xì)庀?流經(jīng)級(jí)間通道后,截面積突擴(kuò),燃?xì)獍l(fā)生分離流動(dòng),在級(jí)間通道下游形成渦流,在級(jí)間根部形成角渦,同時(shí)在渦流結(jié)束處形成氣流再附著點(diǎn),而后流出噴管。

現(xiàn)代高能推進(jìn)劑中添加有金屬粉劑,發(fā)動(dòng)機(jī)內(nèi)流場(chǎng)的特點(diǎn)會(huì)造成顆粒半徑較小的粒子在一脈沖燃燒室前端區(qū)域高度聚集,惡化一脈沖燃燒室的工作環(huán)境,嚴(yán)重時(shí)導(dǎo)致絕熱層防護(hù)體系失效,出現(xiàn)發(fā)動(dòng)機(jī)燒穿故障。因此,在粒子聚集的區(qū)域,絕熱層需進(jìn)一步加厚設(shè)計(jì),保證燃燒室殼體結(jié)構(gòu)完整性。

圖2 發(fā)動(dòng)機(jī)內(nèi)部流線分布圖

3.2 過(guò)載條件下發(fā)動(dòng)機(jī)內(nèi)流場(chǎng)特性 發(fā)動(dòng)機(jī)在不同過(guò)載條件,分別取軸向過(guò)載x=0、橫向過(guò)載y=20g；軸向過(guò)載x=20g,橫向過(guò)載y=0；軸向過(guò)載x=20g,橫向過(guò)載y=20g,對(duì)發(fā)動(dòng)機(jī)內(nèi)流場(chǎng)進(jìn)行數(shù)值模擬。發(fā)動(dòng)機(jī)內(nèi)部流線如圖3所示,與無(wú)過(guò)載情況相比,流場(chǎng)結(jié)構(gòu)相似,一脈沖燃燒室前端回流區(qū)域大小基本相同。

圖3 不同過(guò)載條件下發(fā)動(dòng)機(jī)內(nèi)部流線圖

在不同過(guò)載條件下,發(fā)動(dòng)機(jī)軸線上壓力分布如圖4所示。與無(wú)過(guò)載條件相比,過(guò)載對(duì)軸向上壓力分布影響較大。在有過(guò)載條件下,一、二脈沖燃燒室整體壓強(qiáng)明顯增加,且節(jié)流損失有所增加,造成發(fā)動(dòng)機(jī)燃燒室殼體強(qiáng)度需求有所增加。在不同過(guò)載條件下發(fā)動(dòng)機(jī)軸線壓力分布對(duì)比可看出,一脈沖燃燒室壓強(qiáng)明顯增加,在單獨(dú)橫向過(guò)載條件下最大。雖然不同過(guò)載條件下一脈沖燃燒室前端回流區(qū)域面積基本相同,但過(guò)載的存在,將造成回流區(qū)域燃?xì)饬髦泄虘B(tài)顆粒對(duì)絕熱層的沖刷速率增加,因此在設(shè)計(jì)時(shí),應(yīng)考慮絕熱層沖刷燒蝕速率增加造成的影響,增加絕熱層厚度或采用多層絕熱層結(jié)構(gòu)。從軸線壓力分布還可以看出,過(guò)載的存在也將影響噴管處絕熱材料的燒蝕率,在設(shè)計(jì)時(shí)應(yīng)加以考慮。

圖4 不同過(guò)載條件下發(fā)動(dòng)機(jī)軸線壓力分布

結(jié)論

本文對(duì)雙脈沖固體火箭發(fā)動(dòng)機(jī)在二脈沖燃燒室工作時(shí)的流場(chǎng)進(jìn)行數(shù)值模擬,研究了過(guò)載對(duì)發(fā)動(dòng)機(jī)內(nèi)流場(chǎng)的影響,通過(guò)對(duì)流線分布、壓強(qiáng)分布和節(jié)流損失特性的分析,得出以下結(jié)論:在導(dǎo)彈飛行過(guò)載條件下,應(yīng)考慮絕熱層沖刷燒蝕速率增加造成的影響,增加絕熱層厚度或采用多層絕熱層結(jié)構(gòu)。