何超峰，張俊峰，武義鋒，仰 葉，丁先庚

(1.中國(guó)電子科技集團(tuán)公司第十六研究所，合肥230043;2.安徽萬(wàn)瑞冷電科技有限公司，合肥230088)

1 引言

高溫超導(dǎo)電機(jī)以其體積小、重量輕、效率高和振動(dòng)噪音低等優(yōu)點(diǎn)在全世界范圍內(nèi)引起了廣泛關(guān)注［1－3］。目前已有多個(gè)國(guó)家在高溫超導(dǎo)電機(jī)領(lǐng)域取得了階段性成果，包括美國(guó)、德國(guó)、韓國(guó)、日本等國(guó)家，其中以美國(guó)36.5MW［4］高溫超導(dǎo)電機(jī)最受業(yè)內(nèi)關(guān)注。同時(shí)隨著超導(dǎo)風(fēng)電、超導(dǎo)電纜、超導(dǎo)限流器、超導(dǎo)磁儲(chǔ)能技術(shù)等快速發(fā)展，相信高溫超導(dǎo)的應(yīng)用領(lǐng)域會(huì)越來(lái)越廣。

在所有的超導(dǎo)應(yīng)用中，超導(dǎo)磁體是其中的核心部件。超導(dǎo)磁體一般工作在35 K左右的低溫，為達(dá)到這樣的低溫條件，需要持續(xù)不斷的給磁體輸送冷量。在高溫超導(dǎo)電機(jī)中，冷量是由AL330 GM制冷機(jī)產(chǎn)生，制冷機(jī)冷頭首先冷卻冷媒，將冷媒液化，液態(tài)冷媒進(jìn)入電機(jī)轉(zhuǎn)子，然后蒸發(fā)汽化帶走磁體產(chǎn)生的損耗和漏熱，如圖1－a所示。

在低溫30 K左右下，制冷機(jī)的制冷量與制冷溫度有很大關(guān)系［5］，如圖1－b所示，每變化一度，制冷機(jī)制冷量會(huì)變化較大，因此，為最大限度的利用制冷機(jī)的冷量，需要制冷機(jī)冷頭與冷媒之間的溫差盡可能的小。

2 理論分析

高溫超導(dǎo)電機(jī)低溫系統(tǒng)冷媒采用氖氣，氖氣在標(biāo)準(zhǔn)大氣壓下沸點(diǎn)為27.1 K，氖氣通過(guò)制冷機(jī)冷頭降溫液化為液體。由公式(1)，假設(shè)換熱器溫差為0.8 K，GM制冷機(jī)冷頭溫度為26.3 K，在此溫度下，制冷機(jī)制冷量為62 W。

對(duì)流換熱公式:

式中:Q為制冷機(jī)的制冷量，W;h為對(duì)流換熱系數(shù)，W/(m2·K);A為對(duì)流換熱面積，m2;△T為溫差，K。

由(1)式可知，當(dāng)制冷機(jī)制冷量一定，對(duì)流換熱系數(shù)一定，要減小溫差，須增大對(duì)流換熱面積。

圖1 AL330制冷機(jī)冷頭及性能曲線(xiàn)

GM制冷機(jī)冷量主要集中在冷頭的下端面，因此，需設(shè)計(jì)與此端面相配合的高效換熱器，才能充分利用制冷機(jī)冷量。對(duì)流換熱系數(shù)可由下式［6］計(jì)算:

式中:ρ為液氖在低溫下的密度，kg/m3;λ為液氖的導(dǎo)熱系數(shù)，W/(m·K);l為液氖氣化潛熱，kJ/kg;η為動(dòng)力黏度，Pa·s;ts為飽和蒸汽溫度，K;tw為換熱器壁面的溫度，K;l為豎壁的高度，m。

計(jì)算可得對(duì)流換熱系數(shù)為5 500 W/(m2·K)。

由 Q=62 W，△T=0.8 K，根據(jù)式(1)，(2)可得，對(duì)流換熱面積為 0.014 m2。

3 數(shù)值計(jì)算及分析

根據(jù)理論計(jì)算得出冷頭換熱器需要的換熱面積，設(shè)計(jì)一換熱器結(jié)構(gòu)形式，如圖2所示。材料選用牌號(hào)為T(mén)2［7］的紫銅，考慮換熱器的安裝方便，留出安裝螺栓的空間。

圖2 冷頭換熱器結(jié)構(gòu)形式

根據(jù)此結(jié)構(gòu)形式，建立三維計(jì)算模型，利用ANSYS熱分析模塊進(jìn)行溫度場(chǎng)分布計(jì)算。

3.1 邊界條件設(shè)置

