一體化汽-氣穩(wěn)壓器動態(tài)特性研究

吳 磊，賈海軍，劉 洋，馬喜振

（清華大學(xué) 核能與新能源技術(shù)研究院 先進反應(yīng)堆工程與安全教育部重點實驗室，北京 100084）

一體化壓水反應(yīng)堆將主換熱器或蒸汽發(fā)生器、主循環(huán)泵、穩(wěn)壓器等主回路核心裝置布置在反應(yīng)堆壓力容器內(nèi)。此種布置消除了傳統(tǒng)的冷卻劑回路管道，結(jié)構(gòu)緊湊，避免大LOCA，大幅提高了反應(yīng)堆的安全性。汽-氣穩(wěn)壓器具有自穩(wěn)壓特性，其因結(jié)構(gòu)簡單、省去了加熱和噴淋裝置、防止容積沸騰等特點而在一體化反應(yīng)堆以及其他類型中小型反應(yīng)堆中得到應(yīng)用。清華大學(xué)自主開發(fā)設(shè)計的NHR 系列等即是采用汽-氣穩(wěn)壓器的一體化壓水反應(yīng)堆的典型代表［1］。近年來，NHR 系列的熱工參數(shù)逐步提升，以適用于城市集中供熱、生產(chǎn)工業(yè)蒸汽、海水淡化等非發(fā)電領(lǐng)域。穩(wěn)壓器的穩(wěn)壓能力和動態(tài)特性對整個反應(yīng)堆一回路系統(tǒng)的安全性有重要影響，需要對其進行深入研究。

當(dāng)前，對汽-氣穩(wěn)壓器的動態(tài)特性進行了一些實驗和理論研究。Leonard 等［2］對汽-氣穩(wěn)壓器涌入過程進行實驗研究，重點考察了不同種類非凝結(jié)性氣體，不同的氣體組分以及分層現(xiàn)象等對穩(wěn)壓特性的影響。Hassan等［3］進一步改進和發(fā)展了Relap5程序中相關(guān)計算模型。2006年，Ambrosini等［4］對汽-氣穩(wěn)壓器內(nèi)存在的傳熱傳質(zhì)理論進行總結(jié)并指出，可基于不同的表達式來歸納多種形式的傳熱傳質(zhì)比擬方法。Kim 等［5］通過實驗的方法考察了在高壓下有非凝結(jié)性氣體存在時的冷凝傳熱系數(shù)，并發(fā)展了高壓下的汽-氣穩(wěn)壓器模型［6］。

本文主要是建立汽-氣穩(wěn)壓器計算模型，通過與汽-氣穩(wěn)壓器涌入實驗數(shù)據(jù)進行比較，驗證模型的可靠性，同時基于此模型對采用一體化汽-氣穩(wěn)壓器的NHR-Ⅱ系統(tǒng)的穩(wěn)壓動態(tài)特性進行研究，揭示NHR-Ⅱ系統(tǒng)功率調(diào)節(jié)過程中壓力和溫度跟隨滯后對系統(tǒng)壓力變化的影響機理。

1 理論模型

汽-氣穩(wěn)壓器通過壓縮非凝結(jié)性氣體和蒸汽的混合物進行穩(wěn)壓，一般布置在一體化壓水反應(yīng)堆頂端。汽-氣穩(wěn)壓器按照工質(zhì)特性分為3個區(qū)：氣相混合區(qū)、氣液相界面區(qū)和液相區(qū)，如圖1所示。上部的氣相混合區(qū)主要包含非凝結(jié)性氣體和蒸汽的混合物，由氣體和水蒸氣的分壓決定非凝結(jié)性氣體與蒸汽的含量比值；下部的液相區(qū)主要是對應(yīng)蒸汽分壓的飽和水；氣液界面是氣液混合區(qū)域。在動態(tài)條件下，冷凝、擴散、蒸發(fā)、自然對流等共同決定汽-氣穩(wěn)壓器內(nèi)部物理現(xiàn)象的復(fù)雜性。

