微重力條件下熱管吸熱器瞬態熱分析

桂小紅 宋香娥

(1 中國科學院工程熱物理研究所，北京 100190)(2 北京第二外國語學院，北京 100024)

?

微重力條件下熱管吸熱器瞬態熱分析

桂小紅1宋香娥2

(1 中國科學院工程熱物理研究所，北京 100190)(2 北京第二外國語學院，北京 100024)

基于微重力條件下的導熱控制微分方程，采用焓法對熱管吸熱器相變材料容器進行了二維數值建模與仿真，在同時考慮空穴和相變的情況下，對微重力條件下蓄熱單元相變傳熱進行了模擬計算，分析了空穴率對蓄熱容器內部的溫度場和熱性能的影響，并將計算結果同美國航空航天局(NASA)方案熱管吸熱器蓄熱單元相變傳熱計算結果進行了比較，驗證了文中微重力條件下計算模型的合理性與準確性。研究結果表明：空穴影響著蓄熱單元相變的進程，空穴的存在增加了容器內部的溫度梯度，使得容器的蓄熱能力降低；由于熱管徑向溫差較小，熱管壁溫在相變材料熔點附近變化較小，從而在一定程度上能緩解熱斑和熱松脫現象。

空間太陽能發電系統；熱管吸熱器；相變材料；空穴；相變；航天器

1 引言

空間站在軌道運行時，必然經過太陽陰影期,要保證發電系統持續運轉、連續供電，必須采用蓄熱裝置[1-2]。在太陽能熱動力發電系統中，將吸熱、蓄熱與工質換熱三項功能集合在一起的部件稱為吸熱器。當航天器處于日照區時，反射器收集太陽能并將太陽能聚集到位于拋物面焦點的吸熱器腔口。太陽能進入腔體后，其中一部分能量用來加熱循環工質，另一部分被蓄熱介質吸收儲存起來。當航天器進入陰影區后，蓄熱介質釋放出能量給循環工質，維持其出口溫度在循環所要求的最低溫度值以上，使動力系統在陰影區與日照區一樣運行發電。吸熱器的質量約占發電系統總系統質量的1/3，減少吸熱器的質量對于降低發電系統成本有重要意義[3-4]。

盡管從20世紀80年代末國內外就開始對吸熱器的初步研究，研究的內容主要局限于基本型吸熱器[5-6]，而對于熱管吸熱器等先進吸熱器只有少數概念性的介紹，缺乏深入的實質性研究。文獻[7]中的NASA方案對熱管吸熱器蓄熱單元進行了預設計和初步熱分析，粗略介紹了計算物理模型，其方案中采用一維徑向節點模型，沒有給出具體計算方法和程序，只是分析了單個蓄熱容器徑向節點溫度的變化，沒有分析其內部溫度場。文獻[8]中的NASA方案對熱管吸熱器結構參數進行了優選分析，給出了部分計算結果，所采用計算模型仍為一維徑向節點物理模型。文獻[9]對熱管式吸熱器單元熱管進行了數值分析，其計算模型參照了NASA初步方案，也為一維徑向節點模型，采用焓方法對其進行數值求解，得出了軌道周期內單個相變材料(Phase Change Material，PCM)容器外壁最高溫度、熱管壁溫等主要參數的變化，但沒有給出蓄熱容器內部溫度場，并且沒考慮空穴對蓄熱容器溫度場梯度所帶來的影響。

本文結合美國Allied-signal公司Garrett Airesearch提出的用于空間太陽能熱動力發電系統的布雷頓式熱管吸熱器設計案[7]，以及其他國內外吸熱器的研究進展及成果，對微重力條件下熱管吸熱器蓄熱單元內伴隨有空穴的相變換熱進行建模與仿真，對PCM容器進行了詳盡的熱分析，并進行了具體和深入的研究。

2 單元熱管吸熱器結構設計

圖1 熱管吸熱器 Fig.1 Heat pipe receiver

