EAST阻性換熱器多物理場耦合模擬計算

鄧 威，奚維斌

(中國科學(xué)院 等離子體物理研究所，安徽 合肥 230031)

?

EAST阻性換熱器多物理場耦合模擬計算

鄧 威，奚維斌

(中國科學(xué)院 等離子體物理研究所，安徽 合肥 230031)

阻性換熱器是EAST高溫超導(dǎo)(HTS)電流引線的重要組成部分，目前有三頭螺旋翅片和疊片兩種結(jié)構(gòu)形式，為了比較這兩種阻性換熱器的優(yōu)劣，對它們的熱工水力性能進(jìn)行了多物理場耦合模擬計算，計算結(jié)果表明：兩種阻性換熱器在換熱性能方面基本相當(dāng)，均可滿足快速換熱的要求，但疊片換熱器的流動阻力遠(yuǎn)小于三頭螺旋翅片換熱器的。實(shí)際運(yùn)行過程中，三頭螺旋翅片換熱器中氮冷卻回路的壓力控制較為困難，經(jīng)常需人工調(diào)節(jié)控制閥閥門，而疊片換熱器中氮冷卻回路的壓力控制則較為簡單，不需經(jīng)常調(diào)節(jié)。因此，疊片式結(jié)構(gòu)較三頭螺旋翅片式結(jié)構(gòu)更適合應(yīng)用在EAST阻性換熱器中。

高溫超導(dǎo)電流引線；阻性換熱器；熱工水力性能；多物理場耦合

EAST托卡馬克裝置的核心部件是大型超導(dǎo)磁體系統(tǒng)，由縱場和極向場兩個子系統(tǒng)組成[1-2]。EAST利用13對二元高溫超導(dǎo)電流引線連接室溫電源與低溫超導(dǎo)磁體系統(tǒng)，為超導(dǎo)磁體的運(yùn)行提供勵磁電流通道和熱連接[3-5]。高溫超導(dǎo)電流引線主要由工作在室溫與液氮溫度之間的阻性換熱器段和工作在液氮溫度與液氦溫度之間的高溫超導(dǎo)(HTS)段組成[5-6]。阻性換熱器是HTS段與室溫電源的連接部分，其作用是在液氮溫度以上的溫區(qū)承載勵磁電流。為了帶走阻性換熱器中從室溫端流入的傳導(dǎo)熱和通電時產(chǎn)生的焦耳熱，減少流向HTS段的漏熱，通常的做法是采用冷氦氣對阻性換熱器進(jìn)行強(qiáng)迫對流換熱冷卻。在EAST HTS電流引線中，為了減少液氦的消耗，阻性換熱器采用液氮自然蒸發(fā)對流換熱的方式進(jìn)行冷卻[7-8]。EAST阻性換熱器的冷卻包括液氮自然蒸發(fā)和氮蒸氣流動換熱兩個自洽的過程。當(dāng)?shù)獨(dú)庠谧栊該Q熱器中的流動阻力較大時，氮冷卻回路中壓力的不平衡度較大，一方面會增加冷卻回路的控制難度，另一方面會導(dǎo)致阻性換熱器的運(yùn)行偏離設(shè)計工況，不利于電流引線的長期安全穩(wěn)定運(yùn)行。

了解氮?dú)庠谧栊該Q熱器特殊流道中的流動和換熱狀況將有助于阻性換熱器的設(shè)計和優(yōu)化，提高電流引線的運(yùn)行安全性。EAST阻性換熱器有三頭螺旋翅片和疊片兩種結(jié)構(gòu)形式，本文將對這兩種阻性換熱器的熱工水力性能進(jìn)行多物理場耦合模擬計算，并根據(jù)計算結(jié)果對其熱工水力性能進(jìn)行比較，從而為阻性換熱器的設(shè)計、優(yōu)化和運(yùn)行提供參考。

1 EAST阻性換熱器

