趙月晶 王曉濤 戴 巍* S Vanapalli 陳燕燕 羅二倉

(1中國科學(xué)院理化技術(shù)研究所低溫工程學(xué)重點實驗室 北京 100190)

(2中國科學(xué)院大學(xué) 北京100049)

(3荷蘭特文特大學(xué)能源、材料及系統(tǒng)科學(xué)技術(shù)學(xué)院 荷蘭恩斯赫德7500 AE)

1 引言

脈沖管制冷機是一種近幾十年來廣泛發(fā)展的小型低溫制冷機。因其冷側(cè)不存在運動部件，具有結(jié)構(gòu)簡單、運行可靠性高以及使用壽命長等優(yōu)點［1-2］。隨著現(xiàn)代信息技術(shù)、空間技術(shù)、超導(dǎo)電子學(xué)、紅外探測、低溫生物醫(yī)學(xué)等行業(yè)飛速發(fā)展，尤其是低溫電子學(xué)器件及低溫超導(dǎo)磁體在各領(lǐng)域的推廣應(yīng)用，對小型液氦溫區(qū)制冷機提出更多的需求和更高的要求。目前液氦溫區(qū)小型低溫制冷機主要為GM制冷機或GM型脈沖管制冷機這兩款機型的共同點是需要有閥壓縮機和配氣閥來提供低頻高壓比的壓力波(通常2 Hz以內(nèi)，壓比2—3)。這樣的結(jié)構(gòu)中存在很嚴(yán)重的能量損失:氦氣通過單向閥和配氣閥時由于壓降會造成不可逆損失［3］。

針對上述問題，近年來的一個前沿領(lǐng)域為利用斯特林型脈沖管制冷機獲得液氦溫度，但制冷量以毫瓦計，系統(tǒng)的熱效率也遠低于GM型制冷機［4-8］。根本原因在于，液氦溫區(qū)的回?zé)崞魈盍弦话氵x取磁性材料，形式為顆粒狀或球狀，直徑比較小，一般為百微米級，故而回?zé)崞鞯目紫堵时容^低，采用直線壓縮機進行驅(qū)動，在高頻交變流動狀態(tài)下，流動阻力帶來的損失比較大，進而造成整機的熱效率比較低，而直線壓縮機提供低頻的壓力波又比較困難，故整個技術(shù)發(fā)展還不太成熟。

本課題組提出利用VM型熱壓縮機驅(qū)動低溫脈沖管制冷機這一新結(jié)構(gòu)形式獲得液氦溫度，如圖1所示，①為高效液氮溫區(qū)斯特林型脈沖管制冷機，②為熱壓縮機，③為低溫脈沖管制冷機。其中由高效液氮溫區(qū)斯特林型脈沖管制冷機提供VM型熱壓縮機所需要的冷量，熱壓縮機利用室溫與液氮溫區(qū)之間的溫度梯度產(chǎn)生壓力波動，進而驅(qū)動低溫脈沖管制冷機，最終達到液氦溫區(qū)。本系統(tǒng)中VM熱壓縮機具有運行頻率低、驅(qū)動功小、無噪音、可靠性高等優(yōu)點，藉此可能獲得比已有GM型液氦溫區(qū)制冷機更高的系統(tǒng)效率、更高的功率密度和更高的可靠性及壽命。整個系統(tǒng)結(jié)構(gòu)中，VM型熱壓縮機的性能是能否獲得液氦溫度的關(guān)鍵，并對整機效率以及冷量有很大影響。本文接下來針對無負載熱壓縮機進行了數(shù)值模擬與實驗研究，分析了熱壓縮機排出器位移、工作頻率以及冷端溫度對壓比的影響規(guī)律。

2 實驗裝置介紹

圖1 VM型熱壓縮機驅(qū)動低溫脈沖管制冷機Fig.1 VM thermal compressor driven pulse tube cryocooler

無負載時，熱壓縮機是一個等質(zhì)量和等體積的系統(tǒng)，隨著排出器在高低溫之間的運動，系統(tǒng)內(nèi)部氣體溫度發(fā)生變化，從而產(chǎn)生壓力波動，對應(yīng)的聲場以駐波為主。圖2是測試熱壓縮機工作性能的實驗裝置圖，其中包括預(yù)冷系統(tǒng)和熱壓縮機系統(tǒng)。

