常心悅，王 豪，曹 凡，李 康，閔 浩，唐智禮

(1.中國(guó)航發(fā)四川燃?xì)鉁u輪研究院 民用航空動(dòng)力高空模擬四川省重點(diǎn)實(shí)驗(yàn)室，四川 綿陽(yáng) 621000；2.南京航空航天大學(xué) 航空學(xué)院 非定?？諝鈩?dòng)力學(xué)與流動(dòng)控制工信部重點(diǎn)實(shí)驗(yàn)室，南京 210016)

1 引言

航空發(fā)動(dòng)機(jī)高空模擬試驗(yàn)可以滿足航空發(fā)動(dòng)機(jī)在寬廣的飛行范圍內(nèi)進(jìn)行性能和穩(wěn)定性調(diào)試以及鑒定驗(yàn)收的需要，是航空發(fā)動(dòng)機(jī)裝機(jī)飛行前必不可少的重要環(huán)節(jié)[1]。對(duì)于航空發(fā)動(dòng)機(jī)的氣動(dòng)熱力學(xué)問題、機(jī)械系統(tǒng)問題、匹配性問題以及控制規(guī)律問題，都必須通過高空模擬試驗(yàn)來(lái)進(jìn)行充分的調(diào)試與驗(yàn)證[2]。

高空模擬試驗(yàn)艙(簡(jiǎn)稱高空艙)作為航空發(fā)動(dòng)機(jī)高空模擬試驗(yàn)的重要設(shè)施，其建設(shè)投資巨大。在高空艙建設(shè)前，利用CFD 技術(shù)對(duì)高空艙排氣流場(chǎng)進(jìn)行分析，在預(yù)研階段對(duì)相關(guān)項(xiàng)目進(jìn)行可行性分析與驗(yàn)證，能夠避免高空艙建設(shè)中的風(fēng)險(xiǎn)[3-5]。為此，國(guó)內(nèi)眾多學(xué)者進(jìn)行了大量研究。張建東等[6]利用CFD 手段，分析了內(nèi)外涵道噴管流道形狀對(duì)流場(chǎng)的變化規(guī)律，提出了一種合理的發(fā)動(dòng)機(jī)噴管形狀設(shè)計(jì)方案。李柯等[7]利用PID 控制器，對(duì)高空艙內(nèi)的溫度和壓力進(jìn)行了動(dòng)態(tài)精確控制。但志宏等[8-9]研究了排氣擴(kuò)壓器對(duì)高空艙艙壓的影響，并設(shè)計(jì)了一種基于擴(kuò)張狀態(tài)觀測(cè)器的高空臺(tái)進(jìn)氣環(huán)境模擬主動(dòng)抗擾控制技術(shù)。蘇金友等[10]實(shí)現(xiàn)了高空艙試驗(yàn)中低速流進(jìn)氣條件下的精確測(cè)量。劉志友、侯敏杰等[11-13]分析了高空模擬試驗(yàn)中發(fā)動(dòng)機(jī)推力的確定方法和引起誤差的影響因素。

目前，我國(guó)對(duì)小涵道比發(fā)動(dòng)機(jī)高空艙的試驗(yàn)和分析已經(jīng)具備了一定的經(jīng)驗(yàn)，但針對(duì)大涵道比分開排氣發(fā)動(dòng)機(jī)高空艙的試驗(yàn)測(cè)量和流動(dòng)機(jī)理還缺乏詳細(xì)的分析。為了進(jìn)一步探究大涵道比分開排氣發(fā)動(dòng)機(jī)高空艙內(nèi)復(fù)雜結(jié)構(gòu)對(duì)排氣流場(chǎng)的影響，確定適合的試驗(yàn)及仿真方法，本文針對(duì)某型大涵道比分開排氣發(fā)動(dòng)機(jī)，建立了高空試驗(yàn)艙、排氣擴(kuò)壓器、艙門和發(fā)動(dòng)機(jī)噴管的精細(xì)化聯(lián)合仿真模型。對(duì)不同落壓比下的高空艙排氣流場(chǎng)特性和發(fā)動(dòng)機(jī)推力進(jìn)行了仿真計(jì)算和試驗(yàn)測(cè)量，驗(yàn)證了利用CFD 方法進(jìn)行高空艙排氣流場(chǎng)計(jì)算的可靠性，可為大涵道比分開排氣發(fā)動(dòng)機(jī)高空艙內(nèi)的試驗(yàn)測(cè)量和大型高空艙氣動(dòng)布局設(shè)計(jì)提供參考。

