劉　暢　苗建印　何　江　王玉瑩　寧獻(xiàn)文　呂　?。ū本┛臻g飛行器總體設(shè)計部，空間熱控技術(shù)北京市重點(diǎn)實(shí)驗(yàn)室，北京 100094）

結(jié)構(gòu)參數(shù)對水升華器散熱性能影響的研究

劉暢苗建印何江王玉瑩寧獻(xiàn)文呂巍
（北京空間飛行器總體設(shè)計部，空間熱控技術(shù)北京市重點(diǎn)實(shí)驗(yàn)室，北京 100094）

應(yīng)用一維穩(wěn)態(tài)導(dǎo)熱模型，結(jié)合具有滑移效應(yīng)的達(dá)西定律、氣體通過小孔的流動模型，針對一套水升華器提出了穩(wěn)態(tài)理論分析方法。分析了接觸換熱系數(shù)、排氣通道面積兩個關(guān)鍵結(jié)構(gòu)參數(shù)對水升華器散熱性能的影響，并開展了相關(guān)實(shí)驗(yàn)研究。仿真分析和實(shí)驗(yàn)結(jié)果表明：穩(wěn)態(tài)散熱量隨流體回路冷板與給水腔之間接觸換熱系數(shù)增大而增大，接觸換熱系數(shù)較小時，散熱量和接觸換熱系數(shù)呈強(qiáng)耦合關(guān)系，在接觸換熱系數(shù)較大時，二者呈弱耦合關(guān)系。分析表明：散熱量隨排氣通道面積增大而增大，且存在臨界排氣孔面積，排氣孔面積小于臨界值時，散熱量隨排氣孔面積的增大而急劇升高，排氣孔面積大于臨界值時，散熱量隨排氣孔面積的增大而幾乎不變。文章研究結(jié)果可為空間水升華器的設(shè)計提供參考。

水升華器；接觸換熱系數(shù)；排氣通道；散熱性能

1　引言

水升華器是一種利用水工質(zhì)的升華原理、以水作為消耗型介質(zhì)的相變散熱裝置，適用于航天器熱控系統(tǒng)不能提供足夠散熱面或具有短時大功耗等情況，特別適合于短期宇航任務(wù)。研究表明，與輻射器等其它熱控器件相比，水升華器體積小、質(zhì)量輕，在微重力、月球重力場、加速度等環(huán)境以及熱負(fù)荷變化的條件下均能夠可靠工作。在美國和俄羅斯的航天活動中，水升華器得到了大量應(yīng)用［1-3］。我國水升華器技術(shù)的研究起步較晚［4-6］，目前僅在“飛天”航天服上得到應(yīng)用。隨著我國探月工程等工作的逐步展開以及空間站艙外活動單元的不斷增多，水升華器技術(shù)已成為我國未來航天器熱控中必須突破的關(guān)鍵技術(shù)之一。特別對于嫦娥五號在月晝期間執(zhí)行月壤鉆取和月表采樣任務(wù)，對可高溫啟動、具有大散熱量的水升華器需求迫切。瞬態(tài)啟動特性與穩(wěn)態(tài)散熱特性是研究水升華器熱控技術(shù)必須關(guān)注的兩個方面。對于瞬態(tài)啟動特性，國內(nèi)外研究人員均進(jìn)行過相關(guān)研究［5-8］，主要方向集中在研制更高性能的多孔板以及整體結(jié)構(gòu)優(yōu)化。而水升華器作為散熱裝置，其穩(wěn)態(tài)散熱特性也是必須關(guān)注的重要問題。影響水升華器穩(wěn)態(tài)散熱能力的結(jié)構(gòu)參數(shù)主要有:傳熱熱阻、排氣通道面積、換熱面積、多孔板參數(shù)等等［9］。

關(guān)于水升華器的穩(wěn)態(tài)散熱特性，目前研制過程存在的主要問題有:在水升華器換熱面積一定、多孔板參數(shù)不變的情況下，如何提高相同工況下的穩(wěn)態(tài)散熱量；如何設(shè)計蒸汽排放通道，在盡量簡化布局、節(jié)約資源的同時，能夠滿足穩(wěn)態(tài)散熱量的需求。針對這些問題，本文以一維穩(wěn)態(tài)導(dǎo)熱分析法為基礎(chǔ)，結(jié)合具有滑移效應(yīng)的達(dá)西定律、氣體通過小孔的流動過程建立了水升華器穩(wěn)態(tài)工作的理論模型，針對接觸換熱系數(shù)、排氣通道面積等結(jié)構(gòu)參數(shù)對水升華器穩(wěn)態(tài)散熱性能的影響進(jìn)行了仿真分析，并開展了部分實(shí)驗(yàn)研究。通過仿真與實(shí)驗(yàn)研究相結(jié)合，得出了這些參數(shù)對水升華器穩(wěn)態(tài)散熱性能影響的相關(guān)規(guī)律。