圖2 實驗裝置圖Fig.2 Experiment apparatus picture

預(yù)冷用兩臺制冷機均采取CFIC 2S241 W直線壓縮機驅(qū)動，相當(dāng)于一臺直線壓縮機驅(qū)動兩臺制冷機，壓縮機與兩臺制冷機之間通過紫銅管連接，為了保證安全性，連接管采取水套冷卻的方式進行冷卻。系統(tǒng)運行時真空罩內(nèi)部真空度維持在5.8×10-3Pa左右。

表1為熱壓縮機各部件的主要參數(shù)值，其中回?zé)崞魈畛涞?0目不銹鋼絲網(wǎng)，回?zé)崞鞅诤駷?.5 mm。排出器位移最大值為7 mm，其驅(qū)動電機的最大工作電壓為28 V，工作電流不許超過6 A。

表1 熱壓縮機的各部件參數(shù)值Table 1 Main structure parameters of each component in thermal compressor

3 無負載熱壓縮機性能的實驗研究

兩臺預(yù)冷用制冷機的冷端和熱壓縮機的冷端采用熱橋的形式連接，具體形式如圖3所示，為了緩沖溫度不均勻以及材料的不同產(chǎn)生的熱應(yīng)力，兩臺制冷機的冷端均采取了軟連接的形式。制冷機的冷頭以及熱橋?qū)ΨQ布置了8個溫度計，其中兩臺制冷機的冷頭上標(biāo)號為T2和T7，熱壓縮機冷端標(biāo)號為T3和T5。

圖3 溫度計布置圖Fig.3 Thermometer distribution picture

開展無負載熱壓縮機性能的實驗研究，熱橋最低溫達到了66.2 K，從圖中可以看到兩臺制冷機與熱橋上各點的降溫曲線，大約經(jīng)過3個小時溫度趨于穩(wěn)定。

圖4 熱橋的降溫曲線Fig.4 Cool-down curve of thermal bridge

表2為熱壓縮機非工作狀態(tài)下預(yù)冷系統(tǒng)開啟后熱橋上的溫度分布。熱壓縮機冷端對應(yīng)的最低溫為66 K，兩臺預(yù)冷用制冷機的冷頭分別為61 K和59 K。由于兩臺制冷機接有相對大的熱負載(熱壓縮機)，實際運行時制冷機的冷頭溫度比無負荷溫度高。

表3 熱橋上各點的最低溫Table 3 The lowest temperature of thermal bridge

圖5—圖8給出了無負載下熱壓縮機的運行規(guī)律，圖5、圖6和圖7顯示，在不同頻率運行下，熱壓縮機的壓比隨著排出器位移的增大而增大，基本上呈線性關(guān)系，其中運行在3 Hz，位移為6.0 mm時候，壓比達到了最大值1.3。圖8顯示不同熱壓縮機冷頭溫度對壓比位移的影響。從圖中可以看出隨著熱壓縮機冷端溫度的下降，熱壓縮機的壓比也在升高。另外，從圖9(雖熱壓縮機的冷端溫度不一致，但預(yù)冷用制冷機的工作狀態(tài)一樣)中可以看出，在相同位移下，熱壓縮機工作在3Hz的時候，所獲得的壓比更大些。

圖5 1 Hz時熱壓縮機位移與壓比的關(guān)系曲線Fig.5 Dependence of pressure ratio on displacement when thermal compressor works at 1 Hz

圖6 3 Hz時熱壓縮機位移與壓比的關(guān)系曲線Fig.6 Dependence of pressure ratio on displacement when thermal compressor works at 3 Hz

圖7 5 Hz時熱壓縮機位移與壓比的關(guān)系曲線Fig.7 Dependenceof pressure ratio on displacement when thermal compressor works at 5 Hz

圖8 熱壓縮機冷端溫度與壓比的關(guān)系曲線Fig.8 Dependence of pressure ratio on temperature when thermal compressor works at 3 Hz