2 聯(lián)合仿真模型與數(shù)值計(jì)算方法

為準(zhǔn)確進(jìn)行大涵道比分開排氣發(fā)動(dòng)機(jī)在高空艙內(nèi)的排氣流場(chǎng)的數(shù)值模擬，需要建立合適的聯(lián)合仿真模型，并根據(jù)發(fā)動(dòng)機(jī)與高空艙的物理特性確定合適的計(jì)算邊界條件。

2.1 聯(lián)合仿真模型

大涵道比分開排氣發(fā)動(dòng)機(jī)與小涵道比混合排氣發(fā)動(dòng)機(jī)不同，大涵道比分開排氣發(fā)動(dòng)機(jī)結(jié)構(gòu)尺寸較大，也更加復(fù)雜。針對(duì)其結(jié)構(gòu)特點(diǎn)，對(duì)比國(guó)內(nèi)外排氣系統(tǒng)簡(jiǎn)化相關(guān)研究，參考國(guó)外發(fā)動(dòng)機(jī)分開排氣系統(tǒng)模型進(jìn)行簡(jiǎn)化建模。簡(jiǎn)化后的模型具有內(nèi)、外涵道分開排氣的特點(diǎn)[12-13]，如圖1 所示。

圖1 大涵道比分開排氣發(fā)動(dòng)機(jī)簡(jiǎn)化模型Fig.1 Simplified model of high bypass ratio separated exhaust turbofan engine

排氣擴(kuò)壓器和高空艙外形相對(duì)規(guī)則，但由于艙門的存在使得高空艙呈非軸對(duì)稱結(jié)構(gòu)，只能采用全模型計(jì)算。最終建立了帶有艙門、艙底平臺(tái)、發(fā)動(dòng)機(jī)、高空艙和排氣擴(kuò)壓器的聯(lián)合仿真模型，如圖2 所示。模型比例為1：1，未進(jìn)行縮比，保證了雷諾數(shù)相同，防止了縮比模型尺度效應(yīng)帶來(lái)的模擬誤差。

圖2 高空艙整體三維模型Fig.2 Global 3D model of altitude simulation chamber

由于高空艙內(nèi)部結(jié)構(gòu)復(fù)雜，存在龍門柱、次流入口和艙門等復(fù)雜幾何外形，難以生成全結(jié)構(gòu)網(wǎng)格。為此，在模型規(guī)則區(qū)域劃分結(jié)構(gòu)網(wǎng)格，在模型復(fù)雜部分和拓?fù)浣Y(jié)構(gòu)不規(guī)則部分生成非結(jié)構(gòu)網(wǎng)格，整個(gè)流場(chǎng)為結(jié)構(gòu)加非結(jié)構(gòu)的混合網(wǎng)格。在發(fā)動(dòng)機(jī)內(nèi)外涵道出口附近區(qū)域，由于流動(dòng)情況復(fù)雜，流場(chǎng)中存在膨脹波且壓力梯度大，必須對(duì)其網(wǎng)格進(jìn)行加密，以便更好地模擬流場(chǎng)特性。在遠(yuǎn)離發(fā)動(dòng)機(jī)噴口區(qū)域，由于流動(dòng)速度小，流動(dòng)參數(shù)變化梯度小，可以生成較為稀疏的網(wǎng)格，這樣可以在減少網(wǎng)格數(shù)量的同時(shí)保證計(jì)算精度。此外，對(duì)于方柱區(qū)域、圓柱區(qū)域、艙門區(qū)域以及次流入口等區(qū)域，由于幾何結(jié)構(gòu)復(fù)雜，很難通過構(gòu)造規(guī)則的拓?fù)潢P(guān)系生成結(jié)構(gòu)網(wǎng)格。為此，在這些區(qū)域內(nèi)，對(duì)其生成致密的非結(jié)構(gòu)網(wǎng)格。綜合以上方法生成的混合網(wǎng)格如圖3 所示，網(wǎng)格數(shù)量約為1 400 萬(wàn)，其中結(jié)構(gòu)網(wǎng)格數(shù)量約占70%。

