不同壓力條件下水平平板的表面換熱實驗研究

王 晶，徐 瑩，丁 立，郄殿福（. 北京衛(wèi)星環(huán)境工程研究所 可靠性與環(huán)境工程技術(shù)重點實驗室，北京 00094；. 北京航空航天大學(xué) 生物與醫(yī)學(xué)工程學(xué)院，北京 009）

?

不同壓力條件下水平平板的表面換熱實驗研究

王 晶1，徐 瑩2，丁 立2，郄殿福1
（1. 北京衛(wèi)星環(huán)境工程研究所 可靠性與環(huán)境工程技術(shù)重點實驗室，北京 100094；2. 北京航空航天大學(xué) 生物與醫(yī)學(xué)工程學(xué)院，北京 100191）

摘要：為了解不同壓力下水平平板的氣體對流換熱變化情況，搭建了一個提供不同氣壓和環(huán)境溫度的實驗艙，開展了在不同壓力（0.1Pa、0.1kPa、0.2kPa、0.5kPa、1 kPa、10kPa、50kPa和常壓）與幾種加熱量（75、150、300W/m2）組合條件下的水平平板換熱實驗研究。通過對輻射換熱和自然對流換熱的比較，得到不同壓力下氣體的對流換熱系數(shù)。結(jié)果表明：對流換熱系數(shù)在環(huán)境氣體壓力小于1kPa時非常小，而在1kPa以上時才較大；在大于1kPa時，對流換熱系數(shù)隨壓力的升高呈二次方增加。

關(guān)鍵詞：低氣壓；水平平板；對流換熱系數(shù)；輻射換熱；實驗研究

http://www.bisee.ac.cnE-mail: htqhjgc@126.comTel:(010)68116407, 68116408, 68116544

0　引言

飛行器在升空、停留和返回過程中所經(jīng)歷的環(huán)境參數(shù)變化大，尤其是艙內(nèi)設(shè)備的自然對流換熱系數(shù)變化較大，極易導(dǎo)致飛行器及機載設(shè)備出現(xiàn)“過熱”、“過冷”和“熱分層”等現(xiàn)象[1-2]。Devienne等提出了在Knudsen數(shù)基礎(chǔ)上將低壓熱交換分為4個壓力檔進(jìn)行對流換熱計算的方法，但適用范圍有限[3]；Nelson和Bevans建立了一個基于有限壓力范圍的無因次方程估算對流換熱[4]；Warner等研究了球體和圓柱體在低壓狀態(tài)下的對流換熱[5]；Kyte等研究了在低壓下圓柱體的對流換熱[6]；Churchill 和Chu通過對豎直等溫平板表面自然對流換熱實驗數(shù)據(jù)的整理得到了瑞利數(shù)覆蓋較大的自然對流換熱計算準(zhǔn)則與計算式[7]。這些研究多是建立在某幾個低壓力值條件下的自然對流參數(shù)研究，不具有系統(tǒng)性。

本文將開展不同氣體壓力條件下水平平板的自然對流和輻射換熱實驗研究，旨在系統(tǒng)地描述不同壓力下的對流換熱系數(shù)。

1　實驗方案

1.1實驗原理

對于密封艙內(nèi)水平放置的平板，從理論上講，其通過表面向環(huán)境換熱有2種方式：空氣對流換熱和輻射換熱。即平板換熱量為

其中，Qo、Qc和Qr分別為平板總換熱量、自然對流換熱量和輻射換熱量，單位為W。Qo可采用夾在兩塊平板之間的電加熱薄膜得到，即

其中：U為加熱電壓；R為加熱膜電阻。

從理論上講，輻射換熱量需通過測量平板表面溫度和熱沉溫度得到，即

式中：σ為斯忒藩-玻耳茲曼常量，Ts、Tw分別為平板表面和熱沉的溫度；εs、εw分別為平板表面和熱沉的發(fā)射率；A為平板表面積；Xs,w為平板對熱沉的輻射角系數(shù)；F為熱沉表面積。顯然，式(3)中的參數(shù)很多，要準(zhǔn)確地測量這些參數(shù)難度很大。為此，本文假設(shè)這些參數(shù)與氣體壓力無關(guān)，即認(rèn)為輻射換熱參數(shù)不隨實驗艙內(nèi)的氣體壓力變化，并將這些參數(shù)綜合設(shè)定為C，將式(3)簡化為

隨著艙內(nèi)氣體壓力的減小，空氣含量愈來愈少，則平板通過空氣傳遞的熱量也越來越小，而輻射換熱逐漸增強。當(dāng)艙內(nèi)為真空時，對流換熱減弱至0，只有輻射換熱，即輻射換熱量等于總換熱量。這樣就可以通過式(4)反算出C，再利用確定后的C值及前面的式(1)、(2)、(4)計算出不同壓力下的自然對流換熱量Qc。最后，就可以根據(jù)自然對流換熱公式計算得到自然對流換熱系數(shù)

