張 振，張學(xué)偉

(中國航天空氣動力技術(shù)研究院，北京 100074)

0 引 言

深空探測對于發(fā)展科學(xué)技術(shù)、提高科學(xué)認識、探索宇宙和生命的起源及演化等具有重要意義，所以現(xiàn)在世界航天大國都對深空探測技術(shù)發(fā)展和研究給予了高度關(guān)注和支持?；鹦且驗楹偷厍蛴斜姸嘞嗨浦幎蔀楦鲊羁仗綔y的主要目標并進行了相關(guān)的研究工作[1-3]，但是火星表面的大氣條件和地球大氣條件差別巨大，其對比如表1所示，所以地球上的傳統(tǒng)研究方法不再適用于火星探測研究。從表1可以看出,相比于地球大氣環(huán)境來說，火星大氣表面的大氣壓力、大氣密度、溫度都明顯低于地球表面大氣，屬于低壓低密度的大氣環(huán)境。在對火星進行探索之前，需要對火星探測裝置中應(yīng)用的新技術(shù)和新方法進行性能測試和評估，比如利用火星風洞提供的大氣壓力、速度和溫度，驗證在火星探測裝置中運用新的熱控技術(shù)和新型機翼的升力性能等。在火星探測裝置中火星飛行器因為其覆蓋面廣，偵查精度高等優(yōu)點而受到越來越多研究學(xué)者的重視，但是火星飛行器是在低壓低密度情況下進行工作，所以飛行器的飛行速度必須達到足夠高的速度才能維持足夠的升力(Ma=0.75，Re=104～105)。在此雷諾數(shù)范圍內(nèi)飛行器機翼表面易發(fā)生層流分離現(xiàn)象，并且層流分離渦的形態(tài)結(jié)構(gòu)會隨著雷諾數(shù)和迎角變化出現(xiàn)顯著變化，比如隨著迎角和雷諾數(shù)變化，分離渦先后出現(xiàn)后緣層流分離、分離泡前移、分離再附和分離泡破裂等一系列復(fù)雜現(xiàn)象。同時由于壓縮性效應(yīng)開始顯示[4-5]，流動狀態(tài)的轉(zhuǎn)變和壓縮性的雙重影響會給飛行器的機翼工作狀況帶來嚴重的影響。而現(xiàn)在對于這方面的試驗研究較少，所以很有必要對此臨界狀態(tài)下的飛行特性進行研究進而掌握其基本特點，這就需要為火星探測裝置提供火星表面的大氣環(huán)境，并且能夠提供飛行器需要的高亞聲速條件。而風洞是進行氣動力試驗的重要設(shè)備，可以為各種新技術(shù)和新設(shè)備提供準確可靠的試驗數(shù)據(jù)，是進行飛行器設(shè)計的基礎(chǔ)設(shè)備，所以對低壓低密度風洞設(shè)計及研究就顯得尤其重要。

文獻[7]對美國的MARSWIT火星風洞進行了介紹，該風洞屬于直流式閉口風洞，其中動力部分采用引射器裝置而不是風扇等裝置，這是因為在低壓低密度情況下，風扇等裝置效率急劇降低。工質(zhì)可采用空氣和二氧化碳分別用來模擬地球高空和火星表面的大氣環(huán)境，在500 Pa的低壓環(huán)境下，該風洞試驗段馬赫數(shù)可達到0.8。文獻[8]對能夠模擬火星大氣表面的低密度風洞(Mars wind tunnel,MWT)進行了設(shè)計研究，MWT風洞也是屬于直流式風洞，動力裝置也是采用引射驅(qū)動方式，在來流總壓為1 kPa下，工質(zhì)分別為空氣和二氧化碳的情況下，試驗段馬赫數(shù)可以分別達到0.67和0.74。文獻[9]對火星風洞AWTS進行了研究，AWTS風洞置于真空管中，采用液氮冷卻，采用風扇抽吸的方式可使試驗段的風速達到20 m/s，該風洞主要模擬火星大氣壓力、溫度、氣體組分，可見UV光條件和塵埃環(huán)境。在國內(nèi)方面，關(guān)于火星風洞的研究幾乎處于空白狀態(tài)，根據(jù)已有的資料，國內(nèi)只有天津大學(xué)王彩虹[10]對火星風洞進行了初步的研究設(shè)計，但是也僅僅處于對火星風洞本體的數(shù)值初步計算階段，并沒有對風洞進行深入計算和研究分析。