圖1為EAST二元HTS電流引線的結(jié)構(gòu)組成示意圖，電流引線由室溫銅頭、阻性換熱器、過渡銅塊、HTS段和低溫超導(dǎo)段等組成。阻性換熱器位于HTS段和過渡銅塊的上方，依靠液氮?dú)饣a(chǎn)生的氮蒸氣冷卻：1) 從阻性換熱器流入過渡銅塊的漏熱使液氮筒中的液氮蒸發(fā)產(chǎn)生氮蒸氣；2) 冷氮?dú)庋負(fù)Q熱器軸向流動并和換熱器表面發(fā)生對流換熱，從而帶走阻性換熱器通電產(chǎn)生的焦耳熱和從室溫端流入的熱量。液氮蒸發(fā)速率(氮?dú)赓|(zhì)量流量)由流入過渡銅塊的漏熱量決定，因此冷氮?dú)獾漠a(chǎn)生、流動和對流換熱過程是一自洽的過程。液氮筒中的液氮既是制冷劑，又是熱沉，可起到穩(wěn)定HTS段熱端溫度的作用[7]。

圖1 EAST二元HTS電流引線的結(jié)構(gòu)組成示意圖Fig.1 Scheme of hybrid HTS current leads in EAST

EAST共有13對HTS電流引線，采用的阻性換熱器結(jié)構(gòu)形式有三頭螺旋翅片式和疊片式兩種。三頭螺旋翅片阻性換熱器(fin-HEX)的形狀是在空心銅棒上加工出三頭螺紋作為加強(qiáng)散熱能力的翅片，翅片之間的3道溝槽作為氮?dú)饬鲃拥耐ǖ?。fin-HEX通過中心的銅棒承載電流，通過三頭螺旋槽與冷氮?dú)膺M(jìn)行對流換熱，其優(yōu)點(diǎn)是氮?dú)饬鲃虞^強(qiáng)烈、換熱效果好、結(jié)構(gòu)簡單、制造中無需高溫釬焊[7,9]。疊片阻性換熱器(stack-HEX)由70層薄銅片等間隔層疊而成，兩端采用等厚的薄銅條進(jìn)行間隔固定。stack-HEX通過薄銅片承載電流，通過薄銅片之間的間隙與冷氮?dú)膺M(jìn)行對流換熱，其優(yōu)點(diǎn)是換熱面積大、換熱效果好、流動阻力小。

2 多物理場耦合模擬計算

阻性換熱器與氮?dú)庵g的對流換熱過程較復(fù)雜，除有流動和換熱過程外，還涉及焦耳熱效應(yīng)和熱傳導(dǎo)過程，因此現(xiàn)有的對流換熱經(jīng)驗(yàn)公式不能用來確定阻性換熱器的熱工水力性能。雖然可通過實(shí)驗(yàn)的方法來確定阻性換熱器的熱工水力性能，但是這樣不僅過程十分復(fù)雜，而且時間和經(jīng)濟(jì)成本均會很高。因此相對于可行性不高的實(shí)驗(yàn)方法，可利用三維多物理場耦合模擬計算的方法對氮?dú)庠谧栊該Q熱器特殊流道中的流動和換熱過程進(jìn)行仿真計算，從而獲得阻性換熱器的熱工水力性能。

電流引線通電運(yùn)行時，在阻性換熱器導(dǎo)體中存在焦耳熱效應(yīng)和熱傳導(dǎo)，其數(shù)學(xué)表達(dá)式分別為式(1)和(2)；在換熱器流道中存在流動和對流換熱，其數(shù)學(xué)表達(dá)式為式(3)～(7)。式(1)為焦耳熱效應(yīng)的微分方程，S、J和E分別為體積熱功率、電流密度和電場強(qiáng)度。式(2)為熱傳導(dǎo)過程的微分方程，t、x、y和z分別為時間坐標(biāo)和3個方向坐標(biāo)；Ts為導(dǎo)體的溫度；λs、ρs和cs分別為導(dǎo)體材料的導(dǎo)熱系數(shù)、質(zhì)量密度和比熱容。式(3)～(7)為流體流動和對流換熱的微分控制方程，分別表示質(zhì)量守恒、3個坐標(biāo)方向上的動量守恒和能量守恒。u、v、w、p和Tl分別為流體3個坐標(biāo)方向上的速度、壓力和溫度；ρl、μ、λl和cl分別為流體的密度、動力黏度、導(dǎo)熱系數(shù)和定壓比熱容；Sx、Sy、Sz、Bx、By、Bz和Sl分別為3個坐標(biāo)方向上的動量源項(xiàng)、體積力和熱源項(xiàng)。通過耦合求解上述微分方程，可得到阻性換熱器導(dǎo)體中電流及溫度分布和流道中流體流動和換熱情況。