其中To為空氣溫度。

1.2實驗裝置

本研究的實驗裝置及測試原理如圖1所示。

圖1　實驗裝置及測試原理Fig. 1　Principle of experimental platform

實驗艙（φ500mm×500mm）外接真空泵，用于氣壓調(diào)節(jié)。艙內(nèi)氣壓監(jiān)控主要由ZDY-1和ZDZ-52T-II真空計實現(xiàn)，其中ZDY-1的測量范圍1～105Pa，測量精度10Pa；ZDZ-52T-II的測量范圍10-2～105Pa，測量精度10-2Pa。平板由2塊光滑正方形鋁板（100mm×100mm×1mm）和1片電加熱膜（100mm×100mm，電阻47?）組成，其中電加熱膜被夾在2塊表面經(jīng)噴黑漆處理（ε=0.85）的正方形鋁板中間。電加熱膜的加熱功率由TDGC2交流變壓器和PM6803A-3A功率計（精度為0.01W）進(jìn)行控制。艙內(nèi)環(huán)境溫度控制由加熱和冷卻2部分組成：艙壁外均勻纏繞著電加熱材料，以加熱升高艙壁溫度；艙壁外均勻布置氮氣管路，液氮從杜瓦瓶流入另一個大的壓力罐，氣化后變成低溫氮氣，再與空氣壓縮機吹出的氣體進(jìn)行融合，形成溫度可調(diào)的穩(wěn)定低溫氣體，以冷卻降低艙壁溫度。

針對加熱后的平板上可能出現(xiàn)最大溫度差異的部位，在平板的上下表面共布置了6個測溫點（Pt100熱電偶，其測溫精度為0.1K），具體位置參見圖2所示。取6點溫度的平均值作為平板的表面溫度。艙內(nèi)空氣測溫點布置如圖3所示，用以測出從平板到艙壁的溫度梯度。

圖2　平板的測溫布點Fig. 2　Layout of temperature measurement for the plate

圖3　艙內(nèi)空氣測溫點的布置Fig. 3　Arrangement of environmental temperature measurement points in the cabin

1.3實驗設(shè)計

參考航空航天設(shè)備熱流密度和環(huán)境溫度的情況，實驗采用3種加熱功率，即1.5W（75W/m2），3W（150W/m2），6W（300W/m2）；5個艙壁溫度，即-40，-25，0，25，40℃；8個艙內(nèi)壓力，即常壓，50，10，1，0.5，0.2，0.1kPa及0.1Pa。共120個工況，其中真空測量用于確定輻射換熱的綜合參數(shù)C。

1.4實驗方法

在鋁板上布置相應(yīng)的溫度測量點后，按圖3所示的位置關(guān)系放置在實驗艙內(nèi)。接通電源，電壓調(diào)節(jié)器、真空泵等設(shè)備開始工作，調(diào)節(jié)加熱板的加熱功率使之到達(dá)1.5W，控制艙壁溫度，將艙內(nèi)壓力抽至真空；當(dāng)溫度穩(wěn)定20min后，記錄第1個工況下的平板表面溫度、艙內(nèi)空氣溫度和艙壁溫度。然后逐一進(jìn)行每個工況的測量實驗，直至完成所有實驗工況。所有工況實驗時空氣溫度為20℃。

2　實驗結(jié)果及分析

2.1輻射換熱和對流換熱量的比較

圖4是艙壁溫-40℃時，不同加熱狀態(tài)下，平板的輻射和對流換熱量隨艙內(nèi)壓力變化的情況。可以看出：艙內(nèi)壓力在1kPa以下時，只有輻射換熱，而對流換熱量很小，特別在0.2kPa以下，可以近似認(rèn)為只有輻射換熱。當(dāng)氣壓大于1kPa時，輻射換熱量逐漸減小，對流換熱量逐漸增加。當(dāng)艙壁溫為-25、0、25和40℃時，平板的2種換熱量呈同樣的變化趨勢（圖5～圖8）。但對于1.5W的加熱功率，輻射換熱量與對流換熱量的變化受環(huán)境溫度變化的影響比較大，平板的溫度呈現(xiàn)出與3W和6W加熱功率時不同的變化趨勢（圖9）。

圖4　艙壁溫-40℃時的輻射換熱量與對流換熱量Fig. 4　Radiant and convective heat transfer when the temperature of cell wall is -40℃

圖5　艙壁溫-25℃時的輻射換熱量與對流換熱量Fig. 5　Radiant and convective heat transfer when the temperature of cell wall is -25℃

