王書賢，王 山

(西安航空學(xué)院 飛行器學(xué)院，西安 710077)

可調(diào)收斂-擴張噴管周向內(nèi)型面對流場的影響研究

王書賢，王 山

(西安航空學(xué)院 飛行器學(xué)院，西安 710077)

可調(diào)噴管由調(diào)節(jié)片和密封片組成，噴管周向內(nèi)型面不是嚴格圓形而是多邊形，周向內(nèi)型面對噴管流場的影響會進一步影響噴管性能。以圓形和三種多邊形周向內(nèi)型面收斂-擴張噴管為研究對象，給以相同軸向尺寸和進口條件，利用流體計算分析軟件FLUENT對噴管內(nèi)流場進行了計算分析，得到的影響規(guī)律為可調(diào)噴管結(jié)構(gòu)的優(yōu)化設(shè)計提供參考。研究結(jié)果表明：可調(diào)收斂-擴張噴管調(diào)節(jié)片數(shù)越多，周向內(nèi)型面越接近于圓形，噴管流動損失越小。

可調(diào)收擴噴管；周向內(nèi)型面；流場分析

0 引言

要使具有噴氣發(fā)動機的飛行器在寬廣的飛行條件下維持良好的性能和較好的機動性能，具有多彈道下有效的工作能力，先進的排氣噴管是必不可少的，可調(diào)噴管就是當(dāng)代出現(xiàn)的先進排氣噴管之一?？烧{(diào)噴管由調(diào)節(jié)片和密封片組成，噴管周向內(nèi)型面不同于固定噴管光滑的圓形而是規(guī)則的多邊形。周向內(nèi)型面設(shè)計的好壞以及噴管在工作狀態(tài)下能否按要求保持相應(yīng)的周向內(nèi)型面，影響著尾噴管的工作效率與性能[1-5]。

為了得到周向內(nèi)型面對噴管內(nèi)流場的影響規(guī)律，本文以某固沖發(fā)動機為研究應(yīng)用背景，對同一軸向尺寸，周向內(nèi)型面分別為圓形和多邊形(邊數(shù)取8，12，15)的可調(diào)收斂-擴張噴管進行建模，分析周向內(nèi)型面對噴管內(nèi)溫度、壓強、流速以及流動損失的影響規(guī)律。

1 計算模型

1.1 幾何結(jié)構(gòu)及網(wǎng)格劃分

噴管進出口直徑均為185mm，喉部直徑81.1mm，噴管總長250mm，喉部距離進口90mm。對圓形內(nèi)型面噴管采用軸對稱的二維網(wǎng)格劃分，而對多邊形內(nèi)型面噴管分別取1/4、1/6和1/5進行三維網(wǎng)格劃分。由于壁面附近速度梯度大，相應(yīng)的壓強和溫度也變化很快，因此近壁面網(wǎng)格較密，向軸線方向逐漸過渡變疏。幾何結(jié)構(gòu)及網(wǎng)格劃分情況如圖1所示。

(a) 圓形內(nèi)型面 (b) 八邊形內(nèi)型面

(c) 十二邊形內(nèi)型面 (d) 十五邊形內(nèi)型面圖1 模型網(wǎng)格

1.2 流動性質(zhì)

某固沖發(fā)動機進口靜壓P0=0.6Mpa，進口靜溫T0=2000K，出口背壓Pb=0.05Mpa，氣體質(zhì)量流量Wgi=1kg/s，比熱比k=1.3，氣體常數(shù)R=430.4J/kg·K。進一步求得其它相關(guān)物理參數(shù)，如表1所示。

表1 其它相關(guān)物理參數(shù)

1.3 控制方程及求解方法

綜合考慮研究對象的各項流動條件，本文選取為標準k-ε模型[7-8]?？刂品匠坛|(zhì)量守恒方程、動量守恒方程和能量守恒方程外，還包括k方程、ε方程，控制方程的通用表達式如下：

(1)

具體控制方程如表2所示。本文采用壓力耦合隱式亞松弛迭代求解控制方程組。

表2 控制方程

2 計算結(jié)果分析

由于在計算分析時圓形內(nèi)型面噴管采用的二維建模，多邊形采用的三維建模，為了方便計算結(jié)果的分析，在云圖分析時，對三維模型也同時取軸向截面進行分析。

2.1 壓強分析

通過對不同內(nèi)型面可調(diào)收擴噴管靜壓分布云圖的觀察，如圖2所示(以圓形內(nèi)型面為例)可以看出，噴管從進口到出口靜壓逐漸降低，且不同周向內(nèi)型面噴管整體靜壓分布基本相同，僅在數(shù)值上有所差異。由圖3可以看出圓形內(nèi)型面噴管的靜壓數(shù)值在整個過程都要低于多邊形內(nèi)型面噴管，多邊形內(nèi)型面噴管之間差異不大。再通過對其在出口靜壓數(shù)值上的變化的分析(見表3)，可知圓形內(nèi)型面噴管出口靜壓最小，多邊形內(nèi)型面噴管出口靜壓隨邊數(shù)的增加而減小。

圖2 圓形內(nèi)型面可調(diào)收擴噴管靜壓分布云圖

圖3 不同內(nèi)型面噴管靜壓變化趨勢圖

通過對不同內(nèi)型面可調(diào)收擴噴管總壓分布云圖的觀察，如圖4所示(以圓形內(nèi)型面為例)可以看出，不同內(nèi)型面噴管總壓分布也基本相同，收斂段總壓基本不變，喉部過后總壓變化明顯。

喉部過后總壓變化的原因有三點：第一，壁面附近由于粘性摩擦存在沿程損失，總壓較低，且沿程損失與速度的平方成正比，因此在喉部之后氣流達到超聲速，總壓損失明顯。

第二，湍流脈動也會造成一定的總壓損失，氣流速度越高湍流脈動越強烈，總壓損失越大，因此喉部之后沿軸向總壓降低。

第三，噴管工作在超臨界過膨脹狀態(tài)，氣流通過出口的斜激波將出口壓強增加至外部環(huán)境壓強，斜激波以噴管出口壁面為附著點，對出口壁面附近氣流有一定干擾，使總壓有所損失。

將不同內(nèi)型面可調(diào)收擴噴管沿軸向總壓數(shù)值變化趨勢進行歸納在圖5中對比呈現(xiàn)，從圖中可以看出喉部過后總壓數(shù)值開始下降，圓形內(nèi)型面噴管的總壓下降量最少，八邊形內(nèi)型面噴管下降量最大，十二、十五邊形內(nèi)型面噴管則是比較接近的。從噴管出口總壓的具體數(shù)值(見表3)可以看出，圓形內(nèi)型面噴管出口總壓最大，多邊形內(nèi)型面噴管出口總壓隨邊數(shù)增加而增加。

圖4 圓形內(nèi)型面可調(diào)收擴噴管總壓分布云圖

圖5 不同內(nèi)型面噴管總壓變化趨勢圖

2.2 速度分析

通過對不同內(nèi)型面可調(diào)收擴噴管總壓分布云圖的觀察，如圖6所示(以圓形內(nèi)型面為例)可以看出，不同內(nèi)型面的噴管軸向速度分布基本相同，燃氣在噴管內(nèi)持續(xù)加速，在喉部達聲速，噴管工作在超臨界過膨脹狀態(tài)。將不同內(nèi)型面可調(diào)收擴噴管沿軸向總壓數(shù)值變化趨勢進行歸納在圖7中，可以看出整個過程中圓形內(nèi)型面噴管速度數(shù)值上要大于多邊形內(nèi)型面噴管的數(shù)值。出口氣流速度的具體數(shù)值(見表3)表明：圓形內(nèi)型面噴管出口速度最大，多邊形內(nèi)型面噴管出口速度隨邊數(shù)的增加而增大。

圖6 圓形內(nèi)型面可調(diào)收擴噴管速度分布云圖

圖7 不同內(nèi)型面的噴管速度變化趨勢圖

2.3 溫度分析

總溫在噴管整個流動過程中基本不變。靜溫變化如圖8所示(以圓形內(nèi)型面為例)。不同內(nèi)型面噴管的靜溫分布也基本相同，同樣也是數(shù)值上有所差異，對不同內(nèi)型面可調(diào)收擴噴管沿軸向靜溫數(shù)值變化趨勢進行歸納在圖9中對比呈現(xiàn)，可以看出整個過程中圓形內(nèi)型面噴管靜溫數(shù)值要小于多邊內(nèi)型面噴管的靜溫數(shù)值。噴管出口靜溫(見表3)，與靜壓規(guī)律相同，圓形內(nèi)型面噴管出口靜溫最小，多邊形內(nèi)型面噴管出口靜溫隨邊數(shù)增加而減小。

圖8 圓形內(nèi)型面可調(diào)收擴噴管靜溫分布云圖

圖9 不同內(nèi)型面的噴管靜溫變化趨勢圖

2.4 小結(jié)

將以上計算結(jié)果列于表3進行對比分析，可以得出如下規(guī)律：在進口條件一致的情況下，隨著多邊形周向內(nèi)型面邊數(shù)增加出口總壓增大，壓力損失減??；出口靜壓減小，出口流速增大，動能轉(zhuǎn)化增大；出口總溫不變，出口靜溫減小，熱能損失減小。即對于可調(diào)收斂-擴張噴管調(diào)節(jié)片數(shù)越多，周向內(nèi)型面越接近圓形面，噴管流動損失越小。

表3 不同內(nèi)型面噴管出口參數(shù)

*注：( )為其它情況物理參數(shù)以圓形物理參數(shù)為基準的相對變化百分比

3 結(jié)論

(1)噴管周向內(nèi)型面變化不影響軸向物理參數(shù)的總體變化規(guī)律，即總壓靜壓下降、總溫不變、靜溫下降、速度增大。

(2)噴管周向內(nèi)型面變化影響噴管流動損失，周向內(nèi)型面越接近于圓形，噴管內(nèi)流動損失越小，熱能向動能轉(zhuǎn)化效率越高。

(3)在可調(diào)收斂-擴張噴管的結(jié)構(gòu)設(shè)計中應(yīng)在保證結(jié)構(gòu)可調(diào)及有效熱防護的基礎(chǔ)上，盡量采用較多的調(diào)節(jié)片和密封片，以減少流動損失。

[1] 葛李陽,張國偉,張國棟.基于FLUENT的火箭發(fā)動機噴管流場數(shù)值模擬[J].機電技術(shù),2013,36(3):13-15.

[2] 李存杰.空氣噴氣發(fā)動機可調(diào)尾噴管的幾種方案[J].飛航導(dǎo)彈,1997(6):48-50.

[3] LOPEZ-DIEZ A,RODRIGUEZ M,SEGUNDO M S,et al.Aplications of Variable Area Nozzles for Civil Airplanes[C].34th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit,1998:1-7.

[4] 賴謀榮.沖壓發(fā)動機可調(diào)噴管流場的數(shù)值模擬[D].西安:西北工業(yè)大學(xué),2002.

[5] 周唯陽.串聯(lián)布局TBCC可調(diào)噴管的設(shè)計、仿真與實驗研究[D].南京:南京航空航天大學(xué),2012.

[6] 王新月.氣體動力學(xué)基礎(chǔ)[M].西安:西北工業(yè)大學(xué)出版社,2006:324.

[7] 任玉新,陳海昕.計算流體力學(xué)基礎(chǔ)[M].北京:清華大學(xué)出版社,2006:27-30.

[8] 王福軍.計算流體動力學(xué)分析[M].北京:清華大學(xué)出版社,2008:185.

[責(zé)任編輯、校對：李琳]

Abstract：The variable nozzle is composed of flaps and seals,and the inner contour of the nozzle is not precisely circular but polygonal.The inner contour has influence on flow and nozzle performance.The circular and three kinds of polygonal circumferential inner contour convergent-divergent nozzles are taken as the research object.The same axial dimension and inlet condition are used.The software FLUENT is used to analyze the internal flow field.The influence of different inner contours on the flow field of the variable convergent-divergent nozzle is given,which provides a reference for the optimized design of the variable nozzle structure.The results show that the bigger the number of variable nozzles is,and the more the circumferential inner surface is approximate to circular shape,and the smaller the flow loss from the inlet to the outlet is.

Keywords：variable convergent-divergent nozzle;circumferential inner contour;fluid analysis

StudyofCircumferentialInnerContourInfluenceonVariableConvergent-DivergentNozzleInternalFlowField

WANGShu-xian,WANGShan

(School of Aircraft,Xi′an Aeronautical University,Xi′an 710077,China)

V435

A

1008-9233(2017)05-0003-05

2017-07-05

西安市科技計劃項目(CXY1518(5))

王書賢(1977-)，女，陜西西安人，博士，副教授，主要從事發(fā)動機流動與熱結(jié)構(gòu)研究。