(1)

(2)

(3)

(4)

(5)

(6)

(7)

fin-HEX和stack-HEX均使用在TF HTS電流引線中，設(shè)計電流為14 kA(最大可運(yùn)行電流大于16.5 kA)，材料為無氧銅。以下將對這兩種阻性換熱器中氮?dú)獾牧鲃忧闆r和氮?dú)馀c阻性換熱器的對流換熱情況進(jìn)行模擬計算。計算區(qū)域由固體域和流體域組成，固體域模擬阻性換熱器銅導(dǎo)體部分，通電產(chǎn)生焦耳熱，并傳導(dǎo)熱流；流體域模擬阻性換熱器的流道部分，冷氮?dú)庠诹鞯乐辛鲃?，并與流道壁面發(fā)生換熱。

模擬計算考慮材料的物性隨溫度的變化。氮?dú)獾奈镄詤?shù)來自于GasPak軟件，液氮的氣化潛熱取200 J/g；無氧銅的物性參數(shù)來自于Cryocomp軟件。施加的溫度場、電場和流場邊界條件包括：1) 阻性換熱器的底部浸泡在液氮中，因此固體域低溫端面和氮?dú)馊肟跍囟仍O(shè)置為78 K；2) 控制阻性換熱器室溫端熱流量使端面溫度為300 K；3) 固體域低溫端面施加相應(yīng)的電流密度，室溫端面設(shè)置為接地邊界條件；4) 估算從阻性換熱器流入過渡銅塊的漏熱量，冷氮?dú)馊肟谫|(zhì)量流量的初始值可根據(jù)漏熱估算值進(jìn)行設(shè)定；5) 氮?dú)獾某隹谄骄鶋毫θ?個大氣壓。冷氮?dú)馊肟谫|(zhì)量流量可通過迭代進(jìn)行調(diào)整，以使液氮蒸發(fā)速率與冷氮?dú)馊肟谫|(zhì)量流量相匹配。

3 阻性換熱器熱工水力性能模擬計算

3.1 三頭螺旋翅片阻性換熱器

圖2為fin-HEX的幾何模型，由帶三頭螺旋的無氧銅導(dǎo)體和螺旋槽流道組成。鑒于螺旋流道的復(fù)雜性，湍流模型使用剪切應(yīng)力輸運(yùn)模型，可有效預(yù)測流動分離。fin-HEX幾何形狀較復(fù)雜，模擬計算需使用完整的幾何模型，因此模擬計算所需的網(wǎng)格單元數(shù)目較大。圖3為fin-HEX計算模型的網(wǎng)格無關(guān)性驗(yàn)證曲線，可看出，氮?dú)鈮航惦S網(wǎng)格單元數(shù)目的增加而收斂于7 700 Pa左右。當(dāng)計算模型網(wǎng)格單元數(shù)大于400萬時，相對計算誤差小于5%，此時的計算結(jié)果是可信的。

圖2 fin-HEX的幾何模型Fig.2 Geometry model of fin-HEX

圖3 fin-HEX計算模型的網(wǎng)格無關(guān)性驗(yàn)證曲線Fig.3 Curve of grid independence on calculation model for fin-HEX

