李 臻，陸道綱，曹 瓊,*

(1.華北電力大學(xué) 核科學(xué)與工程學(xué)院，北京 102206；2.非能動核能安全技術(shù)北京市重點實驗室，北京 102206)

隨著未來深空探測以及星表基地供電等任務(wù)需求的日益提高以及太陽能、化學(xué)能供電的局限性，結(jié)構(gòu)緊湊、功率密度高、壽命較長、環(huán)境適應(yīng)性強(qiáng)的空間核反應(yīng)堆電源已成為太空能源應(yīng)用的主要發(fā)展方向。前蘇聯(lián)TOPAZ-Ⅱ是目前投入制造的最先進(jìn)、成熟的空間核反應(yīng)堆電源[1-3]。

輻射屏蔽專門用于減弱輻照敏感裝置周圍輻射場，是空間核反應(yīng)堆電源的關(guān)鍵部件之一。輻射屏蔽的溫度場分析是空間核反應(yīng)堆電源設(shè)計中的必要環(huán)節(jié)。輻射屏蔽因自身溫度的不均勻分布會產(chǎn)生較大熱應(yīng)力，以溫度場分析為基礎(chǔ)的熱力學(xué)分析是反應(yīng)堆熱工設(shè)計中不可缺少的部分。從安全角度來看，在各種工況下，要求輻射屏蔽的溫度峰值不能過高，且溫度分布盡可能均勻以減小熱應(yīng)力，才能保證其長期工作的能力。因此，在技術(shù)可行的前提下，對輻射屏蔽采取適當(dāng)?shù)膬?yōu)化措施，控制熱阻、削弱換熱十分必要。

國內(nèi)外對輻射屏蔽在空間核反應(yīng)堆電源內(nèi)真空環(huán)境下的溫度分布研究資料幾乎沒有。馬曉等[4]利用CFX軟件對中國實驗快堆堆頂固定屏蔽冷卻系統(tǒng)進(jìn)行三維數(shù)值模擬，并對該系統(tǒng)提出了優(yōu)化建議。盛選禹等[5]使用Abaqus 6.7商業(yè)有限元軟件分析了高溫氣冷堆的屏蔽層在額定工況和事故工況下的溫度分布。何杰等[6]針對2 MW液態(tài)釷基熔鹽實驗堆將利用MCNP軟件獲得的反應(yīng)堆功率分布導(dǎo)入Fluent軟件，分析了主屏蔽的三維溫度場。

本文針對TOPAZ-Ⅱ型空間核反應(yīng)堆電源所采用的輻射屏蔽使用Fluent軟件對其在真空環(huán)境下的換熱行為進(jìn)行三維數(shù)值模擬，根據(jù)優(yōu)化措施，對輻射屏蔽方案溫度分布進(jìn)行計算，研究分析這些措施對輻射屏蔽溫度分布特性產(chǎn)生的影響，為相關(guān)實驗與工程設(shè)計提供參考。

1 模擬對象及優(yōu)化措施

TOPAZ-Ⅱ型空間核反應(yīng)堆電源結(jié)構(gòu)如圖1所示，TOPAZ-Ⅱ型空間核反應(yīng)堆電源主要由堆芯及熱離子轉(zhuǎn)換系統(tǒng)、電磁泵、輻射屏蔽、控制驅(qū)動單元、銫蒸氣單元和熱管輻射換熱器組成[7]。TOPAZ-Ⅱ空間核反應(yīng)堆電源輻射屏蔽采用“陰影屏蔽”，會在有限的空間內(nèi)造成輻射水平低的區(qū)域，這個區(qū)域具有截頭圓錐體的形狀。因此輻射屏蔽同樣采用截頭圓錐體的形式，并且是1個由各種不同功能的材料組成的多層結(jié)構(gòu)[8]。在空間核反應(yīng)堆電源系統(tǒng)輻射屏蔽的設(shè)計中，選擇屏蔽材料非常重要。對中子而言，LiH是非常優(yōu)秀的中子屏蔽材料；對γ射線而言，應(yīng)選用原子序數(shù)高的材料，如不銹鋼、鉛、W、B4C等[9]。輻射屏蔽結(jié)構(gòu)主要分為兩種類型:降低中子通量的輕屏蔽，是一充有LiH的鋼外殼;降低伽馬通量的重屏蔽，為厚度幾厘米的容器底[10]。

