壓力容器外部冷卻自然循環(huán)比例分析

劉宇生, 李聰新, 溫麗晶, 譚思超, 張盼

(1.環(huán)境保護(hù)部 核與輻射安全中心，北京 100086；2.哈爾濱工程大學(xué) 核安全與仿真技術(shù)國(guó)防重點(diǎn)學(xué)科實(shí)驗(yàn)室，黑龍江哈爾濱 150001)

壓力容器外部冷卻自然循環(huán)比例分析

劉宇生1, 李聰新1, 溫麗晶1, 譚思超2, 張盼1

(1.環(huán)境保護(hù)部 核與輻射安全中心，北京 100086；2.哈爾濱工程大學(xué) 核安全與仿真技術(shù)國(guó)防重點(diǎn)學(xué)科實(shí)驗(yàn)室，黑龍江哈爾濱 150001)

壓力容器外部壁面和保溫層之間環(huán)形流道內(nèi)的兩相自然循環(huán)是嚴(yán)重事故下實(shí)施壓力容器外部冷卻(ERVC)和熔融物堆內(nèi)滯留的關(guān)鍵過(guò)程，為研究ERVC整體性能試驗(yàn)裝置與核電廠原型間應(yīng)滿足的相似準(zhǔn)則，以守恒方程和一維漂移流模型為基礎(chǔ)，采用分級(jí)雙向比例方法對(duì)ERVC自然循環(huán)過(guò)程進(jìn)行了比例分析，通過(guò)假設(shè)和簡(jiǎn)化，獲得了ERVC自然循環(huán)試驗(yàn)裝置與原型的相似準(zhǔn)則，并給出了自然循環(huán)的穩(wěn)態(tài)解。結(jié)果表明：ERVC自然循環(huán)現(xiàn)象可采用等壓等物性方式進(jìn)行模擬，F(xiàn)roude數(shù)、密度數(shù)、焓升數(shù)和熱源數(shù)等相似準(zhǔn)則均能得到滿足，不存在比例轉(zhuǎn)換帶來(lái)的失真；采用等厚等半徑等熱流密度二維切片式加熱試驗(yàn)段模擬壓力容器下封頭，當(dāng)流道最小間距位置和尺寸與原型一致時(shí)，修正Stanton數(shù)相似比為1，下封頭沸騰換熱和回路自然循環(huán)的相似準(zhǔn)則可同時(shí)得到滿足。

比例分析；嚴(yán)重事故；自然循環(huán)；理論模型；壓力容器；外部冷卻

作為一項(xiàng)重要的嚴(yán)重事故緩解措施，熔融物堆內(nèi)滯留設(shè)計(jì)(IVR)在AP1000、CAP1400、華龍一號(hào)等三代核電技術(shù)中得到了廣泛應(yīng)用[1]。IVR設(shè)計(jì)主要利用堆腔內(nèi)冷卻水的自然循環(huán)實(shí)現(xiàn)壓力容器外部冷卻(ERVC)，從而帶走熔融堆芯的衰變熱，維持壓力容器的完整[2]。現(xiàn)有研究中，ERVC的有效性評(píng)估多通過(guò)程序計(jì)算確定，但因自然循環(huán)具有流量較小、參數(shù)變化敏感性強(qiáng)[3]、容易出現(xiàn)流動(dòng)不穩(wěn)定現(xiàn)象[4]等特點(diǎn)，計(jì)算結(jié)果一般要通過(guò)臺(tái)架進(jìn)行試驗(yàn)驗(yàn)證。實(shí)驗(yàn)臺(tái)架參數(shù)往往需要在原型的基礎(chǔ)上依照模化準(zhǔn)則進(jìn)行縮比，如在AP600和AP1000的安全審評(píng)中發(fā)揮重要作用的APEX試驗(yàn)臺(tái)架，該臺(tái)架就是根據(jù)Reyes等[5]對(duì)非能動(dòng)堆芯冷卻系統(tǒng)(PXS)內(nèi)自然循環(huán)現(xiàn)象模化分析的結(jié)果設(shè)計(jì)建造。

鑒于現(xiàn)階段針對(duì)ERVC自然循環(huán)現(xiàn)象的模化較少，本文以一維漂移流模型為基礎(chǔ)，對(duì)ERVC的自然循環(huán)過(guò)程進(jìn)行了等高度、等物性模化，通過(guò)推導(dǎo)和選擇關(guān)鍵無(wú)量綱準(zhǔn)則數(shù)，實(shí)現(xiàn)試驗(yàn)臺(tái)架與原型過(guò)程間最大程度的相似，為ERVC自然循環(huán)試驗(yàn)臺(tái)架設(shè)計(jì)提供理論依據(jù)。

