基于縮比實(shí)驗(yàn)?zāi)Ｐ偷目於阎魅萜鲀?nèi)氬氣空間傳熱特性數(shù)值研究

羅心越，陸道綱，馮佳琪，于宗玉，張鈺浩,*

基于縮比實(shí)驗(yàn)?zāi)Ｐ偷目於阎魅萜鲀?nèi)氬氣空間傳熱特性數(shù)值研究

羅心越1,2，陸道綱1,2，馮佳琪1,2，于宗玉1,2，張鈺浩1,2,*

（1. 華北電力大學(xué)核科學(xué)與工程學(xué)院，北京 102206；2. 非能動(dòng)核能安全技術(shù)北京市重點(diǎn)實(shí)驗(yàn)室，北京 102206）

池式鈉冷快堆主容器內(nèi)鈉液面上方覆蓋有氬氣，鈉液面通過多種傳熱方式相互耦合，不斷向堆容器上部結(jié)構(gòu)傳遞熱量。為了對(duì)上部結(jié)構(gòu)進(jìn)行溫度載荷下的強(qiáng)度分析，有必要開展氬氣空間傳熱特性研究。一方面，國(guó)內(nèi)外針對(duì)錐頂蓋氬氣空間傳熱特性開展了較多的數(shù)值模擬研究，但大多缺乏實(shí)驗(yàn)驗(yàn)證。另一方面，雖有少量針對(duì)氬氣空間的換熱特性實(shí)驗(yàn)研究，但由于測(cè)量技術(shù)的制約，僅能獲得關(guān)鍵部件的溫度分布，難以獲得其流動(dòng)特性，因此對(duì)氬氣空間復(fù)雜的傳熱特性和傳熱機(jī)理有待進(jìn)一步研究。本文基于快堆鈉液面和氬氣空間換熱特性縮比實(shí)驗(yàn)，開展氬氣空間熱工水力特性數(shù)值模擬，通過與實(shí)驗(yàn)結(jié)果對(duì)比，驗(yàn)證了數(shù)值方法的有效性，并在此基礎(chǔ)上，對(duì)其溫度分布及流動(dòng)特性進(jìn)行詳細(xì)研究。研究結(jié)果表明：在不同工況下，氬氣空間內(nèi)溫度分布總體對(duì)稱，由于徑向熱屏蔽的存在，其上部窄流道區(qū)域溫度明顯低于下部氬氣大空間溫度；徑向熱屏蔽與底部鈉液面間的氬氣大空間內(nèi)，由于瑞利（）數(shù)較高，自然對(duì)流現(xiàn)象強(qiáng)烈，換熱效果較強(qiáng)，在靠近底部鈉液面處形成局部渦流，而徑向熱屏蔽外部窄流道區(qū)域內(nèi)由于數(shù)較小，自然對(duì)流較弱，主要依靠熱輻射和導(dǎo)熱進(jìn)行傳熱。

鈉冷快堆；錐頂蓋；氬氣空間；對(duì)流換熱

池式鈉冷快堆是第四代核電技術(shù)中的重要堆型，其堆本體結(jié)構(gòu)復(fù)雜，包含旋塞、氬氣空間、中間熱交換器（IHX）、獨(dú)立熱交換器（DHX）、泵支承、溫度液位計(jì)等多個(gè)部件。堆本體內(nèi)一回路采用液態(tài)鈉作為冷卻劑，錐頂蓋空間內(nèi)覆蓋氬氣作為保護(hù)氣體。反應(yīng)堆運(yùn)行時(shí)，鈉液面不斷向堆容器的上部結(jié)構(gòu)（徑向熱屏蔽、斜肩等）及氬氣傳遞熱量。錐頂蓋氬氣空間內(nèi)傳熱機(jī)理復(fù)雜，存在多種傳熱方式相互耦合的復(fù)雜行為，錐頂蓋溫度場(chǎng)分布難以準(zhǔn)確獲得。此外，反應(yīng)堆運(yùn)行時(shí)，液鈉蒸氣在氬氣空間內(nèi)形成的氣溶膠對(duì)熱輻射具有吸收作用，錐頂蓋、旋轉(zhuǎn)屏蔽塞和貫穿件支承在徑向和軸向方向上存在局部熱傳導(dǎo)等，這進(jìn)一步增加了錐頂蓋溫度場(chǎng)分布的復(fù)雜性，給錐頂蓋溫度載荷下的強(qiáng)度分析和應(yīng)力評(píng)價(jià)及設(shè)計(jì)帶來困難。

