王宇辰,杜 鵬,解 崢
(北京衛(wèi)星環(huán)境工程研究所,北京 100094)
火星是人類目前為止探測次數(shù)最多的行星。從1960年開始至今,人類已經(jīng)實(shí)施了40余次火星探測任務(wù),其中到達(dá)火星表面的探測器包括“海盜一號”“海盜二號”“勇氣號”“機(jī)遇號”“好奇號”等[1]?;鹦潜砻娲髿鈮阂话銥?00~1000 Pa,主要成分為CO2,表面溫度-110~20 ℃,同時(shí)表面風(fēng)速為0~15 m/s[2],因此火星環(huán)境相比常規(guī)的航天器軌道環(huán)境更為復(fù)雜。歷史上有多個(gè)火星探測器由于對環(huán)境認(rèn)識不充分而在到達(dá)火星表面后失效或局部失效[3]。目前各宇航機(jī)構(gòu)更傾向于在地面對火星探測器開展包括低氣壓、氣體成分、溫度、風(fēng)速的綜合環(huán)境試驗(yàn),以達(dá)到修正航天器熱模型、驗(yàn)證系統(tǒng)在極端環(huán)境下工作性能等目的。
為滿足我國目前正在研制的火星車的試驗(yàn)需求,本文通過分析火星表面復(fù)雜熱環(huán)境和NASA類似試驗(yàn),結(jié)合氣氮調(diào)溫系統(tǒng)進(jìn)行計(jì)算,實(shí)現(xiàn)火星表面瞬態(tài)熱環(huán)境模擬。
火星自轉(zhuǎn)1周用時(shí)約為24 h 37 min,與地球類似,晝夜變化使火星表面大氣溫度相應(yīng)改變。由于火星大氣層很薄,難以通過大氣運(yùn)動傳遞表面的熱量,所以在1個(gè)火星晝夜內(nèi)其表面溫度變化較大?!昂1I號”著陸器對其著陸點(diǎn)的大氣溫度進(jìn)行了測量(見圖1[4]):2個(gè)著陸點(diǎn)的夏季平均溫度為-60 ℃,晝夜溫度變化約為50 ℃;冬季最低溫度達(dá)-120 ℃,日溫度變化達(dá)80 ℃?;鹦莾蓸O的整個(gè)冬季溫度低于-123 ℃,使得火星大氣環(huán)境中的主要成分CO2凍結(jié)成白色沉積物,形成極冠;極冠的季節(jié)性循環(huán),使火星表面的總氣壓波動達(dá)30%[5]。
圖1 “海盜號”著陸點(diǎn)大氣溫度日變化Fig. 1 Diurnal variation of atmospheric temperature at Viking’s landing site
火星巡視器發(fā)射前需要在地面經(jīng)歷充分的驗(yàn)證,這就要求在地面模擬火星表面的大氣環(huán)境,其中火星表面隨著晝夜快速變化的溫度環(huán)境是極為重要的模擬要素。如“火星探路者”計(jì)劃中的“旅居者號”探測器,在地面設(shè)計(jì)了專用試驗(yàn)裝置來模擬晝夜火星環(huán)境,包括火星表面溫度、星表氣體溫度,圖2[6]所示為火星漫游車展開時(shí),在晝夜瞬態(tài)溫度工況下的測試熱響應(yīng)。
圖2 火星環(huán)境溫度模擬曲線Fig. 2 The simulation of Mars ambient temperature
火星表面熱環(huán)境需要模擬較大溫度范圍的低壓環(huán)境,這無法通過常規(guī)的環(huán)模設(shè)備實(shí)現(xiàn),氣氮調(diào)溫系統(tǒng)采用液氮、氣氮調(diào)溫結(jié)合電加熱器的方式實(shí)現(xiàn)溫度調(diào)節(jié),可以滿足試驗(yàn)需求[7]。通過氣氮調(diào)溫系統(tǒng)來實(shí)現(xiàn)火星表面熱環(huán)境模擬的原理如圖3所示,系統(tǒng)的主要部件包括氮?dú)鈮嚎s機(jī)、換熱器、氣體調(diào)節(jié)器、加熱電爐和熱沉等。
壓縮機(jī)為氣氮調(diào)溫系統(tǒng)提供氣體循環(huán)的動力,高溫高壓氮?dú)鈴膲嚎s機(jī)出口進(jìn)入換熱器,與來自熱沉出口的低溫低壓氮?dú)膺M(jìn)行熱交換。換熱器一方面可使回到壓縮機(jī)的低壓氣體溫度得到回升,防止低溫氮?dú)膺M(jìn)入壓縮機(jī)造成設(shè)備損壞;另一方面可回收大量的冷量,以減少系統(tǒng)需要補(bǔ)充的液氮用量。從換熱器出來的高壓低溫氣體進(jìn)入氣體調(diào)節(jié)器,與補(bǔ)充進(jìn)來的少量液氮混合,形成低于目標(biāo)溫度的低溫氣體,再經(jīng)過各路電爐的調(diào)溫進(jìn)入熱沉,可為空間環(huán)境模擬器真空容器提供150~350 K范圍內(nèi)可調(diào)節(jié)的環(huán)境溫度。
圖3 氣氮調(diào)溫系統(tǒng)Fig. 3 The gas nitrogen thermostat system
火星環(huán)境模擬熱試驗(yàn)時(shí),空間環(huán)境模擬器容器內(nèi)部通過壓力控制系統(tǒng)維持在火星極限壓力150 Pa和1400 Pa的低壓環(huán)境,為了防止CO2分子在低溫工況時(shí)凝結(jié)在熱沉上,在試驗(yàn)中用氮?dú)獯鍯O2氣體。氮?dú)鉃榉肿咏Y(jié)構(gòu)對稱的雙原子氣體,無發(fā)射和吸收輻射的能力,為熱輻射的透明體,因此其與熱沉之間的換熱以對流換熱為主。采用熱電偶測量氮?dú)猸h(huán)境溫度,為了減少與熱沉之間的輻射換熱對熱電偶測量精度的影響,在熱電偶熱結(jié)點(diǎn)上貼敷鍍鋁膜,以降低其發(fā)射率,并在熱電偶外圍罩裝聚四氟乙烯圓柱筒,如圖4所示。在試件周圍距離每個(gè)艙板中心0.2 m處各布置1個(gè)這樣的氣體溫度測量傳感器,以準(zhǔn)確測量容器內(nèi)的環(huán)境溫度。
圖4 熱電偶測量氣體溫度Fig. 4 Temperature measurement of gas by thermocouple
氣氮調(diào)溫系統(tǒng)中,氣體調(diào)節(jié)器通過液氮補(bǔ)充閥的開度來實(shí)現(xiàn)其出口的溫度控制,是系統(tǒng)內(nèi)唯一的降溫裝置;電爐通過調(diào)節(jié)加熱功率控制熱沉進(jìn)口處的溫度,可以控制進(jìn)入熱沉的氣體溫度,事實(shí)上形成了如圖5所示通過氣體調(diào)節(jié)器和電爐進(jìn)行二級溫度控制的能量傳遞鏈。
圖5 氣氮調(diào)溫系統(tǒng)能量傳遞鏈Fig. 5 The energy transfer chain of the gaseous nitrogen thermostat system
火星表面溫度隨時(shí)間變化,熱環(huán)境溫度模擬時(shí)期望的目標(biāo)溫度曲線如圖6所示,其中:紅色曲線是150 Pa環(huán)境下的高溫曲線,在日出前到達(dá)最低溫度-92 ℃,隨即快速上升到最高溫度5 ℃,然后溫度下行,日落進(jìn)入黑夜后降溫速度變緩;黑色曲線為1400 Pa環(huán)境下的低溫曲線,與150 Pa環(huán)境下的溫度曲線趨勢一致,按照最低溫度-105 ℃—快速上升—高溫-30 ℃緩和—快速下行—慢速下行的幾個(gè)階段變化。
圖6 火星表面熱環(huán)境模擬目標(biāo)溫度曲線Fig. 6 The target temperature curves of Mars ambient temperature simulation
由于容器內(nèi)氣體和壁面存在對流換熱,試驗(yàn)區(qū)氣體溫度和熱沉之間存在一定的溫度差;同時(shí)為了配合火星車及其他試驗(yàn)設(shè)備,須對試驗(yàn)過程中的特定時(shí)刻溫度進(jìn)行實(shí)時(shí)控制,因此需要消除由于氣體存在熱容而導(dǎo)致的控制時(shí)間差。
在模擬火星晝夜環(huán)境時(shí),首先需要對目標(biāo)氣體溫度所對應(yīng)的熱沉溫度進(jìn)行計(jì)算。試驗(yàn)過程中,系統(tǒng)環(huán)境為150 Pa和1400 Pa的低氣壓,氣體為連續(xù)介質(zhì)流動[8],熱沉與容器壁面之間通過容器內(nèi)的稀薄氣體自然對流和熱輻射進(jìn)行換熱,在熱沉外側(cè)裝有不銹鋼板制成的防輻射屏,熱沉向壁面的輻射傳熱約為
式中:Ts為熱沉溫度,試驗(yàn)過程中最低為-120 ℃;Tb為容器壁邊界溫度,試驗(yàn)過程中基本為常溫20 ℃;Tf為防輻射屏溫度。計(jì)算得出輻射漏熱最大為7.42 W/m2,相比于對流傳熱,可以忽略。
試驗(yàn)區(qū)氣體與熱沉的換熱為
式中Ta為氣體溫度,試驗(yàn)中通過熱電偶進(jìn)行測量。
對于氣體溫度測點(diǎn),有:
式中CV為試驗(yàn)區(qū)氣相空間的平均比熱容。
對流換熱系數(shù)可以使用對流換熱準(zhǔn)則式進(jìn)行計(jì)算,但常規(guī)的準(zhǔn)則式集中針對常壓和高壓流體,同時(shí)氣體溫度測點(diǎn)的等效熱容也難以確定,故只能通過實(shí)際測得的試驗(yàn)數(shù)據(jù)進(jìn)行計(jì)算。
整理式(2)、式(3)可以得到:
試驗(yàn)過程中,可以測得空氣測點(diǎn)和各路熱沉溫度點(diǎn)的溫度。通過前期的調(diào)試,獲得至少1個(gè)周期的試驗(yàn)數(shù)據(jù),取6個(gè)氣體溫度測量點(diǎn)的測溫平均值和所有熱沉溫度的平均值代入式(4)進(jìn)行擬合計(jì)算,得到不同壓力下各項(xiàng)系數(shù)的數(shù)值,如表1所示。
表1 不同傳熱環(huán)境下的傳熱擬合系數(shù)Table 1 Fitting coefficients in different heat transfer environments
從表1中各項(xiàng)數(shù)據(jù)可以看出,不同真空度下輻射傳熱系數(shù)基本一致,主要差異來自于因?yàn)檎婵斩茸兓瘞淼臍怏w自然對流傳熱變化,1400 Pa下自然對流換熱強(qiáng)度大概是150 Pa下的3倍。將換熱系數(shù)及瞬態(tài)工況目標(biāo)曲線離散化后代回式(4),即為關(guān)于每一時(shí)刻的熱沉溫度的4次方程,計(jì)算出不同目標(biāo)氣體溫度對應(yīng)的理想熱沉溫度,可形成定時(shí)溫度控制表單,分階段對氣體調(diào)節(jié)器出口和熱沉進(jìn)口進(jìn)行溫度設(shè)置,并進(jìn)行試驗(yàn)測試,最終試驗(yàn)效果如圖7所示。
圖7 空間環(huán)境模擬器內(nèi)晝夜氣體溫度模擬效果Fig. 7 Diurnal gas temperature simulation by the environmental simulator
由圖7可見,實(shí)測氣體溫度與目標(biāo)溫度曲線吻合較好,進(jìn)一步分析二者的差值發(fā)現(xiàn),如圖8所示,大部分情況下實(shí)測溫度與目標(biāo)溫度間的差值小于5 ℃,最大溫差小于10 ℃,不同周期內(nèi)的溫度偏差較大的時(shí)刻在快速升降溫階段。對比理想熱沉溫度發(fā)現(xiàn),在快速升溫階段,熱沉實(shí)際溫度要比理想溫度有一定的滯后,因?yàn)橛墒?4)計(jì)算出的升降溫轉(zhuǎn)變時(shí)刻理想熱沉溫度會有一個(gè)突變,即離散化求解式(4)時(shí)左側(cè)正負(fù)變化使得最終解構(gòu)成的曲線不連續(xù),需要在溫度快速變化區(qū)進(jìn)一步細(xì)化設(shè)置來消除。
圖8 實(shí)測氣體溫度與目標(biāo)曲線溫差Fig. 8 Temperature difference of gas between measured and target value
為了進(jìn)一步驗(yàn)證曲線的重復(fù)性,對3個(gè)周期內(nèi)對應(yīng)時(shí)刻的溫度進(jìn)行比較,求其標(biāo)準(zhǔn)差
結(jié)果如圖9所示,低溫曲線不同周期內(nèi)對應(yīng)時(shí)刻的標(biāo)準(zhǔn)差小于1 ℃,高溫曲線不同周期內(nèi)對應(yīng)時(shí)刻的標(biāo)準(zhǔn)差小于0.2 ℃,說明溫度曲線具有良好的重復(fù)性,對試驗(yàn)溫度的波動不敏感,具備充分驗(yàn)證溫度環(huán)境的條件。
圖9 不同周期對應(yīng)時(shí)刻溫度標(biāo)準(zhǔn)差Fig. 9 Standard deviation at the corresponding time in different cycles
為了在地面模擬環(huán)境中達(dá)到修正火星探測器熱模型,驗(yàn)證探測器在火星極端晝夜溫度變化環(huán)境下的工作能力等目的,需要對容器內(nèi)低壓環(huán)境下的氣體溫度變化進(jìn)行實(shí)時(shí)監(jiān)控。本文介紹了氣氮調(diào)溫系統(tǒng)的原理,以及實(shí)現(xiàn)火星表面熱環(huán)境溫度模擬的方法。試驗(yàn)結(jié)果表明,該系統(tǒng)可實(shí)現(xiàn)在150 Pa和1400 Pa低氣壓下火星晝夜溫度環(huán)境隨時(shí)間變化的瞬態(tài)工況模擬,試驗(yàn)過程中氣體溫度與目標(biāo)溫度最大相差不超過10 ℃,而且多個(gè)周期的重復(fù)性良好。
隨著火星探測器相關(guān)試驗(yàn)的開展,火星表面晝夜溫度模擬技術(shù)也將得到進(jìn)一步的應(yīng)用,為我國火星探測技術(shù)提供服務(wù)。