Garrett公司采用布雷頓循環(CBC)的熱管吸熱器結構如圖1所示。在吸熱器腔內沿周向排列了20根鈉熱管，每根熱管分為3段，分別是吸熱段、儲熱段和冷凝段(HSHX段)。靠近腔口的一段為吸熱段，該段在熱管上沒有任何附加物。中部為儲熱段，在儲熱段的熱管上套以多個分離的環型截面的相變材料容器，高溫相變材料封裝在容器內。最靠腔底的一段為冷凝段，熱管插入通過工質流體的翅片式換熱器。熱管的冷凝液回流系統由鈮粉燒結的毛細芯和干線組成。毛細芯控制冷凝液的周向分布，而干線則提供冷凝液的軸向回流，整根熱管的外壁面無論從周向還是軸向來看基本是等溫面。

在軌道日照期，吸熱段接受太陽輻射熱流，此時吸熱段成為熱管的蒸發段，提供熱量熔解PCM和加熱流體。在軌道陰影期，吸熱段除有少量通過腔口的輻射熱損失外，基本處于絕熱狀態，儲熱段轉變為蒸發段，此時PCM凝結放熱將熱量提供給冷凝段，使得空間站處于陰影期時仍能連續工作發電。

3 熱管吸熱器數值計算

3.1 計算物理模型

為了簡化計算，便于求解，在建立相變材料容器換熱計算模型時，作了以下假設：

1)熱管蓄熱段壁面溫度一致，即熱管中所有PCM容器的工作情況完全相同。

2)空穴初始容積固定，占容器容積的15%，且空穴分布在外壁處[10]。

圖2 熱管吸熱器相變材料容器計算物理模型 Fig.2 Physical model of PCM canister for heat pipe receiver

3)忽略液態PCM的對流影響。

4)PCM容器與熱管管壁間的接觸熱阻忽略不計，事實上二者間為間隙配合，約有5～10m的間隙。

空穴內的輻射換熱計算基于以下假設[11-12]：

1)所有空穴表面均為漫反射灰體表面。

2)PCM表面吸收所有波長的輻射。

本文采用的熱管吸熱器相變材料容器計算物理模型如圖2所示。其采用二維軸對稱模型，R為徑向，Z為軸向，流體區域為PCM及空穴，固體區域為外壁、側壁及內壁(包含熱管壁)。

3.2 數學模型

(1)PCM區模型

PCM及其容器壁的能量控制方程是能量守恒方程。在二維柱坐標(R，Z)情況下，采用焓法表示的導熱控制能量守恒方程：

式中 以比焓H的形式已把相變的影響考慮進去，適用于整個求解區域。其中焓與溫度的關系式為：

式中H為比焓；ρ為密度；k為導熱系數；T為溫度；t為時間；c為比熱容；Tm為相變溫度；ΔHm為物質單位質量的相變潛熱。

(2)PCM空穴模型

空穴內軸向溫度梯度很小，因此空穴內蒸氣溫度分布按徑向穩態傳熱方程確定：

3.3 計算方法

將汽—液—固區作為統一求解域，采用有限控制容積法對上述方程進行離散，將上述復雜的偏微分方程轉化成代數方程組，然后用壓力耦合方程組的半隱式算法和線—線迭代法進行數值求解。

3.4 初始和邊界條件

(1)邊界條件

式中l為熱管長度；Rco為蓄熱容器外徑；Rho為熱管外徑；khw為熱管壁面導熱系數；Thw為熱管外壁溫度；Q為熱流；A為蓄熱容器內表面面積。

設置蓄熱容器內壁(包含熱管壁)為周期性熱流邊界條件，蓄熱容器外壁絕熱，兩側壁絕熱，PCM區與蓄熱容器壁以及空穴與蓄熱容器壁的傳熱邊界均為導熱和輻射，相變材料為熔鹽80.5%LiF-19.5%CaF2，蓄熱容器內壁、外壁及側壁厚度均為1.5 mm。

(2)初始條件

PCM容器的初始溫度低于PCM熔點的溫度，取為950 K，全部PCM均處于固態。

軌道高度500 km，軌道周期90 min，其中日照期54 min，陰影期36 min。