圖1 一體化汽-氣穩(wěn)壓器區(qū)域劃分Fig.1 Three regions in integrated gas-steam pressurizer

本文數(shù)學(xué)模型是通過求解不同的場方程來獲得整個系統(tǒng)描述需要的壓力（p）、相比內(nèi)能（Ug，Uf）、空泡份額（αg）、相速度（vg，vf）、非凝結(jié)性成分含量（Xg）等變量，獨立變量為時間（t）和距離（x）。提出以下假設(shè)：假定在汽-氣混合區(qū)域內(nèi)非凝結(jié)性氣體和水蒸氣具有同樣的速度和溫度，即vn=vs、Tn=Ts，其中vn和vs分別為非凝結(jié)性氣體和水蒸氣速度，Tn和Ts分別為非凝結(jié)性氣體和水蒸氣溫度，并認為蒸汽和非凝結(jié)性氣體的混合物仍處于熱力學(xué)非平衡態(tài)條件。

基于氣相和液相基本的質(zhì)量、動量和能量守恒方程，在有非凝結(jié)性氣體的存在下，需要添加或修正傳熱傳質(zhì)項，并對相關(guān)物理模型進行建模計算。主要包括以下內(nèi)容：

質(zhì)量方程添加項：當(dāng)有非凝結(jié)性氣體存在時，質(zhì)量守恒方程會隨著非凝結(jié)性氣體的種類N 增加N-1個，形式如下：

其中：ρg 為氣相密度；Xn為氣相中非凝結(jié)性氣體的含量；A 為通道橫截面積。

能量方程添加項：當(dāng)有非凝結(jié)性氣體存在時，能量守恒方程中添加了界面直接加熱項Qgf（式（3））。界面相關(guān)項的定性溫度采用主流的蒸汽分壓所對應(yīng)的飽和溫度，而非局部對應(yīng)的界面蒸汽分壓對應(yīng)的飽和溫度。

動量守恒方程：基本形式未改變，只是在有非凝結(jié)性氣體存在時將所有場方程中蒸汽物性參數(shù)采用蒸汽／非凝結(jié)性氣體混合物物性參數(shù)代替。

界面冷凝的處理：通過Colburn-Hougen迭代方法計算蒸汽和液相表面溫度，應(yīng)用此溫度于熱流密度關(guān)系式；在氣體湍流狀態(tài)或非凝結(jié)性氣體存在時，采用Vierow 和Schrock提出的基于Nusselt層流冷凝系數(shù)的修正因子。

相界面質(zhì)量的交換處理：主要是用氣相分壓參數(shù)代替系統(tǒng)總壓對應(yīng)參數(shù)。相界面交換的質(zhì)量為：其中：Γig為相界面交換的質(zhì)量；ps為水蒸氣分壓；Ts（ps）為水蒸氣分壓對應(yīng)的飽和溫度；Tg和Tf分別為氣相溫度和液相溫度；h＊g和h＊f分別為氣相焓和液相焓；Hif為界面與液相間傳熱系數(shù)。

相界面能量的交換處理：有非凝結(jié)性氣體存在時主流之間的熱量交換為：

其中：pn為非凝結(jié)性氣體分壓；Hgf為氣液相間直接傳熱系數(shù)。

當(dāng)非凝結(jié)性氣體存在時，界面的熱量傳遞和質(zhì)量傳遞都需進行修正，能量方程中用到h＊g和Γig，采用蒸汽分壓代替總壓。

混合相物性參數(shù)處理：有非凝結(jié)性氣體存在時，水蒸氣的參數(shù)直接查水蒸氣參數(shù)表，將非凝結(jié)性氣體當(dāng)作理想氣體處理，采用修正的Gibbs-Dalton混合物計算方法。

其中：Us和Un分別為水蒸氣和非凝結(jié)性氣體比內(nèi)能；υs和υn分別為水蒸氣和非凝結(jié)性氣體比體積；υg為氣相比體積。

擴散系數(shù)的計算采用Fuller方法［7］：

其中：D1，2為 擴 散 系 數(shù)；M1，2=2（1／M1+1／M2）-1，為平均相對分子質(zhì)量；Σv為原子擴散體積。