圖9 不同頻率下熱壓縮機位移與壓比的關(guān)系Fig.9 Dependence of pressure ratio on displacement at different frequencies

4 無負載熱壓縮機性能的實驗結(jié)果與數(shù)值模擬結(jié)果的比較

采用Sage軟件進行數(shù)值模擬分析［9-10］，采用圖形化界面，對熱壓縮機的各個部件分別建立模塊，輸入部件的幾何尺寸、材料物性等參數(shù)，各模塊通過質(zhì)量、壓力、能量的連接形成整機，在給定合適的邊界條件后即可求解。針對各實驗工況，對無負載熱壓縮機性能的實驗結(jié)果與數(shù)值模擬結(jié)果進行比較，見圖10—圖13。從圖中看出數(shù)值模擬計算結(jié)果與實驗結(jié)果趨勢基本一致，尤其是3 Hz、5 Hz時，計算結(jié)果與實驗結(jié)果基本相當(dāng)。本文建立的數(shù)值計算模型可預(yù)測熱壓縮機的工作性能，并且能夠指導(dǎo)熱壓縮機的設(shè)計、實驗以及改進提高。

圖10 1 Hz的實驗結(jié)果與理論計算比較Fig.10 Comparisons between calculated and experimental results at 1 Hz

圖11 3 Hz的實驗結(jié)果與理論計算比較Fig.11 Comparisons between calculated and experimental results at 3 Hz

圖12 5 Hz的實驗結(jié)果與理論計算比較Fig.12 Comparisons between calculated and experimental results at 5 Hz

圖13 溫度對壓比的影響實驗與計算的比較Fig.13 Comparisons between calculated and experimental results at different temperatures

5 總結(jié)

根據(jù)本文實驗結(jié)果以及與數(shù)值模擬結(jié)果的對比可以得出:無負載情況下熱壓縮機所獲得的壓比主要受頻率、排出器位移(即掃氣量)和冷端溫度的影響。排出器位移越大，所獲壓比越高;冷端溫度越低，所獲壓比越高，并且兩者結(jié)果吻合良好。在現(xiàn)有測試情況下，無負載熱壓縮機工作在3Hz時能夠獲得較高壓比，其中無負載情況下熱壓縮機工作在3Hz、位移為6.0mm時，可以獲得最大壓比為1.3。建立的數(shù)值模型可以很好地預(yù)測熱壓縮機的工作性能。

1 Radebaugh R.Status and recent trends in cryocooler research［C］.Proceedings of International Conference on Cryogenics and Refrigeraation，2008:11-22.

2 Radebaugh R.Cryocoolers:the state of the art and recent developments［J］.Journal of Physics:condensed matter，2009，21:164219

3 Gan zhihua，Qiu Limin，et al.A single-stage G-M type pulse tube cryocooler operating at 10.6K［J］.Cryogenics，2009(49):198-201.

4 Olson J，Nast T C，Evtimov B.Development of a 10K pulse tube cryocooler for space application［J］.Cryogenics，2003(12):241-246.

5 Olson J，Moore M，Champagne P，et al.Development of a space-type-4-stage pulse tube cryocooler for very low temperature［C］.Advances in Cryogenic Engineering，New York:AIP Press，2006，51(A):623-631.

6 Nguyen T，Colbert R，Durand D，et al.10K pulse tube cooler［J］.Cryocoolers，2007(14):27-31.

7 Qiu Limin，Cao Qiang，Zhi Xiaoqin，et al.A three-stage Stirling pulse tube cryocooler operating below the critical point of helium-4［J］.Cryogenics，2011(51):609-612.

8 Qiu Limin，Cao Qiang，Zhi Xiaoqin，et al.Operating characteristics of a three-stage Stirling pulse tube cryocooler operating near 5K［J］.Cryogenics，2012(52):382-388.

9 David Gedeon.Sage user’s guide-Stirling pulse-tube and low-T cooler model classes［M］.Gedeon Associates，2011.

10 Zhao Y，Dai W，Chen Y，et al.Numerical analysis of 4K regenerator of a pulse tube cryocooler［C］.Proeedings of the China academy of engineering thermal physics conference，2014.