圖6　艙壁溫0℃時的輻射換熱量與對流換熱量Fig. 6　Radiant and convective heat transfer when the temperature of cell wall is 0℃

圖7　艙壁溫25℃時的輻射換熱量與對流換熱量Fig. 7　Radiant and convective heat transfer when the temperature of cell wall is 25℃

圖8　艙壁溫40℃時的輻射換熱量與對流換熱量Fig. 8　Radiant and convection heat transfer when the temperature of cell wall is 40℃

2.2對流換熱系數(shù)

圖10是艙壁溫-40℃且平板處于加熱狀態(tài)下，自然對流換熱系數(shù)隨艙內(nèi)氣壓變化的情況。可以看出：艙內(nèi)壓力在1kPa以下時，自然對流換熱系數(shù)很小，特別在0.2kPa以下時近似為0。

氣壓大于1kPa時，對流換熱系數(shù)逐漸增大。當(dāng)艙壁溫為-25、0、25、40℃時，平板的自然對流換熱系數(shù)呈同樣的變化趨勢（圖11～圖14）。但對于1.5W的加熱功率，平板的對流換熱系數(shù)的變化受到環(huán)境溫度影響很大，艙壁溫-40℃時，還出現(xiàn)負(fù)數(shù)的錯誤現(xiàn)象；只有當(dāng)艙壁溫度較高時，不同加熱量下的對流換熱系數(shù)變化趨勢才一致；艙壁溫為25 和40℃時，3種加熱量下對流換熱系數(shù)基本相等。

圖11　艙壁溫-25℃時的對流換熱系數(shù)Fig. 11　Convective heat transfer coefficient when the temperature of cell wall is -25℃

圖12　艙壁溫0℃時的對流換熱系數(shù)Fig. 12　Convective heat transfer coefficient when the temperature of cell wall is 0℃

圖13　艙壁溫25℃時的對流換熱系數(shù)Fig. 13　Convective heat transfer coefficient when the temperature of cell wall is 25℃

圖14　艙壁溫40℃時的對流換熱系數(shù)Fig. 14　Convective heat transfer coefficient when the temperature of cell wall is 40℃

2.3實驗分析

2.3.1誤差分析

在各組實驗數(shù)據(jù)中，平板在艙壁溫度為-40℃和-25℃這2個工況的實驗數(shù)據(jù)與其他3個工況的差距較大，主要有3個原因：

1）實驗的控溫系統(tǒng)是通過艙內(nèi)側(cè)壁進(jìn)行溫控，而在艙門部分沒有溫控，使得艙內(nèi)環(huán)境溫度的控制受艙門的影響較大；當(dāng)加熱量小，且艙內(nèi)側(cè)壁與艙門溫差較大時，這種影響就更加顯著。但隨著加熱量的增大，則影響明顯下降。

2）在較低的側(cè)壁溫度環(huán)境下，較小的加熱量使平板表面溫度變化不明顯，設(shè)定的測量精度不能滿足要求，導(dǎo)致數(shù)據(jù)曲線起伏大。

3）-40℃和-25℃的側(cè)壁低溫比較難以準(zhǔn)確控制，實驗過程中要保持壁溫穩(wěn)定調(diào)節(jié)難度很大，因此不穩(wěn)定的溫度可能導(dǎo)致小加熱量工況的數(shù)據(jù)出現(xiàn)較大誤差。

2.3.2輻射換熱量和對流換熱量的比較

如果不考慮誤差的影響，將圖4～圖8中的對流換熱量與輻射換熱量相比（以25℃艙壁溫度為例，見圖15），可以看出：艙內(nèi)氣壓在1kPa以下時，對流換熱量非常小，特別在0.1kPa以下時，對流換熱相對輻射換熱可以忽略不計；當(dāng)氣體壓力大于1kPa時，對流換熱所占的比例按二次方或指數(shù)形式逐步增加；到常壓時，自然對流換熱量和輻射換熱量才基本相當(dāng)。但是，在實際應(yīng)用中還應(yīng)該考慮壁溫、空氣溫度和平板表面溫度之間相互的影響，這之間的差別會導(dǎo)致?lián)Q熱量出現(xiàn)較小的不同。

圖15　艙壁溫25℃時的對流與輻射換熱量之比Fig. 15　Ratio of radiant and convection heat transfer when the temperature of cell wall is 25℃

2.3.3對流換熱系數(shù)

由前面的結(jié)果和分析可以看出：受到實驗艙側(cè)壁溫度較低和加熱量較小的影響，對流換熱系數(shù)在較小換熱時變化較大，導(dǎo)致測量結(jié)果有一定的差異；而在艙壁溫度較高（0、25、40℃）和加熱量較大（3W和6W加熱功率）的情況下，對流換熱系數(shù)受到的影響較小，因而誤差也較小，這與變化趨勢很相似。