2 模型驗證

選擇麻省理工學(xué)院進行的汽-氣穩(wěn)壓器內(nèi)涌入波動實驗進行模型的驗證分析［2］，主要考察汽-氣穩(wěn)壓器壓力響應(yīng)特性。N2作為非凝結(jié)性氣體且質(zhì)量分數(shù)分別為0%、10%、20%，實驗條件列于表1。

表1 驗證實驗條件Table 1 Verification experiment conditions

圖2 驗證性實驗系統(tǒng)節(jié)點劃分Fig.2 Nodalization of experimental system

實驗裝置的節(jié)點劃分如圖2所示。實驗水箱用控制體018表示，劃分為18個小控制體，并通過管道012與時間控制體128連接，邊界條件根據(jù)實驗數(shù)據(jù)以流量為時間函數(shù)形式給出。初始化條件列于表1。由于N2和水蒸氣（H2O）的相對分子質(zhì)量相差不大，初始化時將N2和水蒸氣視為平衡態(tài)的均勻混合性氣體處理。計算結(jié)果如圖3所示。

圖3a為N2質(zhì)量分數(shù)為0%時的壓力變化。過程維持時間約34.8s，計算得到系統(tǒng)壓力從初始的0.53 MPa隨過冷水的涌入逐步上升到峰值0.64 MPa，與實驗結(jié)果基本一致，峰值相對誤差為0.94%。N2質(zhì)量分數(shù)為10%和20%的模擬和實驗結(jié)果如圖3b、c所示。計算得到系統(tǒng)壓力從0.53 MPa分別升至0.79 MPa和0.91 MPa 的峰值，平均上升速率分別為0.007 3和0.011 MPa／s。與實驗結(jié)果對比可知，峰值相對誤差分別為8.78%和12.89%，均在15%以內(nèi)，具有較好的一致性。

圖3 系統(tǒng)壓力變化Fig.3 Pressure response with nitrogen

由圖3可知，計算模型可對非凝結(jié)性氣體存在條件下的穩(wěn)壓特性進行較好的預(yù)測，誤差值在工程應(yīng)用允許范圍內(nèi)。隨N2質(zhì)量分數(shù)的增加，系統(tǒng)壓力峰值和壓力增加速率均增加，這說明非凝結(jié)性氣體的存在抑制熱量的散失或傳遞，且這種抑制作用隨著非凝結(jié)性氣體質(zhì)量分數(shù)的增加而增強。

3 NHR-Ⅱ汽-氣穩(wěn)壓動態(tài)特性

根據(jù)前文研究可知，理論模型能對汽-氣穩(wěn)壓器穩(wěn)壓動態(tài)特性進行較好的預(yù)測?；谏鲜鲇嬎隳Ｐ?，進一步研究采用一體化汽-氣穩(wěn)壓器和全功率自然循環(huán)的NHR-Ⅱ系統(tǒng)穩(wěn)壓動態(tài)特性。

3.1 系統(tǒng)參數(shù)及動態(tài)過程介紹

主要考察NHR-Ⅱ系統(tǒng)功率線性迅速增加等反應(yīng)堆系統(tǒng)極端運行工況。功率調(diào)節(jié)具體過程如下：功率從0s開始保持30%直到1 000s，使得整個NHR-Ⅱ系統(tǒng)穩(wěn)定，然后在60s內(nèi)線性增加到100%，功率保持穩(wěn)定至5 000s，整個過程結(jié)束。整個動態(tài)過程中，一、二回路采用運行控制方案為主回路平均溫度不變的方案，即主回路堆芯進出口溫度平均值保持不變，同時保證二回路壓力不變，而二回路流量增大，進口溫度下降。具體參數(shù)列于表2。

表2 NHR-Ⅱ動態(tài)調(diào)節(jié)過程熱工參數(shù)Table 2 Operation parameters and conditions in NHR-Ⅱ

3.2 系統(tǒng)節(jié)點劃分

