乏燃料運(yùn)輸容器傳熱支撐板對(duì)卸料過(guò)程的影響分析

殷 勇, 李其朋, 馬慶俊

(中廣核研究院有限公司， 廣東深圳 518124)

核電技術(shù)

乏燃料運(yùn)輸容器傳熱支撐板對(duì)卸料過(guò)程的影響分析

殷 勇, 李其朋, 馬慶俊

(中廣核研究院有限公司， 廣東深圳 518124)

針對(duì)自主設(shè)計(jì)的乏燃料運(yùn)輸容器卸料冷卻過(guò)程，開(kāi)展卸料冷卻過(guò)程熱工分析，研究傳熱支撐板對(duì)卸料冷卻過(guò)程的影響。分別建立不包含傳熱支撐板和包含傳熱支撐板的幾何模型，通過(guò)對(duì)兩種幾何模型進(jìn)行簡(jiǎn)化并取1/4為研究對(duì)象，采用k-ε方程湍流模型，利用Fluent完成熱工計(jì)算。分析結(jié)果表明：傳熱支撐板對(duì)卸料冷卻總時(shí)間的影響達(dá)到20%，傳熱支撐板對(duì)乏燃料運(yùn)輸容器卸料熱工分析的影響不能忽略。

乏燃料； 運(yùn)輸容器； 傳熱支撐板； 卸料冷卻； 熱工分析

Abstract： Thermal analysis was conducted on the unloading and cooling process of a self-developed spent fuel transport container with supporting slab, so as to analyze the effects of the slab on heat transfer behavior of the container during unloading and cooling process. Based on models set up respectively for the container with and without supporting slab, taking one-fourth of both the simplified models as an object of study, thermal calculation was carried out usingk-εturbulent model with Fluent software. Results show that the cooling time affected by the supporting slab accounts for 20% of the total during unloading and cooling process, which can not be neglected during thermal analysis of spent fuel transport containers.

Keywords： spent fuel; transport container; heat transfer supporting slab; unloading and cooling; thermal analysis

乏燃料運(yùn)輸容器輔助設(shè)備是用于對(duì)乏燃料運(yùn)輸容器進(jìn)行充水排氣、充氣排水、充氣風(fēng)干的設(shè)備，并通過(guò)對(duì)已裝載乏燃料組件的容器進(jìn)行充水冷卻，實(shí)現(xiàn)乏燃料組件的安全裝卸。其中的卸料冷卻模塊用于排除容器內(nèi)乏燃料的余熱，冷卻過(guò)程包括充水排氣階段和熱水置換階段。充水排氣階段向容器緩慢充入冷卻水，降低乏燃料組件表面溫度。如果此時(shí)的出口水溫仍大于允許水溫，則進(jìn)入熱水置換階段，繼續(xù)向容器內(nèi)連續(xù)注入冷卻水進(jìn)行置換，直至出口水溫低于允許水溫。由于裝載有乏燃料組件的運(yùn)輸容器具有放射性，需要保證卸料冷卻模塊能夠盡可能快速地完成冷卻功能。

筆者采用大型流體計(jì)算軟件Fluent，通過(guò)對(duì)自主設(shè)計(jì)的乏燃料運(yùn)輸容器卸料冷卻過(guò)程進(jìn)行熱工分析，分析容器傳熱支撐板對(duì)卸料冷卻過(guò)程的影響，確定卸料冷卻時(shí)間以及流程關(guān)鍵時(shí)間點(diǎn)的容器內(nèi)部溫度場(chǎng)分布情況，從而及時(shí)、快速和經(jīng)濟(jì)地得到運(yùn)輸容器的卸料冷卻時(shí)間參數(shù)以及容器內(nèi)部的溫度場(chǎng)分布情況，為后續(xù)的結(jié)構(gòu)改進(jìn)、制定運(yùn)行程序和試驗(yàn)策略提供重要的參考依據(jù)。

1 乏燃料運(yùn)輸容器本體結(jié)構(gòu)

乏燃料運(yùn)輸容器為鋼制圓柱形密閉金屬筒，主要由外殼組件、外蓋組件、內(nèi)蓋組件和吊籃組件組成，結(jié)構(gòu)見(jiàn)圖1。吊籃組件主要包括頂板組件、底板組件、貯存套管、傳熱板、支撐板以及螺桿組成。內(nèi)蓋組件分布在筒體上部，進(jìn)(排)氣、進(jìn)(排)水孔分布在內(nèi)蓋組件上，其中進(jìn)(排)水孔會(huì)有管道直通到容器底部。