針對(duì)氬氣空間換熱特性，部分學(xué)者開展了基于簡(jiǎn)化模型的實(shí)驗(yàn)研究。吳強(qiáng)等[1]建立模擬快堆氬氣空間的水-空氣實(shí)驗(yàn)裝置，研究封閉空間內(nèi)氣體溫度分布特征，發(fā)現(xiàn)液面上存在兩種換熱方式-輻射換熱和對(duì)流換熱，且液面的蒸發(fā)現(xiàn)象能夠阻礙兩種換熱效果。李天舒等[2]基于相似原理建立快堆氬氣空間的簡(jiǎn)化實(shí)驗(yàn)裝置，實(shí)驗(yàn)結(jié)果顯示氣體溫度在封閉空間中間大部分區(qū)域維持相對(duì)恒定，同時(shí)提出了氬氣空間高度和主容器冷卻系統(tǒng)設(shè)計(jì)應(yīng)統(tǒng)籌考慮的設(shè)計(jì)思想。

為了探究錐頂蓋及氬氣空間復(fù)雜熱工水力特性，國(guó)內(nèi)一些學(xué)者建立了詳細(xì)的三維模型，張曉露[7]針對(duì)錐頂蓋主要部件結(jié)構(gòu)進(jìn)行精細(xì)建模，分析旋塞與旋塞支承頸氣隙及生物屏蔽柱對(duì)金屬層溫度分布的影響，發(fā)現(xiàn)生物屏蔽柱對(duì)金屬層的輻射效果較強(qiáng)，30 mm氣隙對(duì)錐頂蓋整體的保溫效果最好。余新太等[8]應(yīng)用低雷諾數(shù)湍流模型和Discrete Ordinates（DO）輻射傳熱模型對(duì)正常運(yùn)行工況下示范快堆氬氣空間進(jìn)行了數(shù)值模擬，計(jì)算結(jié)果表明，泵支承、中間熱交換器、獨(dú)立熱交換器周向溫差較大，需注意彎曲變形。林超等[9]對(duì)氬氣空間及相關(guān)結(jié)構(gòu)部件進(jìn)行穩(wěn)態(tài)共軛傳熱計(jì)算，發(fā)現(xiàn)原型設(shè)計(jì)方案泵支承結(jié)構(gòu)位于氬氣空間部分周向溫度分布有較大溫度梯度，泵支承結(jié)構(gòu)靠近主容器支承頸側(cè)溫度明顯高于遠(yuǎn)離支承頸側(cè)。盡管以上數(shù)值模擬揭示了氬氣空間內(nèi)部部分復(fù)雜傳熱行為和機(jī)理特性，但其結(jié)果大多缺乏實(shí)驗(yàn)驗(yàn)證，針對(duì)數(shù)值模擬中的參數(shù)選取、計(jì)算方法也存在較大不確定性。因此，開展主容器頂部空間內(nèi)氬氣的復(fù)雜流動(dòng)特性的實(shí)驗(yàn)研究十分必要。