1 系統(tǒng)描述及模化分析目標(biāo)

在核電廠發(fā)生嚴(yán)重事故時(shí)，在操作員的動(dòng)作下，堆腔淹沒(méi)系統(tǒng)中的冷卻水會(huì)迅速淹沒(méi)反應(yīng)堆堆腔，并在壓力容器下封頭外壁和保溫層之間的環(huán)形流道內(nèi)形成穩(wěn)定的兩相自然循環(huán)，實(shí)現(xiàn)壓力容器外部冷卻，帶走熔融堆芯的衰變熱，如圖1所示。

圖1 ERVC自然循環(huán)示意圖Fig.1 Schematic of natural circulation for ERVC

ERVC自然循環(huán)試驗(yàn)的目標(biāo)是模擬ERVC系統(tǒng)運(yùn)行時(shí)的自然循環(huán)循環(huán)流速，同時(shí)實(shí)現(xiàn)對(duì)RPV沸騰換熱的模擬。因?yàn)镋RVC系統(tǒng)的運(yùn)行狀態(tài)同時(shí)涉及兩相自然循環(huán)和局部換熱過(guò)程，所以本研究采用分級(jí)雙向比例(H2TS)方法進(jìn)行模化分析，H2TS方法是特別適用于綜合性實(shí)驗(yàn)臺(tái)架的一種先進(jìn)的分級(jí)模化方法[6]。模化分析主要分兩個(gè)層次進(jìn)行，流程如圖2所示。

第一層次為ERVC系統(tǒng)模化分析，主要以整個(gè)自然循環(huán)回路為控制體，通過(guò)對(duì)流動(dòng)換熱平衡方程的處理，獲得系統(tǒng)相似準(zhǔn)則；第二層次主要關(guān)注局部現(xiàn)象的相似要求，包括漂移流模型、RPV下封頭沸騰換熱等內(nèi)容，對(duì)系統(tǒng)相似準(zhǔn)則進(jìn)行補(bǔ)充和完善。

圖2 ERVC自然循環(huán)比例分析流程圖Fig.2 Scaling analysis flow diagram for natural circulation of ERVC

2 自然循環(huán)模化分析

2.1 兩相自然循環(huán)控制方程

ERVC系統(tǒng)運(yùn)行時(shí)，以壓力容器下封頭加熱面和流體作為控制體，其平衡方程如下[5]：

混合物質(zhì)量守恒方程

(1)

混合物動(dòng)量守恒方程

(2)

混合物能量方程

(3)

固體能量方程

(4)

考慮為穩(wěn)態(tài)自然循環(huán)，根據(jù)Ishii的研究，可忽略換熱系數(shù)的模擬，可合并固體能量方程和能量方程為新的能量方程[7]：

(5)

通過(guò)分別除以穩(wěn)定狀態(tài)自然循環(huán)下方程中各參數(shù)的值(以下標(biāo)0表示)，對(duì)方程(1)、(2)、(5)進(jìn)行無(wú)量綱化，以無(wú)量綱數(shù)表示的方程為：

混合物質(zhì)量方程

(6)

混合物動(dòng)量方程

(7)

混合物能量方程

(8)

上述方程中，τTP,o為流體混合物的弛豫時(shí)間，定義如下：

(9)

各無(wú)量綱數(shù)定義如下：漂移通量數(shù)：

(10)

Froude數(shù)：

(11)

密度數(shù)：

(12)

Friction數(shù)：

(13)

修正Stanton數(shù)：

(14)

熱源數(shù)：

(15)

焓升數(shù)：

(16)

2.2 穩(wěn)定自然循環(huán)參數(shù)

2.1節(jié)推導(dǎo)得到的無(wú)量綱數(shù)，均以穩(wěn)定狀態(tài)自然循環(huán)下的參數(shù)來(lái)表示，因此求解穩(wěn)定條件下ERVC自然循環(huán)各參數(shù)的值，其中最為重要的參數(shù)是自然循環(huán)流速。

對(duì)ERVC系統(tǒng)內(nèi)穩(wěn)定的自然循環(huán)流動(dòng)進(jìn)行簡(jiǎn)化，假設(shè)兩相間達(dá)到熱平衡和化學(xué)平衡，忽略由于蒸發(fā)和冷凝導(dǎo)致的對(duì)流加速等的影響，僅在形阻損失中考慮粘性效應(yīng)，采用一維均相流模型，穩(wěn)定狀態(tài)的控制方程為