利用計算流體動力學(xué)(CFD)軟件對上述模型進(jìn)行計算，得到fin-HEX的熱工水力性能數(shù)據(jù)。圖4為在不同電流工況下溫度沿fin-HEX長度方向的分布。由圖4可看出：1) 當(dāng)電流較小時，溫度分布曲線呈略微下凹形狀，溫度最高點(diǎn)(熱點(diǎn))位于室溫端面，焦耳熱對冷端漏熱量的貢獻(xiàn)很小，流向過渡銅塊的漏熱量主要是傳導(dǎo)熱；2) 當(dāng)電流大小居中時，溫度分布曲線幾乎呈線性形狀，熱點(diǎn)位于室溫端面，焦耳熱和傳導(dǎo)熱對冷端漏熱量的貢獻(xiàn)相當(dāng)；3) 當(dāng)電流較大時，溫度分布曲線呈上凸形狀，熱點(diǎn)開始從室溫端面向下移動，焦耳熱對冷端漏熱量的貢獻(xiàn)較傳導(dǎo)熱的大。當(dāng)電流為16 kA時，fin-HEX上最高溫度略高于300 K，無燒毀的危險。由圖4可知，fin-HEX可滿足EAST HTS電流引線對溫度分布和散熱的需要。

圖4 不同電流下溫度沿fin-HEX長度方向的分布Fig.4 Temperature vs length of fin-HEX at different currents

圖5 fin-HEX的室溫端熱流、冷端熱流、焦耳熱生成速率和對流換熱速率隨電流的變化Fig.5 Heat flows at room temperature end and low temperature end, Joule heat generation rate and heat exchange rate vs current for fin-HEX

以銅導(dǎo)體為對象，熱流生成和流入為正，流出為負(fù)，則阻性換熱器室溫端熱流、冷端熱流、焦耳熱生成速率和對流換熱速率滿足能量守恒。圖5為fin-HEX的室溫端熱流、冷端熱流、焦耳熱生成速率和對流換熱速率隨電流變化的情況。隨電流增大，室溫端熱流減小、對流換熱速率增大，但焦耳熱生成速率增大的速率更大，因此冷端熱流即流入過渡銅塊的漏熱量是增大的。圖6為fin-HEX的氮?dú)赓|(zhì)量流量和壓降隨電流變化的情況。由于流入過渡銅塊的漏熱量隨電流增大，因此氮?dú)赓|(zhì)量流量和氮?dú)鈮航狄搽S之增大。當(dāng)電流為設(shè)計電流14 kA時，室溫端熱流的絕對值很小，說明fin-HEX的設(shè)計對優(yōu)化點(diǎn)的偏移很小，這表示在設(shè)計工況點(diǎn)阻性換熱器的室溫端面可保持在300 K左右。

圖6 fin-HEX的氮?dú)赓|(zhì)量流量和壓降隨電流的變化Fig.6 Mass flow rate and pressure drop of nitrogen gas vs current for fin-HEX

3.2 疊片阻性換熱器

stack-HEX的實(shí)際形狀為70層薄銅片等間隔層疊，可近似看作69個3層薄板疊加而成。每個3層薄板的中間層為0.4 mm厚的氮?dú)饬鲃訉樱膺厓蓪訛?.2 mm厚的銅片?？紤]到stack-HEX的幾何形狀和載荷具有對稱性，模擬計算時選擇換熱器的一3層薄板作為計算模型，其所有的固體邊界均可作為絕熱邊界處理。圖7為stack-HEX的幾何模型，流體域夾在兩片固體域之間。氮?dú)饬鲃訉拥牧鲃訝顩r為層流流動，不需使用湍流模型，因此計算較為簡單。圖8為stack-HEX計算模型的網(wǎng)格無關(guān)性驗(yàn)證曲線，可看出，氮?dú)鈮航惦S網(wǎng)格單元數(shù)目的增加而收斂于337.5 Pa左右。由圖8可知，當(dāng)計算模型網(wǎng)格單元數(shù)大于140萬時計算結(jié)果是可信的。

圖7 stack-HEX的幾何模型Fig.7 Geometry model of stack-HEX

利用CFD軟件對上述模型進(jìn)行計算，得到stack-HEX的熱工水力性能數(shù)據(jù)。圖9為在不同電流下溫度沿stack-HEX長度方向的分布。stack-HEX的相應(yīng)曲線和fin-HEX的類似，顯示出幾乎相同的形狀和趨勢。由圖9可知，stack-HEX也可滿足EAST HTS電流引線對溫度分布和散熱的需要。

圖8 stack-HEX計算模型的網(wǎng)格無關(guān)性驗(yàn)證曲線Fig.8 Curve of grid independence on calculation model for stack-HEX