系統(tǒng)節(jié)點劃分如圖4、5所示?？刂企w014表示堆芯流通通道，與之相鄰的熱構(gòu)件是堆芯燃料棒；控制體026 表示穩(wěn)壓空間，即一體化汽-氣穩(wěn)壓器，劃分為10個控制體；控制體024為分支，分別連接穩(wěn)壓器、主換熱器、上升段；013、015、017共同表示節(jié)流件，進行系統(tǒng)阻力控制和自然循環(huán)流量調(diào)節(jié)；部件PHE 表示主換熱器；二回路的進出口分別由時間控制體連接，給出二次側(cè)流量、溫度、壓力等參數(shù)；初始化條件列于表2。

圖4 NHR-Ⅱ系統(tǒng)節(jié)點劃分Fig.4 Nodalization of NHR-Ⅱsystem

圖5 穩(wěn)壓器節(jié)點劃分Fig.5 Nodalization of pressurizer

3.3 結(jié)果分析

1）系統(tǒng)壓力變化特性

由于NHR-Ⅱ采用全功率自然循環(huán)和一體化汽-氣穩(wěn)壓器設(shè)計，在動態(tài)運行過程中對壓力跟隨變化特性具有較高的要求。在功率快速提升過程中，重點關(guān)注系統(tǒng)壓力的變化。

（1）壓力變化影響因素

在功率調(diào)節(jié)的動態(tài)過程中，穩(wěn)壓器下部的液體體積和穩(wěn)壓器內(nèi)的溫度變化，主要決定著系統(tǒng)壓力的變化。選定涌入或涌出穩(wěn)壓器內(nèi)流體質(zhì)量和穩(wěn)壓器內(nèi)溫度進行定量分析。

涌入穩(wěn)壓器流體直接壓縮穩(wěn)壓空間內(nèi)的汽-氣混合物，改變其內(nèi)部溫度以及物質(zhì)組分分布，對壓力的變化具有重要的影響。圖6為涌入穩(wěn)壓器流體質(zhì)量隨時間的變化。

圖6 涌入穩(wěn)壓器流體質(zhì)量隨時間的變化Fig.6 Mass of liquid surged into pressurizer vs time

初始階段，功率保持在30%，涌入穩(wěn)壓器的液相質(zhì)量約為0kg。在功率提升初始階段（1 000～1 013s），功 率 從30%線 性 增 加 到45%，涌入穩(wěn)壓器的流體仍基本為0kg；在接下來功率提升的47s（1 014～1 060s），流體涌入穩(wěn)壓器內(nèi)部的質(zhì)量近線性增加，且涌入穩(wěn)壓器的質(zhì)量最大值出現(xiàn)在1 060s，約為1 625kg；在功率穩(wěn)定在100%階段，涌入穩(wěn)壓器的液相質(zhì)量逐漸減少，最終穩(wěn)定在某一數(shù)值。這表明，在此階段開始有部分液相涌出穩(wěn)壓器，主要是由于二回路冷卻能力的增強使主回路平均溫度降低和壓縮反沖作用等原因造成的。

穩(wěn)壓器內(nèi)部氣相和液相區(qū)域溫度直接影響著穩(wěn)壓空間液相的蒸發(fā)和蒸汽的冷凝，從而對系統(tǒng)壓力產(chǎn)生重要影響。圖7為穩(wěn)壓器不同區(qū)域內(nèi)溫度隨時間的變化。

在30%功率穩(wěn)定階段，穩(wěn)壓器內(nèi)溫度的分布保持穩(wěn)定，且存在一定的溫差。在功率從30%提升至100%階段，穩(wěn)壓器下部液相區(qū)域溫度提升的速率明顯大于穩(wěn)壓器混合相區(qū)域，這主要由不同的傳熱機理所造成。穩(wěn)壓器下部液相區(qū)域溫度的變化主要是堆芯出口溫度對流傳熱遷移引起的；穩(wěn)壓器混合相空間溫度變化主要是通過與穩(wěn)壓器下部液相區(qū)域的熱傳導(dǎo)以及混合相空間內(nèi)微弱的自然對流引起的，熱量傳遞效率遠小于前者，從而穩(wěn)壓器下部液相區(qū)域溫度提升的速率明顯大于穩(wěn)壓空間液相區(qū)。在功率保持100%穩(wěn)定階段，穩(wěn)壓器內(nèi)溫度持續(xù)升高，最終穩(wěn)定在最大值。由于考慮了穩(wěn)壓器外壁的散熱，穩(wěn)壓器在豎直方向由下到上溫度逐漸減小，存在溫度梯度。