綜上，由于池式快堆頂部氬氣空間結(jié)構(gòu)復(fù)雜、耦合傳熱機(jī)理復(fù)雜且受到多種外部因素的影響，難以開展等比例實(shí)驗(yàn)研究主容器頂部空間內(nèi)氬氣的復(fù)雜流動(dòng)特性，因此，為了揭示氬氣空間的換熱機(jī)理、驗(yàn)證其換熱模型，基于相似原理建立了快堆鈉液面和氬氣空間換熱縮比實(shí)驗(yàn)臺(tái)架，同時(shí)，鈉液面具有流動(dòng)性，其加熱所產(chǎn)生的鈉蒸氣會(huì)對(duì)氬氣空間內(nèi)部氬氣的熱工水力特性產(chǎn)生影響，為了排除鈉蒸氣氣溶膠對(duì)輻射換熱效果的影響，開展了采用不銹鋼模擬鈉液面運(yùn)行工況及換料工況下的傳熱特性研究實(shí)驗(yàn)，并開展池式快堆主容器頂部氬氣空間熱工水力特性數(shù)值模擬，驗(yàn)證計(jì)算方法的有效性，進(jìn)而探究氬氣空間內(nèi)多種傳熱方式并存的復(fù)雜耦合傳熱機(jī)理，為主容器頂部空間溫度載荷下的強(qiáng)度分析應(yīng)力評(píng)定提供關(guān)鍵參數(shù)。

1　實(shí)驗(yàn)簡(jiǎn)介

在快堆主容器內(nèi)氬氣空間換熱理論分析基礎(chǔ)上，搭建鈉液面與氬氣空間換熱縮比實(shí)驗(yàn)臺(tái)架。實(shí)驗(yàn)裝置主要包括通風(fēng)系統(tǒng)、電加熱系統(tǒng)、實(shí)驗(yàn)裝置以及數(shù)據(jù)采集系統(tǒng)。其中，實(shí)驗(yàn)裝置基本保留了氬氣空間的主要結(jié)構(gòu)，包含錐頂蓋、獨(dú)立熱交換器（DHX）、中間熱交換器（IHX）、泵支承、徑向熱屏蔽、旋塞等。底部為封閉的不銹鋼圓筒，加熱裝置在不銹鋼底板的底面以完成實(shí)驗(yàn)裝置的加熱，實(shí)驗(yàn)裝置結(jié)構(gòu)如圖1所示。因此該實(shí)驗(yàn)裝置可用于模擬有/無鈉液面條件下快堆主容器頂部空間覆蓋氣體的自然循環(huán)特性。

圖1　實(shí)驗(yàn)臺(tái)架整體結(jié)構(gòu)圖

測(cè)量及數(shù)據(jù)采集系統(tǒng)，包括熱電偶、壓力表、流量計(jì)及數(shù)據(jù)采集系統(tǒng)（NI PIXe），用來進(jìn)行數(shù)據(jù)采集、存儲(chǔ)和后處理。其中，通過將熱電偶點(diǎn)焊在錐頂蓋表面以測(cè)量其表面溫度，將熱電偶布置在不同高度的支架上以監(jiān)測(cè)氬氣空間內(nèi)溫度分布，其位置如圖2所示。每個(gè)支架上隨軸向高度布置6個(gè)測(cè)點(diǎn)。

圖2　測(cè)溫點(diǎn)分布示意圖

2　數(shù)值模型及方法

2.1　幾何模型

基于鈉液面與氬氣空間換熱縮比實(shí)驗(yàn)臺(tái)架建立數(shù)值計(jì)算幾何模型，如圖3所示。

圖3　氬氣空間三維模型結(jié)構(gòu)