圖3 混合網(wǎng)格劃分Fig.3 Generation of mixed mesh

2.2 邊界條件

合理的邊界條件設(shè)置既能加快CFD 計(jì)算迭代的收斂，又能保障CFD 計(jì)算模擬的真實(shí)性。本次研究需要根據(jù)相關(guān)試驗(yàn)數(shù)據(jù)對(duì)大涵道比分開排氣發(fā)動(dòng)機(jī)的入口邊界條件進(jìn)行推導(dǎo)。圖4 為大涵道比分開排氣發(fā)動(dòng)機(jī)邊界條件示意圖。

圖4 大涵道比分開排氣發(fā)動(dòng)機(jī)邊界條件Fig.4 Boundary conditions of high bypass ratio separated exhaust turbofan engine

由于試驗(yàn)中并沒有測(cè)得發(fā)動(dòng)機(jī)入口和出口邊界條件的全部數(shù)值，所以現(xiàn)有試驗(yàn)數(shù)據(jù)并不能滿足CFD 計(jì)算所需要的邊界條件數(shù)據(jù)量。為此，本文按照等熵流總參數(shù)和靜參數(shù)關(guān)系式進(jìn)行推算，公式為：

式中：p0為總壓，p為靜壓，ρ0為總密度，ρ為靜密度，T0為總溫，T為靜溫，Ma為氣流馬赫數(shù)，γ為氣流比熱比。

再根據(jù)質(zhì)量流量m=ρAV，就可以得到如下方程：

式中：m為質(zhì)量流量，A為邊界入口面積，V為速度。

解方程(4)就可以得到馬赫數(shù)，再根據(jù)試驗(yàn)測(cè)量的總參數(shù)就可以得到靜參數(shù)，進(jìn)而就有了發(fā)動(dòng)機(jī)壓力入口的相關(guān)邊界條件參數(shù)。

由于大涵道比分開排氣發(fā)動(dòng)機(jī)飛行包線小，內(nèi)外涵道氣流多處于亞臨界狀態(tài)，需要更精確地模擬高空艙內(nèi)環(huán)境壓力——其模擬精度對(duì)發(fā)動(dòng)機(jī)推力評(píng)估有較大影響。因此，需動(dòng)態(tài)調(diào)整出口處的壓強(qiáng)，使得高空艙內(nèi)平均艙壓與模擬高度大氣壓強(qiáng)一致，此時(shí)計(jì)算收斂。

3 仿真與試驗(yàn)對(duì)比驗(yàn)證

利用高空艙試驗(yàn)測(cè)量的方法，對(duì)所研究的大涵道比分開排氣發(fā)動(dòng)機(jī)進(jìn)行了6 種工況下的高空模擬試驗(yàn)。利用試驗(yàn)測(cè)量得到各測(cè)點(diǎn)位置的壓強(qiáng)與發(fā)動(dòng)機(jī)推力的試驗(yàn)值，將試驗(yàn)值與通過數(shù)值模擬方法得到的計(jì)算值進(jìn)行對(duì)比，驗(yàn)證聯(lián)合仿真模型和計(jì)算方法的準(zhǔn)確性。

3.1 測(cè)試布局

