張志浩，許曉勇，田 原，丁兆波，陳旭揚

氫氧推力室噴管超聲速氣膜冷卻參數(shù)研究

張志浩，許曉勇，田 原，丁兆波，陳旭揚

（北京航天動力研究所，北京，100076）

某型氫氧推力室噴管延伸段采用超聲速氣膜冷卻，在面積比35處引入渦輪排氣作為冷卻氣體。通過采用數(shù)值模擬的方法，研究了唇高、吹風比、冷卻劑流量和靜壓比等參數(shù)對氣膜冷卻效果和比沖的影響。結(jié)果表明：隨著唇高的增大，推力室的比沖和推力略微降低，而冷卻效率和壁溫幾乎不變；在主射流壓力匹配且射流量一定的條件下，吹風比增大可以輕微地提高氣膜冷卻效果和發(fā)動機比沖；在壓力匹配且吹風比一定的條件下，射流量增加可以提高氣膜冷卻效果；在射流量一定的條件下，主射流壓力匹配時，氣膜冷卻效果最佳，發(fā)動機比沖最高。

氫氧推力室噴管；超聲速氣膜冷卻；數(shù)值模擬；壓力匹配

0 引言

某型氫氧推力室噴管延伸段采用超聲速氣膜冷卻，冷卻氣體來自渦輪排氣。渦輪排氣相對推力室主流燃氣溫度較低，并以超聲速進入噴管，緊貼壁面流動，形成一層低溫氣體保護膜，這股冷卻介質(zhì)將高溫燃氣與壁面隔開，避免高溫燃氣與壁面直接進行對流換熱，并且，冷卻氣膜還可以帶走高溫燃氣的輻射熱量，以達到保護噴管的效果。文獻[1～6]介紹了噴管延伸段采用超聲速氣膜冷卻方案，噴管延伸段采用單壁高溫合金輻射冷卻噴管，同時引入冷卻氫或渦輪排氣進行氣膜冷卻，內(nèi)壁噴涂隔熱涂層。超聲速氣膜冷卻比亞聲速氣膜有更高的冷卻效率，因為它有更薄的邊界層和更小的湍流尺度，因此有更低的混合率[7]。在對高速氣流的氣膜冷卻研究中，應用修正過的不可壓縮流動的關(guān)聯(lián)式[8]來計算冷卻效率，然而，這些關(guān)聯(lián)式并不能成功預測實驗的測量結(jié)果，這是因為超聲速流動比不可壓縮流動較為復雜，在超聲速流動中，動量方程和能量方程強烈耦合，并伴有激波出現(xiàn)。為了研究某型氫氧推力室噴管超聲速氣膜冷卻的影響因素，本文研究了唇高、吹風比、冷卻劑流量和靜壓比對氣膜冷卻效果的影響，為工程設(shè)計提供了參考。

1 計算方法與模型

1.1 物理模型

某型氫氧發(fā)動機噴管面積比為80，渦輪排氣在噴管面積比35處引入噴管內(nèi)排放，形成氣膜，保護下游噴管，其中，渦輪排氣通過具有收斂擴張段的環(huán)縫式噴嘴引入。推力室噴管氣膜冷卻示意如圖1所示，渦輪排氣引入的局部示意如圖2所示。

圖1 噴管氣膜冷卻示意

圖2 氣膜引入局部示意

1.2 計算模型

假定氫氧推力室噴管內(nèi)主燃氣和渦輪排氣均為水蒸汽和氫氣的混合氣體，表1給出了主燃氣和渦輪排氣的工況參數(shù)，在下文所有算例中，主燃氣參數(shù)不變，而渦輪排氣參數(shù)根據(jù)需要改變總壓或流量，以討論所研究參數(shù)的影響。

表1 工況參數(shù)