2　水升華器工作原理

空間高真空環(huán)境壓力遠(yuǎn)小于水的三相點(diǎn)壓力，液態(tài)水在真空條件下會迅速蒸發(fā)吸熱，由于水的相變潛熱大，水溫迅速降低，溫度降至水的三相點(diǎn)以下便會結(jié)冰，冰繼續(xù)通過升華帶走熱量，形成升華制冷過程。水升華器正是利用這一物理現(xiàn)象，實(shí)現(xiàn)散熱功能。

本文研究的水升華器為平板型，由給水腔和多孔板構(gòu)成，并且通過流體回路冷板與單相流體回路耦合，如圖1所示。水升華器啟動后，給水進(jìn)入給水腔、暴露在真空環(huán)境中，液態(tài)水迅速蒸發(fā)吸熱，水溫下降，直至結(jié)冰，冰繼續(xù)升華吸收熱量。在熱負(fù)荷作用下，冰層消失，給水再次進(jìn)入多孔板蒸發(fā)，直至結(jié)冰，開始新的升華過程。通過多孔板控制水蒸汽排放阻力，即可獲得穩(wěn)定的周期性蒸發(fā)/升華相變過程。流體回路通過流體工質(zhì)的循環(huán)流動收集航天器設(shè)備、元器件產(chǎn)生的廢熱，作為水升華器的熱負(fù)荷。水升華器通過水的相變過程吸收流體回路的熱量，從而形成穩(wěn)定的制冷過程，實(shí)現(xiàn)散熱功能。

圖1　水升華器結(jié)構(gòu)示意圖Fig.1　Sketch map of water sublimator structure

水升華器內(nèi)部的工作機(jī)制涉及多孔介質(zhì)內(nèi)液體流動、真空環(huán)境下工質(zhì)蒸發(fā)/升華等相變傳熱傳質(zhì)耦合過程，是極其復(fù)雜的氣-液-固三相流動換熱問題。受外部環(huán)境及內(nèi)部結(jié)構(gòu)、多孔板參數(shù)等多種因素影響，工作模式復(fù)雜多變，一般的工作模式有蒸發(fā)模式、升華模式、周期模式、混合模式等等［10-11］。

3　水升華器穩(wěn)態(tài)工作理論模型

對于文中所研究的水升華器，可忽略其工作過程中自然對流的影響，將內(nèi)部流體的換熱視為導(dǎo)熱［9］。采用一維穩(wěn)態(tài)導(dǎo)熱等效分析模型（能量守恒），結(jié)合具有滑移效應(yīng)的達(dá)西定律、氣體通過小孔的流動模型（質(zhì)量守恒），建立水升華器穩(wěn)態(tài)工作的理論模型。

3.1能量守恒

水升華器工作過程中，流體回路工質(zhì)熱量的損失與通過水升華器的散熱相平衡。不考慮溫度差異造成的結(jié)構(gòu)熱物性變化，且忽略水升華器與外界的輻射漏熱。水升華器系統(tǒng)散熱量可由下式計算:

式中:Cp為流體回路工質(zhì)比熱容；m·為工質(zhì)質(zhì)量流量；ΔT為流體回路出入口溫差；Tin，Tout分別為流體回路入口、出口處溫度；ρ為工質(zhì)密度；V·為工質(zhì)體積流量。

另外，根據(jù)傳熱關(guān)系，水升華器穩(wěn)態(tài)工作時，散熱量可用冷、熱側(cè)之間工質(zhì)的換熱量表示，即

式中:Rtot為工質(zhì)至冰水界面之間傳熱路徑的總熱阻；ΔTm為冷、熱側(cè)的平均溫差，冷側(cè)冰水界面可認(rèn)為是均勻壁溫，因此平均溫差采用對數(shù)平均溫差，即

式中:Tc為多孔板內(nèi)冰水界面處溫度。