圖1 TOPAZ-Ⅱ型空間核反應(yīng)堆電源

真空環(huán)境內(nèi)稀薄空氣流速緩慢、換熱系數(shù)偏低，對流換熱作用較弱，輻射換熱占據(jù)主導(dǎo)地位。對于輻射屏蔽，輻射換熱主要發(fā)生在冷卻劑管道與輻射屏蔽之間。為削弱它們之間的輻射換熱，通常在冷卻劑管道及陰影屏蔽上冷卻劑管管槽間添加真空遮熱板、在材料表面覆蓋改變發(fā)射率的涂層[11]。

2 數(shù)值模擬

2.1 幾何模型

本文計算區(qū)域為圖1中虛線框部分。根據(jù)對輻射屏蔽傳熱過程的分析，將TOPAZ-Ⅱ空間核反應(yīng)堆電源的輻射屏蔽等效為輕、重屏蔽兩均勻整體，輕屏蔽為一高度0.6 m、上截面半徑0.38 m、下截面半徑0.5 m的截頭圓錐體，重屏蔽為一厚度0.02 m的圓板。并且為模擬輻射屏蔽在真空環(huán)境下的換熱行為，在屏蔽模型外部加上一高度2 m、底面半徑1 m的圓柱體外罩即劃分一模擬真空環(huán)境的腔室；為考慮反應(yīng)堆底部溫度對輻射屏蔽溫度分布的影響，在重屏蔽上部添加一薄板模擬反應(yīng)堆底部。以在冷卻劑管與冷卻劑管管槽間添加一厚度2 mm的真空遮熱板的輻射屏蔽模型為例，簡化處理后的結(jié)果如圖2所示。

圖2 簡化模型

2.2 網(wǎng)格劃分

本研究采用四面體非結(jié)構(gòu)網(wǎng)格，整體網(wǎng)格質(zhì)量較好，滿足計算精度要求。以在冷卻劑管與冷卻劑管管槽間添加一層真空遮熱板的輻射屏蔽模型為例，局部網(wǎng)格分布示于圖3。

3月21日，擔(dān)任銀監(jiān)會主席一年多的郭樹清成為剛剛成立的中國銀行保險業(yè)監(jiān)督管理委員會第一任主席。5天之后，中國人民銀行召開干部大會，中組部相關(guān)負(fù)責(zé)人宣布了他的另一項職務(wù)：郭樹清任央行黨委書記、副行長。

圖3 局部網(wǎng)格分布

2.3 控制方程

考慮到本文模擬真空環(huán)境下輻射屏蔽的換熱行為，選用Surface-to-Surface(S2S)輻射模型，該模型非常適用于無參與性介質(zhì)的封閉空間中的輻射問題[12]。S2S輻射模型的控制方程[12]為:

式中：qout,k、qout,j為表面k、j發(fā)出的輻射熱流密度；εk為物體發(fā)射率(黑度)；σ為斯忒藩-玻耳茲曼常量(黑體輻射常數(shù))，其值為5.67×10-9W/(m2·K4)；Tk為表面溫度；ρk為反射比；Fkj為角系數(shù)。

固體域能量控制方程[12]可表示為：

式中：k為熱導(dǎo)率；T為溫度；Sh為體熱源。

2.4 物性參數(shù)

除輕屏蔽的材料LiH和重屏蔽、反應(yīng)堆底部、真空遮熱板、冷卻劑管的材料不銹鋼外，數(shù)值計算中使用的材料還包括為模擬真空環(huán)境在腔室內(nèi)填充的稀薄空氣，這些材料的物性參數(shù)列于表1。