在軌道周期內，用于熱分析的模型是時間的函數。模型有兩種熱輸入，一種是工質吸收的熱量是熱管壁溫的函數(見圖3)，假定蓄熱段和冷凝段熱管的壁面溫度相同，如前所述，這種假設是合理的；另一種熱輸入是靠吸熱段傳遞給蓄熱段熱管壁的，熱流是時間的函數，在日照期有一定值，在陰影期其值為0。由于太陽熱流對稱，蓄熱段溫度相同，吸熱器內有480個容器，因此每個蓄熱容器所接受的熱流為總熱流的1/480[9]。

可是，輸入熱流的定值并不是已知的，在分析過程中，通過軌道的性能參數可估計不同熱流輸入，熱流輸入在陰影末期(日照初期)使熱管的壁溫達到了1 100 K。熱流通過合適的吸熱器腔口和聚能器尺寸，很容易地轉化為全部的功率輸出。

給定聚能器尺寸，可以用計算機程序 CAV2來預測投射到腔壁內瞬態太陽入射熱流的分布[9]，如圖4所示。

圖3 功率與熱管壁溫之間關系Fig.3 Relationship between receiver power and heat pipe wall temperature

圖4 投射到吸熱器腔壁的太陽入射熱流分布Fig.4 Solar flux of the projection to cavity wall of heat receiver

3.5 計算結果及分析

圖5(a)～(d)分別為不考慮空穴日照期末(54 min)、考慮空穴時日照期末(54 min)、不考慮空穴陰影期末(90 min)及考慮空穴時陰影期末(90 min)PCM容器內部溫度場分布，從中可以看出：

1)在日照期，蓄熱容器內壁(包括熱管壁)溫度高于側壁溫度，側壁溫度高于外壁溫度，外壁溫度高于PCM區溫度，這說明在日照期，熱量是由熱管壁傳向蓄熱容器壁，再傳至PCM相變材料。

2)PCM容器內壁面與熱管壁相接觸，從熱管吸收熱流，容器中部的等溫線較直，而在容器的內外壁等溫線向內呈彎曲狀，整體上從容器的內壁到外壁有層次降低。有空穴時整個溫度場發生了較大的變化，容器中部等溫線更彎曲，容器內部溫度分布較無空穴時溫度分布更不均勻。

3)空穴的存在影響著PCM相變的進程，由于空穴熱阻遠大于容器側壁熱阻，其導熱和輻射使得PCM區溫度梯度顯著提高，可能導致該處熱應力過大，從而降低其使用壽命。日照期內，靠近容器兩側壁的PCM首先熔化，并逐漸沿軸向向內推進，由于空穴的存在使得容器壁(尤其是容器外壁與側壁交界處)內部的溫度梯度顯著，因此容器側壁在整個PCM容器的換熱過程中發揮了重要作用。

4)熱管壁溫在PCM熔點附近變化較小，蓄熱容器溫度變化不大，而基本型吸熱器在日照期末(54.6 min)考慮空穴時PCM溫度遠遠超過PCM熔點，最低超過52 K，最高超過達100 K[6]，PCM容器過熱，因此熱管吸熱器較基本型吸熱器熱性能更穩定。

5)在陰影期，PCM區溫度高于外壁溫度，外壁溫度高于側壁溫度，側壁溫度高于蓄熱容器內壁(包括熱管壁)溫度，這說明在陰影期，熱量是由PCM相變材料傳向蓄熱容器壁，再傳至熱管壁。

6)可以看出進入陰影期之后，PCM釋放潛熱，蓄熱容器的溫度分布都是從中心到外部降低，空穴的存在使得容器的蓄熱能力降低。

7)空穴的存在加大了容器外壁及其壁與交界處在日照期與陰影期的溫差，即增大了這些地方的溫度梯度，從而增大了它們的熱應力。

