一種擺動噴管的流場和傳熱特性研究

劉宇濤,田小濤,鄧 恒,李志浩,孫曉博

(中國兵器工業(yè)第203研究所,西安 710065)

一種擺動噴管的流場和傳熱特性研究

劉宇濤,田小濤,鄧 恒,李志浩,孫曉博

(中國兵器工業(yè)第203研究所,西安 710065)

針對某球窩式擺動噴管,采用兩方程k-ε湍流模型,在擺動角度為0°和6°的情況下,對噴管內(nèi)部流場進行三維數(shù)值計算;根據(jù)計算的流場,對噴管進行了傳熱數(shù)值模擬。分析了兩種不同擺角對噴管內(nèi)流場和熱防護層的傳熱性能產(chǎn)生的影響。研究結(jié)果表明,該型噴管在小范圍擺動時,噴管內(nèi)流場變化較小,噴管入口處的流場發(fā)生偏轉(zhuǎn);噴管外壁面溫度變化趨勢基本保持一致,內(nèi)壁面溫度存在差異,噴管的熱防護性能滿足發(fā)動機的工作需求。

球窩式擺動噴管;數(shù)值計算;流場;傳熱

0 引言

動力系統(tǒng)作為導(dǎo)彈的重要一環(huán),其發(fā)展水平制約著導(dǎo)彈技術(shù)的進步,推力矢量技術(shù)對于提高導(dǎo)彈的控制性能有著重要的影響。目前的推力矢量技術(shù),按系統(tǒng)產(chǎn)生推力偏向的手段不同,可以大致分為擺動噴管致偏、流體二次噴射以及機械導(dǎo)流阻流等三類[1]。擺動噴管是實現(xiàn)彈用發(fā)動機推力矢量控制的可行方案,在現(xiàn)役導(dǎo)彈上已有成熟的應(yīng)用[2],例如美國的標準-3導(dǎo)彈、戰(zhàn)斧巡航導(dǎo)彈、民兵導(dǎo)彈。擺動噴管提供了矢量推力,然而燃燒室內(nèi)部的流動也隨之改變,噴管的熱防護方案也需要重新設(shè)計。

近年來,在擺動噴管流場計算領(lǐng)域有較多研究,劉君[3]等對嵌入式擺動噴管兩相流內(nèi)流場進行了數(shù)值模擬,研究了噴管擺動角度對發(fā)動機喉道附近引起的流動變化以及內(nèi)表面壓力的影響。黃振宇[4]等三維數(shù)值計算了某擺動噴管在不同擺角下的流動和推力矢量,初步得出了推力矢量隨擺角的變化規(guī)律。周紅梅[5]等通過動網(wǎng)格技術(shù)數(shù)值計算了不同擺角下擺動噴管的內(nèi)流場,得出了隨著擺角的增大,氣流的偏轉(zhuǎn)逐漸明顯的結(jié)論。

在噴管的傳熱計算方面,較多采用先計算噴管內(nèi)流場,根據(jù)流場計算結(jié)果確定噴管傳熱的邊界條件然后進行傳熱計算。文獻[6-7]對噴管內(nèi)壁面施加溫度邊界條件,計算噴管的傳熱特性。文獻[8-9]采用巴茲公式確定燃氣及噴管的對流換熱系數(shù),以此作為邊界條件計算噴管的傳熱特性。傅學(xué)金[10]等通過流固耦合的方法,計算了噴管機構(gòu)瞬態(tài)溫度場,分析了溫度場分布狀態(tài)和初期傳播特點,以及喉襯溫度隨時間變化的規(guī)律。

文中針對某球窩擺動噴管,先計算不同擺動角度情況下噴管的內(nèi)流場,根據(jù)流場計算結(jié)果對噴管內(nèi)邊界施加溫度邊界條件的方法對噴管的傳熱進行三維數(shù)值計算,考核噴管的工作性能,從而對該型噴管進行評價。

1 物理模型與計算方法

1.1 物理模型

球窩擺動噴管采用潛入式方案,一方面能夠減小能量損失,另一方面可以縮短發(fā)動機總長。具體結(jié)構(gòu)見圖1,擺動角度分別為0°和6°。主要結(jié)構(gòu)尺寸見表1。

圖1 球窩擺動噴管幾何模型

長度/mm外徑/mm喉徑/mm出口直徑/mm440440130300

球窩噴管的前端深入到燃燒室中,受到高溫高壓燃氣的作用,因此較多的采用耐燒蝕材料。其中喉襯和前端鎖緊環(huán)采用石墨材料,熱防護采用碳酚醛材料,金屬殼體材料為鋼。

1.2 網(wǎng)格劃分

使用ICEM對流體域進行結(jié)構(gòu)化網(wǎng)格劃分。固體域網(wǎng)格在Workbench中進行劃分,其中擴張段金屬殼體進行了局部加密,網(wǎng)格劃分結(jié)果見圖2。

圖2 網(wǎng)格劃分結(jié)果

1.3 計算方法及邊界條件