在ANSYS中建立三維模型，選擇熱分析模塊［8］，設(shè)定各部分物性參數(shù)，劃分網(wǎng)格，邊界條件加載。對(duì)流換熱系數(shù):5 500 W/(m2·K);紫銅低溫的導(dǎo)熱系數(shù)［9］:800 W/(m·K);與冷頭接觸面熱流密度:8 025 W/m2;

3.2 結(jié)果及分析

在ANSYS中求解，對(duì)流換熱系數(shù)為5 500 W/(m2·K)情況下，計(jì)算出換熱器溫度場(chǎng)分布如圖3－a所示。同時(shí)考慮實(shí)際情況，若氖氣冷凝換熱系數(shù)實(shí)際值低于理論計(jì)算值，取保守值3 000 W/(m2·K)，計(jì)算結(jié)果如圖3－b所示。

圖3 冷頭換熱器溫度場(chǎng)分布

從計(jì)算結(jié)果可以看出，對(duì)流換熱系數(shù)為5 500 W/(m2·K)時(shí)，換熱器上最大溫差為0.84 K，高于設(shè)計(jì)值0.8 K，這說(shuō)明換熱器上有一部分換熱面積沒(méi)有起到作用。換熱器的溫差主要集中在收縮段，與冷頭接觸面溫度為26.2 K，低于設(shè)計(jì)值26.3 K，制冷機(jī)在26.2 K 時(shí)制冷量達(dá)不到62 W;對(duì)流換熱系數(shù)為3 000 W/(m2·K)，換熱器上最大溫差變大至1.2 K，且溫差仍集中于收縮段。

從以上分析可以得知，這種換熱器結(jié)構(gòu)形式熱阻較大，造成換熱器上溫差偏大，尤其是當(dāng)對(duì)流換熱系數(shù)較小時(shí)，這對(duì)制冷機(jī)冷頭換熱器效果影響較大。

圖4 改進(jìn)后換熱器形式

4 換熱器優(yōu)化設(shè)計(jì)

根據(jù)上述計(jì)算結(jié)果，換熱器熱阻大，溫差主要集中于收縮段，因此，考慮盡可能減小收縮段，同時(shí)要考慮換熱器安裝方便，給固定螺栓留出安裝空間，即將螺栓孔上方對(duì)應(yīng)的位置挖空，如圖4所示，對(duì)改進(jìn)后的換熱器進(jìn)行仿真分析，邊界條件設(shè)置如前，結(jié)果如圖5所示。

圖5 改進(jìn)后換熱器上溫度場(chǎng)分布

從計(jì)算結(jié)果可以得出，當(dāng)對(duì)流換熱系數(shù)為5 500 W/(m2·K)時(shí)，換熱器上最大溫差為0.37 K，與冷頭接觸面溫度為26.6 K;當(dāng)對(duì)流換熱系數(shù)為3 000 W/(m2·K)時(shí)，換熱器上最大溫差為0.51 K，與冷頭接觸面的溫度為26.5 K，均滿(mǎn)足設(shè)計(jì)要求。從圖5－b還可以看出，當(dāng)對(duì)流換熱系數(shù)小于理論計(jì)算值時(shí)，換熱器上一部分區(qū)域溫度低于27 K，也就是說(shuō)換熱面積不夠，因此，需要適當(dāng)增加換熱面積。

5 結(jié)論

通過(guò)以上計(jì)算分析，得出:

(1)當(dāng)冷頭換熱器形式相同時(shí)，對(duì)流換熱系數(shù)越小，冷頭換熱器溫差越大;

(2)對(duì)于低溫冷凝換熱器，導(dǎo)熱熱阻所占的比例較大，因此，在設(shè)計(jì)過(guò)程中要盡量減小導(dǎo)熱段的長(zhǎng)度，并適當(dāng)增加導(dǎo)熱面積。

［1］BIAO Chen，GU Guobiao，ZHANG Guoqiang，et al.Analysis and Design of Cooling System in High Temperature Superconducting Synchronous Machines［J］.IEEE Transaction On Applied Superconductivity，2003，17(2):1577－1560.

［2］MICHAEL Frank，PETER Van Habelt，PETER Kummeth，High-Temperature Superconducting Rotating Machines for Ship Applications［J］.IEEE Transaction On Applied Superconductivity，June 2006，16(2):1465－1468.

［3］SWARM S Kalsi.Development Status of Rotating Machines Employing Superconducting Field Windings［J］.Proceedings of The IEEE，Octo ber 2004，92(10):1688－1704.

［4］Swarn S Kalsi.HTS Ship Propulsion Motors for Podded Application［J］.American Superconductor Corporation，2006，6(3):52－54.

［5］Cryomech Inc.AL330 CP970 User Manual，2007，46.

［6］楊世銘.傳熱學(xué)［M］.北京:高等教育出版社，1998，186－188.

［7］陳國(guó)邦.低溫工程材料［M］.杭州:浙江大學(xué)出版社，1998，56.