圖9 不同電流下溫度沿stack-HEX長度方向的分布Fig.9 Temperature vs length of stack-HEX at different currents

圖10 stack-HEX的室溫端熱流、冷端熱流、焦耳熱生成速率和對流換熱速率隨電流的變化Fig.10 Heat flows at room temperature end and low temperature end, Joule heat generation rate and heat exchange rate vs current for stack-HEX

圖10為stack-HEX的室溫端熱流、冷端熱流、焦耳熱生成速率和對流換熱速率隨電流變化的情況。隨電流增大，室溫端熱流減小、對流換熱速率增大，但焦耳熱生成速率增大的速率更大，因此冷端熱流即流入過渡銅塊的漏熱量是增大的。圖11為stack-HEX的氮?dú)赓|(zhì)量流量和氮?dú)鈮航惦S電流變化的情況。由于流入過渡銅塊的漏熱量隨電流增大，因此氮?dú)赓|(zhì)量流量和氮?dú)鈮航狄搽S之增大。當(dāng)電流為設(shè)計電流14 kA時，室溫端熱流的絕對值很小，說明stack-HEX的設(shè)計對優(yōu)化點(diǎn)的偏移很小，這表示在設(shè)計工況點(diǎn)阻性換熱器室溫端面可保持在300 K左右。

圖11 stack-HEX的氮?dú)赓|(zhì)量流量和壓降隨電流的變化Fig.11 Mass flow rate and pressure drop of nitrogen gas vs current for stack-HEX

4 兩種阻性換熱器對比及討論

4.1 換熱性能對比

圖12 兩種阻性換熱器的氮?dú)赓|(zhì)量流量隨電流的變化Fig.12 Mass flow rate of nitrogen gas for two HEXs vs current

圖12為fin-HEX和stack-HEX的氮?dú)赓|(zhì)量流量隨電流變化的情況，其中氮?dú)赓|(zhì)量流量(液氮自蒸發(fā)速率)表征從阻性換熱器流入過渡銅塊的漏熱量。由圖12可知，兩種阻性換熱器流入過渡銅塊的漏熱量基本相同。圖13為fin-HEX和stack-HEX與氮?dú)獾膶α鲹Q熱速率隨電流變化的情況。由圖13可知，這兩種阻性換熱器的對流換熱速率基本相同，均可滿足快速換熱的要求。因此，從對低溫系統(tǒng)的熱負(fù)荷和對流換熱速率來看，兩種阻性換熱器在換熱性能方面基本相同。

圖13 兩種阻性換熱器與氮?dú)獾膶α鲹Q熱速率隨電流的變化Fig.13 Heat exchange rate between two HEXs and nitrogen vapor vs current

比較兩種換熱器的幾何結(jié)構(gòu)可知，fin-HEX和氮?dú)庵g的有效換熱面積雖較小，但氮?dú)饬鲃訛橥牧髁鲃?，且氮?dú)庋亓鞯缆菪搅鲃?，流動和熱邊界層被破壞，因此換熱能力較強(qiáng)；而stack-HEX中氮?dú)饬鲃与m為層流流動，但由于換熱器和氮?dú)庵g的有效換熱面積較大，因此換熱能力也較強(qiáng)。

4.2 流動性能對比

圖14為fin-HEX和stack-HEX中氮?dú)鈮航惦S電流變化的情況，可看出，stack-HEX的流動阻力較fin-HEX的至少小1個數(shù)量級。圖15為fin-HEX和stack-HEX中氮?dú)庥辛烤V摩擦系數(shù)隨電流變化的情況，可看出，stack-HEX的有量綱摩擦系數(shù)同樣較fin-HEX的至少小1個數(shù)量級，fin-HEX的有量綱摩擦系數(shù)隨電流上升，但stack-HEX的有量綱摩擦系數(shù)隨電流下降?？赏茢啵?dāng)電流繼續(xù)增大時，fin-HEX的流動阻力會迅速增大，而stack-HEX的流動阻力的增大較緩慢。從氮?dú)饬鲃幼枇湍Σ料禂?shù)來看，stack-HEX在流動性能方面明顯好于fin-HEX。