相比于國內(nèi)，國外的火星風洞具有很大的技術(shù)優(yōu)勢，并且在進行火星風洞設(shè)計時，由于在低壓環(huán)境下風扇或者渦輪的工作效率不高，所以要達到高亞聲速的要求，一般采用引射器作為動力系統(tǒng)。相比于風扇或者渦輪，引射器還具有結(jié)構(gòu)簡單，制造容易，無需外力等優(yōu)點。所以對引射器的性能研究成為低壓低密度風洞設(shè)計研究的重點，很多研究學(xué)者也都對引射器進行了研究。

文獻[11]對超聲速引射器進行了一維計算，與國內(nèi)其他機構(gòu)的結(jié)果進行了對比，發(fā)現(xiàn)計算結(jié)果與試驗結(jié)果符合較好，并得出了引射器各參數(shù)對引射器性能的影響。文獻[12]對多噴管引射器進行了試驗研究和數(shù)值模擬，研究表明使用FLUENT的標準k-e湍流模型和壓力-速度耦合的SIMPLE算法，可以很好地對五噴管引射器模型進行精確結(jié)算，試驗曲線結(jié)果和計算曲線結(jié)果的平均偏差在4.42%。文獻[13]對等截面引射器啟動性能進行了數(shù)值研究，研究了等截面引射器長度、引射器與噴管的間隙對啟動性能的影響，并得到了其影響規(guī)律，發(fā)現(xiàn)采用k-w 過渡剪切應(yīng)力傳輸模型(Transition shear-stress transport model,SST)可以很好地模擬引射器的工作情況。由上述研究可知，數(shù)值計算準確性得到越來越多研究學(xué)者的認可，其中k-w SST湍流模型的適用性得到了普遍的驗證。

本次研究基于上述低壓低密度風洞的研究現(xiàn)狀及相應(yīng)特點，對火星風洞進行了設(shè)計并對其性能進行了初步數(shù)值研究。本次設(shè)計的風洞基本形式為直流式風洞，其基本型面如圖1所示，風洞本體分為六個部分，分別是穩(wěn)定段、收縮段、試驗段、第一擴散段、動力段和第二擴散段，其中穩(wěn)定段上接唇口段。風洞本體放置在一個低壓環(huán)境艙中，環(huán)境艙直徑1800 mm，總長度5400 mm。通過控制環(huán)境艙底部的抽氣口的抽氣量來實時保證環(huán)境艙內(nèi)的總壓維持在1000 Pa左右，其中低壓環(huán)境艙內(nèi)的壓力檢測由裝在風洞穩(wěn)定段的壓力傳感器得到，傳輸?shù)娇刂葡到y(tǒng)進而對艙內(nèi)壓力進行控制。

圖1 風洞結(jié)構(gòu)圖

1 風洞組成部分設(shè)計

風洞試驗段的設(shè)計主要考慮以后對弦長為12 mm機翼表面在低密度低雷諾數(shù)下的層流分離特性的研究，同時也對火星探測飛行器或者其他裝置提供試驗平臺，雖然國內(nèi)也有相關(guān)研究學(xué)者對機翼層流分離對機翼特性進行了相關(guān)研究[14-16]，但是在低壓低密度情況下的機翼特性未進行過研究，認識還不夠深入。本次設(shè)計試驗段的截面形狀設(shè)計為矩形，進口尺寸為100 mm×150 mm，長度為2.5倍試驗段進口截面當量直徑，即300 mm，同時需要考慮到在低密度情況下氣流邊界層對試驗段氣流品質(zhì)的影響，對試驗段的上、下壁面各采用1°的擴散角。

