谷 俊
(中航工業(yè)沈陽發(fā)動機設計研究所,沈陽 110015)
管殼式燃滑油熱交換器在航空發(fā)動機上的應用非常普遍,無論是軍、民用發(fā)動機(如CFM56系列發(fā)動機)都采用這種結構簡單可靠的裝置,為滑油系統(tǒng)提供冷卻或為燃油提供加溫需求。燃滑油散熱器換熱特性的初步計算是滑油系統(tǒng)熱分析及循環(huán)量設計的關鍵環(huán)節(jié),換熱性能決定了滑油系統(tǒng)在包線范圍內的熱負荷耐受能力。因此,掌握管殼式燃滑油散熱器的換熱特性計算方法對于滑油系統(tǒng)設計尤為重要。
管殼式燃滑油熱交換器的管內流動形式簡單,采用經(jīng)典管內流動和換熱計算方法即可達到足夠精度,而管外殼程的流動復雜,需要考慮的結構因素較多,換熱特性計算均為試驗擬合出的經(jīng)驗方法,目前常用的殼程換熱計算方法主要有2種:kern法及Bell-Delaware法[1],本文在此基礎上提出了1種按流動特性分段計算殼程換熱的方法。
本文通過以上3種殼程換熱計算方法對某型燃滑油散熱器進行了換熱特性的計算,并與2009~2010年所做的數(shù)次航空燃滑油散熱器換熱特性試驗進行對比驗證,得到具有足夠精度的換熱特性計算方法,應用于航空發(fā)動機潤滑系統(tǒng)的熱分析計算中。
某型燃滑油散熱器的結構如圖1所示。散熱器為典型的管殼式結構,內部管程流過燃油RP-3為冷卻介質,殼程通4050型滑油為被冷卻介質。散熱芯體的直徑為94 mm,管程數(shù)為5,殼程數(shù)為6,采用圓缺型折流板,缺口面積約20%,散熱管為915根Φ2 mm×0.2 mm的不銹鋼管,長度為200 mm,采用蜂窩式排列,管中心距為3 mm。
圖1 某型燃滑油散熱器結構
為對比散熱器散熱性能,需要計算并繪制不同燃油流量下的散熱器散熱量變化曲線,即換熱特性曲線。為與試驗數(shù)據(jù)進行對比,計算采用的散熱器進口燃、滑油溫度和流量均與試驗參數(shù)保持一致。
根據(jù)穩(wěn)態(tài)傳熱公式,燃、滑油之間換熱量應為
式中:Q為換熱量;K為總換熱系數(shù);A為換熱面積;Δtm為換熱的2種流體之間的平均溫差。
該散熱器采用冷熱逆流,燃油流量遠大于滑油,因此較小的進口溫差為滑油出口溫度和燃油進口溫度之差Δt1=Th,O-Tc,j而較大的溫差為滑油進口和燃油出口溫度之差Δt2=Th,O-Tc,j。由文獻[2]可知,當兩側的進口溫差,采用算術平均溫差即可達到足夠精度,因此該散熱器計算的Δtm采用算術平均溫差
式中:Ft為溫差修正系數(shù),查圖求取,此處Ft=1。由放熱和吸熱的熱量守恒可知
式中:mc為冷卻介質質量流量;Cp,c為冷卻介質定壓比熱容;Tc,o為冷卻介質出口溫度;Tc,i為冷卻介質進口溫度;mh為被冷卻介質質量流量;Cp,h為被冷卻介質定壓比熱容;ηc為熱損失系數(shù),通常取ηc=0.97~0.98。
當 K、A、Cp,c、Cp,h、mc、mh、ηc、Th,i、Tc,i已知時,聯(lián)立式(1)~(3)組成方程組可求得換熱量Q。該方程組中除K和A外的其他參數(shù)與散熱器結構無關,只與進口條件有關。因此,求得K值是計算散熱器換熱量的關鍵。
根據(jù)散熱器結構可知,該散熱器為管殼式換熱器,散熱管是冷熱流體交換的界面,其內為冷流體,其外為熱流體,在不考慮輻射換熱的前提下,該處的熱傳遞路線如下:熱流體與散熱管外壁之間產生對流換熱,散熱管外壁與內壁之間產生熱傳導,散熱管內壁與冷流體之間產生對流換熱,因此根據(jù)對流換熱和熱傳導的穩(wěn)態(tài)傳熱公式可求得以管外側面積為基準的換熱系數(shù)K
式中:do為散熱管外徑;di為散熱管內表面換熱面積;λw為散熱管材料的導熱系數(shù);ro,ri分別為管外壁、內壁污垢熱阻;α0、αi分別為殼程、管程對流換熱系數(shù)。
由于散熱器內流動的復雜性,對于管內對流換熱系數(shù)(又稱管程換熱系數(shù))和管外對流換熱系數(shù)(又稱殼程換熱系數(shù))有很多基于傳熱學理論和試驗所得的計算方法。計算前首先確定冷熱流體的物性,冷卻流體為RP-3燃油,熱流體為4050滑油,二者物性[3-4]如下:
RP-3 燃油:密度 ρ=0.7796-0.7606×(t-20)×10-3,比熱容Cp=Cp(01+αt),導熱系數(shù))4/3·λ20,黏度 lg lg(v+0.73)=8.908-3.8265lg T。
4050 滑油:ρ=0.9729-0.000035t,Cp=1.76+0.0031t,λ=0.1591-0.000134t,ln ln(v+0.6)=21.52-3.54ln T。
由上述試驗公式可知,物性隨介質溫度的變化而改變,對于溫度變化的流體,工程上取其平均溫度為定性溫度進行計算。在計算中以試驗數(shù)據(jù)來確定燃油和滑油的定性溫度:燃油定性溫度為90℃,滑油定性溫度為125℃。通過以上公式并查取手冊[3-4]中相關數(shù)據(jù),計算得出:
RP-3燃油在90℃下的物性為:
4050滑油在125℃下的物性為:
計算管內換熱系數(shù)αi,根據(jù)文獻[5],對于Re>10000的湍流,管內對流換熱系數(shù)根據(jù)(Sieder-Tate)實驗關聯(lián)式[6]計算
式中:μi為按流體均溫為定性溫度的動力黏度;μw為按散熱管內壁溫為定性溫度的動力黏度。
殼程的流動非常復雜,特別是在折流板存在的情況下,因此對于殼程換熱系數(shù)的計算方法均為經(jīng)驗公式,主要采用 kern 法[2]和 Bell-Delaware 法[1,7]。
2.2.1 Kern法
式中:μi為殼程傳熱系數(shù);μw為管外液體熱傳導率;de為特征管徑尺寸。
式中:Pt為管中心距;d0為管外徑。
從以上公式可見,kern法計算過程相對簡單,將殼程流動簡化為1種管外橫流的換熱。
2.2.2 Bell-Delaware法
Bell-Delaware法根據(jù)實際的流動形式,對殼程流動的橫流區(qū)、折流板、側泄分別進行考慮,由試驗得到不同結構對理想傳熱系數(shù)的校正因子,具體計算過程如下。
2.2.2.1 理想殼程傳熱系數(shù)
根據(jù)Bell-Delaware法,理想傳熱系數(shù)公式為
2.2.2.2 已知條件及輔助計算
根據(jù)Bell-Delaware法,為計算換熱系數(shù)的修正因子,首先針對該結構散熱器結構和流動參數(shù)進行如下計算。
(2)殼程橫流面積 Sm=Lbc[Lbb+n×(Ltp-d)]=26.5×[7.5+30×(3-2)]=993.75 mm2。
(5)流動方向有效橫流管排數(shù)Ntcc=17和缺口區(qū)管排數(shù)Ntcw=8。
(6)管束與殼體間的C流路和F流路的面積參數(shù)Sb=Lbc(DS-Dotl)=26.9×(98-90.5)=198.75 mm2,F(xiàn)sbp=Sb/Sm=20.3%。
2.2.2.3 殼程傳熱修正因子
考慮結構上的折流板、旁路、泄漏等因素,Bell-Delaware法對理想的殼程橫流管束傳熱系數(shù)通過每一路的修正因子加以修正,故殼程傳熱系數(shù)
(1)弓形折流板缺口修正因子Jc=0.55+0.72Fc=0.989。
(3)管束旁路傳熱修正因子Jb=exp{-CbhFsbp·[1-(2rss)1/3]},根據(jù)流動類型選擇Cbh的經(jīng)驗值,對于Re>100 的情況為 1.35,rss=Nss/Ntcc,由于該散熱器無旁路密封擋板所以為 0,Jb=e-1.35×0.203=0.76。
(4)層流下逆向溫差修正因子Jr在Re<100時才考慮,而對于該散熱器不適用,可不考慮。
2.2.3 分段法計算殼程換熱系數(shù)
根據(jù)文獻[8],殼程的流動可以分解為如圖2所示的各條支流,在不考慮泄漏、旁路等因素的情況下,殼程理想流路分為2部分B、F和C流動,B、F部分是在折流板缺口區(qū),近似為順管路方向的管外對流,C部分在折流板中間區(qū)近似于管外純橫流,在這2部分流動的基礎上加以修正系數(shù),可以得到較滿意的結果,因此下文對殼程分段進行了傳熱系數(shù)的計算得到理想殼程傳熱系數(shù)。
根據(jù)文獻[9],流經(jīng)管排的純橫流傳熱系數(shù)為
在缺口區(qū)將流動看成管外順流,傳熱系數(shù)按非圓形截面的管內流動換熱計算
由以上推算計算出管殼程的傳熱系數(shù),其中管程換熱計算均采用式(5),而對于較復雜的殼程流動采用了3種計算方法,根據(jù)式(4)分別計算其總換熱系數(shù)K。
為對比計算和試驗的結果,采用了與試驗條件相同的進口條件進行散熱器的換熱特性計算(滑油進口溫度140℃,燃油進口溫度80℃),該散熱器基于換熱管外側的總換熱面積0.91094 m2,由于試驗過程散熱器有保溫措施,因此計算中不考慮殼程與外界環(huán)境換熱,聯(lián)立式(1)~(3),計算Q值及換熱特性η=,并繪制換熱特性隨燃油流量的變化曲線。
在2009~2010年,在動力傳輸航空科技重點實驗室的A611試驗臺上進行了3套該型散熱器的換熱特性試驗。散熱器放置在保溫箱中,有單獨的滑油和燃油循環(huán)系統(tǒng)保證燃滑油的循環(huán)供應,并通過加溫、冷卻和溫控系統(tǒng)保持散熱器進口燃滑油溫度的穩(wěn)定。試驗過程中采用4050滑油,流量穩(wěn)定在50 L/min,進口溫度為140℃;采用RP-3燃油,流量在1000~7000 L/h范圍內可調,進口溫度保持在80℃。進行了不同燃油流量下的換熱特性試驗,并繪制試驗曲線,試驗現(xiàn)場照片如圖3所示。
圖3 試驗現(xiàn)場
根據(jù)以上介紹的3種殼程換熱系數(shù)計算方法,分別進行某型燃滑油散熱器在不同燃油流量下的換熱特性計算,繪制換熱特性曲線,與試驗數(shù)據(jù)進行對比,如圖4所示。
圖4 某型燃滑油散熱器換熱特性曲線對比
通過圖4中換熱特性曲線的對比可得:
(1)殼程采用kern法計算所得的散熱器換熱特性曲線整體較高,隨著燃油流量的增加,與試驗值的偏差加大,最大偏差量為25%,該方法計算的結果為理想狀態(tài)下的換熱特性;
(2)殼程采用Bell-Delaware法計算的燃滑油散熱器換熱特性曲線整體偏低,隨著燃油流量的增加,其偏差量有所增加,在散熱器最大燃油流量下單位換熱量低于試驗值約12%;
(3)殼程采用分段計算并加以修正因子的方法所得的換熱特性曲線更接近試驗所得曲線,最大偏差約3%,與試驗曲線有良好的整體相似度。
(1)對于管殼式燃滑油散熱器的計算主要包括管內流動換熱系數(shù)和殼程換熱系數(shù)2部分,由于散熱管內的流動簡單,其對流換熱系數(shù)采用(Sieder-Tate)實驗關聯(lián)式即可以達到良好的工程精度,而管外殼程由于受折流板、旁路和散熱管束的結構因素影響,流動極為復雜,很難準確分析計算,需要通過試驗數(shù)據(jù)的擬合修正才能得到相對準確的結果。對于結構類似的管殼式散熱器,采用經(jīng)試驗修正的計算方法可以滿足航空發(fā)動機潤滑系統(tǒng)熱平衡計算的需求。
(2)以某型燃滑油散熱器為例,采用3種不同的殼程換熱計算方法進行計算,并與試驗結果進行對比可知: 采用Kern法計算殼程參數(shù)方法簡單、應用方便,但計算結果為理想狀態(tài)下的換熱特性,高于實際換熱能力; 采用Bell-Delaware法對結構影響因素考慮較為全面,其計算所得的換熱特性略低于試驗數(shù)據(jù),但計算過程較為繁瑣,且部分經(jīng)驗公式的適用范圍有限;采用分段計算殼程換熱,并代入流動、側泄、非等跨等結構修正因子的方法所得的換熱特性曲線更接近試驗曲線,過程相對簡單,可在采用管殼式燃滑油散熱器的航空發(fā)動機潤滑系統(tǒng)熱分析計算中應用。
由于試驗所用散熱器均為新件,本次計算未考慮污垢熱阻的影響。在散熱器使用中,隨工作時間增加,污垢熱阻將造成散熱器換熱特性衰減,其形成過程復雜且在不同介質和工作環(huán)境下差異較大,應進行后續(xù)試驗研究。
[1]Bell K J.Final report of the cooperative research program on shell and tube heat exchanger[J].ASHRAE Transactions,1991:79-89.
[2]錢頌文.換熱器設計手冊[M].北京:化學工業(yè)出版社,2002:55-108.QIAN Songwen.Heat exchanger design manual[M].Beijing:Chemistry Industry Publisher,2002:55-108.(in Chinese)
[3]劉濟灜.中國噴氣燃料 [M].北京:中國石化出版社,1991:45-70.LIU Jiying.Jet fuel of China[M].Beijing:China Petrochemical Industry Publisher,1991:45-70.(in Chinese)
[4]楊九高.航空潤滑油 [M]//林基恕.航空發(fā)動機設計手冊:第12冊.北京:航空工業(yè)出版社,2002:447-456.YANG Jiugao.Aero lubrication[M]//LIN Jisu.Aero engine design manual, Twelfth volume. Beijing:Aero Industry Publisher,2002:447-456.(in Chinese)
[5]楊世銘.傳熱學[M].北京:高等教育出版社,1987:41-101.YANG Shiming.Heat transfer subject[M].Beijing:Higher Education Publisher,1987:41-101.(in Chinese)
[6]毛希瀾.換熱器設計 [M].上海:上??萍汲霭嫔纾?988:51-55.MAO Xilan.Heat exchanger design[M].Shanghai:Shanghai Science and Technology Publisher,1988:51-55.(in Chinese)
[7]閻皓峰,甘永平.新型換熱器與傳熱強化[M].北京:宇航出版社,2001:32-100.YAN Haofeng,GAN Yongping.New type heat exchanger and heat transfer indentify[M].Beijing:Yu Hang Publisher,2001:32-100.(in Chinese)
[8]Bell K.J Delaware.Method for shell side design[J].Heat Exchanger Thermal Hydraulic Fundamentals and Design.Washington D.C:Hemisphere/McCraw-Hill,1981:581-618.
[9]Webb R L.Principles of enhanced heat transfer[M].New York:John wiley and sons,Inc,1994 :373-426.
[10]蘭州石油研究所.換熱器[M].北京:中國石化出版社,1988:51-75.Lanzhou Petroleum Research Institute.Heat exchanger[M].Beijing:China Petrochemical Industry Publisher,1988:51-75.(in Chinese)