(17)

(18)

式中：Lth為冷熱中心高度差，下標(biāo)in表示RPV流道的入口，i表示通道內(nèi)任意位置，l表示液相。

分別定義熱平衡含汽率和兩相平均密度為

(19)

(20)

式中：γ為汽化潛熱，下標(biāo)ls，gs分別表示飽和水和飽和蒸汽。

利用方程(18)和方程(19)，并令hsub=hls-hin，可得

(21)

方程(21)代入方程(20)可得

(22)

將方程(22)代入方程(17)，整理后可得

(23)

式中∶φ1、φ2、φ3為方程的系數(shù)，由熱功率、冷卻水物性、摩阻損失及形阻損失確定，其計(jì)算公式為

(24)

(25)

(26)

(27)

方程(23)即為ERVC系統(tǒng)自然循環(huán)流速的三次多項(xiàng)式，可用于求解循環(huán)流速。可知，ERVC自然循環(huán)流速的相似主要與入口過(guò)冷度、流體密度、下封頭加熱功率、幾何位置等因素有關(guān)。

方程(25)、(26)和(27)中，F(xiàn)為與回路阻力相關(guān)的參數(shù)，ERVC自然循環(huán)回路的阻力主要考慮入口局部損失、加熱段局部損失，上升段摩阻損失、出口局部損失及下降段摩阻損失及局部損失。

2.3 自然循環(huán)比例分析

Reyes證明，對(duì)于方程(23)，存在常數(shù)因子β使方程的系數(shù)嚴(yán)格滿足比例，保證模型與原型中的流體速度相似[8]，即：

(28)

(29)

(30)

(31)

(32)

式中∶下標(biāo)m和p分別代表模型值和原型值。

將方程(24)、(25)、(26)和(27)代入方程(29)～(32)中，當(dāng)ERVC穩(wěn)定運(yùn)行時(shí)，堆腔保溫層底部入口的水溫接近飽和，即過(guò)冷焓為0，于是得到：

(33)

(34)

(35)

由于ERVC系統(tǒng)的工作壓力接近常壓，所以試驗(yàn)臺(tái)架的工作壓力及溫度變化范圍可選定為與原型相同，并選定相同的工質(zhì)，即有

(36)

此外在試驗(yàn)臺(tái)架中，可使用孔板等阻力件，保證自然循環(huán)回路中各個(gè)部件的阻力特性相等，即有

(37)

式(33)、(35)、(36)和式(37)結(jié)合，可得

(38)

將式(38)代入(33)中，得

(39)

方程(21)與(39)結(jié)合，可得質(zhì)量含汽率比為

(40)

體積含汽率比為

(41)

兩相平均密度比為

(42)

可知在飽和情況下，當(dāng)運(yùn)行壓力和自然循環(huán)回路各部分阻力滿足相似時(shí)，質(zhì)量含汽率、體積含汽率和兩相平均密度都可滿足相似要求。

系統(tǒng)的時(shí)間比(即弛豫時(shí)間)變?yōu)?/p>

(43)

3 局部現(xiàn)象比例分析

3.1 兩相漂移速度

為保證漂移通量數(shù)相等，需評(píng)估漂移速度是否滿足要求，這里使用文獻(xiàn)[9-10]提出的方程，即

(44)

可知，漂移速度比主要與物性參數(shù)和速度有關(guān)，當(dāng)采用等高度等物性模擬時(shí)，漂移通量數(shù)準(zhǔn)則可通過(guò)速度相似和流體物性相似自動(dòng)滿足。

3.2 下封頭沸騰換熱

根據(jù)文獻(xiàn)[11]的研究，對(duì)RPV下封頭沸騰換熱過(guò)程的模擬，需要保證壓力容器外表面上的核態(tài)沸騰過(guò)程、兩相邊界層中的動(dòng)量傳遞過(guò)程以及經(jīng)最小間距的蒸汽排放過(guò)程三個(gè)傳遞過(guò)程相似。

對(duì)于使用與ERVC系統(tǒng)原型中相同的的工作流體，表征熱力學(xué)狀態(tài)的溫度范圍和壓力范圍也基本一致的試驗(yàn)臺(tái)架，表征沸騰過(guò)程的無(wú)量綱數(shù)Πb為

(45)

兩相邊界層中流動(dòng)的動(dòng)量傳遞無(wú)量綱數(shù)Πm為

(46)