本研究使用STAR-CCM+軟件對(duì)主容器錐頂蓋進(jìn)行建模。錐頂蓋內(nèi)包含獨(dú)立熱交換器（DHX）、中間熱交換器（IHX）、泵支承、徑向熱屏蔽、旋塞等多個(gè)部件，同時(shí)存在許多窄縫流道。為了能夠?qū)崿F(xiàn)模型關(guān)鍵部件的主要運(yùn)行功能，模擬鈉冷快堆頂部氬氣空間的自然循環(huán)現(xiàn)象，同時(shí)考慮到整體網(wǎng)格質(zhì)量、生成難度以及保證計(jì)算精度，根據(jù)結(jié)構(gòu)特點(diǎn)和流動(dòng)特性，對(duì)各部件單獨(dú)進(jìn)行網(wǎng)格劃分，將氬氣空間分為貫穿件、錐頂蓋金屬層保溫層、旋塞、凸臺(tái)等不同模塊。同時(shí)，在建模過程中盡可能保留錐頂蓋的主要幾何特征，將螺栓、熱電偶、液位計(jì)等尺寸較小的部件省略，將貫穿件及其周圍屏蔽簡(jiǎn)化為實(shí)體或空心圓柱，將旋塞屏蔽層簡(jiǎn)化為實(shí)體結(jié)構(gòu)，以便有效模擬錐頂蓋氬氣空間內(nèi)的關(guān)鍵熱工特性。

基于所建立的幾何模型開展非結(jié)構(gòu)化網(wǎng)格劃分，并進(jìn)行網(wǎng)格敏感性分析，由于鈉冷快堆頂部氬氣空間內(nèi)部結(jié)構(gòu)復(fù)雜，建立2 500萬、4 900萬、8 000萬的網(wǎng)格開展網(wǎng)格無關(guān)性分析計(jì)算，對(duì)比錐頂蓋平均溫度，結(jié)果如圖4所示。

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

由圖4可見，2 500萬網(wǎng)格和4 900萬網(wǎng)格模擬結(jié)果相差較大，4 900萬網(wǎng)格與8 000萬網(wǎng)格的模擬結(jié)果差異小于3%，因此，能夠認(rèn)為，4 900萬網(wǎng)格計(jì)算結(jié)果滿足計(jì)算要求。最終選取約4 900萬網(wǎng)格開展驗(yàn)證分析。

2.2　數(shù)值模型及邊界條件

2.2.1數(shù)值模型

在近壁面區(qū)域，湍流應(yīng)力不起作用，因此，在近壁面區(qū)域，采用壁面函數(shù)法，無需將壁面區(qū)域加密，用一組半經(jīng)驗(yàn)公式將壁面上的物理量與湍流核心區(qū)內(nèi)相應(yīng)物理量相聯(lián)系，直接得到與壁面相鄰控制體積的節(jié)點(diǎn)變量值[11]，在本文中，采用STAR-CCM+默認(rèn)的壁面模型Ally+wall treatment計(jì)算。

DO模型能夠較好地解決絕大多數(shù)光學(xué)深度區(qū)間的輻射問題，能夠考慮到氣體介質(zhì)與顆粒間的輻射交換，能量不平衡率最小，在存在局部熱源的情況下，能夠較好地反應(yīng)模型輻射換熱情況。

本文選擇DO模型以求解有限數(shù)量的離散立體角的輻射運(yùn)輸方程，其原理是通過計(jì)算每個(gè)離散縱坐標(biāo)上的偏微分方程來求解強(qiáng)度[9]，輻射傳熱方程形式如下：

——射線路徑長(zhǎng)度；

——當(dāng)?shù)販囟取?/p>

DO模型是將球坐標(biāo)系下的立體角轉(zhuǎn)化為直角坐標(biāo)系下的方向，把沿方向的輻射傳熱方程視為某個(gè)場(chǎng)方程，于是式（3）變形為：

2.2.2計(jì)算工況及邊界條件

為了驗(yàn)證不同工況下主容器頂部空間的傳熱特性，在實(shí)驗(yàn)臺(tái)架中，針對(duì)不銹鋼底面模擬鈉液面溫度開展運(yùn)行工況及換料工況兩組關(guān)鍵工況的實(shí)驗(yàn)，在模擬計(jì)算中，不同工況下模擬鈉液面溫度的設(shè)置與實(shí)驗(yàn)臺(tái)架保持一致。分別為：

（1）運(yùn)行溫度工況：不銹鋼底面（模擬鈉液面）溫度為540 ℃，通過分析該工況下各堆內(nèi)構(gòu)件的關(guān)鍵三維熱動(dòng)水力參數(shù)，評(píng)估在反應(yīng)堆正常運(yùn)行下頂部空間的壓力邊界安全性。