圖14 兩種阻性換熱器中氮?dú)鈮航惦S電流的變化Fig.14 Pressure drop of nitrogen gas for two HEXs vs current

圖15 兩種阻性換熱器中有量綱摩擦系數(shù)隨電流的變化Fig.15 Dimensional friction coefficient of two HEXs vs current

比較兩種換熱器的幾何結(jié)構(gòu)可知，fin-HEX中流道為螺旋形式，氮?dú)獾牧鲃影殡S著較為強(qiáng)烈的渦旋和湍流，因此流動阻力較大；而stack-HEX中流道為薄板形式，氮?dú)饬鲃訛閷恿髁鲃?，因此流動阻力較小。

4.3 運(yùn)行對比

EAST的13對HTS電流引線分為5對和8對兩組，分別安裝在兩個真空容器(電流引線罐)中。圖16為二元電流引線在電流引線罐中布置位置和冷卻結(jié)構(gòu)示意圖，電流引線的冷卻流程包括超導(dǎo)部分的氦冷卻回路和阻性換熱器的氮冷卻回路。換熱器下部的液氮筒與電流引線罐中的大液氮槽相連，液氮從大液氮槽中進(jìn)入液氮筒，液氮?dú)饣a(chǎn)生的氮蒸氣流經(jīng)阻性換熱器表面后從回氣口流出，形成氮冷卻回路。電流引線通電運(yùn)行過程中，若液氮筒中壓力較低，氮?dú)饬鲃訒霈F(xiàn)阻塞或部分阻塞的情況，從而導(dǎo)致阻性換熱器上部溫度急劇上升，產(chǎn)生過熱甚至燒毀的嚴(yán)重后果；若液氮筒中壓力過大，則會出現(xiàn)液氮隨氣流噴出的情況，從而導(dǎo)致電流引線室溫端的損壞和液氮消耗過快。為保證電流引線的正常運(yùn)行，有必要平衡和控制氮冷卻回路中的壓力。

圖16 電流引線在電流引線罐中布置位置和冷卻結(jié)構(gòu)示意圖Fig.16 Schematic diagram of location and cooling loop for current leads in current lead tank

阻性換熱器中過大的流動阻力和劇烈的流動阻力變化均會給壓力的平衡和控制帶來不利影響：過大的流動阻力使氮冷卻回路中產(chǎn)生較大的壓力不平衡，需調(diào)節(jié)的壓力范圍大；劇烈的流動阻力變化使氮冷卻回路中壓力調(diào)節(jié)的難度大。stack-HEX中氮?dú)獾牧鲃幼枇白兓^fin-HEX的小，因此其氮冷卻回路中壓力的平衡和控制較為容易。

從運(yùn)行結(jié)果看，使用stack-HEX的電流引線室溫端的平均溫度為280.0 K，較靠近正常室溫；而使用fin-HEX的電流引線室溫端的平均溫度為261.7 K，偏離正常室溫較遠(yuǎn)；stack-HEX中氮冷卻回路中壓力的平衡和控制較fin-HEX中的簡單，前者比后者更易調(diào)節(jié)到設(shè)計工況。因此，對于EAST阻性換熱器的設(shè)計和優(yōu)化，不僅要求其有足夠的換熱能力，還要盡量降低冷卻氣體的流動阻力。

5 結(jié)論

由EAST阻性換熱器的熱工水力性能的模擬計算結(jié)果可知，兩種阻性換熱器在換熱性能方面基本相當(dāng)，均可滿足快速換熱的要求；stack-HEX在流動性能方面明顯好于fin-HEX，即stack-HEX中的流動阻力遠(yuǎn)小于fin-HEX中的。fin-HEX中氮?dú)獾牧鲃訛閺?qiáng)烈的湍流流動，流動阻力大，這使氮冷卻回路中壓力平衡較難實(shí)現(xiàn)，從而使電流引線的運(yùn)行偏離設(shè)計工況較遠(yuǎn)，不利于電流引線的長期安全運(yùn)行。stack-HEX中氮?dú)獾牧鲃訛閷恿髁鲃?，流動阻力小，氮冷卻回路中壓力平衡易實(shí)現(xiàn)，因而電流引線運(yùn)行時偏離設(shè)計工況程度較小。從EAST阻性換熱器的冷卻回路壓力平衡控制和運(yùn)行安全性考慮，疊片式結(jié)構(gòu)較三頭螺旋翅片式結(jié)構(gòu)更適合應(yīng)用在EAST阻性換熱器中。

