江晟昊 袁艷平 曹曉玲

GDC電極冷卻流道系統(tǒng)結(jié)構(gòu)優(yōu)化及適用范圍分析

江晟昊 袁艷平 曹曉玲

（西南交通大學(xué)機(jī)械工程學(xué)院 成都 610031）

利用數(shù)值模擬的方法研究了各區(qū)域流道結(jié)構(gòu)尺寸對GDC電極冷卻流道系統(tǒng)傳熱效果的影響。分析結(jié)果表明減小第一區(qū)域冷卻流道厚度及寬度，增大電極體最外壁頂邊的倒圓角半徑及第二、三、四區(qū)域流道管徑均可降低GDC電極的最大溫度及平均溫度。綜合這些影響因素對冷卻流道系統(tǒng)結(jié)構(gòu)進(jìn)行了優(yōu)化設(shè)計(jì)，在保證足夠安全裕量的前提下有效降低了GDC電極的最大溫度及平均溫度。針對優(yōu)化后的GDC電極冷卻流道設(shè)計(jì)，研究了GDC電極在不同工況下所受熱負(fù)荷及冷卻水入口溫度大小對電極體各限制參數(shù)的影響，并確定了GDC電極冷卻流道系統(tǒng)的適用范圍。

GDC電極；溫度場；數(shù)值模擬；壓降

0 引言

輝光放電（GDC）電極作為托卡馬克裝置項(xiàng)目的子項(xiàng)目，旨在清除第一壁表面的核燃料粒子和其它雜質(zhì)粒子，為高約束性能等離子體提供清潔的和低再循環(huán)的第一壁[1]。GDC電極運(yùn)行期間會經(jīng)歷兩個時(shí)期，等離子體操作期及烘烤期。在等離子體操作期，GDC電極最外層壁面會作為托卡馬克裝置內(nèi)壁面的一部分，并接收到大量的熱負(fù)荷。

而烘烤期，GDC電極會伸入托卡馬克裝置的內(nèi)腔進(jìn)行輝光放電清除雜質(zhì)，此時(shí)接收到的熱負(fù)荷較小。由于GDC電極在運(yùn)行期間承受了較大的熱負(fù)荷，為了保障電極體的穩(wěn)定運(yùn)行，在電極內(nèi)需要設(shè)計(jì)冷卻流道并通過循環(huán)冷卻水來達(dá)到降溫的目的[2]。Lijun Cai、Hongbing Xu等對GDC電極進(jìn)行過冷卻流道設(shè)計(jì)與數(shù)值計(jì)算[3,4]，但設(shè)計(jì)的主要目的是滿足托卡馬克裝置對各設(shè)計(jì)參數(shù)的要求，缺少結(jié)構(gòu)參數(shù)對冷卻流道系統(tǒng)散熱效果影響情況的研究，無法保障GDC電極在部分極端工況下的穩(wěn)定運(yùn)行。

本文在GDC電極傳熱數(shù)值研究的基礎(chǔ)上，建立了合適的幾何、數(shù)學(xué)模型，使用Fluent軟件進(jìn)行數(shù)值模擬。通過對比較大范圍冷卻水入口流量條件下不同冷卻流道結(jié)構(gòu)參數(shù)對GDC電極最大溫度、進(jìn)出口壓降、平均流速等參數(shù)的影響，找到最優(yōu)的設(shè)計(jì)參數(shù)。針對優(yōu)化后的流道設(shè)計(jì)，通過增大GDC電極在等離子體操作期承受的熱負(fù)荷，研究冷卻流道系統(tǒng)可承受的熱負(fù)荷范圍，進(jìn)一步保障GDC電極在部分極端工況下的穩(wěn)定性。

1 模型及參數(shù)確定

1.1 幾何模型

GDC電極冷卻流道系統(tǒng)模型如圖1所示，采用了根據(jù)熱負(fù)荷大小不同冷卻流道形式與密度相應(yīng)變化的多區(qū)域流道設(shè)計(jì)。

圖1 冷卻流道系統(tǒng)結(jié)構(gòu)示意圖

擬變動的結(jié)構(gòu)參數(shù)為第一區(qū)域流道厚度及寬度，電極體最外壁頂邊的倒圓角半徑及第二、三、四區(qū)域流道管徑，具體取值見表1所示。