使用Fluent計算噴管內(nèi)的流場,計算區(qū)域為噴管段和燃燒室入口段的燃氣流動區(qū)域。由于發(fā)動機在實際工作中建立穩(wěn)定流場的時間很短,初始階段燃氣對噴管的傳熱效應(yīng)可以忽略,噴管內(nèi)部的流場可以視為定常流動。采用RANS方法進行穩(wěn)態(tài)求解,選擇RNGk-ε兩方程模型模擬湍流粘性的影響,計算采用基于壓力的SIMPLE算法進行求解。流體域選擇燃氣作為工質(zhì),進行流場計算,入口選擇壓力入口,給定總壓值為8 MPa,總溫3 200 K,出口選擇壓力出口,設(shè)為一個標準大氣壓,初始溫度300 K,壁面為無滑移絕熱壁面。RNGk-ε模型的k方程和ε方程分別為[11]:

(1)

(2)

流場計算完成之后,將流場外邊界溫度作為噴管內(nèi)表面溫度邊界條件進行非穩(wěn)態(tài)傳熱計算。在噴管內(nèi)部的導(dǎo)熱控制方程為[12]:

(3)

式中:T為溫度;C為材料的比熱容;k為材料的導(dǎo)熱系數(shù);Q′為單位體積的內(nèi)熱源。

計算中,不計燃氣與壁面的輻射換熱,忽略噴管金屬殼體與空氣之間的對流換熱,設(shè)置噴管固體材料初始溫度為300 K,材料的導(dǎo)熱系數(shù)見表2。

表2 材料的導(dǎo)熱系數(shù)

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

2.1 流場分析

選擇噴管流體域的中間截面作為特征面對噴管中的流動進行分析,壓力云圖和馬赫數(shù)云圖分別見圖3、圖4,可以看出,雖然擺動角度不同,但噴管中的流動是相似的。兩個狀態(tài)的出口總壓均約為7.7 MPa,差異小于萬分之一,一方面說明潛入式噴管的流動損失小,另一方面擺動角度為6°時不會對發(fā)動機的能量發(fā)揮帶來較為明顯的影響,初步驗證了該結(jié)構(gòu)方案的可行性。

圖5為噴管入口處外邊界速度云圖,圖6為入口處中間截面速度云圖,從圖中可以看出,擺動角度為0°,噴管入口處氣流沒有偏轉(zhuǎn),擺動角度為6°時,噴管入口處氣流發(fā)生偏轉(zhuǎn),變?yōu)榉菍ΨQ的流場,一側(cè)的流速大于另一側(cè)的流速,并且入口段的最大流速稍大于0°時的最大流速,說明隨著噴管的擺動,噴管入口段流場會發(fā)生偏轉(zhuǎn)和局部增速,這對噴管入口段的熱防護可能會造成一定影響,在設(shè)計時應(yīng)當予以考慮。

圖3 壓力云圖

圖4 馬赫數(shù)云圖

圖5 噴管入口處外邊界速度云圖

圖6 噴管入口處中間截面速度云圖

2.2 傳熱分析

該型噴管熱防護層應(yīng)保證在10 s工作時間內(nèi)具有良好的熱防護能力,考核指標為金屬外壁面溫度不大于300℃。選定噴管工作時的兩個極限狀態(tài),即擺動角度為0°和6°時,對熱防護層的防護能力進行考核。

圖7 溫度分布云圖

不同時刻噴管結(jié)構(gòu)溫度云圖見圖7,可以看出,在10 s的工作時間內(nèi),噴管外壁面溫度沒有大幅度升高,均小于80℃,說明噴管的熱防護層滿足工作要求。在開始工作之后,噴管內(nèi)部的喉襯溫度迅速升高,這是因為石墨喉襯的導(dǎo)熱率較高,且處在接觸燃氣溫度最高的區(qū)域,因此傳熱劇烈導(dǎo)致迅速升溫。在10 s工作時間結(jié)束之后,0°和6°兩種擺動角度所對應(yīng)的喉襯最低溫度分別為300℃和430℃,且6°工況的溫度呈現(xiàn)非對稱式的分布,這是因為噴管擺動之后,喉襯一側(cè)與燃氣的接觸面積增大,另一側(cè)接觸面積減小,面積增加的一側(cè)熱傳導(dǎo)加劇,從而導(dǎo)致喉襯溫度比0°時有所提高。

繼續(xù)計算至30 s,監(jiān)測金屬殼體外壁面最高溫度,具體數(shù)值見表3。由表3可以看出,兩種擺動角度條件下外壁面溫度隨時間變化基本一致,工作時間超過10 s后,噴管最高溫度基本呈線性趨勢迅速增加。擺動角度為0°和6°的兩種情況,對噴管的熱防護結(jié)構(gòu)整體造成的影響較小。從表中也可以看到,當工作時間增加到25～30 s時,外壁面最高溫度達到300℃,最高溫度出現(xiàn)在噴管出口區(qū)域,前端熱防護層的設(shè)計仍有較大的余量,但是并不應(yīng)該減薄前端熱防護層,主要是因為擺動噴管的分離線位于燃燒室內(nèi),緊鄰高溫高壓燃氣,密封困難,如果在發(fā)動機工作過程中,高溫燃氣從分離線進入配合間隙中,將會對熱防護帶來嚴重的問題,因此,前端的熱防護設(shè)計應(yīng)有較大的余量,以保證擺動噴管的正常工作。