（2）換料溫度工況：不銹鋼底面（模擬鈉液面）溫度為250 ℃，基于該工況下氬氣空間的溫度分布特性，驗(yàn)證在裝卸料情況下，反應(yīng)堆頂部空間的壓力邊界完整性。

不同工況下邊界條件如表1所示。

表1　各位置邊界條件

基于上述關(guān)鍵邊界條件的設(shè)置，開展鈉液面與氬氣空間換熱實(shí)驗(yàn)數(shù)值模擬計(jì)算。計(jì)算過程中所涉及的不銹鋼、保溫層及氬氣的密度、導(dǎo)熱系數(shù)、比熱容等物性參數(shù)均來源于《載熱質(zhì)熱物性計(jì)算程序及數(shù)據(jù)手冊(cè)》[10]。

3　氬氣空間傳熱特性分析

3.1　數(shù)值方法驗(yàn)證

在氬氣空間數(shù)值模擬中，運(yùn)行工況、換料工況下錐頂蓋溫度分布如圖5所示。

圖5　錐頂蓋溫度分布

如圖5所示，錐頂蓋上溫度總體呈對(duì)稱分布，運(yùn)行工況下，錐頂蓋溫度基本穩(wěn)定在260～400 ℃，換料工況下，錐頂蓋溫度穩(wěn)定在80～120 ℃。整體溫度分布較均勻，但在貫穿件附近，溫度出現(xiàn)明顯上升，推測(cè)因?yàn)闅鍤馀c貫穿件的對(duì)流換熱加強(qiáng)，導(dǎo)致靠近貫穿件側(cè)的部分溫度升高。

對(duì)照實(shí)驗(yàn)臺(tái)架布置的測(cè)點(diǎn)位置，在數(shù)值模擬計(jì)算中選擇相同位置處（見圖3）的溫度進(jìn)行對(duì)比分析。選擇Line 2作為監(jiān)測(cè)線，對(duì)比如圖6所示，分別分析實(shí)驗(yàn)運(yùn)行工況、模擬運(yùn)行工況、實(shí)驗(yàn)換料工況、模擬換料工況的溫度變化。

圖6　Line 2溫度變化曲線

通過溫度變化對(duì)比曲線可以看出，數(shù)值模擬計(jì)算結(jié)果的溫度變化趨勢(shì)與實(shí)驗(yàn)結(jié)果基本一致，與實(shí)驗(yàn)結(jié)果相比，數(shù)值模擬得到溫度平均相對(duì)誤差在 5%以內(nèi)?？梢则?yàn)證，數(shù)值模擬所采用的計(jì)算模型、計(jì)算方法合理，能夠較好地描述主容器頂部氬氣空間內(nèi)覆蓋氣體的熱工水力特性。

3.2　氬氣空間傳熱特性分析

3.2.1氬氣空間溫度分布

在氬氣空間內(nèi)選擇監(jiān)測(cè)面截面（1）、截面（2），得到氬氣空間的溫度分布如圖7所示。可以發(fā)現(xiàn)，氬氣空間的溫度分布基本呈對(duì)稱狀態(tài)，運(yùn)行工況下最高溫度為 540 ℃，最低溫度約為411 ℃，溫差約 131 ℃，換料工況最高溫度為250 ℃，最低溫度約為 151 ℃，溫差約 99 ℃，垂直高度上溫度自下而上逐漸降低，最高溫度出現(xiàn)在鈉液面附近。由于徑向熱屏蔽的存在，上部窄流道區(qū)域溫度明顯低于下部氬氣大空間溫度，其存在非常有效的阻隔了溫度的傳導(dǎo)，使金屬層和保溫層的溫度始終處于主容器的設(shè)計(jì)限值內(nèi)。

圖7　監(jiān)測(cè)面溫度分布