大涵道比分開排氣發(fā)動(dòng)機(jī)高空模擬試驗(yàn)在某型高空艙內(nèi)進(jìn)行。試驗(yàn)設(shè)備由供氣機(jī)組、抽氣機(jī)組、高空艙、加降溫系統(tǒng)、排氣冷卻系統(tǒng)、測(cè)控系統(tǒng)等組成[14]。發(fā)動(dòng)機(jī)安裝在高空艙內(nèi)推力測(cè)量臺(tái)架上，采用了多種先進(jìn)的試驗(yàn)測(cè)量方法，對(duì)發(fā)動(dòng)機(jī)進(jìn)氣流量和推力進(jìn)行測(cè)量。試驗(yàn)原理如圖5 所示。

圖5 高空模擬試驗(yàn)原理圖Fig.5 Schematic diagram of high altitude simulation test

試驗(yàn)選取高度為7.35 km 的工況，高空模擬試驗(yàn)采用總壓、靜壓復(fù)合測(cè)量耙進(jìn)行測(cè)量。在高空艙內(nèi)外涵道噴管出口處放置1 支12 點(diǎn)總壓、靜壓復(fù)合測(cè)量耙。測(cè)點(diǎn)數(shù)量為12 個(gè)，其中第1～第5 測(cè)點(diǎn)間距為50 mm，第5 個(gè)測(cè)點(diǎn)距高空艙內(nèi)龍門柱前沿距離為200 mm，距離外涵道噴管出口距離為433 mm，距離軸線距離為900 mm；第5 點(diǎn)之后各測(cè)點(diǎn)之間的間距為133 mm，第6 個(gè)測(cè)點(diǎn)距離高空艙軸線距離為1 033 mm。測(cè)量耙上測(cè)量裝置為皮托管，皮托管長(zhǎng)度為200 mm，固定皮托管的圓管伸出高空艙內(nèi)龍門架長(zhǎng)度為350 mm，圓管直徑為20 mm。測(cè)量裝置在流場(chǎng)中的位置如圖6 所示。

圖6 壓力測(cè)量耙及測(cè)點(diǎn)布局示意圖Fig.6 Layout of pressure measuring rake and measuring points

大涵道比分開排氣發(fā)動(dòng)機(jī)推力作為高空模擬試驗(yàn)最重要的測(cè)量參數(shù)，通過推力測(cè)量系統(tǒng)測(cè)量。推力測(cè)量系統(tǒng)由動(dòng)架、推力預(yù)載、推力現(xiàn)場(chǎng)校準(zhǔn)、推力測(cè)量和附加阻力測(cè)量等多個(gè)子系統(tǒng)組成，其簡(jiǎn)化結(jié)構(gòu)如圖7 所示。

圖7 推力測(cè)量系統(tǒng)簡(jiǎn)化結(jié)構(gòu)示意圖Fig.7 Simplified structure diagram of thrust measurement system

試驗(yàn)時(shí)，發(fā)動(dòng)機(jī)推力通過發(fā)動(dòng)機(jī)安裝架傳遞到動(dòng)架上，由測(cè)力天平臺(tái)架測(cè)得。由于存在大量的管線及相關(guān)測(cè)量裝置，發(fā)動(dòng)機(jī)測(cè)量推力與真實(shí)飛行狀態(tài)下推力會(huì)有一定的誤差。

試驗(yàn)測(cè)得了同一高空艙模擬高度下6 組不同發(fā)動(dòng)機(jī)工況下的測(cè)量數(shù)據(jù)，各工況下發(fā)動(dòng)機(jī)轉(zhuǎn)速不同導(dǎo)致內(nèi)外涵道噴管截面的總壓、總溫不同。6 組工況下發(fā)動(dòng)機(jī)的各項(xiàng)參數(shù)如表1 所示，其中推力的計(jì)算結(jié)果依據(jù)線性歸一化方法統(tǒng)一進(jìn)行了無(wú)量綱的歸一化處理。線性歸一化能夠?qū)⒃紨?shù)據(jù)進(jìn)行線性變化，投影到[0，1]的區(qū)間內(nèi)，其公式如下：