采用ΙCEM軟件對推力室噴管劃分了結(jié)構(gòu)網(wǎng)格，對壁面附近和流動梯度較大的位置進行網(wǎng)格加密處理。使用Fluent軟件基于密度的求解器進行求解計算，采用SST kω?湍流模型，推力室入口與渦輪排氣入口均設(shè)為壓力入口邊界條件，給定入口的總壓、總溫、各氣體組分的濃度、水力直徑和湍流強度等參數(shù)，所有壁面設(shè)為絕熱無滑移條件?；旌蠚饷芏劝凑绽硐霘怏w考慮，各組分氣體粘性系數(shù)、定壓比熱和熱傳導系數(shù)由物理性質(zhì)查詢軟件NΙST_refprop 8和NASA報告[9]得到。

超聲速氣膜冷卻效率定義為[10]

式中 Taw，Tri分別為絕熱壁面溫度和冷卻氣膜入射處的氣膜恢復溫度；rT∞為主流的恢復溫度。

恢復溫度是指氣流在絕熱的固體表面上被滯止到零速度時的溫度。一般認為在絕熱條件下，主流恢復溫度rT∞與無膜冷卻的絕熱壁面溫度相等。

1.3 湍流模型

計算采用SSTkω?兩方程模型，集合了kε?模型和kω?模型的特點，在邊界層內(nèi)采用kω?模型，在邊界層邊緣和邊界層以外的主流燃氣區(qū)采用kε?模型（kω?形式），邊界層和主流燃氣區(qū)通過一個混合函數(shù)來過渡[11]。SSTkω?模型的輸運方程與標準的kω?模型形式相似，即：

式中 ρ為密度；k為湍動能；ui為分速度；G~k為由于平均速度梯度引起的湍動能生成項；Gω為ω的生成項；Γk，Γω分別為k和ω的擴散效應；Yk，Yω分別為由于湍流引起的k和ω的耗散；Dω為交叉擴散項；Sk，Sω分別為源項。

湍流粘性定義為

式中 S為應變率大?。籉1為混合函數(shù)。

湍流普朗特數(shù)定義為

式中 F2為混合函數(shù)。

模型中的部分常數(shù)設(shè)為：σk,1=1.176，σω,1=2，σk,2=1，σω,2=1.168；a1=0.31，對于高雷諾數(shù)流動為α*=1。

2 計算結(jié)果分析

2.1 唇高的影響

唇高h是一個不可避免的結(jié)構(gòu)參數(shù)，本文分別選取了h=3 mm，h=8 mm，h=11 mm進行計算，3種算例下的渦輪排氣射流量均為20.88 kg/s，噴嘴喉部高度和面積比不變，分別為13.6 mm和1.75。圖3給出了沿軸向壁面絕熱溫度，圖4為3種唇高h的冷卻效率。

圖3 不同唇高h下沿軸向絕熱壁面溫度曲線

圖4 不同唇高h下氣膜冷卻效率曲線

由圖3、圖4可知，在x/s=120之前，壁面溫度在600～620 K之間，壁溫緩慢下降，冷卻效率緩慢上升，并大于1，主要是由于氣膜噴出噴嘴后，氣膜在噴管內(nèi)膨脹，溫度進一步降低，使得壁溫比出口處氣膜的恢復溫度還低，所以使得冷卻效率大于1；此外，在氣膜出口附近，即x/s=75附近，絕熱壁面有一點溫度迅速下降又迅速回升，冷卻效率在這一點也迅速升高又迅速降低，分析此點為唇高部位所發(fā)出的膨脹波在壁面上的反射點，由于膨脹波的緣故，使得附近壁面溫度降低。

3種算例下絕熱壁面溫度的變化趨勢基本一致。在x/s=120之前，三者的壁面溫度幾乎相等。在噴管尾端，唇高h=3 mm的壁面溫度比h=8 mm和h=11 mm的壁面溫度低20 K左右，h=8 mm和h=11 mm的壁面溫度基本一致，冷卻效率變化規(guī)律與溫度曲線一致，三者的冷卻效率也相等。經(jīng)過此點后，唇高h=3 mm的冷卻效率略高于唇高h=8 mm和h=11 mm的冷卻效率，h=8 mm和h=11mm的冷卻效率很接近。