通過冰層的導(dǎo)熱量為

式中:ki為冰層與多孔材料等效導(dǎo)熱系數(shù)；；Tf為冰水界面溫度；Ts為升華界面溫度；δice為冰層厚度。

升華界面通過相變吸收的熱量為

式中:Qc為升華界面相變吸收的熱量；M為水蒸汽質(zhì)量流量；ΔHs為水的升華相變潛熱。

根據(jù)能量守恒，有

式中:ΔHf為水的凝固相變潛熱。

3.2質(zhì)量守恒

水升華器工作過程中，給水的質(zhì)量流量與水蒸氣的質(zhì)量流量相等。水蒸汽質(zhì)量流量可從式（5）和式（6）計算得出:

由于水升華器穩(wěn)態(tài)工作時升華表面蒸汽壓低，蒸汽分子平均自由程較大，而多孔板屬于微孔結(jié)構(gòu)，克努曾數(shù)Kn＞1.5，需要考慮因氣體稀薄性而引起的滑動現(xiàn)象［12］。因此，水蒸汽流經(jīng)多孔板的過程采用具有滑移效應(yīng)的達(dá)西定律描述。根據(jù)達(dá)西定律，水蒸汽質(zhì)量流量為

式中:Kr為考慮氣流滑移的表觀滲透率；A為多孔板表面積；ρv為水蒸汽密度；Ps為升華表面蒸汽壓；Pvc為蓋板腔內(nèi)水蒸汽壓強(qiáng)；ηv為水蒸汽黏性系數(shù)；δpv為升華表面至多孔板上表面距離；K為材料滲透率；C0=4.018，為常數(shù)。

對于蓋板腔體內(nèi)水蒸汽通過排氣孔的流動［13］，由于小孔尺寸較大，λ＜＜dh（λ為氣體分子平均自由程，dh為小孔尺寸），慣性力與摩擦力起主要作用。對于排氣孔出口處，氣體流出小孔后立即波浪式擴(kuò)張。在這種流動下，流量通常隨環(huán)境壓力P0的降低而增加，但當(dāng)P0/Pvc=r值小到一定臨界值rc后，流量就不再隨P0的降低而變化。根據(jù)流體力學(xué)，rc值以及氣體質(zhì)量流量為

式中:μ為水的摩爾質(zhì)量；γ為水蒸汽的比熱容比；Tvc為蒸汽溫度；Ah為排氣通道開孔面積。

水升華器穩(wěn)態(tài)工作時傳熱路徑熱阻包括流體工質(zhì)與冷板表面?zhèn)鳠釤嶙?、冷板和給水腔之間接觸熱阻、給水腔整體熱阻、多孔板下表面至冰水界面導(dǎo)熱熱阻。流體回路冷板工質(zhì)與冷板表面等效換熱系數(shù)由研制單位提供；給水腔與多孔板為多孔結(jié)構(gòu)，等效導(dǎo)熱系數(shù)按照多孔材料的計算方法進(jìn)行評估［14］，即

式中:εp為材料的孔隙度；λs和λf分別為固體材料和流體材料的導(dǎo)熱系數(shù)。實(shí)際工作過程中冰水界面位置會隨不同工況而變化，熱阻也會發(fā)生相應(yīng)改變。

根據(jù)以上各式，以流體回路工質(zhì)流量以及冷板入口溫度作為輸入，首先假設(shè)出口溫度的值，由式（1）及式（7）計算出散熱量、水蒸汽質(zhì)量流量，進(jìn)一步由式（2）、式（11）可得到冰水界面位置和蓋板內(nèi)蒸汽壓，根據(jù)式（4）以及式（8）可以求得升華界面位置以及相應(yīng)的表面蒸汽壓。由以下限制條件:

可得流體回路冷板出口溫度的可行域。

4　理論分析與實(shí)驗(yàn)研究

4.1模型準(zhǔn)確性驗(yàn)證

4.1.1實(shí)驗(yàn)系統(tǒng)