表1 冷卻流道系統(tǒng)各參數(shù)值

1.2 網(wǎng)格及湍流模型

利用Workbench中的Design Modeler對GDC電極冷卻流道系統(tǒng)進(jìn)行模型建立，并利用Meshing進(jìn)行網(wǎng)格劃分，采用四面體網(wǎng)格[5]，網(wǎng)格總數(shù)控制在450～550萬。

使用Fluent模擬計(jì)算時(shí)，設(shè)置進(jìn)口邊界條件為質(zhì)流量及冷卻劑入口溫度；出口條件為定壓；設(shè)置壁面條件為非滑移壁面。湍流模型為模型。

1.3 運(yùn)行工況及設(shè)計(jì)參數(shù)

GDC電極主要會經(jīng)歷兩個主要時(shí)期，分別是等離子體操作期和GDC烘烤期。托卡馬克裝置文件Summary of GDC Heat Load[6]給出了兩個時(shí)期GDC電極的熱負(fù)荷參數(shù)。

對于等離子體操作期，由于GDC電極頭等離子體表面面向托卡馬克裝置內(nèi)腔，所以會接收到輻射和電荷交換熱以及中子熱。首先，輻射和電荷交換熱為0.35MW/m2，熱負(fù)荷加載在電極頭的等離子體表面上。其次，中子熱是延徑向的體積熱負(fù)荷的指數(shù)衰減，具體公式如下：

根據(jù)托卡馬克裝置要求[7]，為了保證電極體的安全穩(wěn)定運(yùn)行，托卡馬克裝置要求電極體最高溫度不能超過400℃，冷卻液平均流速為0.5～5m/s。在等離子體操作期，冷卻水進(jìn)出口壓力降小于1.35MPa，冷卻水進(jìn)出口溫差小于126±10℃。在GDC烘烤期，出于安全考慮冷卻水溫度小于4.4MPa下水的沸騰溫度250.4℃。兩個不同工作時(shí)期冷卻水運(yùn)行工況如表2所示。

表2 不同時(shí)期冷卻液運(yùn)行工況

2 結(jié)構(gòu)參數(shù)對冷卻流道系統(tǒng)性能的影響

為了保證冷卻流道系統(tǒng)其余參數(shù)均滿足托卡馬克裝置要求并具有足夠的安全裕量，取托卡馬克裝置對冷卻水進(jìn)出口壓降、進(jìn)出口溫差、平均流速等參數(shù)限制條件的80%作為參數(shù)上限，則冷卻水進(jìn)出口壓降需不超過1.08MPa，平均流速需不超過4m/s。

2.1 第一區(qū)域流道厚度

如圖2所示，隨著第一區(qū)域流道厚度減小，電極體最大溫度逐漸降低，并且下降幅度較為均勻，流道厚度每減小1mm電極體最大溫度約平均下降2℃。這是由于減小流道厚度則減小了流道橫截面積，在流量不變的情況下增大了流道內(nèi)冷卻水流速，增大了換熱系數(shù)。為了保證電極體最大溫度不超過400℃，當(dāng)流道厚度不小于13mm時(shí)，需保證入口流量不小于0.5kg/s，當(dāng)流道厚度小于13mm時(shí)，在全入口流量范圍內(nèi)都能滿足托卡馬克裝置要求。

圖2 不同厚度下最大溫度隨入口流量變化

圖3 不同厚度下進(jìn)出口壓降隨入口流量的變化

如圖3所示，隨著第一區(qū)域流道厚度減小，冷卻水進(jìn)出口壓降逐漸增大，流道厚度每減小1mm冷卻水進(jìn)出口壓降平均上升約0.1MPa。為了保證冷卻水進(jìn)出口壓降不超過1.08MPa，當(dāng)?shù)谝粎^(qū)域流道厚度為14mm時(shí)，需保證入口流量不大于1.2kg/s，當(dāng)?shù)谝粎^(qū)域流道厚度為12mm～13mm時(shí)，需保證入口流量不大于1.1kg/s，當(dāng)?shù)谝粎^(qū)域流道厚度為11mm時(shí)，需保證入口流量不大于1.0kg/s。