對于穩(wěn)定段的設(shè)計，收縮比一般取為7～10，并且大的收縮比有利于減小試驗段的湍流強度，但是過大的收縮比會增加收縮段的加工制造難度，綜合考慮選擇收縮比為9，則穩(wěn)定段的截面尺寸為(高×寬)300 mm×450 mm，考慮到穩(wěn)定段需要布置蜂窩器和阻尼網(wǎng)進行整流，并考慮到靜流段需要足夠的長度使經(jīng)過蜂窩器和阻尼網(wǎng)的氣流得到充分的均勻穩(wěn)定，取穩(wěn)定段長度720 mm。在本次研究過程中，由于蜂窩器具有一定的軸向厚度，采用多孔介質(zhì)模型進行模擬，而阻尼網(wǎng)軸向厚度很小，所以采用多孔階躍模型進行模擬[17]。為了進一步使進入到風洞穩(wěn)定段的氣流具有較高的均勻度，在穩(wěn)定段進口設(shè)置了半徑為50 mm的唇口。

對于收縮段的收縮曲線采用雙三次曲線，其相對于韋氏曲線和五次曲線能夠產(chǎn)生比較均勻的氣流品質(zhì)[18]，雙三次曲線的方程如式(1)所示，需要注意的是本次所設(shè)計風洞的上下壁面和左右壁面曲線不同，需分別計算求得，收縮段的長度通常取進口當量直徑的0.5～1倍，在本次設(shè)計中取收縮段長度為300 mm。對于第一擴散段和第二擴散段進出口面積比為2，當量擴散角選為5°。

雙三次曲線方程：

(1)

式中：xm是兩曲線前后連接點，在這里取xm=0.5,D2是收縮段出口截面半徑，D1是收縮段進口截面半徑，D是軸向距離為x處的截面直徑，原點在收縮段進口截面中心。

2 引射器設(shè)計及性能參數(shù)

引射器是風洞的動力裝置，其性能的好壞直接關(guān)系到風洞性能能否達到設(shè)計要求，所以必須對引射器的設(shè)計進行特別關(guān)注。在引射器設(shè)計過程中，必須對引射器的布置形式和關(guān)鍵參數(shù)進行謹慎選擇設(shè)計。

1)引射器型式選擇

2)引射器參數(shù)選擇

引射器的計算理論模型中，假定引射氣流與被引射氣流兩股氣流在噴嘴出口截面上是均勻的，并且在混合室出口，兩股氣流已經(jīng)充分混合，可以運用氣體一維理論進行分析。

由流量守恒定理，并假定高、低壓氣流總溫相等(在氣流總溫溫度比不大于2時對結(jié)果影響很小，可以忽略不計)，

P01q(λ1)A1+P′0q(λ′)A′=P″0q(λ″)A″

(2)

(3)

由動量守恒定理可得：

κz(λ1)+z(λ′)=(1+κ)z(λ″)

(4)

根據(jù)引射系數(shù)定義：

(5)

同時由于引射氣流和被引射氣流在摻混過程中會產(chǎn)生壓縮和擴散現(xiàn)象，在開始混合階段假設(shè)引射氣流和被引射氣流還未發(fā)生明顯摻混，可以視為兩股流體分別受到壓縮和擴散，由于引射氣流壓力較高會發(fā)生擴散，而被引射氣流則會受到壓縮，當被引射氣流達到馬赫數(shù)為1時，達到臨界狀態(tài)，引射系數(shù)不再增加，所以要對引射器是否達到臨界狀態(tài)進行檢驗，運用式(6)進行檢驗，式(6)可由質(zhì)量方程和動量方程得出。

(6)