不同唇高下噴嘴附近的壓力分布如圖5所示。的波系的強弱也將不同，h=11 mm時所產(chǎn)生的波系比h=3 mm時強烈得多，唇高越大，有更多的氣膜和主流燃氣改變流動方向以填補主射流之間的間隙，主射流之間的摻混就越多，從而導致冷卻效率降低。

表2給出了3種唇高對推力室推力、比沖的影響。

由表2可知，推力室的比沖和推力均隨著唇高的增大而降低，說明唇高越大，唇部產(chǎn)生的波系越強，對噴管流場的擾動就越大，降低了噴管性能。然而，唇高從h=3 mm增加到h=11 mm時，推力室比沖減小不到0.06%。因此，綜合考慮唇高對氣膜冷卻效果和推力性能的影響，唇部尺寸的選取可以在保證結(jié)構(gòu)的基礎(chǔ)上，盡量減小唇高。

2.2 壓力匹配條件下吹風比的影響

吹風比λ又叫密流比，是氣膜冷卻研究中一個重要的參數(shù)，其定義為射流與主流單位面積質(zhì)量流量之比的大小，由氣動函數(shù)[12]推出：

式中 下標c，∞分別為氣膜和主流；k為比熱比；M為分子量。

由式（7）可知，吹風比主要與主/射流的靜壓比pc/p∞、總溫比Tc*/T∞*及主/射流馬赫數(shù)、分子量、比熱比有關(guān)。

為了探究吹風比的影響，保持射流量不變，通過改變射流總壓和噴嘴結(jié)構(gòu)以改變吹風比，并保證主射流壓力匹配，此時，射流馬赫數(shù)也隨之變化。表3給出了壓力匹配條件下不同吹風比的計算參數(shù)，其中也列出了每種情況下小噴嘴的喉部高度和面積比。

由圖5可知，由于唇高不同，因此唇高處所發(fā)出

表3 壓力匹配不同吹風比主射流計算參數(shù)

圖6給出了壓力匹配條件下不同吹風比的壁面溫度分布情況。

圖6 壓力匹配下不同吹風比的軸向壁面溫度分布曲線

由圖6可知，隨著吹風比和射流馬赫數(shù)的增大，壁面溫度有減小的趨勢，說明在流量一定并保證壓力匹配的情況下，吹風比增大可以輕微地提高冷卻效果，其原因是隨著吹風比的增大，射流馬赫數(shù)也隨之增大，氣膜靜溫隨之降低，因此冷卻效果提高。

表4給出了在壓力匹配情況下不同吹風比的發(fā)動機推力與比沖的關(guān)系。

表4 壓力匹配下不同吹風比的發(fā)動機推力及比沖

由表4可知，在壓力匹配的情況下，隨著吹風比的提高，發(fā)動機推力及比沖略有提高，其原因是在射流量一定的情況下，吹風比和射流馬赫數(shù)同時變大，氣膜速度也變大，氣膜部分提供的推力由此變大，所以發(fā)動機推力和比沖也越大。

2.3 壓力匹配條件下冷卻劑流量的影響

為保持入射總壓不變，改變射流量，通過改變噴嘴結(jié)構(gòu)以保證主射流壓力匹配，以此實現(xiàn)吹風比不變，分析冷卻效果隨射流量的變化情況。表5給出了不同射流量下的計算參數(shù)。

表5 不同射流量下主射流計算參數(shù)

圖7給出了3種射流量下的壁面溫度分布情況。由圖7可知，最高壁溫隨射流量增大而減小，說明冷卻劑流量增加可以提高冷卻效果。

圖7 不同射流量下的軸向壁面溫度

圖8 給出了3種射流量下軸向壁面氫氣摩爾分數(shù)分布情況。