圖1 乏燃料運(yùn)輸容器結(jié)構(gòu)示意圖

在卸料冷卻工藝過(guò)程中，外蓋組件不安裝，冷卻水(25 ℃)通過(guò)進(jìn)(排)水孔以及容器內(nèi)的管道進(jìn)入容器底部，氣體(或者置換水)從進(jìn)(排)氣孔排出。

2 計(jì)算程序及方法

筆者采用專(zhuān)業(yè)前處理軟件ICEM程序創(chuàng)建網(wǎng)格模型，國(guó)際通用的計(jì)算流體計(jì)算軟件Fluent求解結(jié)果和后處理。首先，通過(guò)ICEM對(duì)容器、傳熱支撐板以及流體區(qū)域建立了規(guī)則的六面體網(wǎng)格，再將網(wǎng)格導(dǎo)入Fluent程序進(jìn)行瞬態(tài)計(jì)算。

整個(gè)卸料冷卻過(guò)程中的衰變熱主要通過(guò)強(qiáng)制對(duì)流和輻射效應(yīng)來(lái)進(jìn)行散熱，筆者忽略輻射效應(yīng)的影響，將乏燃料組件等效為具有均勻等效導(dǎo)熱系數(shù)的體熱源[1-2]。在數(shù)值計(jì)算中，選擇標(biāo)準(zhǔn)k-ε方程湍流模型，增強(qiáng)壁面函數(shù)，采用Eulerian多相流的沸騰模型，并考慮了重力的影響，采用coupled算法進(jìn)行三維模擬瞬態(tài)分析，以確定卸料冷卻所需時(shí)間以及關(guān)鍵時(shí)間點(diǎn)的溫度場(chǎng)分布。

由于容器內(nèi)部包括大量的傳熱板和支撐板，使得容器的幾何模型和內(nèi)部流域非常復(fù)雜，分析起來(lái)較為困難。報(bào)告中采用兩種幾何模型來(lái)進(jìn)行分析：一是忽略傳熱板和支撐板；二是建立包括傳熱板和支撐板的模型。通過(guò)比較兩種模型的結(jié)果，分析有無(wú)傳熱支撐板對(duì)卸料冷卻計(jì)算時(shí)間造成的差異，以確定傳熱支撐板對(duì)卸料冷卻過(guò)程的影響。

3 數(shù)值模型

3.1 計(jì)算模型

乏燃料運(yùn)輸容器的實(shí)際結(jié)構(gòu)比較復(fù)雜，建立有限元模型時(shí)，需要對(duì)實(shí)際的幾何結(jié)構(gòu)進(jìn)行一定的簡(jiǎn)化和保守設(shè)定，具體簡(jiǎn)化內(nèi)容如下：

(1) 因螺栓緊固件等不影響整體傳熱，忽略所有的緊固件以及細(xì)小零部件。

(2) 為保守計(jì)算，忽略吊籃組件中螺紋桿的影響。

(3) 因進(jìn)、排氣(水)口位置對(duì)傳熱影響極低，為簡(jiǎn)化模型，假設(shè)進(jìn)、排氣(水)口位于容器軸線上。

(4) 容器堆焊層占比容器壁厚較低，且不銹鋼堆焊層與低碳鋼相比傳熱系數(shù)更低，忽略容器內(nèi)層的不銹鋼堆焊層。

(5) 筒體、內(nèi)蓋組件以及乏燃料單元組件(包括貯存套管)用均勻固體區(qū)域來(lái)等效。

(6) 為保守計(jì)算，忽略輻射效應(yīng)對(duì)散熱的貢獻(xiàn)。

(7)采用等效體熱源來(lái)處理乏燃料組件區(qū)域，忽略乏燃料組件區(qū)域的對(duì)流傳熱。

(8) 充水排氣階段冷卻水緩慢流入容器過(guò)程，簡(jiǎn)化成容器內(nèi)部一直充滿水靜置3 500 s。

(9)熱水置換階段將以充滿水的容器靜置3 500 s后的溫度場(chǎng)作為初始條件。

3.2 網(wǎng)格劃分

冷卻水從下向上流過(guò)乏燃料組件和容器內(nèi)壁之間的間隙，帶走衰變熱。衰變熱傳遞過(guò)程涉及多層結(jié)構(gòu)材料的共軛傳熱(固體壁面的對(duì)流換熱和結(jié)構(gòu)材料內(nèi)部導(dǎo)熱)。不考慮容器外部空氣域，對(duì)整個(gè)運(yùn)輸容器的固體區(qū)域和流體區(qū)域建模。