在對引射器性能進行計算時，上述式(2)～(5)方程組可簡化為式(7)進行引射器設(shè)計。

(7)

在引射器設(shè)計時需要對設(shè)計參數(shù)中的部分參數(shù)根據(jù)經(jīng)驗進行選定。在本次設(shè)計中，由于來流壓力較低，選擇引射馬赫數(shù)為3，面積比選擇為0.05，以得到較高的引射系數(shù)。引射噴嘴放置在風洞擴壓器出口、混合室進口附近，引射器噴嘴與混合室軸線夾角選為0。并且由于引射器的氣源比較充足，在這里暫時不考慮氣源和運行時間問題。在增壓比方面，由于跨聲速風洞的增壓比一般在1.05～1.25之間，在這次設(shè)計中選擇增壓比為1.1。

在引射馬赫數(shù)、引射面積比、被引射氣流馬赫數(shù)、被引射氣流總壓、引射增壓比確定后，就可以利用式(7)得出引射噴嘴的膨脹比為12。

3 多孔介質(zhì)和多孔階躍模型及其阻力參數(shù)

多孔介質(zhì)模型可用于篩網(wǎng)，過濾器等多孔物質(zhì)，將多孔區(qū)域簡化為增加了阻力源項的流體區(qū)域，從而省去建立具體多孔模型的麻煩。而多孔階躍模型多應(yīng)用于薄孔板，而且此模型的計算穩(wěn)定性好，應(yīng)在不影響流場的情況下盡量采用。本次風洞設(shè)計的蜂窩器選擇為正六邊形蜂窩器，邊長為4 mm，長度取15倍邊長，即60 mm。阻尼網(wǎng)設(shè)計為5層，網(wǎng)絲直徑為0.1 mm，阻尼網(wǎng)之間的距離為60 mm，根據(jù)蜂窩器和阻尼網(wǎng)的結(jié)構(gòu)參數(shù)可以得出不同來流情況下各自的損失系數(shù)[17]，進而得出在不同來流速度條件下的壓力損失，可進一步得出相應(yīng)的多孔阻力參數(shù)。

多孔模型的一般動量方程源項[18]：

(8)

式中：Si是不同方向上的動量壓力源項，D和C是指定的矩陣，μ是黏性系數(shù)，右側(cè)第一項為黏性損失項，第二項為慣性損失項。

當多孔介質(zhì)為各向同性時，上述方程變?yōu)閇17]：

(9)

式中：α是滲透率，其倒數(shù)為黏性阻力系數(shù)，C2是慣性阻力系數(shù)。

由以上的參數(shù)可以得出蜂窩器的黏性阻力系數(shù)和慣性阻力系數(shù)分別為0和6.7，每層阻尼網(wǎng)的黏性阻力參數(shù)和慣性阻力參數(shù)分別為0和1384。

4 仿真校驗

本次計算中，運用商業(yè)軟件FLUENT進行計算，該軟件能夠?qū)α鲌銮闆r進行準確預(yù)測，在流體計算中得到廣泛的應(yīng)用。在流場求解方面，采用隱式ROE格式和二階迎風格式進行流場計算。在邊界條件方面，噴管的入口為壓力入口邊界條件，抽氣口為壓力出口邊界條件。因為SST湍流模型在近壁區(qū)和射流等復(fù)雜流動方面比k-w SST湍流模型更加精確，而且也能對層流流動進行比較準確的預(yù)測，所以湍流模型采用SST四方程模型。壁面采用無摩擦絕熱邊界條件。

為了計算方便，將每個噴嘴簡化為一個收擴噴管，不考慮噴管上游進氣管路對流場的影響，并對風洞進行分區(qū)劃分網(wǎng)格，以便于快速計算收斂，考慮到本次計算的對象是軸對稱模型，所以對風洞的四分之一模型進行計算和劃分網(wǎng)格，劃分的網(wǎng)格模型如圖2～圖3所示。