表1 材料的物性參數(shù)

2.5 邊界條件

邊界條件的設(shè)置包括：1) 根據(jù)文獻(xiàn)[10]，3根從堆芯帶出冷卻劑的冷卻劑管道表面溫度843 K，3根向堆芯帶入冷卻劑的冷卻劑管道表面溫度743 K；2) 保守地將模擬真空環(huán)境的腔室頂面設(shè)置絕熱邊界條件；3) 考慮到輻射屏蔽底部的堆內(nèi)構(gòu)件的溫度限值，將模擬真空環(huán)境的腔室底面溫度設(shè)置為473 K；4) 考慮到堆芯底部的氦氣腔傳熱較差，將反應(yīng)堆底部溫度設(shè)置為550 K；5) 輻射屏蔽俘獲中子、吸收γ射線后，會有釋熱，因此對屏蔽添加體熱源70 W/m3；6) 考慮到腔室為模擬真空環(huán)境，因此將真空腔室外罩頂面、底面及側(cè)壁發(fā)射率設(shè)為0.9；7) 輕屏蔽、重屏蔽及反應(yīng)堆底部表面常噴涂低發(fā)射率材料，設(shè)為0.3。

3 計算結(jié)果與分析

為得到輻射屏蔽的溫度分布狀況，在輻射屏蔽上選取一截面A，并在截面A上取直線B、C，分析直線B、C上溫度計算結(jié)果，截面A及直線B、C的位置如圖4所示。對于直線B、C，取箭頭方向為正方向，分別為x、y方向，其中直線B正方向為向堆芯帶入冷卻劑的冷卻劑管側(cè)向，并分別取截面中心(0，0)點、(0，0.04)點為原點。

圖4 截面A及直線B、C位置

3.1 網(wǎng)格無關(guān)性分析

圖5 網(wǎng)格無關(guān)性分析

3.2 環(huán)境溫度敏感性分析

輻射屏蔽溫度分布的數(shù)值計算依賴于其所處的環(huán)境是穩(wěn)定的，因此進(jìn)行環(huán)境溫度敏感性分析十分必要。為模擬在不同環(huán)境溫度下輻射屏蔽的換熱行為，將真空腔室的側(cè)壁溫度設(shè)置為303、313、323、333、343 K等5種工況，環(huán)境溫度敏感性分析結(jié)果如圖6所示。由圖6可見，當(dāng)輻射屏蔽所處的環(huán)境溫度在改變時，輻射屏蔽的溫度分布也會隨之改變，輻射屏蔽的溫度變化值與其所處環(huán)境的溫度變化值呈正比，且每種工況下溫度在x、y方向上的分布趨勢不變。

圖6 環(huán)境溫度敏感性分析

3.3 真空遮熱板影響

1) 真空遮熱板層數(shù)影響

為減少冷卻劑管道向輻射屏蔽的輻射傳熱，在冷卻劑管管槽的外側(cè)設(shè)置厚度為2 mm的真空遮熱板。在其他邊界條件相同的情況下，對無真空遮熱板、單層真空遮熱板、3層真空遮熱板3種工況進(jìn)行數(shù)值計算，結(jié)果如圖7、8所示。

圖7 不同真空遮熱板層數(shù)工況下輻射屏蔽溫度場

真空遮熱板層數(shù)分別為0、1、3時輻射屏蔽的溫度場如圖7所示。由圖7可見，隨著真空遮熱板層數(shù)的增加，輻射屏蔽的溫度下降明顯。并且3種工況下截面A上最低溫度均出現(xiàn)在調(diào)節(jié)機(jī)構(gòu)傳動管側(cè)，而最高溫度出現(xiàn)位置則隨真空遮熱板層數(shù)的增加逐漸由從堆芯帶出冷卻劑的冷卻劑管側(cè)向截面中心移動。

