航空發動機燃燒室噴嘴內燃油傳熱特性的數值研究

苗　輝，王慧汝，賈　真，扈鵬飛（.中航空天發動機研究院有限公司，北京0008；.沈陽發動機設計研究所，沈陽005）

航空發動機燃燒室噴嘴內燃油傳熱特性的數值研究

苗輝1，王慧汝1，賈真1，扈鵬飛2
（1.中航空天發動機研究院有限公司，北京100028；2.沈陽發動機設計研究所，沈陽110015）

摘要：采用數值模擬方法，研究航空發動機燃燒室噴嘴內燃油的傳熱特性：使用三維建模軟件對噴嘴進行精細的幾何建模，噴嘴模型包括旋流器葉片和隔熱套管等細節特征；采用商業CFD軟件對比研究不同噴嘴內燃油的傳熱特性，確定噴嘴隔熱套管的隔熱效果。結果表明，在所給工況條件下，隔熱套管能減少噴嘴桿部溫升約50%，對主油路溫升的抑制作用隨著主油路流量的增加而逐漸減小，對副油路燃油溫升的抑制作用比較明顯且基本不隨主油路流量變化。工程算法和數值模擬計算的噴嘴出口燃油溫度之間相差明顯。

關鍵詞：航空發動機；燃燒室；燃油噴嘴；熱防護；傳熱；結焦；數值模擬

1　引言

先進航空發動機的總壓比越來越高，壓氣機出口溫度也顯著增加。如推重比10一級發動機，壓氣機出口總溫可達1 000 K［1］。而壓氣機出口溫度的增加，導致噴嘴內燃油溫度隨之上升，增加了燃油沉積和結焦的風險［2-4］，為發動機帶來極大的安全隱患。很多學者［5-7］從燃油結焦特性方面對航空發動機燃燒室燃油噴嘴進行了研究，結果顯示燃油管道沉積率等參數與燃油溫度直接相關。

為控制燃油結焦風險，在燃燒室設計階段就要對噴嘴內的燃油溫度進行計算。工程上，設計者們往往采用傳熱學基本公式［8］進行簡單估算：把燃油噴嘴簡單分為桿部和頭部兩部分，分別選用實驗關聯式進行計算。對桿部換熱采用流體橫掠圓管公式，對頭部換熱則采用流體外掠平板公式［9］。

工程算法雖然簡單方便，但是由于其簡化模型與實際噴嘴結構相差較大，且噴嘴周圍流動又非常復雜，所以只能定性計算燃油溫度，無法獲得詳細的流場、溫度場特征和為噴嘴熱防護提供依據。采用詳細的數值模擬計算，可彌補工程算法的不足。

本文對實際噴嘴進行三維數值模擬，考察流場和溫度場特征，并通過計算結果對比，研究噴嘴隔熱套管的隔熱效果。

2　噴嘴模型和計算工況

數值模擬所用到的噴嘴模型為氣動噴嘴，具有主油路和副油路兩個燃油通道，包含兩級旋流器，各有6個葉片，如圖1所示。圖中加裝隔熱套管的噴嘴，其隔熱套管和噴嘴桿部本體之間有空氣間隙，以產生隔熱效果。數值計算中的噴嘴模型，包含內錐、分油環和葉片等精細化的結構特征。

圖1　燃油噴嘴結構模型Fig.1 Schematic views of fuel injector

計算條件模擬推重比10一級發動機典型燃燒室的工況參數，主流氣流速度為50 m/s，溫度約為810 K。假設主、副油路燃油進口溫度均為373.15 K，副油路流量恒定為60 L/min，主油路流量范圍為100～600 L/min。

3　數值計算方法

數值模擬采用三維計算域（圖2），包含燃油域、固體結構域和空氣域三部分，流體和固體結構的交界面為耦合面，采用流/熱耦合方式計算傳熱性能。由于主要研究對象為燃油在噴嘴內的傳熱性能，所以在數值計算中，不考慮噴嘴出口下游燃油的霧化、蒸發、摻混、燃燒過程對噴嘴內燃油傳熱的影響，將燃油域從噴嘴出口絕熱延長至計算域出口［9］。為節省計算資源，在噴嘴下游的計算域長度較短。