如圖4所示，隨著第一區(qū)域流道厚度減小，冷卻水平均流速逐漸增大，流道厚度每減小1mm冷卻水平均流速增大約0.1m/s。為了保證冷卻水平均流速不超過4m/s，當(dāng)流道厚度為14mm時(shí)，需保證入口流量不大于1kg/s，當(dāng)流道厚度為11mm～13mm時(shí)，需保證入口流量不大于0.9kg/s。

圖4 不同厚度下冷卻水平均流速隨入口流量變化

綜合上述分析結(jié)果，在可運(yùn)行范圍內(nèi)，第一區(qū)域流道厚度為11mm時(shí)電極體最大溫度具有最小值363.5℃，該結(jié)構(gòu)參數(shù)下冷卻流道系統(tǒng)設(shè)計(jì)具有更好的散熱效果。

2.2 第一區(qū)域流道寬度

圖5 不同寬度下最大溫度隨入口流量變化

如圖5所示，隨著第一區(qū)域流道寬度減小，電極體最大溫度逐漸降低，流道寬度由17mm減小至16mm時(shí)，溫度下降幅度較小，流寬度由16mm減小至14mm時(shí)，寬度每減小1mm電極體最大溫度下降幅度有所增大。為了保證電極體最大溫度不超過400℃，當(dāng)?shù)谝粎^(qū)域流道寬度不小于16mm時(shí)，需保證入口流量不小于0.5kg/s，當(dāng)寬度小于16mm時(shí)，在0.4kg/s～1.2kg/s范圍內(nèi)均可滿足托卡馬克裝置要求。

如圖6所示，隨著第一區(qū)域流道寬度減小，冷卻水進(jìn)出口壓降逐漸增大。為了保證冷卻水進(jìn)出口壓降小于1.08MPa，當(dāng)流道寬度為14mm時(shí)，需保證入口流量不大于1.1kg/s，當(dāng)流道寬度大于14mm時(shí)，在0.4kg/s～1.2kg/s范圍內(nèi)均可滿足限制條件。

圖6 不同寬度下進(jìn)出口壓降隨入口流量的變化

如圖7所示，隨著第一區(qū)域流道寬度減小，冷卻水平均流速逐漸增大。為了保證冷卻水平均流速不超過4.0m/s，當(dāng)流道寬度不小于16mm時(shí)，需保證入口流量不大于1.0kg/s，當(dāng)流道寬度小于16mm時(shí)，需保證入口流量不大于0.9kg/s。

圖7 不同寬度下冷卻水平均流速隨入口流量變化

綜合上述分析，在可運(yùn)行范圍內(nèi)，第一區(qū)域流道寬度為14mm時(shí)電極體最大溫度具有最小值360℃，該結(jié)構(gòu)參數(shù)下冷卻流道系統(tǒng)設(shè)計(jì)具有更好的散熱效果。

2.3 最外層壁面頂邊倒角半徑

如圖8所示，隨著電極頭最外層壁面頂邊倒角半徑的增大，電極體最大溫度顯著降低。倒角半徑由6mm增大至7mm時(shí)電極體最大溫度降低幅度相對較小，半徑由7mm增大至9mm時(shí)，每增大1mm半徑電極體最大溫度降低幅度有所增大。為了保證電極體最大溫度小于400℃，當(dāng)?shù)菇前霃讲淮笥?mm時(shí)，需保證入口流量不小于0.5kg/s，當(dāng)?shù)菇前霃酱笥?mm時(shí)，在0.4kg/s～1.2kg/s流量范圍內(nèi)均可滿足要求。

圖8 不同半徑下最大溫度隨入口流量的變化

綜合上述分析，在可運(yùn)行范圍內(nèi)，當(dāng)?shù)菇前霃綖?mm時(shí)，電極體最大溫度具有最小值302.5℃，該結(jié)構(gòu)參數(shù)下冷卻流道系統(tǒng)設(shè)計(jì)具有更好的散熱效果。

2.4 第二、三、四區(qū)域流道管徑

圖9 不同管徑下最大溫度隨入口流量的變化

如圖9所示，隨著第二、三、四區(qū)域流道管徑減小，電極體最大溫度逐漸減小，但下降幅度較小，在0.4kg/s～1.2kg/s入口流量范圍內(nèi)，不同管徑的冷卻流道間電極體最大溫度最大差值僅為2.2℃。為了保證電極體最大溫度小于400℃，四種管徑下冷卻流道系統(tǒng)入口流量需保證不小于0.5kg/s。