此外，從截面溫度云圖中能夠明顯看出，在氬氣空間底部存在兩個(gè)對(duì)稱位置的局部渦流區(qū)域，在該區(qū)域內(nèi)，氬氣溫度分布呈明顯的不規(guī)律性，形成兩個(gè)小的漩渦。提取換料溫度工況下截面速度矢量分布如圖8所示?？拷ǖ榔帘蔚囊粋?cè)渦流區(qū)域，由于氬氣在該位置同時(shí)與通道屏蔽和金屬層進(jìn)行對(duì)流傳熱和輻射換熱，溫度明顯較低，密度增大，向下回流，從而形成局部渦流。

圖8　氬氣空間整體流場(chǎng)

在數(shù)值模擬中選擇與實(shí)驗(yàn)裝置相同位置處（見圖 3）的監(jiān)測(cè)線進(jìn)行溫度對(duì)比分析，如圖9所示。

基于線上各點(diǎn)的溫度分布可知，氬氣空間內(nèi)溫度隨垂直高度的增加而降低，變化基本趨勢(shì)不變，Line 2、Line 3在氬氣空間內(nèi)呈對(duì)稱狀態(tài)，其溫度變化趨勢(shì)基本一致，進(jìn)一步驗(yàn)證在整個(gè)氬氣空間內(nèi)，溫度分布基本呈對(duì)稱狀態(tài)，分析Line 4可以看到，在貫穿件附近，由于通道屏蔽和凸臺(tái)的存在，氬氣的自然循環(huán)高度降低，因此，相較Line 1，Line 4的變化趨勢(shì)更平緩、溫度梯度更小。而由于氬氣與貫穿件的熱交換作用明顯，氬氣溫度會(huì)出現(xiàn)小范圍波動(dòng)。

3.2.2換熱機(jī)理分析

在傳熱過程中，氬氣在空間內(nèi)流動(dòng)，將熱量傳遞給徑向熱屏蔽，經(jīng)熱傳導(dǎo)傳遞到徑向熱屏蔽外表面，最后通過外表面?zhèn)鬟f到金屬層。自然循環(huán)過程中，自然對(duì)流傳熱是氬氣空間內(nèi)引起氬氣與物體表面之間的熱量傳遞方式。

圖9　不同工況下監(jiān)測(cè)線溫度曲線圖

瑞利（）數(shù)是與自然對(duì)流相關(guān)的無量綱數(shù)，反映自然對(duì)流效果，其公式為：

式中，格拉曉夫（）數(shù)描述了流體的浮力和粘度之間的關(guān)系，普朗特?cái)?shù)（）描述了動(dòng)量擴(kuò)散系數(shù)和熱擴(kuò)散系數(shù)之間的關(guān)系。普朗特?cái)?shù)（）反映了動(dòng)量擴(kuò)散厚度和能量擴(kuò)散厚度的相對(duì)大小。

——重力加速度為9.8 m/s2；

為了分析氬氣空間內(nèi)基于自然對(duì)流機(jī)理下的氬氣換熱特性，在氬氣空間內(nèi)選擇垂直于徑向熱屏蔽的兩條分段平行的溫度線，如圖2所示。提取各監(jiān)測(cè)線溫度，計(jì)算對(duì)應(yīng)瑞利數(shù)如表2所示。

表2　瑞利(Ra) 數(shù)表

圖10為不同工況下溫度線在無量綱高度上的溫度分布。從圖中可知，溫度隨無量綱高度增加逐漸降低，無量綱高度為0～0.2，溫度變化劇烈且變化梯度大，無量綱高度在0.2～1，溫度下降趨勢(shì)變緩，結(jié)合瑞利（）數(shù)分析可得，在徑向熱屏蔽下氬氣空間大區(qū)域內(nèi)，數(shù)達(dá)到105～107量級(jí)，空間越大，瑞利（）數(shù)越大，自然對(duì)流特征越明顯，驅(qū)動(dòng)力越大，因此會(huì)產(chǎn)生局部渦旋（見圖7）。而在徑向熱屏蔽外部窄流道區(qū)域內(nèi)，數(shù)較小，自然對(duì)流影響較小，溫度分布更均勻，主要依靠輻射和導(dǎo)熱進(jìn)行換熱。