由于容器幾何結(jié)構(gòu)對(duì)稱(chēng)，取1/4區(qū)域?yàn)檠芯繉?duì)象，采用ICEM程序創(chuàng)建六面體網(wǎng)格模型(見(jiàn)圖2)。

圖2 乏燃料運(yùn)輸容器計(jì)算域網(wǎng)格

3.3 求解設(shè)置

利用Fluent對(duì)網(wǎng)格模型進(jìn)行計(jì)算。

流體區(qū)域雷諾數(shù)：

(1)

式中：ρ為冷卻劑密度，1 000 kg/m3；v為冷卻劑流速，3.73 m/s；μ為運(yùn)動(dòng)黏度，0.001 m2/s；d為流道直徑，25 mm。

其中，流體區(qū)域是不可壓縮的湍流流動(dòng)，采用k-ε方程湍流模型。固體區(qū)域內(nèi)部主要是熱傳導(dǎo)，流體內(nèi)部傳熱通過(guò)建立沸騰模型來(lái)實(shí)現(xiàn)，流體和固體之間設(shè)置交界面。

邊界條件：環(huán)境溫度為30 ℃，入口水溫為25 ℃，出口允許水溫為36 ℃。

具體條件設(shè)置如下：

(1) 對(duì)稱(chēng)邊界。

(2) 速度入口。

(3) 壓力出口。

(4) 容器外壁面設(shè)備環(huán)境溫度30 ℃。

(5) 乏燃料組件區(qū)域設(shè)置體熱源，q=6 235.8 W/m3。

(6) 乏燃料組件、筒體及內(nèi)蓋區(qū)域設(shè)置等效導(dǎo)熱系數(shù)，不同區(qū)域的等效導(dǎo)熱系數(shù)見(jiàn)表1。

表1 不同區(qū)域的等效導(dǎo)熱系數(shù)W/(m·K)

4 計(jì)算結(jié)果及分析

通過(guò)建立包含與不包含傳熱支撐板的兩種模型，并分別進(jìn)行充水和熱水置換兩個(gè)階段的瞬態(tài)計(jì)算。按照簡(jiǎn)化假設(shè)，充水階段簡(jiǎn)化成初始條件為容器滿水，之后靜置3 500 s；而熱水置換的初始條件則為容器充滿水后3 500 s的溫度場(chǎng)，直到出口溫度降到36 ℃停止計(jì)算。

4.1 溫度隨時(shí)間的變化曲線

圖3為兩種模型在3 500～10 000 s的出口溫度隨時(shí)間的變化曲線，圖4為兩種模型在0 ～10 000 s的流域溫度隨時(shí)間的變化曲線。分析溫度變化曲線發(fā)現(xiàn)：模型一(不包含傳熱支撐板)在經(jīng)過(guò)3 500 s的充水階段后，出口最高水溫達(dá)到57 ℃左右，進(jìn)入熱水置換階段后5 850 s后出口水溫小于允許水溫(36 ℃)；模型二(包含傳熱支撐板)在經(jīng)過(guò)3 500 s的充水排氣階段后，出口最高水溫能達(dá)到100 ℃左右，進(jìn)入熱水置換階段后7 300 s后出口水溫小于允許水溫(36 ℃)。其中，模型二出口最高水溫比較高的原因是受出口回流的影響，實(shí)際出口最高水溫小于373 K(100 ℃)。為了減小出口回流的影響，監(jiān)測(cè)溫度的出口面選在距離實(shí)際出口往容器內(nèi)部300 mm的位置，相對(duì)于模型一，模型二的計(jì)算模型更接近真實(shí)。

圖3 出口水溫隨時(shí)間的變化情況(3 500 s～10 000 s)

圖4 流域水溫隨時(shí)間的變化情況(0～10 000 s)

4.2 關(guān)鍵時(shí)間點(diǎn)的溫度場(chǎng)分布

圖5、圖6為兩種模型在關(guān)鍵時(shí)間點(diǎn)計(jì)算域內(nèi)部溫度分布云圖，左邊為充水排氣階段末的溫度分布云圖，右邊為熱水置換階段末的溫度分布云圖。

圖5 模型一 關(guān)鍵時(shí)間點(diǎn)計(jì)算域內(nèi)部溫度分布云圖

圖6 模型二 關(guān)鍵時(shí)間點(diǎn)計(jì)算域內(nèi)部溫度分布云圖