如圖10所示，隨著第二、三、四區(qū)域流道管徑減小，冷卻水進(jìn)出口壓降逐漸增大。為了保證冷卻水進(jìn)出口壓降小于1.08MPa，當(dāng)管徑為17mm時(shí)，需保證入口流量不超過1.0kg/s，當(dāng)管徑大于17mm時(shí)，在0.4kg/s～1.2kg/s范圍內(nèi)均可滿足要求。

圖10 不同管徑下進(jìn)出口壓降隨入口流量的變化

如圖11所示，隨著第二、三、四區(qū)域流道管徑減小，冷卻水平均流速逐漸增大。為了保證冷卻水平均流速小于4m/s，當(dāng)管徑為20mm時(shí)，需保證入口流量不大于1.1kg/s，當(dāng)管徑為18mm～19mm時(shí)，需保證入口流速不大于1.0kg/s，當(dāng)管徑為17mm時(shí)，需保證入口流速不大于0.9kg/s。

圖11 不同管徑冷卻水平均流速隨入口流量變化

綜合綜合上述分析，在可運(yùn)行范圍內(nèi)，當(dāng)?shù)沟诙?、三、四排流道管徑?0mm時(shí)，電極體最大溫度具有最小值363.8℃。該結(jié)構(gòu)參數(shù)下冷卻流道系統(tǒng)設(shè)計(jì)具有更好的散熱效果。

2.5 優(yōu)化設(shè)計(jì)

通過前文的影響因素分析可知，減小第一區(qū)域流道厚度、寬度，第二、三、四區(qū)域流道管徑，增大電極體最外層壁面頂邊倒圓角半徑均可改善電極體最大溫度及平均溫度。根據(jù)前文的分析結(jié)果，選取第一區(qū)域流道厚度為11mm，寬度為14mm，倒角半徑為9mm，第二、三、四區(qū)域流道管徑為18mm，對冷卻流道系統(tǒng)結(jié)構(gòu)進(jìn)行優(yōu)化。在可運(yùn)行范圍內(nèi)，取入口流量為0.9kg/s時(shí)，在滿足所有設(shè)計(jì)參數(shù)要求的前提下GDC電極最大溫度具有最小值299℃。

3 適用范圍分析

3.1 中子熱對冷卻流道系統(tǒng)的影響

圖12所示為不同大小的中子熱影響下，GDC電極最大溫度的變化情況。隨著熱負(fù)荷的增大，GDC電極的最大溫度在熱流密度增大到12MW/m3前增長較緩，每增大1MW/m3熱流密度GDC電極最大溫度僅升高約3℃，在熱流密度超過12MW/m3以后，電極體最大溫度的升高幅度增加，每增大1MW/m3熱流密度GDC電極最大溫度升高約19℃。當(dāng)熱流密度超過12MW/m3時(shí)電極桿區(qū)域的GDC電極最大溫度便超過了最外層壁面，最終當(dāng)熱流密度達(dá)到17MW/m3時(shí)達(dá)到398.8℃，非常接近托卡馬克裝置對于電極體最大溫度不超過400℃的要求。

圖12 不同大小中子熱對GDC電極最大溫度影響

3.2 輻射和電荷交換熱對冷卻流道系統(tǒng)的影響

圖13所示為不同大小的輻射和電荷交換熱影響下，GDC電極最大溫度的變化情況。隨著輻射和電荷交換熱的增大，GDC電極最大溫度呈線性增大，熱流密度每增加0.035MW/m2，電極體最大溫度增加約21℃。當(dāng)輻射和電荷交換熱的熱流密度達(dá)到0.525MW/m2時(shí)，GDC電極最大溫度達(dá)到406℃，超過托卡馬克裝置對于電極體最大溫度的要求。

圖13 不同大小的輻射和電荷交換熱對GDC電極最大溫度的影響

4 結(jié)論