圖2 風洞整體網(wǎng)格

圖3 噴管網(wǎng)格

為保證計算的準確性，分別使用四套網(wǎng)格進行了網(wǎng)格無關(guān)性驗證，分別是A,B,C,D四套網(wǎng)格，四套網(wǎng)格的基本信息如表2所示。

表2 網(wǎng)格基本信息

在出口馬赫數(shù)為3，膨脹比為12時試驗段中心線上的馬赫數(shù)分布如圖4所示，從計算結(jié)果可以看出，網(wǎng)格A明顯與其他三套網(wǎng)格的計算結(jié)果存在差別，B，C和D網(wǎng)格的計算結(jié)果相對來說比較接近。網(wǎng)格B的計算結(jié)果相對于C和D網(wǎng)格始終偏小，存在一定的偏差，C和D網(wǎng)格的計算結(jié)果比較接近。試驗段截面(X=0.45 m，Y=0位置處)中心線(Z=-0.05 m～0.0 m)馬赫數(shù)變化如圖5所示，圖4展現(xiàn)的結(jié)果與圖5相似，網(wǎng)格A的計算結(jié)果與其他三個網(wǎng)格存在明顯不同，網(wǎng)格B,C,D的計算結(jié)果基本一致。綜合考慮到計算資源、計算時間和精確度，選擇網(wǎng)格C作為此次研究計算的網(wǎng)格。

圖4 試驗段中心軸線馬赫數(shù)變化

圖5 試驗段中間截面(X=0.45 m，Y=0)，沿著Z方向馬赫數(shù)變化

5 計算結(jié)果分析

1)模式一(膨脹比12)

(1)引射器性能分析

由于本次研究采用引射器作為風洞的動力系統(tǒng)，所以很有必要對引射器的性能進行分析，根據(jù)理論計算，要達到試驗段馬赫數(shù)為0.7，引射器的膨脹比應(yīng)取為12(來流總壓為1000 Pa)。引射器按照理論設(shè)計時的性能計算結(jié)果如表3所示，此次設(shè)計的引射器引射系數(shù)為3，增壓比為1.085，都達到了比較高的水平[18]。

表3 引射器性能參數(shù)(模式一)

(2)風洞性能分析

本次研究首先對介質(zhì)為空氣的風洞進行了研究，風洞ZX截面(即寬度截面)的壓力分布和馬赫數(shù)分布如圖6和7所示，XY截面(即高度截面)的馬赫數(shù)分布如圖8所示。由圖6可知，低壓艙內(nèi)的壓力穩(wěn)定在989 Pa，與先前設(shè)計值1000 Pa壓力比較符合，達到了設(shè)計目的。從圖7和圖8的馬赫數(shù)分布可以看出，低壓艙內(nèi)馬赫數(shù)較高的區(qū)域在風洞的出口，馬赫數(shù)達到了0.3，而在風洞進口的馬赫數(shù)非常低，只有0.04。并且從風洞進口的馬赫數(shù)和壓力等值線可以看出，風洞進口的速度和壓力都相當均勻，保證了風洞進口流體高質(zhì)量的要求。