圖10　監(jiān)測(cè)位置溫度分布變化

4　結(jié)論

本文基于快堆鈉液面和氬氣空間換熱臺(tái)架運(yùn)行工況及換料工況下的傳熱特性研究實(shí)驗(yàn)，對(duì)鈉冷快堆頂部氬氣空間開展整體三維數(shù)值模擬分析。獲得主要結(jié)論如下：

（2）氬氣空間溫度分布整體基本呈對(duì)稱狀態(tài)，在垂直高度上自下而上出現(xiàn)明顯溫度分層，平均溫差達(dá)到150 ℃。由于徑向熱屏蔽有效阻隔了熱量的傳導(dǎo)，上部窄流道區(qū)域溫度明顯低于下部氬氣大空間溫度。

（3）在徑向熱屏蔽下的氬氣空間大區(qū)域內(nèi)，數(shù)能夠達(dá)到105～107量級(jí)，自然對(duì)流換熱特征明顯，形成局部渦流，但在徑向熱屏蔽外部窄流道區(qū)域，數(shù)較小，自然對(duì)流換熱影響較小，主要依靠輻射和導(dǎo)熱進(jìn)行換熱。

［1］ 吳強(qiáng)，張東輝，吳純良，等. CEFR主容器氬氣空間溫度分布模擬實(shí)驗(yàn)研究［G］. 中國(guó)原子能科學(xué)研究院年報(bào)，2009：9-10.

［2］ 吳強(qiáng)，張東輝，李天舒，等. 基于相似理論的快堆氬氣空間溫度場(chǎng)實(shí)驗(yàn)研究［J］. 核科學(xué)與工程，2016，36 （04）：465-469.

［3］ 吳強(qiáng)，張東輝，劉云焰. 快堆主容器內(nèi)氬氣空間傳熱實(shí)驗(yàn)原理及分析［J］. 核科學(xué)與工程，2009，29（01）：27-32.

［4］ 許義軍. CEFR大小旋塞氬氣空間流固耦合傳熱計(jì)算的數(shù)值分析［C］. 計(jì)算流體力學(xué)研究進(jìn)展——第十二屆全國(guó)計(jì)算流體力學(xué)會(huì)議論文集，2004：726-731.

［5］ 張維忠，賈斗南，秋穗正. 高溫氣冷堆一回路艙室的輻射換熱計(jì)算模型［J］. 核動(dòng)力工程，2002（02）：10-16.

［6］ 郭磊. 對(duì)FLUENT輻射模型的數(shù)值計(jì)算與分析［J］. 制冷與空調(diào)（四川），2014，28（03）：358-360.

［7］ 張曉露. 快堆主容器頂部復(fù)雜空間中的覆蓋氣體熱工水力行為研究［D］. 北京：華北電力大學(xué)（北京），2019.

［8］ 余新太，高鑫釗，馬曉，等. 示范快堆主容器內(nèi)氬氣空間數(shù)值模擬［J］. 核科學(xué)與工程，2021，41（04）：695-702.

［9］ 林超，楊紅義，周志偉. CFR600泵支承氬氣空間部分三維傳熱數(shù)值模擬［J］. 原子能科學(xué)技術(shù)，2021，55（04）：654-659.

［10］居懷明，徐元輝，李懷萱. 載熱質(zhì)熱物性計(jì)算程序及數(shù)據(jù)手冊(cè)［M］. 北京：原子能出版社，1990.

［11］謝東. 地下水電站廠房氣流組織CFD數(shù)值模擬方法研究［D］. 重慶：重慶大學(xué)，2015.