（1）通過改變GDC電極冷卻流道結(jié)構(gòu)參數(shù)，對GDC電極進(jìn)行模擬分析發(fā)現(xiàn)，減小第一區(qū)域冷卻流道厚度及寬度，增大電極體最外壁頂邊的倒圓角半徑，第二、三、四區(qū)域流道管徑均可降低GDC電極的最大溫度及平均溫度。減小第一區(qū)域冷卻流道厚度、寬度及后排區(qū)域管徑會增大冷卻液進(jìn)出口壓力降及平均流速。而增大電極體最外壁頂邊倒圓角半徑對這幾類參數(shù)影響均很小。以對GDC電極最大溫度的改善效果為評價(jià)指標(biāo)，選取第一區(qū)域流道厚度為11mm，寬度為14mm，頂邊倒角半徑為9mm，后排流道管徑為18mm作為研究范圍內(nèi)的最優(yōu)流道結(jié)構(gòu)尺寸。優(yōu)化后冷卻流道設(shè)計(jì)具有更好的冷卻效果，在可運(yùn)行范圍內(nèi)取入口流量為0.9kg/s，電極體最大溫度為299℃。

（2）通過增大等離子體操作期GDC電極承受的熱負(fù)荷，對優(yōu)化后的GDC電極進(jìn)行模擬分析發(fā)現(xiàn)，在保持輻射和電荷交換熱的熱負(fù)荷大小不變的情況下，GDC電極冷卻流道系統(tǒng)最大可承受中子熱的熱流密度為17MW/m3，電極體最大溫而在保持中子熱熱負(fù)荷大小不變的情況下，GDC電極冷卻流道系統(tǒng)最大可承受輻射和電荷交換熱的熱流密度為0.525MW/m2。

[1] 畢海林,胡建生,余耀偉,等.HT-7托卡馬克全金屬壁及鋰化條件下輝光放電清洗的研究[J].真空科學(xué)與技術(shù)學(xué)報(bào),2014,34(7):731-736.

[2] 曾建爾,鄒桂清,冉紅.ITER輝光放電電極的冷卻計(jì)算[J].核聚變與等離子體物理,2007,(3):203-207.

[3] Lijun Cai, Tao Lin, Yingqiao Wang, et al. Preliminary fluid channel design and thermal-hydraulic analysis of glow discharge cleaning permanent electrode[J]. Fusion Engineering and Design, 2016:109-111.

[4] Hongbing Xu, Yinglong Yuan, Yong Lu, et al. The preliminary design of cooling channel for ITER GDC[J]. Fusion Engineering and Design, 2017:121.

[5] 張秀杰,諶繼明,康偉山,等.ITER屏蔽包層屏蔽塊熱工水力分析[J].核聚變與等離子體物理,2010,30(3):267- 271.

[6] Y Yang. Heat load of GDC electrode during different phases and operation states[M]. ITER, 2016.

[7] V Udintsev. System Design Description Document (DDD) Diagnostic First Wall DFW[R]. ITER IDM Report, 2014.

Structural Optimization and Application Range Analysis of GDC Electrode Cooling Runner System

Jiang Chenghao Yuan Yanping Cao Xiaoling

( College of Mechanical Engineering, South west Southwest Jiaotong University of china, Chengdu, 610031 )

In this paper, the effect of structural size of each channel on the heat transfer effect of GDC electrode cooling runner system is studied by numerical simulation. The analysis results show that reducing the thickness and width of the cooling channel in the first region, increasing the rounding radius of the top edge of the outermost wall of the electrode body and the diameter of the second, third and fourth region flow channels can reduce the maximum temperature and average temperature of the GDC electrode. Based on these influencing factors, the structure of the cooling runner system is optimized, and the maximum temperature and average temperature of the GDC electrode are effectively reduced under the premise of ensuring sufficient safety margin. According to the optimized GDC electrode cooling flow channel design, the effects of the heat load and cooling water inlet temperature on the limiting parameters of the GDC electrode under different working conditions were studied, and the applicable range of the GDC electrode cooling channel system was determined.

GDC electrode; Temperature field; Numerical simulation; Pressure drop

1671-6612（2019）06-629-06

TL639

A

江晟昊（1993-），男，碩士研究生，E-mail：chenghaojiang@my.swjtu.edu.cn

袁艷平（1973-），男，博士后，博士生導(dǎo)師，E-mail：ypyuan@home.swjtu.edu.cn

2019-05-15