圖2　計算域示意圖（隔熱噴嘴）Fig.2 Computational domain

采用ICEM生成非結構化網格，對燃油域、固體結構域和空氣近壁區域進行加密，總網格數量約為1 000萬。采用商業軟件Fluent 14.5進行數值計算，空氣和燃油的物性參數為溫度的函數。湍流模型采用標準k-ε模型，壁面采用標準壁面函數處理?？諝庥嬎阌虿捎盟俣热肟?、壓力出口的邊界條件。

采用SIMPLEC壓力速度耦合算法，二階迎風差分格式，收斂精度為10-4。在網格無關性驗證中，計算采用的網格數目分別為1 000萬和1 800萬，對燃油出口平均溫度的計算結果相差不到1%，故下面選取節點數為1 000萬左右的網格進行計算。

4　結果與討論

4.1典型工況對比分析

圖3所示為主油路流量500 L/h、副油路流量60 L/h條件下，原始噴嘴和隔熱噴嘴固體結構區域的溫度云圖?？梢?，旋流器外環和隔熱套管的溫度接近氣流溫度，旋流器葉片溫度有明顯的過渡。傳熱過程中，旋流器葉片起到的類似翅片的作用，增強了主油路在頭部的換熱。

圖4示出了上面算例中主、副油路燃油溫升情況，圖中將噴嘴簡單分為桿部和頭部兩部分分別考察。由圖中可知，原始噴嘴中，副油路在桿部的溫升大于在頭部的溫升，而主油路在頭部的溫升大于在桿部的溫升。這是因為主油路在頭部的傳熱面積明顯比在桿部的大，而副油路在頭部的傳熱面積比在桿部的?。ǜ庇吐吩陬^部的橫截面浸潤周長與桿部的相當，但沿燃油流向頭部的長度比桿部的?。?。

隔熱套管的主要作用，是在從熱空氣流到燃油的傳熱過程中增加了一段空氣夾層，借助空氣低導熱系數的特點，減小了整個傳熱過程的傳熱系數。在算例中，隔熱套管使副油路的桿部溫升下降了48.0%，使主油路的桿部溫升下降了54.5%。故隔熱套管使桿部的總傳熱系數大約減小了50.0%，對溫升的抑制作用明顯，但對頭部換熱基本沒有影響。

4.2主油路燃油流量的影響

航空發動機在一個工作周期內，副油路燃油流量變化較小，而主油路燃油流量調節范圍較大。主油路的燃油流量變化對噴嘴溫度場和燃油溫升等有重要影響，如圖5所示?？梢?，隨著主油路流量的增加，副油路出口溫度基本恒定，主油路出口溫度顯著降低。因此，隔熱套管可明顯抑制副油路的溫升，約25%，且基本保持恒定。而隨著主油路流量的增加，隔熱套管對主油路溫升的抑制作用明顯降低。

4.3工程算法與數值模擬計算結果對比

圖6對比了工程算法與數值模擬方法的燃油出口溫度計算結果。工程算法的詳細推導過程見文獻［9］。由圖中可知，工程算法與數值模擬方法計算的燃油出口溫度具有相同的趨勢，隨著主油路流量的增加，主油路溫度降低且降低幅度越來越小，但數值模擬方法計算的出口燃油溫度明顯比工程算法的高。其原因可能是工程算法忽略了很多影響傳熱效果的細節特征，如旋流器翅片等。在副油路流量保持恒定時，主油路流量越小，工程算法與數值模擬方法計算結果之間的偏差越明顯。如當主油路流量為100 L/h時，工程算法計算的原始噴嘴主油路溫升，

圖3　典型工況下噴嘴結構溫度場Fig.3 The temperature contours of injector at typical conditions