［12］李光正，馬洪林. 封閉腔內(nèi)高瑞利數(shù)層流自然對(duì)流數(shù)值模擬［J］. 華中科技大學(xué)學(xué)報(bào)（城市科學(xué)版），2004（03）：14-17.

［13］邱志民. AP1000非能動(dòng)余熱排出系統(tǒng)一、二次側(cè)耦合傳熱特性實(shí)驗(yàn)研究［D］. 北京：華北電力大學(xué)（北京），2021.

［14］梁江濤，陸道綱，趙海琦，等. 池式鈉冷快堆雙環(huán)路12%差異非對(duì)稱功率運(yùn)行及流量調(diào)節(jié)緩解工況的三維數(shù)值模擬［J］. 科學(xué)技術(shù)與工程，2021，21（25）：10703-10709.

［15］ Dinesh Nath，Mahendra K.Verma.Numerical simulation of convection of argon gas in fast breeder reactor［G］. Annals of Nuclear Energy，2014：63.

［16］ Xiaoxue Huang，Shuisheng He，Numerical modelling of cover gas thermal hydraulics in Sodium-cooled Fast Reactors［J］. Nuclear Engineering and Design，2019，355：110347.

［17］ Dinesh Nath，Mahendra K.Verma.Numerical simulation of convection of argon gas in fast breeder reactor［G］. Annals of Nuclear Energy，2014：63.

［18］ Guleren K M，Pinarbasi A. Numerical simulation of the stalled flow within a vaned centrifugal pump. Proceedings of the Institution of Mechanical Engineers，Part C［J］. Journal of Mechanical Engineering Science，2004；218（4）：425-435.

Numerical Simulation of the Argon Space for the Heat Transfer Characteristics in the Main Vessel Based on the Scaled-down Test Model

LUO Xinyue1,2，LU Daogang1,2，F(xiàn)ENG Jiaqi1,2，YU Zongyu1,2，ZHANG Yuhao1,2,*

（1. School of Nuclear Science and Engineering，North China Electric Power University，Beijing 102206，China；2. Beijing Key Laboratory of Passive Safety Technology for Nuclear Energy，Beijing 102206，China）

The argon gas is covered above the sodium level in the main vessel of pool-type sodium-cooled fast reactor (SFR). The sodium level continuously transfers heat to upper structure of the roof slab through the coupled heat transfer mechanism. In order to analyze the strength and evaluate the stress of roof slab, it is necessary to study the heat transfer characteristics of the argon space. On the one hand, many numerical simulation studies have been carried out on the space heat transfer characteristics of the argon gas in the roof slab. Most of the numerical simulation studies on the roof slab and the argon space lack of experimental verification. In addition, due to the limitation of measurement technique, it is easy to obtain the temperature characteristics, but it is difficult to obtain the flow characteristics in the existing study of numerical simulation. Therefore, the complex heat transfer characteristics and mechanism of the argon space still needs study. In this paper, a numerical simulation of the thermal-hydraulic characteristics in the argon space was carried out based on the heat-transfer experimental facility for the argon space at the top of the main vessel (HEFA), whose results were compared with the experimental data. The results show that the numerical method is effective, and the calculated results are reasonable. By analyzing the three-dimensional temperature distribution of the domain, it is found that under different working conditions, the temperature distribution in the argon space is generally symmetric. Due to the existence of radial heat shielding, the temperature in the upper narrow region is obviously lower than that in the external narrow channel area. In addition, there is a vortex forming in the local argon space. The heat transfer effect of natural convection is obvious, while thenumber of argon near the radial thermal shield is small in the external narrow channel area. The main heat transfer processes are heat conduction and radiation heat transfer.

Sodium cooled fast reactor; Roof slab; Argon space; Convective heat transfer

TL33

A

0258-0918（2023）05-0961-10

2022-10-04

羅心越（1996—），女，甘肅天水人，碩士研究生，現(xiàn)主要從事反應(yīng)堆熱工水力學(xué)方面的研究

張鈺浩，E-mail：zhangyuhao@ncepu.edu.cn