本文計算熱管壁面和PCM容器外壁面在0～20個軌道周期內的溫度變化及在一個穩定周期內的溫度變化分別如圖6(a)～(b)所示，本文計算結果與NASA方案計算結果比較[9]如圖7所示。從中可以看出：本文數值計算所得一個周期內熱管壁面、蓄熱容器外壁面溫度變化范圍和趨勢與NASA方案計算所得到的溫度變化范圍和趨勢比較接近，對應處各點的溫度值和NASA方案計算值相差較小。由于本文在計算時采用的邊界條件和NASA方案完全相同，所不同的是計算模型中將相變材料LiF的相變溫度設置為1 117～1 125 K的溫度區間，更加符合LiF的實際熔化過程，而NASA方案計算模型中采用了1 121 K的相變溫度點，另外本文計算模型中既考慮了導熱，也考慮了輻射，這更加全面地考慮了實際換熱過程，而NASA方案計算模型中僅僅考慮了導熱，因此本文計算所得的溫度響應速度比NASA方案計算結果要慢，并且對應各處溫度更加準確。

圖5 PCM容器溫度等值線分布Fig.5 Temperature field distribution in PCM canister

圖6 熱管吸熱器蓄熱單元熱循環計算結果Fig.6 Calculation results of heat pipe receiver during orbit cycles

圖7 本文計算結果與NASA方案計算結果比較分析 Fig.7 Numerical results comparison between ours and NASA′s

通過分析后可以發現：熱管壁溫和PCM容器外壁溫度無論是日照期還是陰影期，在PCM熔點附近波動較小，熱管吸熱器性能較穩定。這是由于熱管有很強的導熱性和良好的等溫性，熱管軸向溫差很小，這就使得熱管不同位置上蓄熱容器內的PCM在日照期時都能均勻、同步、一致地從熱管吸收熱量而蓄熱，在陰影期時都能均勻、同步、一致地釋放潛熱給熱管，從而對冷凝段的循環工質實現均勻加熱，于是循環工質出口溫度較小，保證系統的穩定可靠運行。

4 結束語

在同時考慮空穴和相變的情況下，對微重力條件下熱管吸熱器相變材料容器進行了熱分析，并與NASA方案計算結果進行了比較。

1)通過比較后發現：本文數值計算所得到的一個周期內熱管壁面、蓄熱容器外壁面溫度變化范圍和趨勢與NASA方案計算所得到的溫度變化范圍和趨勢比較接近，對應處各點的溫度值和NASA方案計算值相差不大，從而驗證了本文所建立的微重力條件下熱管吸熱器蓄熱單元計算模型的合理性和準確性，能夠較好地模擬熱管吸熱器蓄熱單元蓄/放熱過程，為今后熱管吸熱器的進一步研究分析打下了良好的基礎。

2)空穴的存在影響著PCM相變的進程，空穴熱阻和PCM區熱阻遠大于容器側壁熱阻，容器側壁在整個PCM容器的換熱過程中發揮了重要作用。日照期內，靠近容器兩側壁的PCM首先熔化，并逐漸沿軸向向內推進，容器中部的等溫線較直，而在容器的內外壁等溫線向內呈彎曲狀，整體上從容器的內壁到外壁有層次降低；陰影期內，容器中的溫度分布都是從中心到外部降低，空穴的存在使得容器的蓄熱能力降低，有空穴時部分PCM處于糊態。

3)日照期末，絕大部分PCM熔化，發生相變，由固態變成液態，并儲存潛熱；陰影期末，絕大部分PCM凝固，發生相變并釋放潛熱。由于熱管徑向溫差較小，熱管壁溫在PCM熔點附近變化較小，從而在一定程度上能緩解熱斑和熱松脫現象，熱管吸熱器熱性能較穩定。

[1] 崔海亭，袁修干，邢玉明，等.空間站太陽能熱動力發電系統研究進展[J]. 中國空間科學技術，2002，22(6)：34-42.

CUI HAITING, YUAN XIUGAN, XING YUMING, et al. The research and development of space solar dynamic power system [J]. Chinese Space Science and Technology, 2002,22(6)：34-42.

[2] NICOLAS PETER. Space power and its implications—the case of europe[J].Acta Astronautica, 2010, 66(3-4):348-354.