圖7 穩(wěn)壓器不同區(qū)域溫度Fig.7 Temperature in different regions of pressurizer

（2）系統(tǒng)壓力變化特性及壓力跟隨滯后現(xiàn)象

圖8為功率調(diào)節(jié)過程中系統(tǒng)壓力的變化。當(dāng)系統(tǒng)功率增大到100%后，系統(tǒng)壓力仍會不斷上升，但增加速率減緩，直到壓力穩(wěn)定，且最終的穩(wěn)定值即為系統(tǒng)壓力的最大值。

圖8 整個動態(tài)過程系統(tǒng)壓力變化Fig.8 Pressure change in whole progress

系統(tǒng)壓力在初始的1 000s內(nèi)達到穩(wěn)定值，約為5.56MPa；從1 000s開始，系統(tǒng)壓力隨功率的提升不斷增大，最終穩(wěn)定值（亦為最大值）約為6.85MPa。整個過程壓力增加1.29 MPa，壓力增加速率為0.000 19 MPa／s。在功率提升階段（1 000～1 060s），壓力從5.56 MPa上升至6.08MPa，占整個調(diào)節(jié)過程壓力增加值的40.4%；壓力增長的平均速率為0.009 9MPa／s，此時壓力增加速率遠大于整個調(diào)節(jié)過程的變化速率。這表明壓力的增加大部分是在恒定功率（100%）下完成，并非在壓力提升過程（1 000～1 060s）完成，壓力的變化整體滯后于功率的提升。

為進一步研究系統(tǒng)壓力動態(tài)變化特性，選定功率提升過程進行分析。選定考察時間范圍為980～1 100s，此過程包含功率提升階段（1 000～1 060s）。在選定階段，壓力變化的斜率逐漸增加，并定義：在功率提升階段，小于平均壓力變化斜率20%的階段為壓力變化滯后階段。

根據(jù)壓力變化滯后的定義，壓力變化滯后現(xiàn)象出現(xiàn)在功率提升的初始13s（1 000～1 013s），此階段功率從30%增加到45%，系統(tǒng)壓力變化速率約為整個過程平均速率的20%，系統(tǒng)壓力從5.56 MPa增加到5.58 MPa，壓力增加值僅占整個功率調(diào)節(jié)階段壓力增加值的3.25%；緊隨壓力變化滯后階段，功率增加到100%，壓力變化斜率增大，壓力會隨著功率的增加而近線性增大，此時系統(tǒng)壓力變化的跟隨性較好，認為是壓力近線性跟隨階段（1 014～1 060s）；下一階段功率保持100%不變，壓力不斷增加但增加率逐漸減小，此階段被定義為功率提升后壓力跟隨階段，如圖9所示。

圖9 選定階段系統(tǒng)壓力和功率變化Fig.9 Change of pressure and power in selected stage

根據(jù)對系統(tǒng)壓力變化各影響因素分析，將系統(tǒng)壓力變化的三階段（壓力變化滯后階段、壓力近線性跟隨階段、功率調(diào)節(jié)后壓力跟隨階段）劃分為4個過程：

第1 過程t1：功率從30%提升至45%，NHR-Ⅱ系統(tǒng)穩(wěn)壓器下部液相體積變化很小，同時加熱段出口溫度尚未遷移到穩(wěn)壓空間底部液相區(qū)域，穩(wěn)壓空間內(nèi)溫度基本保持不變，從而系統(tǒng)壓力基本保持不變。此過程持續(xù)13s，對應(yīng)壓力變化滯后階段（1 000～1 013s）。