從圖7和圖8可以看出，試驗段馬赫數(shù)沿著軸向基本保持不變。從圖7中的混合段的馬赫數(shù)分布可以看出，由于引射的作用，在引射器位置附近產(chǎn)生了高速區(qū)，對來流氣體產(chǎn)生了抽吸作用。這點也可以從圖6的壓力分布中看出，由于引射器的作用，在引射器出口附近區(qū)域產(chǎn)生了低速區(qū)，對來流產(chǎn)生了明顯的抽吸作用。同時也可以看到沿著軸線方向不同截面方向上的馬赫數(shù)分布規(guī)律并不相同：在XY截面上，邊界層增長較慢，厚度較薄，而在XZ截面上，邊界層較厚，并且等馬赫數(shù)線出現(xiàn)了突尖的現(xiàn)象。造成這個現(xiàn)象主要是由兩方面導(dǎo)致的，一方面由于引射器是在混合室Y方向和Z方向按照混合室截面不同方向的尺寸大小均勻分布的，在XY截面上噴管之間的距離是75 mm，而在XZ截面上噴管之間的距離是47 mm，噴管與壁面不同的距離對試驗段的速度分布產(chǎn)生了影響，噴管距離噴管較近時影響主要體現(xiàn)在邊界層方面，所以不會出現(xiàn)等馬赫數(shù)線突尖的現(xiàn)象，而距離較大時，噴管則會影響邊界層之外的區(qū)域，從而形成了突尖。另一方面由于試驗段上下壁面具有1°的向外傾角來平衡壁面邊界層增長的影響，所以導(dǎo)致在XZ方向會出現(xiàn)邊界層較厚而XY邊界層較薄的現(xiàn)象。從ZX截面和XY截面穩(wěn)定段的壓力分布可以看出，經(jīng)過穩(wěn)定段的流速和壓力很均勻，這說明采用多孔介質(zhì)和多孔跳躍模型能夠很好地模擬防湍網(wǎng)和蜂窩器的作用。

圖6 風洞ZX截面壓力分布

圖7 風洞ZX截面馬赫數(shù)分布

圖8 風洞XY截面馬赫數(shù)分布

試驗段軸線方向風洞的靜壓分布如圖9所示，從靜壓分布中可得沿著軸線的靜壓變化梯度為0.025 Pa/mm，這說明試驗段上下壁面的擴張角大小能夠正好抵消試驗段內(nèi)邊界層的增加，保證了試驗段內(nèi)壓力和速度分布的均勻性，使試驗段保持了良好的品質(zhì)，有效地提高試驗測量的精確度。同時也能看到在試驗段進口段沿著軸向壓力降低，這是由于采用雙三次曲線的收縮段會產(chǎn)生輕微的過沖現(xiàn)象[19]，與文獻[26]中的描述相符。

圖9 試驗段軸線靜壓變化

風洞的湍流強度大小直接影響到風洞內(nèi)流體的品質(zhì)，所以需對風洞內(nèi)流體的湍流強度變化情況進行研究，湍流強度定義為氣流三個方向的脈動氣流速度的方根平均值與主流平均速度之比。風洞不同軸向截面的湍流強度大小如圖10所示，其中曲線上標出的位置依次是穩(wěn)定段出口、蜂窩器進出口、五層湍流網(wǎng)和穩(wěn)定段出口截面的位置。從圖10可以看出，在穩(wěn)定段進口和蜂窩器前湍流強度不斷加強，而在經(jīng)過蜂窩器和湍流網(wǎng)后湍流強度不斷減小，說明蜂窩器和湍流網(wǎng)能夠很好地降低氣體湍流度，起到梳理流體的作用，也從另一方面反映了多孔模型能夠模擬蜂窩器和湍流網(wǎng)的降湍作用。從第五層湍流網(wǎng)到穩(wěn)定段出口反映了經(jīng)過蜂窩器和湍流網(wǎng)后湍流強度進一步衰減，到穩(wěn)定段出口衰減到最小。

圖10 穩(wěn)定段湍流強度變化

2)模式二(膨脹比15)

為了進一步探索此次設(shè)計引射器的性能特點及驗證能否運用引射器作為動力系統(tǒng)使試驗段馬赫數(shù)達到0.7，本次研究采用膨脹比為15的噴管作為動力進行研究，其他設(shè)計參數(shù)和膨脹比為12的引射器相同。

(1)引射器性能分析

模式二的計算性能參數(shù)如表4所示，模式一和模式二的引射器性能參數(shù)如表5所示，從表5可以看出，相比于模式一的膨脹比為12的引射系統(tǒng)，膨脹比為15的引射器的引射系數(shù)略微降低，降低幅度為4.76%，但是增壓比卻得到了比較大的提高，相比于模式一提高了11.33%。