綜合上述實驗結(jié)果，可以擬合出一個與加熱量、艙壁溫度無關(guān)的平板對流換熱系數(shù)α和環(huán)境氣壓p間的關(guān)系式

用于計算1kPa壓力以上的氣體對流換熱系數(shù)，有很好的精度。

3　結(jié)論

通過對實驗結(jié)果的分析，對于低加熱量、有限空間內(nèi)水平平板的自然對流換熱得到如下結(jié)論：

1）環(huán)境氣壓小于1kPa時，對流換熱量非常小，只占輻射換熱量的不到10%。

2）環(huán)境氣壓大于1kPa時，平板的對流換熱系數(shù)隨環(huán)境氣壓的升高呈二次方增加。

參考文獻(xiàn)（References）

[1] 姚偉, 李勇, 王文雋, 等. 平流層飛艇熱力學(xué)模型和上升過程仿真分析[J]. 宇航學(xué)報, 2007, 28(3): 603-607 Yao Wei, Li Yong, Wang Wenjun, et al. Thermodynamic model and numerical simulation of a stratospheric airship take-off process[J]. Journal of Astronautics, 2007, 28(3): 603-607

[2] 方賢德, 王偉志, 李小建. 平流層飛艇熱仿真初步探討[J]. 航天返回與遙感, 2007, 28(2): 5-9 Fang Xiande, Wang Weizhi, Li Xiaojian. A study of thermal simulation of stratospheric airships[J]. Spacecraft Recovery & Remote Sensing, 2007, 28(2): 5-9

[3] Devienne F M. Low density heat transfer[M]//Harnett J P, Irvine T F. Advances in heat transfer. New York: Academic Press Inc., 1965: 271-356

[4] Nelson K E, Bevans J T. Errors of the calorimetric method of total emittance measurement[M]//Richmond J C. Measurement of thermal radiation properties of solids. Washington D. C.: NASA, 1963:55-65

[5] Warner C Y, Arpaci V S. An experimental investigation of turbulent natural convection in are at low pressure along a vertical heated flat plate[J]. International Journal of Heat and Mass Transfer, 1968, 11: 397-406

[6] Kyte J R, Madden A J, Piret E L. Natural-convection heat transfer at reduced pressure[J]. Chemical Engineering Progress, 1953, 49(12): 653-662

[7] Churchill S W, Chu H H S. Correlating equations for laminar and turbulent free convection from a vertical plate[J]. Int J Heat Mass Transfer, 1975, 18: 1323-1329

（編輯：肖福根）

Experimental investigation of the thermal performance of horizontal flat’s surface under different pressures

Wang Jing1, Xu Ying2, Ding Li2, Qie Dianfu1

(1. Science and Technology on Reliability and Environmental Engineering Laboratory, Beijing Institute of Spacecraft Environment Engineering, Beijing 100094, China;

2. School of Biological Science and Medical Engineering, Beihang University, Beijing 100191, China)

Abstract:In order to understand the convective heat transfer of the horizontal flat at different pressures, a sealed cabin that can provide different pressures and temperatures is designed, and the thermal performance of the horizontal flat in several fixed heating input(75W/m2, 150W/m2, 300W/m2)under different pressures(0.1Pa, 0.1kPa, 0.2kPa, 0.5kPa,1kPa, 10kPa, 50kPa, ambient pressure)is experimentally investigated. By comparing the radiant heat exchange with the natural convection heat transfer, the convective heat transfer coefficient under different pressures is obtained. The experimental results show that the convective heat transfer coefficient is very small when the absolute environmental pressure is less than 1kPa; it might be meaningful only when the absolute pressure is more than 1kPa; and it increases quadratically along with the increase of the environmental pressure.

Key words:low air pressure; horizontal flat; convective heat transfer coefficient; radiant heat exchange; experimental study

作者簡介：王 晶（1980—），男，博士學(xué)位，高級工程師，主要從事航天器熱試驗技術(shù)研究；E-mail: wangjing19800510@ buaa.edu.cn。通信作者：丁 立（1971—），男，博士學(xué)位，副教授，從事生物與醫(yī)學(xué)工程研究；E-mail: ding1971316@buaa.edu.cn。

基金項目：可靠性與環(huán)境工程技術(shù)重點實驗室開放基金

收稿日期：2015-10-24；修回日期：2016-03-14

DOI:10.3969/j.issn.1673-1379.2016.02.010

中圖分類號：O551.1

文獻(xiàn)標(biāo)志碼：A

文章編號：1673-1379(2016)02-0167-06