表1 不同工況下的試驗(yàn)測(cè)量數(shù)據(jù)Table 1 Test measurement data under different work conditions

3.2 總壓與靜壓的對(duì)比驗(yàn)證

高空模擬試驗(yàn)所選取的12 個(gè)測(cè)點(diǎn)沿著發(fā)動(dòng)機(jī)外涵道出口依次排開。圖8 示出了工況1 的壓力校核曲線?？梢?，CFD 計(jì)算的總壓和靜壓值與高空模擬試驗(yàn)的測(cè)量值相近，最大誤差出現(xiàn)在測(cè)點(diǎn)1 到測(cè)點(diǎn)4 之間，超過了10%。這是因?yàn)檫@幾個(gè)測(cè)點(diǎn)處于發(fā)動(dòng)機(jī)外涵道噴管噴出的高速主流范圍之內(nèi)，總壓、靜壓復(fù)合測(cè)量耙的皮托管在高速氣流的作用力下容易發(fā)生結(jié)構(gòu)變形，導(dǎo)致測(cè)量值誤差較大；再加上靜壓耙測(cè)量裝置對(duì)流場(chǎng)的干擾，試驗(yàn)的測(cè)量值也會(huì)存在一定誤差。其余測(cè)點(diǎn)的測(cè)量值與計(jì)算值之間的相對(duì)誤差較小，在5%以內(nèi)，計(jì)算準(zhǔn)確度較高。

圖8 工況1 壓力校核曲線Fig.8 Pressure check curve under condition 1

3.3 發(fā)動(dòng)機(jī)推力的對(duì)比驗(yàn)證

推力的精確計(jì)算在發(fā)動(dòng)機(jī)性能預(yù)測(cè)和穩(wěn)態(tài)性能仿真中具有重要作用，為發(fā)動(dòng)機(jī)的設(shè)計(jì)提供了依據(jù)[15]。發(fā)動(dòng)機(jī)總推力等于動(dòng)量推力加上壓力推力，動(dòng)量推力通過質(zhì)量流量乘以速度得到，壓力推力通過壓強(qiáng)乘以面積得到，公式如下：

式中：Fg為發(fā)動(dòng)機(jī)總推力，W9為噴管出口質(zhì)量流量，c9為噴管出口射流速度，p9為出口處壓強(qiáng)，pe為環(huán)境壓強(qiáng)，A9為噴管出口面積。

大涵道比分開排氣發(fā)動(dòng)機(jī)內(nèi)外涵道分開排氣，相當(dāng)于有2 個(gè)噴管在排氣。在發(fā)動(dòng)機(jī)高空模擬試驗(yàn)時(shí)，不考慮發(fā)動(dòng)機(jī)短艙結(jié)構(gòu)，因此取推力計(jì)算控制體如圖9 所示。根據(jù)模型的對(duì)稱性，可以先計(jì)算半?？刂企w的推力，然后由半模的推力乘以2 得到整個(gè)發(fā)動(dòng)機(jī)的推力。圖中，1 為外涵道外壁面，2 為內(nèi)外涵道間隔壁面，3 為發(fā)動(dòng)機(jī)內(nèi)涵道外壁面，A91為外涵道出口面積，A92為內(nèi)外涵道間隔出口面積，A93為內(nèi)涵道出口面積。

圖9 發(fā)動(dòng)機(jī)推力計(jì)算控制體Fig.9 Engine thrust calculation control body

將計(jì)算域離散，在每個(gè)網(wǎng)格上對(duì)動(dòng)量推力和壓力推力分別進(jìn)行計(jì)算。計(jì)算公式分別為：

式中：W9i、W9j、W9k分別為外涵道出口，內(nèi)外涵道間隔出口和內(nèi)涵道出口每個(gè)網(wǎng)格處的質(zhì)量流量；c9i、c9j、c9k分別為外涵道出口，內(nèi)外涵道間隔出口和內(nèi)涵道出口每個(gè)網(wǎng)格處的速度；n1、n2、n3分別為外涵道出口，內(nèi)外涵道間隔出口和內(nèi)涵道出口面網(wǎng)格數(shù)量。