圖8 不同射流量下軸向壁面氫氣摩爾分數(shù)

由圖8可知，氣膜有效冷卻長度隨射流量的增加而增加，即說明射流量越大，氣膜厚度越厚，主流燃氣越難以完全穿透氣膜，緊貼壁面的氣膜維持的距離越長。

2.4 壓力不匹配的影響

為保持入射總壓和射流量不變，改變主射流靜壓比，觀察壓力不匹配對冷卻效果和發(fā)動機性能的影響，表6給出了壓力不匹配條件下的計算參數(shù)。圖9給出了壓力不匹配條件下的壁面溫度分布情況。

表6 壓力不匹配條件下的計算參數(shù)

圖9 壓力不匹配條件下軸向的壁面溫度

靜壓比為0.98時可認為主射流達到壓力匹配，由圖9可見，3種情況下，有效冷卻長度基本相等，而超過有效冷卻長度后，靜壓比為0.98比靜壓比為0.6的壁面溫度最大低約35 K，靜壓比為0.98比靜壓比為1.6時的壁面溫度最大低20 K左右。說明不管氣膜入射壓力高于主流還是低于主流的當?shù)貕毫?，壓力不匹配均會略微降低氣膜冷卻效果。

圖10為速度矢量與氫氣摩爾分數(shù)。

圖10 速度矢量與氫氣摩爾分數(shù)

圖10 中箭頭代表速度矢量，氣膜區(qū)與主流區(qū)之間的兩條實線劃分出了氣膜與主流的摻混區(qū)。由圖10可見，靜壓比為1.6時，由于氣膜靜壓大于當?shù)刂髁黛o壓，氣膜進入噴管后向主流膨脹，加劇了主流和氣膜的擴散摻混；靜壓比為0.98時，氣膜與主流壓力基本相等，氣膜進入噴管后，基本與主流沿切向平行流動，它們之間的混合層比較??；靜壓比為0.6時，氣膜靜壓低于主流靜壓，氣膜進入噴管后被主流壓縮，這同樣加快了主流與氣膜的摻混。所以，壓力不匹配導致冷卻效果降低是由于加劇了主射流間的摻混。

表7給出了主射流壓力不匹配下的發(fā)動機推力及比沖。由表7可見，發(fā)動機推力和比沖由大到小，依次是靜壓比0.98（壓力匹配）、靜壓比0.6（氣膜過膨脹）和靜壓比1.6（氣膜欠膨脹）。說明主射流壓力匹配的情況下，發(fā)動機性能最優(yōu)；氣膜過膨脹和欠膨脹都會降低發(fā)動機性能，但欠膨脹情況下對發(fā)動機性能最不利，其原因是由于氣膜靜壓大于主流壓力時，會對主流造成一定程度的壓縮，主流加速膨脹受到影響，使發(fā)動機性能有所降低，而當氣膜過膨脹時，主流壓縮氣膜，影響氣膜的加速膨脹，但推力室推力主要來自主流，所以性能降低不明顯。

表7 壓力不匹配下的發(fā)動機推力及比沖

3 結(jié) 論

本文采用數(shù)值模擬的方法，應用SST湍流模型研究了某型氫氧推力室噴管超聲速氣膜冷卻的影響因素，得出如下結(jié)論：

a）隨著唇高的增大，壓力波動越強，推力室的比沖和推力隨之降低，而冷卻效率和壁溫幾乎不變，但唇高從h=3 mm增加到h=11 mm，推力室比沖減小不到0.06%，所以唇部尺寸的選取可以在保證結(jié)構(gòu)的基礎(chǔ)上，盡量減小唇高；

b）在主射流壓力匹配且射流量一定的條件下，吹風比增大可以輕微地提高氣膜冷卻效果和發(fā)動機比沖；

c）在壓力匹配且吹風比一定的條件下，射流量增加可以提高氣膜冷卻效果；