表4 引射器性能參數(shù)(模式二)

表5 模式一和模式二性能比較表

(2)風洞流場性能分析

膨脹比為15時風洞ZX截面的壓力分布和馬赫數(shù)分布如圖11和圖12所示，從壓力分布可以看出低壓艙的壓力維持在988 Pa不變，和預(yù)期值1000 Pa很接近，說明本次設(shè)計能夠滿足風洞來流總壓為1000 Pa的要求。并且從風洞進口的等壓分布可以看出風洞進口的壓力分布很均勻，說明此次風洞進口唇口設(shè)計的合理性。

圖11 風洞ZX截面壓力分布

圖12 風洞ZX截面馬赫數(shù)分布

從圖11和圖13可以看出，試驗段的馬赫數(shù)均可達到0.77，并且沿著試驗段的軸線方向馬赫數(shù)幾乎不發(fā)生變化。從圖11可以看出風洞不同組成部分的壓力保持穩(wěn)定。試驗段中心軸線的壓力變化如圖14所示，壓力的梯度變化為0.037 Pa/mm，可以看出相對于膨脹比為12的試驗段壓力變化來說，壓力梯度變化增大，同樣也可以看到試驗段同樣存在輕微的過沖現(xiàn)象。模式一和模式二試驗段X=0.45 m和Y=0 m截面位置從Z=0 m到Z=0.075 m的馬赫數(shù)分布如圖15所示，從兩種模式的馬赫數(shù)分布來說，模式二的邊界層較薄，這主要是因為邊界層厚度與雷諾數(shù)呈負相關(guān)。并且，在XZ截面上試驗段中等馬赫數(shù)線未出現(xiàn)模式一中出現(xiàn)的突尖現(xiàn)象，這主要是因為引射總壓提高之后引射系數(shù)降低，對上游的影響減弱，但是在擴壓段里仍然可以看到有類似突尖的現(xiàn)象出現(xiàn)。進一步，可以看出模式二的試驗段馬赫數(shù)明顯增大，這說明增加膨脹比能夠使試驗段的馬赫數(shù)得到大的提高，是提高試驗段雷諾數(shù)的一個有效措施。

圖13 風洞XY截面馬赫數(shù)分布

圖14 試驗段中心軸線壓力變化

圖15 在X=0.45和Z=0時在Y軸方向上的馬赫數(shù)分布

6 結(jié) 論

本次研究設(shè)計了一個低密度高亞聲速引射風洞，并對所設(shè)計的風洞進行了數(shù)值研究。針對不同膨脹比下引射器的性能特點和風洞試驗段的變化進行了對比分析，得出以下結(jié)論：

1)運用引射器作為動力系統(tǒng)能夠使試驗段的速度達到比較高的馬赫數(shù)，并且通過設(shè)置合理的試驗段傾角，可使試驗段沿著軸線的靜壓基本不變。但是膨脹比增加，試驗段中心軸線的靜壓梯度也增加。當膨脹比從12變化到15時，試驗段中心軸線的靜壓梯度由0.025 Pa/mm增加到0.037 Pa/mm。

2)提高引射膨脹比能夠有效提高試驗段馬赫數(shù)，提高試驗雷諾數(shù)，增加被引射流體的流量，提高增壓比，但是會降低引射系數(shù)，降低引射效率。當引射器膨脹比從12增加到15時，試驗段馬赫數(shù)從0.64提高到0.76，增壓比從1.085提高到1.208，但是引射系數(shù)從3.0降低到2.86。

3)通過此次研究，說明采用引射器作為低壓低密度風洞動力系統(tǒng)能夠使風洞達到高亞聲速的要求，可以作為將來低壓低密度風洞設(shè)計和建設(shè)的動力系統(tǒng)，為將來低壓低密度風洞的建設(shè)提供有力支持。