式中：右邊第1 項(xiàng)為控制體2 壁面的壓強(qiáng)推力，其中m2為控制體外壁面上的網(wǎng)格數(shù)量，pi為每個(gè)網(wǎng)格處的壓強(qiáng)，psch為高空艙艙壓，為每個(gè)網(wǎng)格的單位法向量，為x方向的單位向量，與飛行器實(shí)際飛行方向相反。類似地，第2 項(xiàng)為控制體3 壁面的壓力推力，第3 項(xiàng)為發(fā)動(dòng)機(jī)內(nèi)涵尾椎的壓力推力，第4 項(xiàng)為外涵道出口截面的壓力推力，第5 項(xiàng)為內(nèi)外涵道間隔出口的壓力推力，第6 項(xiàng)為發(fā)動(dòng)機(jī)內(nèi)涵出口的壓力推力。

圖10示出了不同落壓比下發(fā)動(dòng)機(jī)推力試驗(yàn)結(jié)果與CFD計(jì)算結(jié)果對(duì)比。圖中，NPR為落壓比?？煽闯觯S著落壓比增大，發(fā)動(dòng)機(jī)推力也增大，且呈線性變化；CFD 計(jì)算推力與試驗(yàn)測(cè)量推力之間有一定的誤差，但誤差在5%以下，在可接受范圍內(nèi)。誤差產(chǎn)生的原因是，試驗(yàn)中高空艙內(nèi)有管線，支撐發(fā)動(dòng)機(jī)臺(tái)架的前、后龍門柱產(chǎn)生的阻力也會(huì)傳遞到發(fā)動(dòng)機(jī)測(cè)力天平上，這些因素都會(huì)影響發(fā)動(dòng)機(jī)推力測(cè)量。

4 大涵道比分開排氣發(fā)動(dòng)機(jī)排氣流場(chǎng)特性分析

依照聯(lián)合仿真模型，對(duì)大涵道比分開排氣發(fā)動(dòng)機(jī)排氣流場(chǎng)特性進(jìn)行數(shù)值模擬分析。大涵道比分開排氣發(fā)動(dòng)機(jī)內(nèi)外涵道噴出的氣流形成主流射流，其膨脹過程會(huì)造成氣流流道的變化，處于接近完全膨脹狀態(tài)。外涵射流進(jìn)入高空艙后，受到次流的影響，在內(nèi)涵噴口截面之前符合等溫自由淹沒射流模型，當(dāng)其與內(nèi)涵氣流接觸混合后符合圓形噴口伴隨射流模型[16-17]。次流流量明顯小于主流流量，受發(fā)動(dòng)機(jī)主流引射作用影響，次流速度逐漸增大。此外，由于艙門外形和龍門柱的影響，艙門與龍門柱附近區(qū)域形成較多的旋渦結(jié)構(gòu)，不利于氣體的排出。在整個(gè)排氣流場(chǎng)中，發(fā)動(dòng)機(jī)高速射流與次流都由排氣擴(kuò)壓器排出，形成了一種內(nèi)壓式引射器，如圖11 所示。

圖11 高空艙內(nèi)速度流線Fig.11 Velocity streamlines in altitude simulation chamber