本文實(shí)驗(yàn)系統(tǒng)如圖2（a）所示，由4部分組成:①真空系統(tǒng)，包括真空室及真空機(jī)組。②測量系統(tǒng)，對實(shí)驗(yàn)中溫度、壓力、流量等數(shù)據(jù)進(jìn)行實(shí)時測量采集。溫度采用銅/康銅熱電偶進(jìn)行測量，熱電偶經(jīng)過標(biāo)定，測溫誤差在±0.2℃內(nèi)。真空壓力測量采用ZJ-1P型壓阻真空規(guī)管，測量精度為±1.0%。流量測量采用DMF-1-1-A型質(zhì)量流量計，測量精度為±0.5%。③流體回路系統(tǒng)，作為熱負(fù)荷，具有工質(zhì)溫度、流量連續(xù)可調(diào)功能。④水升華器實(shí)驗(yàn)件（見圖2（b））及其供水系統(tǒng)。實(shí)驗(yàn)件與流體回路通過板式換熱器（以下簡稱冷板）進(jìn)行熱耦合，實(shí)驗(yàn)件的給水腔與冷板之間涂抹導(dǎo)熱硅脂。

圖2　水升華器實(shí)驗(yàn)系統(tǒng)Fig.2　Experimental system of water sublimator

4.1.2模型仿真結(jié)果與實(shí)驗(yàn)結(jié)果對比

對水蒸汽排放通道面積分別為980mm2和1960mm2兩種狀態(tài)的水升華器，流體回路工質(zhì)流量保持在（150±2）L·h—1，在不同冷板入口溫度情況下進(jìn)行散熱性能實(shí)驗(yàn)。本文所有實(shí)驗(yàn)工況，散熱量的測量相對最大誤差均在5%以內(nèi)。根據(jù)第3章所建立的分析理論模型，對不同實(shí)驗(yàn)工況進(jìn)行仿真。計算過程中，水升華器傳熱路徑各部分熱阻取值見表1（冷板熱阻為流體回路流量150 L·h—1情況下的數(shù)據(jù)）。

實(shí)驗(yàn)結(jié)果與仿真結(jié)果對比情況見圖3。結(jié)果表明，在不同工況下，仿真結(jié)果與實(shí)驗(yàn)測量值的相對誤差均不超過4.5%，均在實(shí)驗(yàn)測量誤差范圍內(nèi)。因此，可以根據(jù)本文模型對此種狀態(tài)水升華器在不同輸入條件（流體回路工質(zhì)流量及冷板入口溫度）下的散熱量進(jìn)行預(yù)測。

表1　水升華器傳熱路徑各部分熱阻Table1　Thermal resistance of water sublimator heat transfer path

圖3　不同工況下實(shí)驗(yàn)結(jié)果與模型仿真結(jié)果對比Fig.3　Comparison results in different working condition between experiment and simulation

4.2接觸換熱系數(shù)對散熱量影響分析與實(shí)驗(yàn)研究

實(shí)際上，對于不同的水升華器，冷板和給水腔之間可以通過干接觸、涂抹導(dǎo)熱脂、焊接等不同方式進(jìn)行連接，即使連接方式相同，接觸換熱系數(shù)也會因工藝尺寸、表面狀態(tài)等存在差別而具有離散性。另外，對于4.1節(jié)的水升華器實(shí)驗(yàn)件狀態(tài)，冷板和給水腔之間的接觸熱阻占水升華器傳熱路徑中總熱阻73.46%（見表1），所占比重最大。因此，須要研究冷板和給水腔之間接觸換熱系數(shù)對水升華器穩(wěn)態(tài)散熱量的影響。

排氣通道面積為1960mm2，水升華器其他工藝參數(shù)不變，僅改變冷板和給水腔之間接觸換熱系數(shù)，利用第3章所建立的分析模型進(jìn)行仿真。仿真所得流體回路工質(zhì)體積流量150 L·h—1、冷板入口溫度15℃時的穩(wěn)態(tài)散熱量結(jié)果如圖4所示。

仿真結(jié)果表明，冷板和給水腔之間接觸換熱系數(shù)對穩(wěn)態(tài)散熱量的影響是顯著的，穩(wěn)態(tài)散熱量隨接觸換熱系數(shù)增大而增大，且分為強(qiáng)耦合作用與弱耦合作用區(qū)域。在強(qiáng)耦合作用區(qū)域，接觸換熱系數(shù)增大，散熱量提升效果比較明顯；在弱耦合作用區(qū)域，接觸換熱系數(shù)增大，散熱量的提升效果并不明顯。對于本文研究的水升華器，兩個區(qū)域分界處接觸換熱系數(shù)為1328 W·m—2·K—1。接觸換熱系數(shù)增加時，接觸熱阻占傳熱路徑總熱阻比例逐漸降低，其影響也逐漸變小。接觸換熱系數(shù)為1328 W·m—2· K—1時，接觸熱阻占傳熱路徑總熱阻的一半。

