張忠利，周立新，張蒙正

（西安航天動力研究所，陜西 西安710100）

0 引言

為了提高性能，高空發(fā)動機通常采用大面積比噴管，這給發(fā)動機地面試驗帶來了很大困難。高空發(fā)動機在地面整機熱試車過程中，由于外界環(huán)境壓力遠高于噴管出口壓力，會在噴管內(nèi)產(chǎn)生激波，使總壓降低、傳熱惡化并伴隨振動、氣流分離或火焰偏擺等現(xiàn)象，對準確評定其在真空環(huán)境下的推力及比沖性能以及工作可靠性帶來嚴重影響。為此，高空發(fā)動機在地面試驗時需用引射裝置以達到噴管滿流的目的。高性能引射器設(shè)計成為高空發(fā)動機研制中的關(guān)鍵技術(shù)之一。

引射器正常工作時，經(jīng)受著高溫燃氣與壁面的強迫對流換熱和輻射換熱，因此引射器的冷卻方案設(shè)計[1-3]也十分重要。引射器冷卻結(jié)構(gòu)設(shè)計時通?？紤]燃氣流量、燃燒室燃氣壓力、燃氣溫度及引射器的氣動結(jié)構(gòu)等參數(shù)。本文介紹引射器冷卻結(jié)構(gòu)設(shè)計時所進行的傳熱特性分析方法，同時以地面試驗用引射器為算例[4]對其進行傳熱分析，分析其地面試車時引射器不被燒蝕所需的冷卻水流量范圍。

1 傳熱分析模型

本文的研究對象是再生冷卻式引射器，冷卻通道數(shù)量依據(jù)流阻損失和能否可靠冷卻來確定，因為引射器為地面設(shè)備，為了節(jié)約試驗費通常情況下冷卻劑選取常溫自來水。對于均勻分布的冷卻通道，通常認為周向溫度基本一致，因此可以不考慮周向的傳熱。在傳熱分析時通常取一條通道作為研究對象，分析其溫度變化。

對于圖1的冷卻通道，一邊為冷卻水進口、另一邊為冷卻水出口。為了方便傳熱分析，可以不考慮集液環(huán)內(nèi)流體帶走的能量，原因有二：其一是這樣得到的壁溫及冷卻水溫升略高實際測量值，如果計算分析結(jié)果可以可靠冷卻，那么實際工作時引射器肯定能被可靠冷卻；另外一個原因是集液環(huán)內(nèi)液體帶走的能量經(jīng)過精確傳熱分析后僅占總傳熱量的1%，因此可以忽略集液環(huán)的影響，把冷卻通道入口作為冷卻劑的入口可以滿足工程設(shè)計要求。

圖1 物理模型Fig.1 Physical model

對圖1的物理模型應(yīng)用二維傳熱計算方法[5-7]對模型簡化處理后，可以得到問題的近似解，但應(yīng)基于一定的假設(shè)條件。首先，計算冷卻劑與冷卻通道之間的對流換熱時把肋簡化為一維散熱片處理，這種處理方法會給計算結(jié)果帶來一些偏差，但作為工程設(shè)計計算完全可以接受；其次，不考慮冷卻劑與外殼體之間的換熱，因為該換熱量與冷卻劑和內(nèi)壁之間的換熱量相比很小，可以忽略；再次，認為內(nèi)壁與冷卻通道是線型變化的，考慮通道流通面積與換熱面積大小。另外，計算時應(yīng)將冷卻水的最高溫度控制在當?shù)貕毫Φ姆悬c以下，全部按照液體性質(zhì)進行計算，傳熱計算分析模型見圖2。

圖2 傳熱分析模型Fig.2 Analytical model of heat transfer

如圖2所示，將從燃氣到冷卻劑之間的換熱看作是一個串聯(lián)換熱問題，燃氣與內(nèi)壁面之間的熱量傳遞由對流換熱和輻射換熱兩部分組成；引射器壁內(nèi)部為導熱；冷卻劑與冷卻通道之間為對流換熱，這種熱量傳遞關(guān)系可以表示為：

式中：qtot為總換熱的熱流密度；qr為輻射熱流密度；qcv為對流熱流密度，其中對流換熱是主要形式；hg為燃氣與引射器壁的對流換熱系數(shù)；hl,eq為冷卻劑與冷卻通道的等效換熱系數(shù)；Tad為絕熱壁溫；δ為引射器的內(nèi)壁厚；λ為壁材料的導熱系數(shù)。

計算冷卻劑與冷卻通道對流換熱時，引入肋片效率概念，其等效對流換熱系數(shù)為：