排氣擴(kuò)壓器效率是衡量排氣擴(kuò)壓器性能的一項(xiàng)重要參數(shù)，是在給定幾何條件、引射條件和噴管膨脹比條件下，排氣擴(kuò)壓器出口混合氣流總壓與噴管進(jìn)口氣流總壓之比。相關(guān)研究表明，排氣擴(kuò)壓器的效率對(duì)排氣流場(chǎng)的影響很大，同時(shí)還會(huì)影響高空艙環(huán)境壓力模擬的準(zhǔn)確性[18-19]。不同工況下排氣擴(kuò)壓器的效率如圖12 所示?？梢园l(fā)現(xiàn)，大涵道比分開排氣擴(kuò)壓器在圖示工況范圍內(nèi)，隨著落壓比的增大，排氣擴(kuò)壓器效率逐漸降低。小落壓比情況下氣流更容易被引射進(jìn)排氣擴(kuò)壓器。圖13 示出了高空艙中心軸線馬赫數(shù)分布。由圖可知，落壓比越大，排氣擴(kuò)壓器中氣流速度越大，整體呈先上升后下降的趨勢(shì)，在排氣擴(kuò)壓器出口處趨于平緩。氣流在排氣擴(kuò)壓器中總的能量損失相同，排氣擴(kuò)壓器出口處速度越小，排氣擴(kuò)壓器效率越高。因此，需合理設(shè)計(jì)排氣擴(kuò)壓器，防止排氣擴(kuò)壓器效率降低后無(wú)法及時(shí)排出高溫氣流。

圖12 排氣擴(kuò)壓器效率隨落壓比的變化Fig.12 Curve of exhaust diffuser efficiency changing with nozzle pressure ratio

圖13 高空艙中心軸線馬赫數(shù)分布Fig.13 Mach number distribution along the central axis of altitude simulation chamber

不同落壓比下高空艙內(nèi)發(fā)動(dòng)機(jī)噴口之后的馬赫數(shù)云圖如圖14 所示?？梢?，馬赫數(shù)隨著落壓比的增大而增大，在排氣擴(kuò)壓器入口處產(chǎn)生了一個(gè)低速區(qū)域，發(fā)動(dòng)機(jī)主流射流進(jìn)入排氣擴(kuò)壓器之后流動(dòng)速度不斷衰減。當(dāng)落壓比為1.67 時(shí)，發(fā)動(dòng)機(jī)尾噴口的馬赫數(shù)還未達(dá)到1，處于亞聲速狀態(tài)。當(dāng)落壓比為2.14 時(shí)，外涵道出口馬赫數(shù)開始逐漸等于或大于1，已經(jīng)達(dá)到臨界或超臨界狀態(tài)。氣流噴出發(fā)動(dòng)機(jī)外涵道尾噴管出口以后繼續(xù)膨脹加速，達(dá)到超聲速狀態(tài)，這與發(fā)動(dòng)機(jī)噴管理論相符。

圖14 不同落壓比下馬赫數(shù)云圖Fig.14 Mach number contour under NPR=1.67 and NPR=2.14

圖15 示出了不同噴管落壓比下發(fā)動(dòng)機(jī)主流的影響范圍。由圖可知，隨著噴管落壓比的增大，主流流動(dòng)范圍不斷外擴(kuò)。在與發(fā)動(dòng)機(jī)軸線距離小于0.85 m 時(shí)，氣流都處于發(fā)動(dòng)機(jī)外涵道的主流之中，速度很大；此后，隨著與發(fā)動(dòng)機(jī)軸線距離的增大，速度不斷下降，在距離大于1.00 m 以后速度趨近于平緩，這時(shí)氣流已經(jīng)完全不在發(fā)動(dòng)機(jī)主流區(qū)域范圍內(nèi)，氣流的速度受次流速度和射流引射作用影響。

圖15 不同噴管落壓比下發(fā)動(dòng)機(jī)主流的影響范圍Fig.15 Influence range of main flow of engine under different nozzle pressure ratio

落壓比為1.67 時(shí)，主流范圍為873.2 mm，此時(shí)主流影響范圍最小。由于高空艙龍門柱加強(qiáng)筋與發(fā)動(dòng)機(jī)中心軸線的距離為841.0 mm，所以龍門柱位于發(fā)動(dòng)機(jī)主流范圍內(nèi)，主流會(huì)直接流向發(fā)動(dòng)機(jī)外涵道出口處的龍門柱與發(fā)動(dòng)機(jī)臺(tái)架裝置，龍門柱產(chǎn)生的反作用力傳遞到發(fā)動(dòng)機(jī)測(cè)力天平上，影響發(fā)動(dòng)機(jī)推力的測(cè)量。在后續(xù)的大涵道比發(fā)動(dòng)機(jī)高空模擬試驗(yàn)中，可以通過改進(jìn)龍門柱安裝位置，減小其對(duì)主流的影響，從而降低測(cè)量誤差。