[1] WU S, EAST Machine Team. EAST machine assembly and its measurement system[J]. Fusion Engineering and Design, 2005, 73: 135-140.

[2] WU S, EAST Team. An overview of the EAST project[J]. Fusion Engineering and Design, 2007, 82: 463-471.

[3] 畢延芳，馬登奎. 高溫超導(dǎo)電流引線在EAST裝置工程調(diào)試中投運(yùn)[J]. 低溫物理學(xué)報，2007，23(1)：131-135.

BI Yanfang, MA Dengkui. Operation of HTS current leads for EAST superconducting Tokamak commissioning[J]. Chinese Journal of Low Temperature Physics, 2007, 23(1): 131-135(in Chinese).

[4] 畢延芳. EAST超導(dǎo)托克馬克磁體的大型供電供冷外饋線[J]. 低溫物理學(xué)報，2009，31(4):333-337.

BI Yanfang. Large-scale outer feeders for SC magnets of EAST Tokamak[J]. Chinese Journal of Low Temperature Physics, 2009, 31(4): 333-337(in Chinese).

[5] DING K, BI Y, FENG H, et al. Development of ITER HTS current lead at ASIPP[J]. Physics Procedia, 2012, 36: 931-936.

[6] 周挺志. ITER特大電流高溫超導(dǎo)電流引線設(shè)計與技術(shù)研究[D]. 合肥：中國科學(xué)院合肥物質(zhì)科學(xué)研究院，2011.

[7] 畢延芳. EAST裝置15 kA高溫超導(dǎo)電流引線研發(fā)[J]. 低溫物理學(xué)報，2005，27(5)：1 074-1 079.

BI Yanfang. R&D on 15 kA HTS current leads for EAST Tokamak[J]. Chinese Journal of Low Temperature Physics, 2005, 27(5): 1 074-1 079(in Chinese).

[8] 朱銀鋒，武松濤，畢延芳，等. EAST縱場磁體系統(tǒng)16 kA高溫超導(dǎo)電流引線的研究[J]. 核聚變與等離子體物理，2007，27(2)：101-105.

ZHU Yinfeng, WU Songtao, BI Yanfang, et al. Study of 16 kA high temperature superconducting current leads used in EAST’s Toroidal field magnets system[J]. Chinese Journal of Nuclear Fusion and Plasma Physics, 2007, 27(2): 101-105(in Chinese).

[9] MAEHATA K, KAWASAKI S, ISHIBASHI K, et al. Operational performance of spiral-fin current leads[J]. Cryogenics, 1993, 33(7): 680-685.

Multi-physics Coupling Simulation of Resistive Heat Exchanger in EAST

DENG Wei, XI Wei-bin

(InstituteofPlasmaPhysics,ChineseAcademyofSciences,Hefei230031,China)

Two types of resistive heat exchangers which are cooled with nitrogen vapor are designed for the HTS current leads in EAST. The study of the thermal-hydraulic performance with multi-physics coupling simulation for the two types of resistive heat exchangers was presented in this paper. The simulation results show that the two resistive heat exchangers are almost the same in heat transfer performance, but the foil stack heat exchanger is much better than three-helical fin heat exchanger in flow performance with a low flow resistance. Considering the large pressure drop of nitrogen vapor in the three-helical fin heat exchanger, which can cause troubles for the control of pressure in the nitrogen cooling loop, the foil stack is a better choice than three-helical fin for the resistive heat exchanger in EAST HTS current leads.

HTS current lead; resistive heat exchanger; thermal-hydraulic performance; multi-physics coupling

2014-03-18;

2014-08-25

鄧 威(1989—)，男，湖北孝感人，碩士研究生，核能科學(xué)與工程專業(yè)

TM89

A

1000-6931(2015)07-1330-08

10.7538/yzk.2015.49.07.1330