式中：A1為肋片表面積；A2為肋基未裝肋的光滑表面積；A0為未裝肋時壁表面積；hl為對流換熱系數(shù)，根據(jù)管內(nèi)湍流試驗關(guān)聯(lián)式確定；η為肋片效率，對于等截面直肋，η計算公式為：

在引射器壁的冷卻通道內(nèi)均勻劃分計算節(jié)點，冷卻通道內(nèi)相鄰節(jié)點之間通過流體流動傳遞的熱量聯(lián)系起來，即認為冷卻通道某一節(jié)點內(nèi)，引射器壁傳給冷卻劑的熱量全部傳給下一節(jié)點內(nèi)的冷卻劑，但在引射器壁內(nèi)部不考慮熱量沿軸向傳遞，因此計算程序嚴格地講應(yīng)該準二維。計算時，邊區(qū)燃氣總溫依據(jù)流場分析結(jié)果。

1.1 燃氣對流熱流分析

燃氣對流換熱[9]是引射器內(nèi)燃氣向壁傳熱的主要形式。在引射器內(nèi)，對流熱流在燃氣進口處常占總熱流的80%以上，在引射器出口附近可達98%以上。因此確定燃氣對流熱流，是分析引射器傳熱狀況和采取正確冷卻措施的首要工作。

巴茲在計算近壁層燃氣與壁的換熱系數(shù)時，用下列方法得到的壁溫與實際測量值比較接近。

1.2 燃氣輻射熱流分析

引射器內(nèi)壁受燃氣中二氧化碳和水蒸氣的氣體輻射，以及燃燒過程產(chǎn)生的自由碳粒子的熱輻射作用。氣體輻射對波長有明顯的選擇性，二氧化碳的主要輻射譜帶在波長2.7 μm，4.3 μm和15 μm附近；水蒸氣的主要輻射譜帶在波長1.9 μm，2.8 μm，6.7 μm 和 21 μm 附近。氣體輻射的波長均在紅外輻射范圍內(nèi)，所以有時也稱非發(fā)光輻射。固態(tài)自由碳粒的輻射具有連續(xù)的光譜，因此稱為發(fā)光輻射。燃氣對引射器壁的輻射換熱熱流密度可由公式（9） 計算：

式中：σ為斯蒂芬-波爾茲曼常數(shù)；εw為引射器壁的發(fā)射率（黑度）；εg為燃氣的發(fā)射率（黑度）；Tg和Twg分別為燃氣和引射器內(nèi)壁的溫度，K。

2 引射器傳熱分析

2.1 圓柱型引射器傳熱分析

對于地面試車用引射器，考慮到研發(fā)成本，采用的引射器為夾套式水冷圓柱型結(jié)構(gòu)，如圖3所示，流量參數(shù)見表1。

圖3 圓柱型引射器冷卻結(jié)構(gòu)簡圖Fig.3 Cooling structure of cylindrical ejector

表1 引射器流量參數(shù)Tab.1 Flux parameters of ejector

從圖3可看出，該引射器采用水冷夾套結(jié)構(gòu)，在夾套內(nèi)焊接60條筋板，筋板是內(nèi)外壁聯(lián)接的橋梁，兼顧支撐內(nèi)外壁和分流冷卻水作用，使冷卻水在夾套內(nèi)不致產(chǎn)生死區(qū)，從而產(chǎn)生引射器局部燒蝕。依據(jù)1節(jié)的分析方法，應(yīng)用圖3和表1的數(shù)據(jù)，計算得到圓柱型引射器氣壁溫Twg、液壁溫Twl及冷卻水溫度Tl，對流熱流qcv、輻射熱流qr及總熱流qtot沿軸線方向的變化過程，分析結(jié)果見圖4～5和表2，起點為冷卻水的入口。

圖4 引射器沿軸線溫度變化曲線Fig.4 Temperature variation of ejector along axis

圖5 引射器沿軸線熱流變化曲線Fig.5 Heat flux variation of ejector along axis

從圖4～5得到，引射器在冷卻水的作用下，隨著冷卻水流量的增加，氣壁溫Twg在不斷減小，液壁溫Twl也在不斷減小，冷卻水出口溫度也在不斷減?。豢偀崃髟诓粩嗌?。從以上分析結(jié)果來看，該引射器的冷卻水流量在300～600 kg/s范圍內(nèi)工作時是安全的。

表2 圓柱型引射器傳熱分析結(jié)果Tab.2 Analysis results of heat transfer of cylindrical ejector

2.2 二次喉道型引射器傳熱分析