第2 過 程t2：功 率 從45% 線 性 增 加 到100%，主回路液相持續(xù)被加熱，平均溫度升高，涌入穩(wěn)壓空間的液相體積逐步變大并達到最大值。穩(wěn)壓空間底部液相區(qū)域溫度從537.5K 增加到544.6K，穩(wěn)壓空間的溫度也增加4K。此過程系統(tǒng)壓力近線性增加，對應(yīng)壓力近線性跟隨階段（1 014～1 060s）。

第3過程t3：此時功率穩(wěn)定在100%，涌入穩(wěn)壓空間的液相總量逐步減小，即有液相從穩(wěn)壓器內(nèi)涌出，主要歸因于二回路的冷卻增強和壓縮反沖作用。同時，穩(wěn)壓空間的溫度快速增加。此過程系統(tǒng)壓力持續(xù)增加，但增加速率明顯低于第2過程。

第4過程t4：此過程功率恒定，整個系統(tǒng)穩(wěn)壓液相區(qū)域的溫度分布穩(wěn)定，無液相涌入或涌出穩(wěn)壓空間。同時，在穩(wěn)壓器中存在熱傳導(dǎo)和微弱的自然對流，穩(wěn)壓空間的溫度緩慢增加。此過程整個系統(tǒng)的壓力仍增加并最終趨于穩(wěn)定，但增加速率較第3過程緩慢。

綜上可知，壓力變化滯后發(fā)生在第1過程，壓力近線性跟隨階段發(fā)生在第2過程，且第1、2過程屬于功率調(diào)節(jié)階段；第3、4 過程屬于功率調(diào)節(jié)后壓力跟隨階段。

2）系統(tǒng)溫度變化特性

（1）溫度滯后現(xiàn)象

在功率調(diào)節(jié)過程中，堆芯出口溫度最先變化，且能實時跟隨功率提升。堆芯出口溫度變化需經(jīng)上升段傳遞到穩(wěn)壓空間下部液相區(qū)域，此過程稱為溫度遷移滯后（圖10）。

圖10中，兩條曲線分別代表堆芯出口溫度變化和穩(wěn)壓空間下部液相區(qū)域溫度變化，在時間軸上呈現(xiàn)的是“曲線平移”。平移時間thys可通過下式估算得到：

其中：Lu為上升段長度；vu為上升段平均流速。本例中平移時間約為35s。

系統(tǒng)溫度變化除溫度遷移滯后外，還存在溫度傳遞滯后，主要發(fā)生在穩(wěn)壓器內(nèi)。在一體化汽-氣穩(wěn)壓器中，當(dāng)熱量被傳遞到穩(wěn)壓器底部液相區(qū)域后，需通過熱傳導(dǎo)和穩(wěn)壓空間微弱的自然對流將熱量傳遞到整個穩(wěn)壓器內(nèi)，稱為溫度傳遞滯后。由于熱量的傳導(dǎo)速率遠小于熱量的對流導(dǎo)熱，從而在穩(wěn)壓空間內(nèi)溫度的變化遠滯后于系統(tǒng)溫度的變化。正是由于這兩種溫度滯后現(xiàn)象共同組成系統(tǒng)溫度滯后，影響系統(tǒng)壓力的變化。

圖10 溫度遷移滯后Fig.10 Temperature transport hysteresis

（2）堆芯出口過冷度變化

由于壓力滯后階段的存在，影響堆芯出口過冷度變化并威脅系統(tǒng)安全。系統(tǒng)在一定堆芯出口過冷度下運行，堆芯出口溫度隨功率迅速提升而增加，在壓力變化滯后階段，系統(tǒng)總壓保持不變，對應(yīng)系統(tǒng)飽和溫度保持不變，從而造成堆芯出口過冷度大幅減小甚至出現(xiàn)局部沸騰，威脅系統(tǒng)的安全性。圖11 為NHR-Ⅱ系統(tǒng)功率快速提升情況下堆芯出口過冷度的變化。堆芯過冷度ΔTc算式為：