真空遮熱板層數(shù)對溫度分布的影響如圖8所示。由圖8可見，添加單層真空遮熱板與無真空遮熱板相比，輻射屏蔽溫度下降明顯且溫度分布更均勻；添加3層真空遮熱板與添加單層真空遮熱板相比，輻射屏蔽溫度依然下降但下降幅度極大減小，并且x方向上溫度分布均勻性上升程度極大減弱、y方向上溫度分布均勻性反而惡化，說明當(dāng)真空遮熱板層數(shù)超過1層，繼續(xù)增加真空遮熱板層數(shù)對優(yōu)化輻射屏蔽溫度分布已十分有限。因此，綜合考慮應(yīng)用場景尺寸限制、溫度分布優(yōu)化程度等因素，單層真空遮熱板為最佳選擇。

圖8 真空遮熱板層數(shù)對溫度分布的影響

2) 真空遮熱板表面發(fā)射率影響

真空遮熱板表面發(fā)射率對溫度分布的影響如圖9所示。圖9中，ε為真空遮熱板表面發(fā)射率。由圖9可見，輻射屏蔽的溫度隨著真空遮熱板表面發(fā)射率下降而降低，真空遮熱板表面發(fā)射率取0.01時，輻射屏蔽溫度最低，而且輻射屏蔽溫度的下降幅度隨著真空遮熱板表面發(fā)射率的下降而增大。并且隨著真空遮熱板表面發(fā)射率下降，輻射屏蔽溫度分布更加均勻，但在發(fā)射率為0.01時，考慮到此時反應(yīng)堆底部溫度對輻射屏蔽溫度分布的影響超過冷卻劑管道溫度的影響，x方向上溫度分布均勻性上升程度大為減弱、y方向上溫度分布均勻性反而惡化。因此在工程實踐中，應(yīng)盡可能降低真空遮熱板表面的表面發(fā)射率但不能太低，發(fā)射率為0.04～0.07最為適宜。

圖9 真空遮熱板表面發(fā)射率對溫度分布的影響

3.4 表面涂層影響

通過在冷卻劑管道及冷卻劑管管槽表面噴射不同發(fā)射率的涂層，可改變其表面發(fā)射率，起到改變輻射換熱效果的作用。

1) 冷卻劑管道表面發(fā)射率

冷卻劑管道表面發(fā)射率分別為0.1、0.5、0.9時輻射屏蔽的溫度場如圖10所示。由圖10可見，隨著冷卻劑管道表面發(fā)射率的降低，輻射屏蔽溫度下降明顯。截面A上最低溫度均出現(xiàn)在調(diào)動機(jī)構(gòu)傳動管側(cè)，而最高溫度出現(xiàn)位置則隨著冷卻劑管道表面發(fā)射率的降低逐漸由從堆芯帶出冷卻劑的冷卻劑管側(cè)向截面中心移動。冷卻劑管道表面發(fā)射率對溫度分布的影響如圖11所示。圖11中，ε為冷卻劑管道表面發(fā)射率。由圖10、11可見，輻射屏蔽溫度隨著冷卻劑管道表面發(fā)射率下降而降低，但輻射屏蔽溫度下降幅度隨著冷卻劑管道表面發(fā)射率下降而增大。且隨著冷卻劑管道表面發(fā)射率下降，輻射屏蔽溫度分布更均勻。但在發(fā)射率為0.3時，輻射屏蔽溫度分布的均勻性反而稍優(yōu)于發(fā)射率為0.1時的。因此，經(jīng)綜合考慮，冷卻劑管道表面發(fā)射率為0.1～0.3最適宜。

圖10 不同冷卻劑管道表面發(fā)射率工況下輻射屏蔽溫度場

圖11 冷卻劑管道表面發(fā)射率對溫度分布的影響

2) 冷卻劑管管槽表面發(fā)射率