二次喉道型引射器也采用夾套式水冷圓柱形結(jié)構(gòu)，如圖6所示，流量參數(shù)見表1。

圖6 二次喉道型引射器冷卻結(jié)構(gòu)簡圖Fig.6 Cooling structure of secondary-throat ejector

依據(jù)1節(jié)的分析方法，應(yīng)用圖6和表1的數(shù)據(jù)，計算得到二次喉道型引射器氣壁溫Twg、液壁溫Twl及冷卻水溫度Tl，對流熱流qcv、輻射熱流qr及總熱流qtot沿軸線方向的變化過程，分析結(jié)果見圖7～8和表3，起點位置為冷卻水的入口。

從圖7～8得到，引射器在冷卻水的作用下，隨著冷卻水流量的增加，氣壁溫Twg在不斷減小，液壁溫Twl也在不斷減小，冷卻水出口溫度也在不斷減小；總熱流在不斷升高。從以上分析結(jié)果來看，該引射器的冷卻水流量在300～600 kg/s范圍內(nèi)工作時是安全的。

圖7 二次喉道型引射器沿軸線溫度變化曲線Fig.7 Temperature variation of secondary-throat ejector along axis

圖8 引二次喉道引射器沿軸線熱流變化曲線Fig.8 Heat flux variation of secondary-throat ejector along axis

表3 二次喉道型引射器傳熱分析結(jié)果Tab.3 Analytical results of heat transfer of secondary-throat ejector

3 試驗研究

某型液體火箭發(fā)動機是正在研制的二級發(fā)動機，該發(fā)動機的室壓高、噴管面積比大，出口壓力遠低于大氣壓，因此在地面試車時需使用引射器。為了簡化設(shè)計、節(jié)約研制成本，試車時采用了圓柱型引射器。冷卻結(jié)構(gòu)與圓柱型結(jié)構(gòu)相同，引射器安裝簡圖見圖9。

圖9 引射器安裝簡圖Fig.9 Installation diagram of ejector

該型引射器采用自來水作為冷卻劑，冷卻水從引射器的后部進入，前端流出。冷卻水流量約為600 kg/s，冷卻水入口溫度約為20℃。冷卻水測試數(shù)據(jù)見表4，其中圖10為某三次試車冷卻水溫的測試值變化曲線。

從表4和圖10可看出，對于該型號發(fā)動機試應(yīng)用引射器試車時，當冷卻水流量為600 kg/s時冷卻水溫升約為10℃，與第2節(jié)分析結(jié)果一致。

表4 冷卻水測試數(shù)據(jù)Tab.4 Tested values of cooling water

圖10 試車時冷卻水溫測試值隨試車時間變化曲線Fig.10 Variation of cooling water temperature with hot fire time

4 結(jié)論

1） 建立了引射器傳熱分析模型，應(yīng)用數(shù)值計算方法得到了圓柱型和二次喉道型引射器在發(fā)動機正常工作時冷卻水流量范圍為300～600 kg/s時的冷卻液溫升、液壁溫、氣壁溫、熱流等參數(shù)，從這些參數(shù)的變化情況來看引射器結(jié)構(gòu)是安全的。

2） 本文分析的引射器結(jié)構(gòu)在發(fā)動機正常工作時冷卻液流量為600 kg/s，試車時測量的冷卻液溫升值約為10℃，仿真分析得到的冷卻液溫升值約為9.3℃，二者一致。

[1]徐萬武,王振國.環(huán)型超聲速空氣引射器零二次流流場數(shù)值研究[J].推進技術(shù),2003,(2):36-39.

[2]吳繼平,陳健,王振國.帶二次流的多噴管超聲速引射器性能實驗研究[J].強激光與粒子束,2007,19(9):1439-1443.

[3]吳繼平,王振國.第二喉道超聲速引射器啟動性能理論研究[J].航空動力學報,2008,23(5):803-803.

[4]張忠利,超音速二次喉道擴壓器氣動特性研究[J].火箭推進,2001,27(3):14-22.

[5]LU Kuen tzong,KOU Hong Sen,LAN Ting hsien,Geometrically and thermally non.optimum ejector heat pump analysis[J].Energy Converts MGMT,1993,34(12):1287-1297.

[6]HSU C T.Investigation of an ejector heat pump byanalytical methods[R].USA:Oak Ridge National Lab,1984.

[7]HAMMER R M.An investigation of an ejector compression refrigeration cycle and its application to heating,cooling and energy conservation[D].USA:The University of Alabam,1978.

[8]KENAN J H.An investigation of ejector design by analysis and experiment[J].J.App1.Mech Trans AMSE,1950(5):299-309.