圖4　冷板和給水腔之間接觸換熱系數(shù)與穩(wěn)態(tài)散熱量關(guān)系仿真結(jié)果Fig.4　Simulative result of steady-state heat dissipation vs heat transfer coefficient

其他條件相同，針對強(qiáng)耦合作用區(qū)域的兩種接觸換熱系數(shù)情況下的水升華器進(jìn)行實(shí)驗(yàn)，實(shí)驗(yàn)結(jié)果見表2（流體回路工質(zhì)體積流量保持在150±2 L· h—1）。結(jié)果表明，在強(qiáng)耦合作用區(qū)域，接觸換熱系數(shù)從250 W·m—2·K—1增大至500 W·m—2·K—1，水升華器穩(wěn)態(tài)散熱量增加35.07%。與理論分析相符。弱耦合區(qū)接觸熱阻的影響尚缺乏實(shí)驗(yàn)數(shù)據(jù)，后續(xù)會繼續(xù)開展實(shí)驗(yàn)。

表2　兩種接觸換熱系數(shù)情況下實(shí)驗(yàn)結(jié)果Table2　Experimental results in different heat transfer coefficients

根據(jù)以上分析，在強(qiáng)耦合作用區(qū)域，通過控制接觸表面平面度、控制硅脂涂抹工藝等方式增大接觸換熱系數(shù)提升穩(wěn)態(tài)散熱量的效果比較明顯；在弱耦合作用區(qū)域，為提高散熱量，通過更加嚴(yán)格地控制接觸表面平面度、改進(jìn)硅脂涂抹工藝等方式進(jìn)一步增大接觸換熱系數(shù)的做法意義不大，應(yīng)該考慮改變結(jié)構(gòu)設(shè)計來大幅降低熱阻，比如流體回路冷板與給水腔直接通過焊接連接或者采用3D打印技術(shù)進(jìn)行一體化加工。

4.3排氣通道面積對散熱量影響分析與實(shí)驗(yàn)研究

由于總體布局、使用環(huán)境等限制，水升華器的蒸汽排放通道狀態(tài)會受到影響。因此有必要進(jìn)行排氣通道面積對水升華器散熱性能影響的研究。對于4.1節(jié)的水升華器實(shí)驗(yàn)件狀態(tài)，其他工藝參數(shù)不變，僅改變排氣通道面積，利用第3章所建立的分析模型進(jìn)行仿真。通過仿真所得流體回路工質(zhì)體積流量150 L·h—1、冷板入口溫度15℃時的穩(wěn)態(tài)散熱量結(jié)果如圖5所示。

圖5　排氣孔面積與穩(wěn)態(tài)散熱量關(guān)系的仿真結(jié)果Fig.5　Simulative result of steady-state heat dissipation vs exhaust port area

仿真結(jié)果表明，排氣孔面積對水升華器的穩(wěn)態(tài)散熱量的影響十分顯著，且存在臨界值（對于本文研究的水升華器，排氣孔面積臨界值約為78.4mm2）。排氣孔面積小于臨界值時，水升華器的穩(wěn)態(tài)散熱量隨排氣孔面積的增大急劇升高；排氣口面積大于臨界值時，水升華器的穩(wěn)態(tài)散熱量隨排氣孔面積的增大繼續(xù)增大，但增幅極小，幾乎不增加。

針對水升華器實(shí)驗(yàn)件排氣通道面積分別為980mm2和1960mm2的兩種狀態(tài)，其他結(jié)構(gòu)參數(shù)不變，利用4.1節(jié)所述實(shí)驗(yàn)系統(tǒng)開展實(shí)驗(yàn)。流體回路工質(zhì)流量保持在（150±2）L·h—1，不同冷板入口溫度條件下，水升華器穩(wěn)態(tài)散熱量實(shí)驗(yàn)結(jié)果見圖3。結(jié)果表明，排氣孔面積大于臨界值時，增大排氣通道面積，在不同冷板入口溫度條件下，散熱量存在一定程度升高，但變化幅度很小，幾乎沒有變化。與仿真結(jié)果相符。排氣通道面積存在臨界值的結(jié)論尚缺乏實(shí)驗(yàn)數(shù)據(jù)支持，后續(xù)將繼續(xù)開展相關(guān)實(shí)驗(yàn)。