圖16 示出了不同狀態(tài)下發(fā)動(dòng)機(jī)水平面的壓力等值線?？梢?，隨著落壓比的增大，發(fā)動(dòng)機(jī)內(nèi)外涵道進(jìn)口總壓也增大。當(dāng)落壓比小于2.14 時(shí)，內(nèi)外涵道出口的靜壓與高空艙內(nèi)靜壓基本相等，氣體處于亞臨界狀態(tài)，發(fā)動(dòng)機(jī)內(nèi)外涵道內(nèi)的氣流完全膨脹。當(dāng)落壓比大于2.14 后，外涵道噴管壓力逐漸大于高空艙壓力，表明氣體在噴管內(nèi)并未完全膨脹，氣體流出外涵道噴管以后繼續(xù)膨脹加速，氣流動(dòng)壓不斷增大，靜壓不斷減小，主流與高空艙內(nèi)次流互相耦合，最終達(dá)到壓強(qiáng)相等。而所有計(jì)算狀態(tài)下，發(fā)動(dòng)機(jī)內(nèi)涵道尾噴口始終處于亞臨界狀態(tài)，氣體在涵道中完全膨脹，但由于外涵道流出的高速氣流與內(nèi)外涵道間隔處氣流干擾影響，使得內(nèi)涵道出口處壓強(qiáng)發(fā)生變化。

圖16 不同落壓比下發(fā)動(dòng)機(jī)水平面壓力等值線Fig.16 Pressure contour in engine horizontal surface under different nozzle pressure ratio

5 結(jié)論

以大涵道比分開排氣發(fā)動(dòng)機(jī)為研究對(duì)象，建立了高空艙、發(fā)動(dòng)機(jī)、艙門和排氣擴(kuò)壓器的精細(xì)化聯(lián)合仿真模型。分析了大涵道比分開排氣發(fā)動(dòng)機(jī)不同噴管落壓比下的高空艙內(nèi)排氣流場(chǎng)特性，計(jì)算了發(fā)動(dòng)機(jī)推力，并與試驗(yàn)數(shù)據(jù)進(jìn)行了對(duì)比驗(yàn)證。主要結(jié)論為：

(1) 高空艙內(nèi)總壓和靜壓計(jì)算結(jié)果與試驗(yàn)結(jié)果在測(cè)點(diǎn)1 到測(cè)點(diǎn)4 之間誤差較大(最大誤差超過10%)，在測(cè)點(diǎn)5 到測(cè)點(diǎn)12 之間誤差小于5%，說明發(fā)動(dòng)機(jī)噴管高速射流影響測(cè)量精度，造成了一定的誤差；發(fā)動(dòng)機(jī)推力整體上隨著落壓比增大而線性增大，整體誤差保持在5%以內(nèi)，數(shù)值模擬計(jì)算精度較高。

(2) 隨著落壓比的增大，外涵道噴管射流馬赫數(shù)不斷增大，并逐漸達(dá)到超臨界狀態(tài)，主射流區(qū)域的影響范圍也不斷增大，排氣擴(kuò)壓器效率逐漸降低。為防止氣流阻塞，需合理設(shè)計(jì)排氣擴(kuò)壓器的直徑與長(zhǎng)度，保證氣體高效排出。

(3) 測(cè)力臺(tái)架的龍門柱對(duì)外涵道射流產(chǎn)生了一定的阻礙作用，影響了測(cè)量精度，且在龍門柱與艙門區(qū)域容易形成回流。應(yīng)對(duì)高空艙的排氣氣動(dòng)布局進(jìn)行適當(dāng)改進(jìn)以減少回流，提高試驗(yàn)測(cè)量精度。