[3] 楊陽，段寶巖，黃進,等.OMEGA型空間太陽能電站聚光系統設計[J].中國空間科學技術，2014，34(5)：18-23.

YANG YANG, DUAN BAOYAN, HUANG JIN,et al.SSPS-OMEGA: A new concentrator system for SSPS [J]. Chinese Space Science and Technology, 2014,34(5)：18-23.

[4] 王希季.中國載人航天工程的外部設計[J]. 中國空間科學技術，2002，22(5)：1-8.

WANG XIJI. External design of chinese manned spaceflight engineering [J]. Chinese Space Science and Technology, 2002,22(5)：1-8.

[5] MURAT M KENISARIN. High-temperature phase change materials for thermal energy storage [J]. Renewable and Sustainable Energy Reviews, 2009: 11-16.

[6] ABDALLA ALRASHDAN, AHMAD TURKI MAYYAS, SAID AL-HALLAJ. Thermo-mechanical behaviors of the expanded graphite-phase change material matrix used for thermal management of Li-Ion battery packs [J]. Journal of Materials Processing Technology, 2010, 210(1): 174-179.

[7] STRUMPF H J, COOMBS M G. Advanced heat receiver conceptual design study [R]. N88-25977, Lewis Research Center, 1988.

[8] ROGER A CRANE. Thermal evaluation of advanced solar dynamic heat pipe receiver performance[R]. NASA Contractor Report,185117,1989.

[9] 崔海亭.空間站太陽能熱動力發電系統吸熱/蓄熱器優化研究[D]. 北京：北京航空航天大學，2003.

CUI HAITING. Optimization research of heat receiver for space solar thermal power system [D]. Beijing: Beijing University of Aeronautics and Astronautics, 2003.

[10] KESLAKE T W, IBRAHIM M B. Two-dimensional model of space station freedom thermal energy storage canister[J]. Journal of Solar Energy Engineering, 1992, 114(5): 114-121.

[11] 董克用,袁修干.微重力條件下相變材料容器的二維瞬態熱分析[J].太陽能學報，1998，19(1):41-47.

DONG KEYONG, YUAN XIUGAN. Two dimensional transient heat transfer analysis of a phase change material container [J]. Journal of Solar Energy, 1998, 19(1):41-47.

[12] 董克用．高溫固液相變蓄熱容器的熱分析[D].北京：北京航空航天大學,1999.

DONG KEYONG. Thermal analysis of high temperature solid-liquid phase change thermal energy storage container [D]. Beijing: Beijing University of Aeronautics and Astronautics, 1999.

桂小紅 1976年生，2009年獲北京航空航天大學人-機-環境工程專業博士學位，副研究員。研究方向為空間太陽能熱動發電系統設計及分析。

(編輯：王曉宇)

Transient Thermal Analysis of Heat Pipe Receiver under Microgravity

GUI Xiaohong1SONG Xiang′e2

(1 Institute of Engineering Thermophysics, Chinese Academy of Sciences, Beijing 100190) (2 Beijing International Studies University, Beijing 100024)

High temperature heat pipe receivers are commonly used as a core component of space solar dynamic power system. On the base of heat conduction governing differential equation under microgravity, two-dimensional numerical analysis of PCM container was conducted by a enthalpy method. With both void cavity and phase change considered, thermal performance of heat pipe receiver was analyzed. Numerical results were compared with NASA results.The accuracy of calculation model under gravity was verified. The results indicate that void cavity influences the process of phase change. The void cavity reduces the utility and thermal storage ability of PCM. The temperature gradient of PCM zone is very significant because of the void cavity. The PCM contained in the integrated heat pipe performs the averaging function of the heat loads. Normal working of wick ensures the uniformity of heat pipe, thus heat pipe receiver alleviates thermal spot and thermal ratcheting.

Space solar power system; Heat pipe receiver; Phase change material; Void cavity; Phase change;Spacecraft

國家自然科學基金(51476172)資助項目

2015-03-20。收修改稿日期：2015-04-30

10.3780/j.issn.1000-758X.2015.04.007