根據(jù)上述分析，進(jìn)行水升華器結(jié)構(gòu)設(shè)計時，應(yīng)當(dāng)注意排氣孔面積臨界值的存在。排氣孔面積大于臨界值時，可以通過適當(dāng)減少排氣通道數(shù)目對總體布局進(jìn)行簡化。

5　結(jié)束語

本文以一維穩(wěn)態(tài)導(dǎo)熱分析法為基礎(chǔ)，結(jié)合具有滑移效應(yīng)的達(dá)西定律和氣體通過小孔的流動過程，建立了水升華器穩(wěn)態(tài)工作的理論模型。針對接觸熱阻、排氣通道面積等因素對水升華器穩(wěn)態(tài)散熱性能影響進(jìn)行分析，并開展了部分實(shí)驗(yàn)研究。分析表明:

（1）實(shí)驗(yàn)結(jié)果與模型仿真結(jié)果相符，可以根據(jù)本文模型對此種狀態(tài)水升華器在不同輸入條件（流體回路工質(zhì)流量及冷板入口溫度）下的穩(wěn)態(tài)散熱量進(jìn)行預(yù)測。

（2）冷板和給水腔之間接觸熱阻對穩(wěn)態(tài)換熱量的影響比較顯著，穩(wěn)態(tài)散熱量隨接觸換熱系數(shù)增大而增大，且分為強(qiáng)耦合作用與弱耦合作用區(qū)域。在強(qiáng)耦合作用區(qū)域，接觸換熱系數(shù)增大，散熱量提升效果比較明顯；在弱耦合作用區(qū)域，接觸換熱系數(shù)增大，散熱量的提升效果并不明顯。

（3）排氣孔面積對水升華器的穩(wěn)態(tài)散熱量的影響十分顯著，且存在臨界值。排氣孔面積小于臨界值時，穩(wěn)態(tài)散熱量隨排氣孔面積的增大急劇升高。排氣孔面積大于臨界值時，穩(wěn)態(tài)散熱量隨排氣孔面積的增大繼續(xù)增大，但增幅極小，幾乎不增加。因此，排氣孔面積大于臨界值時，可以適當(dāng)減少排氣孔面積或減少排氣通道數(shù)目以簡化布局。

由于條件有限，本文實(shí)驗(yàn)工況較少，僅開展了強(qiáng)耦合區(qū)接觸熱阻的影響實(shí)驗(yàn)，且排氣通道面積存在臨界值的結(jié)論缺乏實(shí)驗(yàn)數(shù)據(jù)支持，后續(xù)將繼續(xù)開展相關(guān)實(shí)驗(yàn)。另外，由于水升華器內(nèi)部工作機(jī)制十分復(fù)雜，本文所建模型進(jìn)行了簡化，是一種等效的分析方法，后續(xù)還須進(jìn)行深入的理論分析，并通過引入可視化等手段開展進(jìn)一步研究。

（References）

［1］J L Wilson，BM Lawson.Investigation into venting and non-venting technologies for the space station freedom extravehicular activity life support system，SAE 901319 ［R］.Houston，Texas:NASA，1990

［2］H Standard.Phase 1 engineering and technical data report for the thermal control extravehicular life support system，N75-24360［R］.Washington D.C.:NASA，1975

［3］A I Skoog，I P Abramov，A Y Stoklitsky，et al.The Soviet-Russian space suits a historical overview of the 1960's［J］.Acta Astronautic，2002，51（1-9）:113-131

［4］吳志強(qiáng)，袁修干，沈力平.水升華器散熱系統(tǒng)分析［J］.航空學(xué)報，1999，20（z）:16-18 Wu Zhiqiang，Yuan Xiugan，Shen Liping.Analysis of water sublimator heat rejection systems［J］.Acta Aeronautica et Astronautica Sinica，1999，20（z）:16-18（in Chinese）

［5］吳志強(qiáng)，沈力平，韓力軍，等.憎水涂層多孔板對水升華器散熱性能影響的實(shí)驗(yàn)研究［J］.航天醫(yī)學(xué)與醫(yī)學(xué)工程，2003，16（4）:287-291 Wu Zhiqiang，Shen Liping，Han Lijun，et al.Experimental research on the effects of hydrophobic coating porous plates on heat dissipation of water sublimator［J］.Space Medicine&Medical Engineering，2003，16（4）:287-291 （in Chinese）