其中：Ts（p）為系統(tǒng)壓力下對應(yīng)飽和溫度；Tco為堆芯出口溫度。

圖11 NHR-Ⅱ系統(tǒng)功率快速提升過程堆芯出口過冷度變化Fig.11 Change of core outlet subcooled temperature during power regulation stage

在NHR-Ⅱ系統(tǒng)功率快速提升過程中，堆芯出口過冷度隨著功率增加立即減小，在壓力滯后階段從7.3K 減小到最小值2.7K，隨后逐步增加。在NHR-Ⅱ系統(tǒng)功率迅速提升過程中，堆芯出口過冷度最小值為2.7 K，表明NHR-Ⅱ系統(tǒng)是安全的。

4 結(jié)論

本文研究了一體化汽-氣穩(wěn)壓器動態(tài)穩(wěn)壓特性，主要結(jié)論如下：

1）計算模型能較好地預(yù)測一體化汽-氣穩(wěn)壓器動態(tài)特性，誤差控制在可接受范圍內(nèi)，可滿足工程應(yīng)用要求。

2）穩(wěn)壓空間混合相的壓縮和穩(wěn)壓器內(nèi)熱量的傳遞共同決定NHR-Ⅱ系統(tǒng)中壓力的變化。這兩個因素先后在功率提升過程占主導(dǎo)作用，致使整個調(diào)節(jié)過程經(jīng)歷壓力滯后、壓力近線性變化、壓力跟隨等3個階段。壓力滯后出現(xiàn)在功率調(diào)節(jié)初始階段，且此階段堆芯過冷度會急劇減小。

3）系統(tǒng)溫度變化滯后主要是由于系統(tǒng)不同區(qū)域傳熱機理差異引起的，并對系統(tǒng)壓力的變化有重要影響。溫度遷移滯后主要發(fā)生在從堆芯出口到穩(wěn)壓器下部液相區(qū)域，由對流傳熱決定；溫度傳遞滯后主要發(fā)生在穩(wěn)壓器內(nèi)，由穩(wěn)壓器下部液相區(qū)域?qū)岷头€(wěn)壓混合相空間微弱的自然對流決定。

4）在功率提升過程中的壓力滯后階段，堆芯出口過冷度迅速下降并達到最小值，可能會威脅整個系統(tǒng)的安全性。在NHR-Ⅱ系統(tǒng)功率突升過程中，堆芯出口過冷度在安全范圍內(nèi)，系統(tǒng)安全。

［1］ 張亞軍，王秀珍.200 MW 低溫供熱堆研究進展及產(chǎn)業(yè)化發(fā)展前景［J］.核動力工程，2003，24（2）：180-188.ZHANG Yajun，WANG Xiuzhen.Application of the 200MW low temperature nuclear heating rector［J］.Nuclear Power Engineering，2003，24（2）：180-188（in Chinese）.

［2］ LEONARD M T.The effect of a non-condensable gas on pressurizer in-surge transientients［R］.US： Massachusetts Institute of Technology，1983.

［3］ HASSAN Y A，BANERJEC S.Implementation of a non-condensable mode in Relap5／mod3［J］.Nuclear Engineering and Design，1996，162：281-300.

［4］ AMBROSINI W，F(xiàn)ORGIONE N，MANFREDINI A，et al.On various forms of the heat and mass transfer analogy：Diffusion and application to condensation experiments［J］.Nuclear Engineering and Design，2006，236：1 013-1 027.

［5］ KIM J W，LEE Y G，AHN H K，et al.Condensation heat transfer characteristic in the presence of non-condensable gas on natural convection at high pressure［J］.Nuclear Engineering and Design，2009，239：688-698.

［6］ KIM T W，KIM J W，PARK G C.Development of non-equilibrium pressurizer model with noncondensable gas［J］.Nuclear Engineering and Design，2006，236：375-384.

［7］ POLING B E，PRAUSNITZ J M，O’CONNELL J P.The properties of gases and liquids：Ch.11［M］.5th ed.New York：WCB／McGraw-Hill，2001.