因試驗(yàn)使用相同工質(zhì)，上述條件能夠滿足。

對(duì)于經(jīng)流動(dòng)通道最小間距的蒸汽排放過(guò)程，當(dāng)工質(zhì)相同，壁面熱流密度相似時(shí)，該過(guò)程的特征時(shí)間比Πv為

(47)

(48)

3.3 兩相流型轉(zhuǎn)變

結(jié)合ERVC系統(tǒng)各部分的特點(diǎn)，可將自然循環(huán)回路簡(jiǎn)化為豎直下降段、下封頭兩相段和豎直上升段三部分，其中，豎直下降段可視為單相絕熱流動(dòng)過(guò)程，下封頭兩相段視為兩相加熱流動(dòng)過(guò)程，豎直上升段可視為兩相絕熱流動(dòng)過(guò)程。文獻(xiàn)[12-13]的研究表明，ERVC中的兩相沸騰主要為過(guò)冷沸騰，在下封頭兩相段，兩相流動(dòng)主要為蒸汽相偏向下封頭外壁一側(cè)的不對(duì)稱流型。在豎直上升段中，主要以泡狀流為主，當(dāng)堆芯衰變熱功率較高時(shí)，可能出現(xiàn)彈狀流。因此ERVC系統(tǒng)兩相流型的轉(zhuǎn)變主要發(fā)生在豎直環(huán)狀間隙通道中，為不均勻泡狀流向均勻泡狀流或彈狀流的轉(zhuǎn)變。根據(jù)Reyes的研究[5]，流型轉(zhuǎn)變相似準(zhǔn)則主要與物性參數(shù)相關(guān)，等高等物性模擬時(shí)，可自動(dòng)滿足。

3.4 熱源模擬體相似

RPV下封頭模擬體通常選用相近或相同材料的加熱模擬體，即有

(49)

根據(jù)最小間距和加熱面寬度的比例，可以得到流道面積與加熱面積的比為

(50)

定義ds為加熱段厚度，將方程(38)、(49)和(50)代入Stanton數(shù)相似準(zhǔn)則有

(51)

式(51)表明，試驗(yàn)臺(tái)架的加熱段固體厚度和溫度應(yīng)與原型相同，加熱段固體厚度相似條件可通過(guò)幾何相似滿足；對(duì)于加熱段固體溫度相似條件，考慮到ERVC穩(wěn)態(tài)自然循環(huán)過(guò)程中流體溫度變化緩慢，當(dāng)RPV換熱過(guò)程相似時(shí)，換熱系數(shù)也比較穩(wěn)定，可認(rèn)為固體溫度相似可自動(dòng)滿足。

4 ERVC自然循環(huán)試驗(yàn)臺(tái)架比例

如前文分析所述，鑒于ERVC系統(tǒng)的運(yùn)行特點(diǎn)，ERVC自然循環(huán)臺(tái)架可采用等壓等溫同工質(zhì)的模擬方式，這樣可保證臺(tái)架與原型間工質(zhì)物性的一致。

一般來(lái)說(shuō)，試驗(yàn)臺(tái)架的幾何參數(shù)越接近原型的尺寸越能反映實(shí)際工況。對(duì)于自然循環(huán)現(xiàn)象，模型和原型的幾何比例過(guò)小可能會(huì)造成局部換熱現(xiàn)象的失真，因此本臺(tái)架將高度比選定為1∶1；對(duì)壓力容器下封頭的模擬，選擇工程上更容易實(shí)現(xiàn)的切片式二維試驗(yàn)段，半徑比為1∶1，流動(dòng)通道間距比為1∶1，綜合考慮用電功率、經(jīng)濟(jì)性等要求，切片寬度與原型的比wR為1∶100。

本文比例分析中涉及到的兩相自然循環(huán)、兩相流型轉(zhuǎn)變、下封頭換熱等過(guò)程均能完全模擬。對(duì)IVR工程應(yīng)用中的臨界熱流密度，因不是本文研究的重點(diǎn)，該比例分析只能對(duì)CHF現(xiàn)象中的核態(tài)沸騰過(guò)程進(jìn)行模擬。

綜上，根據(jù)前文結(jié)論，可得到等物性模擬時(shí)，試驗(yàn)裝置主要的縮比比例，如表1所示。

表1 ERVC自然循環(huán)試驗(yàn)?zāi)；Y(jié)果表

5 結(jié)論