［6］王玉瑩，鐘奇，寧獻(xiàn)文，等.水升華器空間應(yīng)用研究［J］.航天器工程，2013，22（3）:105-112 Wang Yuying，Zhong Qi，Ning Xianwen，et al.Overview of space application and development of water sublimator［J］. Spacecraft Engineering，2013，22（3）:105-112（in Chinese）

［7］F L Worley，WI Honeywell，A H Blair.Causes and explanation of‘breakthrough phenomenon'when LEM cooling system sublimator is fed with chlorinated feedwater，NAS-CB-45735［R］.Houston，Texas:NASA Manned Spacecraft Center，1967

［8］Wang Yuying，Zhong Qi，Ning Xianwen，et al.Transient study about the heat transfer of sublimator combined with fluid loop，IAC-13-C2.7.12［C］//64thInternational Astronautical Congress.Beijing:International Astronautical Federation，2013

［9］王玉瑩.空間水升華器相變傳熱傳質(zhì)動態(tài)特性及穩(wěn)定性研究［D］.北京:中國空間技術(shù)研究院，2014 Wang Yuying，Transient work performance and stability study on the phase change heat and mass transfer of space water sublimator［D］.Beijing，China Academy of Space Technology，2014（in Chinese）

［10］A J Chapman.A fundamental study of sublimation through a porous surface，NASA-CR-115153［R］.Houston，Texas:NASA Manned Spacecraft Center，1971

［11］D W Graumann，G R Woods.Research study on instrument unit thermal conditioning heat sink concepts，NAS-8-11291［R］.Huntsville，Alabama: George C.Marshall Space Flight Center，1968

［12］吳偉，雷樹業(yè)，杜建華，等.氣體通過顆粒層突破壓力和達(dá)西滲透率的關(guān)系［J］.化工學(xué)報200354（6）:738-742 Wu Wei，Lei Shuye，Du Jianhua，et al.Relationship of threshold diameter and Darcean permeability for gas flow through unconsolidated porous medium［J］. Journal of Chemical Industry and Engineering，2003，54（6）:738-742（in Chinese）

［13］王欲知，陳旭.真空技術(shù)［M］.北京:北京航空航天大學(xué)出版社，2007 Wang Yuzhi，Chen Xu.Vacuum technology［M］.Beijing:Beihang University Press，2007（in Chinese）

［14］劉偉，范愛武，黃曉明.多孔介質(zhì)傳熱傳質(zhì)理論與應(yīng)用［M］.北京:科學(xué)出版社，2006 Liu Wei，F(xiàn)an Aiwu，Huang Xiaoming.Theory and application of heat and mass transfer in porous medium ［M］.Beijing:Science Press，2006（in Chinese）

（編輯：張小琳）

Research on Effects of Structure Parameters on Heat Dissipation of Water Sublimator

LIU Chang MIAO Jianyin HE Jiang WANG Yuying NING Xianwen LYU Wei
（Beijing Key Laboratory of Space Thermal Control Technology，Beijing Institute of Spacecraft System Engineering，Beijing 100094，China）

Based on the model of one-dimension steady-state thermal conductance，Darcy's Law with slip effect and the model of gas flowing through a hole，a steady-state theoretical analysis method of the water sublimator is proposed.The effects of contact heat transfer coefficient and area of exhaust channel on heat dissipation are analyzed，and relevant experimental study is carried out.Both the numerical and experimental results reveal that the heat dissipation increases with increasing contact heat transfer coefficient between the coldplate of fluid loop and the feedwater chamber.Moreover，the coupling between heat dissipation and contact heat transfer coefficient is strong when the contact heat transfer coefficient is small，while the coupling is weak when the contact heat transfer coefficient is big；there exists a critical area of exhaust channel.With increase of the area of exhaust channel，the heat dissipation increases rapidly when the area of exhaust channel is smaller than the critical area，otherwise，the heat dissipation hardly changes. The results of this paper can offer some reference for the design of space water sublimator.

water sublimator；contact heat transfer coefficient；exhaust channel；heat dissipation

V476.2

A

10.3969/j.issn.1673-8748.2016.03.009

2016-02-25；

2016-04-12

國家自然科學(xué)基金（11472040）

劉暢，男，碩士，研究方向?yàn)楹教炱鳠峥刂?、水升華器熱控技術(shù)。Email：liuchang0918201@163.com。