圖5顯示模型一在充水排氣階段末乏燃料熱源區(qū)域溫度較高、壁面溫度滿足303.15 K(30 ℃)的恒溫條件，入口溫度是303.15 K(30 ℃)，出口溫度達(dá)到330 K(57 ℃)左右溫度隨著乏燃料區(qū)域往外遞減，基本滿足實(shí)際情況，分析結(jié)果基本合理。熱水置換階段末，計(jì)算域的整體溫度和最高溫度下降，最高溫度位于乏燃料熱源區(qū)域中心，最高溫度約為353.4 K(80 ℃)，出口溫度由3 500 s時(shí)的330 K(57 ℃)左右下降到306.6 K(34 ℃)左右，滿足出口允許水溫要求。

圖6顯示模型二在充水排氣階段末，流域最高溫度位于乏燃料熱源附近區(qū)域(不考慮出口回流高溫的影響)，最高溫度約為342.3 K(69 ℃)，支撐板的導(dǎo)熱系數(shù)相對(duì)較低，兩端水域溫度較低，基本滿足實(shí)際情況，分析結(jié)果基本合理。熱水置換階段末，計(jì)算域的整體溫度和最高溫度下降，最高溫度還是位于乏燃料熱源附近區(qū)域中心，最高溫度約為338.8 K(66 ℃)，出口溫度下降到307.9 K(35 ℃)左右，滿足出口允許水溫要求。

5 結(jié)語(yǔ)

筆者嘗試?yán)肍luent計(jì)算軟件，通過(guò)合理的模型假設(shè)，基于自主設(shè)計(jì)的乏燃料運(yùn)輸容器，進(jìn)而對(duì)容器的卸料冷卻過(guò)程進(jìn)行熱工仿真分析，獲得卸料冷卻過(guò)程中容器內(nèi)部溫度場(chǎng)分布情況以及出口溫度變化曲線，從而確定此設(shè)備的卸料冷卻過(guò)程總時(shí)長(zhǎng)、充水排氣時(shí)間以及熱水置換時(shí)間。

計(jì)算結(jié)果表明：模型一的卸料冷卻過(guò)程的總時(shí)長(zhǎng)為2.5 h，模型二的卸料冷卻過(guò)程的總時(shí)長(zhǎng)為3 h。傳熱支撐板對(duì)乏燃料運(yùn)輸容器的卸料冷卻時(shí)間計(jì)算影響較大，在相關(guān)的熱工分析過(guò)程中不能忽略。

通過(guò)本次熱工仿真分析結(jié)果，可以為后續(xù)結(jié)構(gòu)改進(jìn)，制定運(yùn)行程序和試驗(yàn)策略提供重要的參考依據(jù)。盡可能減少運(yùn)行人員暴露在放射性環(huán)境下的時(shí)長(zhǎng)，降低人員操作過(guò)程以及試驗(yàn)過(guò)程中的劑量率，并且通過(guò)識(shí)別關(guān)鍵時(shí)間節(jié)點(diǎn)及時(shí)有效地監(jiān)控整個(gè)卸料冷卻過(guò)程。

[1] 周有新, 宋磊, 陳軍. 乏燃料干式貯存容器衰變熱導(dǎo)出分析[C]//第十四屆全國(guó)反應(yīng)堆熱工流體學(xué)術(shù)會(huì)議暨中核核反應(yīng)堆熱工水力技術(shù)重點(diǎn)實(shí)驗(yàn)室2015年度學(xué)術(shù)年會(huì)論文集. 北京: 清華大學(xué)先進(jìn)反應(yīng)堆工程與安全教育部重點(diǎn)實(shí)驗(yàn)室, 2015.

[2] 葉水祥. 乏燃料水池內(nèi)流動(dòng)與傳熱數(shù)值分析[D]. 哈爾濱: 哈爾濱工程大學(xué), 2012.

InfluenceofSupportingSlabonHeatTransferBehaviorofaSpentFuelTransportContainerduringUnloadingProcess

Yin Yong, Li Qipeng, Ma Qingjun

(China Nuclear Power Technology Research Institute, Shenzhen 518124, Guangdong Province, China)

2016-08-23；

2016-12-27

殷 勇(1970—)，男，研究員級(jí)高級(jí)工程師，主要從事核電設(shè)備研發(fā)及國(guó)產(chǎn)化工作。

E-mail: yin_yong@cgnpc.com.cn

TL248

A

1671-086X(2017)05-0336-04