冷卻劑管管槽表面發(fā)射率分別為0.1、0.5、0.9時輻射屏蔽的溫度場如圖12所示。由圖12可見，隨著冷卻劑管管槽表面發(fā)射率改變，截面A上最低溫度、最高溫度的位置變化趨勢與冷卻劑管道表面發(fā)射率改變時的位置變化趨勢相同。但對比圖10、12可見：分別降低冷卻劑管道、冷卻劑管管槽表面發(fā)射率至同一值的情況下，當(dāng)發(fā)射率較低時，前者最高溫度出現(xiàn)位置更加接近截面中心；當(dāng)發(fā)射率較高時，后者最高溫度出現(xiàn)位置更加接近截面中心。

圖12 不同冷卻劑管管槽表面發(fā)射率工況下輻射屏蔽溫度場

冷卻劑管管槽表面發(fā)射率對溫度分布的影響如圖13所示。圖中，ε為冷卻劑管管槽表面發(fā)射率。由圖13可見，隨著冷卻劑管管槽表面發(fā)射率下降，輻射屏蔽溫度下降且分布均勻性上升。因此冷卻劑管管槽表面發(fā)射率為0.1左右最為適宜。對比圖13、15可見：分別降低冷卻劑管道、冷卻劑管管槽表面發(fā)射率至同一值的情況下，當(dāng)發(fā)射率較低時，前者溫度更低且溫度分布均勻性更高；當(dāng)發(fā)射率較高時，后者溫度更低、溫度分布更加均勻?？傊档屠鋮s劑管道表面發(fā)射率較降低冷卻劑管管槽表面發(fā)射率能更好優(yōu)化輻射屏蔽溫度分布。

圖13 冷卻劑管管槽表面發(fā)射率對溫度分布的影響

4 結(jié)論

本文針對TOPAZ-Ⅱ型空間核反應(yīng)堆電源的輻射屏蔽在真空環(huán)境下的輻射換熱行為，對優(yōu)化輻射屏蔽的溫度分布所采取了相應(yīng)措施：在冷卻劑管及其管槽間設(shè)置真空遮熱板、在冷卻劑管道表面及冷卻劑管管槽表面噴涂合適發(fā)射率的涂層。利用Fluent軟件進(jìn)行模擬計算，分析了對輻射屏蔽溫度分布產(chǎn)生的影響，得出如下結(jié)論。

1) 在真空環(huán)境條件下，來自冷卻劑管道的輻射傳熱是影響輻射屏蔽溫度分布的主要因素；添加真空遮熱板和噴涂改變發(fā)射率的涂層均可有效降低輻射屏蔽溫度峰值，并提升其溫度分布均勻性。

2) 添加真空遮熱板優(yōu)化輻射屏蔽溫度分布的效果主要受真空遮熱板層數(shù)和其表面發(fā)射率影響，隔熱效果與真空遮熱板層數(shù)呈正比，與表面發(fā)射率呈反比；在工程實踐中，單層真空遮熱板、噴涂表面發(fā)射率為0.04～0.07的涂層為最佳選擇。

3) 冷卻劑管道表面、冷卻劑管管槽表面噴涂改變發(fā)射率的涂層對輻射屏蔽溫度分布的優(yōu)化效果與發(fā)射率呈反比，且冷卻劑管道表面發(fā)射率對輻射屏蔽溫度分布影響更顯著；在工程實踐中，應(yīng)盡可能降低這兩個表面的表面發(fā)射率，特別是冷卻劑管道表面，冷卻劑管管槽表面發(fā)射率為0.1左右、冷卻劑管道表面發(fā)射率為0.1～0.3最適宜。

本文實現(xiàn)了TOPAZ-Ⅱ型空間核反應(yīng)堆電源的輻射屏蔽結(jié)構(gòu)的數(shù)值模擬，為相關(guān)實驗和工程實踐中輻射屏蔽結(jié)構(gòu)的優(yōu)化提供了參考。