圖4　典型工況下燃油溫升情況Fig.4 Temperature rise of fuel at typical conditions

圖5　主油路燃油流量對燃油出口溫度的影響Fig.5 Effect of flow rate of primary fuel channel on fuel outlet temperature

圖6　工程算法與數值模擬方法計算結果對比Fig.6 The comparison of engineering method and numerical simulation results

僅為數值模擬方法計算結果的28%左右。

5　結論

（1）隔熱套管對主、副油路在桿部溫升的抑制效果明顯。所給算例中，可使桿部溫升下降約50%。

（2）副油路流量不變時，隨著主油路流量的增加，隔熱套管對副油路溫升的抑制作用保持一定，約為25%；對主油路溫升的抑制作用逐漸降低。

（3）工程算法和數值模擬計算的噴嘴出口燃油溫度之間相差明顯。當主油路流量為100 L/h時，工程算法計算的原始噴嘴主油路溫升，僅為數值算法結果的28%左右。

參考文獻：

［1］ 彭澤琰，劉剛，桂幸民，等.航空燃氣輪機原理［M］.北京：國防工業出版社，2008.

［2］ 王洪銘.加熱后RP-3煤油性能指標變化的分析報告［D］.北京：北京航空航天大學，1998.

［3］ Spadaccini L J，Sobel D R，Huang H.Deposit formation and mitigation in aircraft fuels［J］.Journal of Engineering for Gas Turbines and Power，1999，123（4）：741—746.

［4］ Brown S P，Frederick Jr R A.Laboratory-scale thermal stability experiments on RP-1 and RP-2［J］.Journal of Propulsion and Power，2008，24（2）：206—212.

［5］ 賈春燕，王洪銘.航空發動機燃燒室噴嘴內部燃油結焦研究［J］.航空發動機，2011，37（5）：1—44.

［6］ Chin J S，Lefebvre A H.Influence of flow conditions on depositsfromheatedhydrocarbonfuels［R］.ASME 1992-GT-114，1992.

［7］ Katta V R，Jones E G，Roquenore W M.Modeling of depositions process in liquid fuels［J］.Combustion Science and Technology，1998，139（1）：75—111.

［8］ 楊世銘，陶文銓.傳熱學［M］.北京：高等教育出版社，2006.

［9］ 黃義勇.航空發動機主燃油噴嘴熱防護設計技術研究［D］.成都：電子科技大學，2011.

中圖分類號：V231.1

文獻標識碼：A

文章編號：1672-2620（2015）03-0030-03

收稿日期：2014-06-04；修回日期：2014-08-11

作者簡介：苗輝（1984-），男，河南滑縣人，高級工程師，博士，主要從事航空發動機傳熱與燃燒研究。

Numerical investigation of heat transfer characteristics of aero-engine combustor injector

MIAO Hui1，WANG Hui-ru1，JIA Zhen1，HU Peng-fei2
（1.AVIC Academy of Aeronautic Propulsion Technology，Beijing 100028，China；2.Shenyang Aero-engine Research Institute，Shenyang 110015，China）

Abstract：Numerical study on the heat transfer performance of aero-engine combustor injector was conducted.Elaborate geometry models with swirler blades and heat insulation casing（HIC）were presented. Comparative heat transfer researches were performed to estimate the effect of HIC on temperature rise of aviation kerosene through the injector.It is found that HIC decrease the temperature rise in bar department of injector by 50%in typical cases.Additionally，HIC has a more obvious effect on thermal protection for secondary fuel channel than for primary fuel channel.Furthermore，with increase in flow rate of primary fuel channel，the thermal protection influence of HIC on primary fuel channel become weaker，and no significant change of that influence on secondary fuel channel were found.Finally，engineering method would get a remarkable lower temperature rise than numerical simulation.

Key words：aero-engine；combustor；fuel injector；heat insulation injector；heat transfer；fuel coke；numerical simulation