本文以一維漂移流模型為基礎(chǔ)，對(duì)ERVC系統(tǒng)的自然循環(huán)過(guò)程進(jìn)行了比例分析，通過(guò)方程無(wú)量綱化，得到了試驗(yàn)臺(tái)架與原型ERVC兩相自然循環(huán)的相似準(zhǔn)則，主要結(jié)論如下：

1)漂移通量數(shù)、Froude數(shù)、密度數(shù)、Friction數(shù)、修正Stanton數(shù)、焓升數(shù)、熱源數(shù)、核態(tài)沸騰數(shù)、動(dòng)量傳遞數(shù)、間隙蒸汽排放過(guò)程特征數(shù)是表征ERVC整體試驗(yàn)臺(tái)架與核電廠原型間自然循環(huán)現(xiàn)象的相似準(zhǔn)則數(shù)，ERVC整體試驗(yàn)裝置的設(shè)計(jì)應(yīng)根據(jù)相似準(zhǔn)則數(shù)確定縮比比例和參數(shù)；

2)結(jié)合ERVC系統(tǒng)自然循環(huán)的特點(diǎn)，當(dāng)試驗(yàn)采用等壓等物性模擬時(shí)，F(xiàn)roude數(shù)、密度數(shù)、焓升數(shù)和熱源數(shù)相似準(zhǔn)則均能得到滿足，且所有臺(tái)架參數(shù)的比值均可以寫為冷熱芯高差比值函數(shù)的形式；

3)當(dāng)下封頭模擬體的半徑、熱流密度分布、加熱段厚、流道最小間距的位置和尺寸均與原型一致，且熱源的加熱面積與熱源處流通面積縮比比例一致時(shí)，修正Stanton數(shù)和表征RPV局部換熱過(guò)程的準(zhǔn)則數(shù)可得到滿足。

[1]徐紅, 周志偉. 大型先進(jìn)壓水堆熔融物堆內(nèi)滯留初步研究[J]. 原子能科學(xué)技術(shù), 2013, 47(6)∶ 969-974. XU Hong, ZHOU Zhiwei. Preliminary study on in-vessel retention in large-scale advanced PWR[J]. Atomic energy science and technology, 2013, 47(6)∶ 969-974.

[2]金越, 鮑晗, 劉曉晶, 等. 大功率先進(jìn)壓水堆IVR有效性評(píng)價(jià)分析[J]. 核動(dòng)力工程, 2015, 36(3)∶ 135-141. JIN Yue, BAO Han, LIU Xiaojing, et al. Assessment of in-vessel retention for advanced large size PWRs[J]. Nuclear power engineering, 2015, 36(3)∶ 135-141.

[3]周濤, 劉亮, 李宇, 等. 矩形窄縫自然循環(huán)密度波不穩(wěn)定特性實(shí)驗(yàn)研究[J]. 原子能科學(xué)技術(shù), 2015, 49(5)∶ 848-851. ZHOU Tao, LIU Liang, LI Yu, et al. Experiment study of density wave instability of natural circulation in narrow rectangle channel[J]. Atomic energy science and technology, 2015, 49(5)∶ 848-851.

[4]周濤, 李精精, 琚忠云, 等. 非能動(dòng)自然循環(huán)技術(shù)的發(fā)展與研究[J]. 核安全, 2013, 12(3)∶ 32-36. ZHOU Tao, LI Jingjing, JU Zhongyun, et al. The development and study on passive natural circulation[J]. Nuclear safety, 2013, 12(3)∶ 32-36.

[5]REYES Jr J N, HOCHREITER L. Scaling analysis for the OSU AP600 test facility (APEX)[J]. Nuclear engineering and design, 1998, 186(1/2)∶ 53-109.

[6]ZUBER N, WILSON G E, ISHII M, et al. An integrated structure and scaling methodology for severe accident technical issue resolution∶ Development of methodology[J]. Nuclear engineering and design, 1998, 186(1/2)∶ 1-21.

[7]盧冬華, 肖澤軍, 陳炳德. 壓水堆自然循環(huán)比例模化基本方程及相似準(zhǔn)則數(shù)的研究[J]. 核動(dòng)力工程, 2009, 30(3)∶ 72-84, 94. LU Donghua, XIAO Zejun, CHEN Bingde. Investigation on basic equations and scaling criteria of PWR natural circulation[J]. Nuclear power engineering, 2009, 30(3)∶ 72-84, 94.

[8]REYES Jr J N. Scaling single-state variable catastrophe functions∶ an application to two-phase natural circulation[J]. Nuclear engineering and design, 1994, 151(1)∶ 41-48.