d）在射流量一定的條件下，主射流壓力匹配時，氣膜冷卻效果最佳，發(fā)動機比沖最高。

[1] Caporicci M, Eriksson L E, Onofri M, et al. Advanced nozzle technologies for Ariane 5 Vulcain engine[R]. AΙAA-94-3263, 1994.

[2] Grosdemange H, Haggander J, Mattstedt T. Advanced nozzle extension demonstration for Vulcain Mark ΙΙ application[R]. AΙAA-95-2537, 1995.

[3] 中國航天工業(yè)總公司. 世界導彈與航天發(fā)動機大全[M]. 北京: 軍事科學出版社, 1999.

[4] Morris C Ι, Ruf J H. Validation of supersonic film cooling modeling for liquid rocket nozzle applications[R]. Submitted to The 46th AΙAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, M10-0883, 2010.

[5] Watanabe Y, Sakazume N, Tsuboi M. LE-7A engine nozzle problems during the transient operations[R]. AΙAA 2002-3841, 2002.

[6] Ogawara A, Onga T, Manako T. LE-X Japanese next liquid booster engine[R]. AΙAA 2008-4665, 2008.

[7] Aupoix B, Mignosi A, Viala S. Experimental and numerical study of supersonic film cooling[R]. AΙAA, 1998, 36(6): 915-923.

[8] Juhany K A, Hunt M L, Sivo J M. Ιnfluence of injectant Mach number and temperature on supersonic film cooling. Journal of thermophysics and heat transfer[R]. AΙAA, 1994, 8(1): 59-67.

[9] Mcbride B J, Gordon S, Reno M A. Cofficients for calculating thermodynamic and transport properties of individual species[R]. NASA, N94-17624, 1993.

[10] Konopka M, Meinke M, Schroder W. Large-eddy simulation of supersonic film cooling at laminar and turbulent injection[R]. AΙAA 2011-2250, 2011.

[11] 閻超. 計算流體力學方法及應用[M]. 北京: 北京航空航天大學出版社, 2006.

[12] 任漢芬, 劉國球, 等. 液體火箭發(fā)動機原理[M]. 北京: 中國宇航出版社, 1993.

Parameters Study of Supersonic Film Cooling of LOX/LH2Thrust Chamber Nozzle

Zhang Zhi-hao, Xu Xiao-yong, Tian Yuan, Ding Zhao-bo, Chen Xu-yang
(Beijing Aerospace Propulsion Ιnstitute, Beijing, 100076)

Supersonic gas film cooling is used at nozzle extension of LOX/LH2thrust chamber, the turbine exhaust gas is imported into the extension section where the area ratio is 35 to form cooling gas film. Numerical simulation is used to study the effect of some parameters (the height of lip, blowing ratio, coolant flow rate, static pressure ratio) to specific impulse and the efficiency of gas film cooling. The research shows that the increasing of the height of lip can result in reducing of specific impulse and thrust, while the efficiency of cooling and the temperature of wall are almost no changing. Specific impulse and thrust decrease less than 0.06% when the height of lip increases from 3mm to 11mm. Raising the blowing ratio can improve the efficiency of cooling and the specific impulse of thrust slightly, when the pressures of main-flow and jet-flow are matching and the jet-flow mass is constant. Ιncreasing the jet-flow mass can improve the effect of gas film cooling when the pressures are matching and the blowing ratio is constant. The efficiency of gas film cooling can achieve the best results and the engine can get the highest specific impulse when the pressures of main-flow and jet-flow are matching and the jet-flow mass is constant.

LOX/LH2thrust chamber nozzle; Supersonic gas film cooling; Numerical simulation; Pressure matching

TP65

A

1004-7182(2017)04-0048-05 DOΙ：10.7654/j.issn.1004-7182.20170412

2016-03-30；

2017-06-11

張志浩（1990-），男，助理工程師，主要研究方向為液體火箭發(fā)動機系統(tǒng)設(shè)計、噴管冷卻技術(shù)