表3 不同時刻外壁面的最高溫度

3 結(jié)論

1)噴管在0°到6°小擺動角度的情況下流場相似,對發(fā)動機能量的發(fā)揮帶來的差異較小;但是噴管入口的流場會發(fā)生偏轉(zhuǎn)和速度增大,所以在噴管設(shè)計時,應(yīng)當考慮入口處流場的偏轉(zhuǎn)和流速的增大對噴管產(chǎn)生的影響。

2)在10 s工作時間內(nèi),兩種擺角下噴管的外壁面溫度沒有超過規(guī)定溫度,內(nèi)壁面溫度存在差異,初步驗證了噴管的熱防護設(shè)計能滿足發(fā)動機工作的需求。但是對于較長時間工作的情況,噴管的熱防護需要加強。

3)文中傳熱計算中選擇的溫度邊界條件,可以很大程度上降低求解的難度和計算的時間,精度上也可以滿足工程上噴管初步設(shè)計的要求,在下一步的工作中,對于噴管傳熱計算將采用精度更高的流固耦合的方法,對現(xiàn)有噴管的熱防護進行優(yōu)化。

[1] 陳汝訓(xùn). 固體火箭發(fā)動機設(shè)計與研究: 下 [M]. 北京: 中國宇航出版社, 1991: 230-232.

[2] 高峰, 唐勝景, 師嬌. 推力矢量技術(shù)在導(dǎo)彈上的應(yīng)用 [J]. 飛航導(dǎo)彈, 2010(12): 52-59.

[3] 劉君, 黃琳, 夏智勛. 嵌入式擺動噴管兩相流內(nèi)流場數(shù)值模擬 [J]. 固體火箭技術(shù), 2003, 26(4): 18-21.

[4] 黃振宇, 徐文燦, 毛鴻羽. 噴管流場及其推力矢量的數(shù)值計算 [J]. 兵工學(xué)報, 2000, 21(1): 31-34.

[5] 周紅梅, 于勝春, 高劼, 等. 動網(wǎng)格技術(shù)在具有擺動噴管發(fā)動機流場數(shù)值模擬中的應(yīng)用 [J]. 海軍工程學(xué)院學(xué)報, 2009, 24(1): 17-19.

[6] 沈勇. 碳/碳固體火箭發(fā)動機柔性噴管熱機熱結(jié)構(gòu)耦合作用研究 [D]. 武漢: 華中科技大學(xué), 2004.

[7] EATON Andrew M, MATHIAS Edward C. Simulating heat transfer to a solid rocket motor nozzle-case joint thermal barrier: AIAA 2000-3807[R].2000.

[8] 高峰, 王建輝, 肖紅雨. 高過載下長尾噴管整體結(jié)構(gòu)傳熱數(shù)值模擬 [J]. 彈箭與制導(dǎo)學(xué)報, 2014, 34(4): 111-114.

[9] 韓珺禮, 蒙上陽, 楊軍輝, 等. 固體火箭發(fā)動機噴管熱防護設(shè)計 [J]. 計算機輔助工程, 2013(1): 160-163.

[10] 傅學(xué)金, 鄧恒, 寇軍強, 等. 復(fù)合噴管結(jié)構(gòu)溫度場流固耦合仿真分析 [J]. 固體火箭技術(shù), 2014, 37(6): 781-786.

[11] 張兆順, 崔桂香, 許春曉. 湍流理論與模擬 [M]. 北京: 清華大學(xué)出版社, 2005: 209-214.

[12] 胡曉平. 傳熱傳質(zhì)分析 [M]. 長沙: 國防科技大學(xué)出版社, 2012: 8-9.

Research on the Flow Field and Heat Transfer Characteristics of a Gimbaled Nozzle

LIU Yutao,TIAN Xiaotao,DENG Heng,LI Zhihao,SUN Xiaobo

(No.203 Research Institute of China Ordnance Industries, Xi’an 710065, China)

For a ball-and-socket type gimbaled nozzle, when the swinging angle were 0 degree and 6 degree, the three-dimensional numerical calculation of nozzle interior flow field was carried out by using two equationk-εturbulence model. According to the calculated flow field, the heat transfer numerical simulation of nozzle was implemented. The effect of two different swinging angles on the heat transfer performance of nozzle interior flow field and thermal protection shield was analyzed. The results showed that the nozzle interior flow field had small variation when the swinging angel range was small, but the flow field at the entrance of the nozzle was deflected. The variation trend of the temperature on the outer wall of the nozzle was basically the same, while the temperature on the inner wall was different. The thermal protection performance of the nozzle met the engine work requirements.

ball-and-socket type gimbaled nozzle；numerical simulation；flow field；heat transfer

2016-11-29

劉宇濤(1986-),男,陜西乾縣人,助理工程師,碩士研究生,研究方向:固體火箭發(fā)動機設(shè)計。

V430

A