[9]HIBIKI T, ISHII M. One-dimensional drift-flux model and constitutive equations for relative motion between phases in various two-phase flow regimes[J]. International journal of heat and mass transfer, 2003, 46(25)∶ 4935-4948.

[10]閆超星, 閻昌琪, 孫立成, 等. 基于漂移流模型的矩形通道內(nèi)泡狀流-彈狀流轉(zhuǎn)變準(zhǔn)則[J]. 原子能科學(xué)技術(shù), 2013, 47(10)∶ 1718-1722. YAN Chaoxing, YAN Changqi, SUN Licheng, et al. Transitional criterion from bubbly flow to slug flow in vertical rectangular channels based on drift-flux model[J]. Atomic energy science and technology, 2013, 47(10)∶ 1718-1722.

[11]CHEUNG F, LIU Y. Critical heat flux experiments to support in-vessel retention feasibility study for an evolutionary advanced light water reactor, Technical Results 1003101[R]. Palo Alto, CA∶ EPRI, 2001.

[12]CHU T Y, BENTZ J H, SLEZAK S E, et al. Ex-vessel boiling experiments∶ laboratory- and reactor-scale testing of the flooded cavity concept for in-vessel core retention Part II∶ Reactor-scale boiling experiments of the flooded cavity concept for in-vessel core retention[J]. Nuclear engineering and design, 1997, 169(1/2/3)∶ 89-99.

[13]李永春, 劉曉晶, 程旭. ERVC流道內(nèi)兩相局部分布實(shí)驗(yàn)研究[J]. 核動(dòng)力工程, 2014, 35(5)∶ 153-156. LI Yongchun, LIU Xiaojing, CHENG Xu. Experimental study on local two phase distributions in ERVC flow channel[J]. Nuclear power engineering, 2014, 35(5)∶ 153-156.

Scaling analysis of natural circulation for external cooling of a reactor vessel

LIU Yusheng1，LI Congxin1，WEN Lijing1，TAN Sichao2，ZHANG Pan1

(1. Nuclear and Radiation Safety Center, Ministry Environmental Protection, Beijing 100084, China; 2. Key Discipline Laboratory of Nuclear Safety and Simulation Technology, Harbin Engineering University, Harbin 150001, China)

The two-phase natural circulation phenomenon occurring in the gap between the external wall of a pressure vessel and the insulation is a key process for external reactor vessel cooling (ERVC) and core corium in-vessel retention in case of a severe accident. To obtain the similarity criteria between ERVC integral performance test facility and a prototype nuclear power plant, the hierarchical two-tiered scaling (H2TS) analysis method was used to conduct a scaling analysis for an ERVC natural circulation process. The H2TS analysis was based on the conservation equation and the one-dimensional drift flow model. The similarity criteria between the ERVC natural circulation test facility and the prototype was attained. In addition, the steady-state solution for natural circulation was given. The results show that the natural circulation phenomenon in ERVC can be simulated in an integral facility of the same operating pressure and the same working fluid, and the similarity criteria of the Froude number, density number, enthalpy rise quantity, heat source volume, etc., can all be met. There is no obvious distortion due to the transition of scaling ratios. The two-dimensional slice type of the heat section with the same thickness, radius, and heat flux distribution was applied to simulate the lower end socket of the pressure vessel. When the minimum interval, location and amount of flow pass are consistent with the prototype, the ratio of similitude of the revised Stanton number is 1, and the similarity criteria for boiling heat transfer and natural circulation can be met simultaneously.

scaling analysis; severe accident; natural circulation; theoretical model; external cooling; reactor vessel

2015-11-24.

日期：2016-11-14.

國(guó)家科技重大專項(xiàng)(2015ZX06002007).

劉宇生(1986-), 男, 工程師； 譚思超(1979-), 男, 教授, 博士生導(dǎo)師.

譚思超，E-mail∶ 13091441949@163.com.

http∶//www.cnki.net/kcms/detail/23.1390.u.20161114.1036.012.html

TL333

A

1006-7043(2017)02-0318-06

劉宇生, 李聰新, 溫麗晶，等. 壓力容器外部冷卻自然循環(huán)比例分析[J]. 哈爾濱工程大學(xué)學(xué)報(bào)， 2017, 38(2)∶ 318-323. LIU Yusheng，LI Congxin，WEN Lijing， et al. Scaling analysis of natural circulation for external cooling of a reactor vessel[J]. Journal of Harbin Engineering University, 2017, 38(2)∶ 318-323.

DOI∶10.11990/jheu.201511057