汪子迪，曹 臻，劉曉晶，程 旭

(上海交通大學 核科學與工程學院，上海 200240)

?

超臨界水實驗回路熱工水力分析

汪子迪，曹 臻，劉曉晶，程 旭

(上海交通大學 核科學與工程學院，上海 200240)

作為超臨界水堆失水事故分析的關(guān)鍵現(xiàn)象，跨臨界過程(即超臨界水堆的壓力從超臨界狀態(tài)降到次臨界狀態(tài)22.1 MPa以下)受到國內(nèi)外的關(guān)注。上海交通大學的超臨界流體多功能實驗回路(SWAMUP)計劃對這一泄壓過程進行實驗研究。為確保該實驗裝置在實驗過程中的安全性能，采用系統(tǒng)程序ATHLET-SC對該實驗回路進行預計算分析，主要針對該系統(tǒng)在泄壓跨臨界過程中的熱工水力參數(shù)，包括系統(tǒng)壓力、冷卻劑流量、加熱棒壁面溫度等展開計算，并討論一些重要參數(shù)如泄壓速度、加熱棒加熱功率等對計算結(jié)果的影響。計算結(jié)果表明，修改后的ATHLET-SC程序可模擬跨臨界瞬態(tài)過程，在實驗過程中，加熱棒壁面溫度不會超過設計上限溫度，然而，回路中換熱器的內(nèi)外最高壓差將會達6 MPa，這一點需在實驗中特別考慮。

SWAMUP實驗裝置；跨臨界瞬態(tài)；ATHLET-SC程序；預計算

超臨界水堆(SCWR)為6種第4代核反應堆中的唯一水冷反應堆，各方面的特性在國內(nèi)外得到廣泛研究[1-2]。根據(jù)第4代核能系統(tǒng)國際論壇(GIF)超臨界水堆系統(tǒng)籌劃委員會制定的技術(shù)路線和研發(fā)計劃，將于2020年建設一座小型(≤150 MWt)超臨界水冷原型堆[3-4]。歐盟第七框架協(xié)議資助項目“超臨界水冷堆燃料性能驗證實驗(SCWR-FQT)”的主要研究內(nèi)容為超臨界水冷回路及安全和輔助系統(tǒng)的小型燃料組件的安全許可申請所需的設計分析和驗證，并以此作為建設原型堆的基礎(chǔ)[5]?？缗R界泄壓過程是超臨界水堆失水事故分析的關(guān)鍵。目前很多系統(tǒng)分析程序可模擬這一泄壓過程，但依然缺少對這些跨臨界模型的實驗驗證。作為SCWR-FQT項目的參與單位，上海交通大學的超臨界流體多功能實驗回路(SWAMUP)[6]承擔了對該回路的安全許可申請的驗證工作之一，即對跨臨界泄壓過程進行實驗研究。本工作應用修改后的系統(tǒng)分析程序ATHLET-SC[7-8]對該回路進行預計算，對該系統(tǒng)泄壓過程中的熱工水力參數(shù)，包括系統(tǒng)壓力、冷卻劑流量、加熱棒壁面溫度等進行分析。此外，本文還將討論一些重要的參數(shù)，包括泄壓速度、加熱棒加熱功率以及臨界熱流密度(CHF)等對結(jié)果的影響。

1 系統(tǒng)介紹

SWAMUP實驗裝置如圖1所示，該系統(tǒng)為一開式回路，包括1個主實驗回路和2個冷卻二回路。主回路主要包括柱塞泵、換熱器、預熱器、實驗段和混合水箱。柱塞泵提供穩(wěn)定的流量，使整個回路的壓力保持在25 MPa左右。水箱中的水經(jīng)由柱塞泵出來的冷卻劑分為兩條主支路，其中1條支路被旁排而不經(jīng)實驗段加熱，以用來調(diào)整另一流經(jīng)實驗段支路的流量達到設計要求。經(jīng)由實驗段加熱的支路中的流體先經(jīng)過換熱器(圖2)，被已經(jīng)過預熱器和實驗段加熱的高溫流體預熱。該換熱器為管式換熱器，被預熱的冷卻劑在管內(nèi)流動(二次側(cè))，由實驗段流出的高溫冷卻劑在管內(nèi)的9根小管(一次側(cè))內(nèi)反向流動。之后這部分流體再被預熱器和實驗段加熱。預熱器為電加熱，其作用是使實驗段進口的溫度達到設計值。實驗段為一直流均勻加熱棒束。實驗段出口處的流體溫度可達450 ℃。從實驗段出口出來的這部分高溫流體再次經(jīng)由換熱器預熱二次側(cè)流體后，與之前被旁排的那支流體在混合箱中混合，溫度有所降低，最終被兩個平行的冷卻二回路繼續(xù)冷卻至60 ℃。主回路流體在最終進入開式水箱前，接有1個穩(wěn)壓器，之后會經(jīng)過2個閥門K1和K2。穩(wěn)態(tài)時，閥門K2打開，閥門K1關(guān)閉。瞬態(tài)開始后，在保持閥門K2打開的同時，閥門K1會打開，這時泄壓開始，回路壓力開始降低，進入跨臨界過程。

圖1 SWAMUP實驗裝置Fig.1 SWAMUP test facility

圖2 換熱器Fig.2 Heat exchanger

2 建模

SWAMUP實驗裝置的ATHLET模型如圖3所示，該計算模型包含185個控制體、193個連接體、37根管、71個熱構(gòu)件。其中建模的關(guān)鍵參數(shù)已在圖中標出，柱塞泵被60 ℃、10 m3/h流量流體的邊界條件所替代，流經(jīng)實驗段支路的流量為3.6 m3/h，預熱器功率為637.4 kW，實驗段功率為63.6 kW，實驗段出口流體溫度為450 ℃，穩(wěn)壓器初始壓力為25 MPa，水箱的邊界條件為60 ℃、0.1 MPa。

圖3 SWAMUP實驗裝置的ATHLET模型Fig.3 ATHLET model of SWAMUP test facility

3 穩(wěn)態(tài)分析

在進行瞬態(tài)分析前，對程序進行穩(wěn)態(tài)計算，程序運行2 000 s后，基本達到穩(wěn)態(tài)，系統(tǒng)此時的關(guān)鍵參數(shù)接近回路設計值，關(guān)鍵熱工水力參數(shù)設計值與計算值的對比列于表1。

圖4為ATHLET-SC程序計算得到的實驗段的運行參數(shù)隨時間的變化。圖5為穩(wěn)態(tài)時主實驗回路的溫度分布。預熱器進口的流體溫度為385 ℃，在擬臨界溫度附近。實驗段進口溫度為439 ℃，出口溫度達450 ℃。

表1 熱工水力參數(shù)設計值與計算值的對比Table 1 Comparison of design and calculation values for thermal-hydraulic parameter

圖4 實驗段的運行參數(shù)隨時間的變化Fig.4 Operation parameters of test section vs. time

圖5 主實驗回路的溫度分布Fig.5 Temperature distribution of main test loop

4 瞬態(tài)分析

本工作的主要目的是分析回路裝置在實驗工況下各安全指標是否滿足設計要求，所以，泄壓過程中加熱棒溫度和系統(tǒng)壓力的變化是本文分析的重點。

2 000 s后，瞬態(tài)開始，閥門K2保持打開的同時，閥門K1在0.1 s內(nèi)打開，而實驗段功率在泄壓過程中保持63.6 kW不變。閥門K1打開后，實驗段的壓力開始下降，在約2 004.3 s時降到臨界壓力(22.1 MPa，圖6a)，之后，實驗段內(nèi)流體的空泡份額從1.00開始下降，最終穩(wěn)定在0.85左右(圖6b)。

圖6c為實驗段棒束與冷卻劑之間的換熱模型隨時間的變化。當模型代號為55時，流體處于超臨界狀態(tài)。而代號為49、48、45表示流體處于膜態(tài)沸騰區(qū)，代號35和25分別表示流體處于過渡沸騰區(qū)和核態(tài)沸騰區(qū)。另外，模型代號49表示此時流體換熱形式為為氣相換熱，而25至48之間的數(shù)字代號則表示此時流體換熱形式為兩相換熱。

當實驗段壓力降至臨界壓力以下時，實驗段內(nèi)的冷卻劑流動由超臨界流動變?yōu)閮上嗔?，? 004.4 s時，冷卻劑的換熱模型變?yōu)槟B(tài)沸騰，在2 060 s時，變?yōu)檫^渡沸騰。圖6d為換熱器兩端的壓差曲線。當泄壓發(fā)生后，首先受影響的是流動的下游方向，即換熱器一次側(cè)壓力顯著下降，而換熱器二次側(cè)的壓力下降則較為緩慢，換熱器的壓差開始增加，在擬臨界點附近達到極值。而后隨著實驗段流體變?yōu)闅庖簝上啵瑑上嗄Σ磷枇ο禂?shù)較大，這導致實驗段和預熱器內(nèi)的壓差上升。但在經(jīng)過臨界點之后的一段時間內(nèi)，實驗段和預熱器內(nèi)的冷卻劑由超臨界狀態(tài)變?yōu)閮上嗔鳎@個相變過程使得壓差有個暫時的下降。最大的壓差出現(xiàn)在冷流體的進口和熱流體的出口，即一次側(cè)的出口和二次側(cè)的進口，且最大壓差將達6 MPa。

圖6 瞬態(tài)時的系統(tǒng)變化Fig.6 System behavior in transient state

瞬態(tài)時系統(tǒng)各參數(shù)隨時間的變化示于圖7。從圖7a可看出，泄壓開始后，壓差增大，實驗段流量上升，在擬臨界點附近達到峰值1.6 kg/s。而后由于實驗段內(nèi)流體發(fā)生了相變，流動阻力暫時增大，而流量取決于整個回路的壓差和阻力，這導致實驗段流量迅速降低，隨泄壓的繼續(xù)，壓差逐漸加大使實驗段流量逐步回升，并保持在較初始流量更高的狀態(tài)(圖6d)。從圖7b可看出，在達到臨界點之前，隨壓力的下降，比定壓熱容顯著上升，伴隨流速的增加，使換熱性能大幅提高；而系統(tǒng)壓力達到臨界壓力后，流體的流速、比定壓熱容和汽化潛熱顯著下降，導致?lián)Q熱性能大幅降低。在2 030 s前，實驗段內(nèi)流體處于膜態(tài)沸騰區(qū)，氣相的換熱系數(shù)幾乎等于整體的換熱系數(shù)，但2 030 s后，空泡份額降低，氣相的換熱系數(shù)依然維持在較低水平，但液相的換熱系數(shù)增加并占主導，使整體的換熱系數(shù)有一上升，直至2 060 s時，流體進入過渡沸騰區(qū)后保持基本穩(wěn)定。

加熱棒壁面溫度和冷卻劑流體溫度的變化是實驗段流量和換熱系數(shù)共同作用的結(jié)果。圖7c為加熱棒壁面溫度隨時間的變化曲線。由上述換熱系數(shù)的分析可知，隨泄壓的進行，在達到臨界壓力前，換熱系數(shù)增加，使加熱棒壁面溫度一直下降，直至達到擬臨界點附近，換熱系數(shù)的降低較流量下降的速度更快，壁面溫度開始上升，出現(xiàn)峰值。而后空泡份額開始變化，如前所述，實驗段流體的換熱由氣相膜態(tài)沸騰轉(zhuǎn)變?yōu)閮上嗄B(tài)沸騰，液相的換熱系數(shù)開始增加，使得實驗段流體總的換熱系數(shù)開始緩慢上升，而此時實驗段流量仍很低，這使得實驗段出口處壁面溫度在短時間內(nèi)穩(wěn)定在一定幅度內(nèi)，而2 030 s后隨實驗段流體換熱系數(shù)和流量的同時快速上升，壁面溫度開始下降直至緩慢穩(wěn)定。從圖7d可看出，實驗段流體溫度的變化基本與壁面溫度的變化趨勢一致。

圖7 瞬態(tài)時系統(tǒng)參數(shù)隨時間的變化Fig.7 System parameter vs. time in transient state

5 敏感性分析

5.1 泄壓速率的影響

為模擬不同泄壓速率對回路的影響，通過改變不同的閥門開度來實現(xiàn)，即在100%閥門開度的基礎(chǔ)上分別減小到50%和增加到200%。泄壓速率的敏感性分析列于表2。不同工況下系統(tǒng)參數(shù)的變化示于圖8。從圖8可看出，在50%和200%的閥門開度情況下，分別形成了兩個峰值，并隨著閥門破口面積的增大，實驗段的降壓速率增加，實驗段流量的峰值也增加，這使得在閥門開度為200%的工況下，加熱棒壁面溫度在擬臨界點附近的峰值更低。之后，隨著流體流量和換熱系數(shù)的降低，實驗段流體溫度又出現(xiàn)一個上升，在2 030 s左右，達到另一個峰值。與工況1相比，在工況3下，實驗段流體的換熱系數(shù)更小，而此時流量與其他工況基本相同，故而工況3的第2個峰值會更高。

表2 泄壓速率的敏感性分析Table 2 Sensitivity analysis on depressurization rate

5.2 加熱棒加熱功率的影響

為了研究加熱棒加熱功率對結(jié)果的影響，在保證實驗段出口溫度(450 ℃)不變，即保證預熱器和實驗段加熱總功率不變的前提下，改變預熱器和實驗段之間的功率分配。功率分配方式為：方式1，實驗段加熱功率63.6 kW，預熱器加熱功率637.4 kW；方式2，實驗段加熱功率100 kW，預熱器加熱功率601 kW；方式3，實驗段加熱功率150 kW，預熱器加熱功率551 kW。不同功率分配方式下系統(tǒng)參數(shù)的變化示于圖9。由圖9可知，以方式1作為基準，方式2、3的實驗段加熱功率逐漸增加，實驗段的壁面溫度和實驗段流體溫度在泄壓過程中也會更高。如圖9a所示，改變加熱功率的分配會影響實驗段流體的換熱模型，對于方式3，實驗段內(nèi)流體一致保持在膜態(tài)沸騰區(qū)，未進入過渡沸騰區(qū)。而對于方式2，其實驗段內(nèi)流體進入過渡沸騰區(qū)的時間點也較方式1有所推遲，這影響到冷卻劑的換熱系數(shù)，從而使實驗段的壁面溫度和流體溫度有高低之差。

圖8 不同工況下系統(tǒng)參數(shù)的變化Fig.8 System parameters of different cases

圖9 不同功率分配方式下系統(tǒng)參數(shù)的變化Fig.9 System parameters of different power distributions

5.3 臨界熱流密度系數(shù)的影響

在本節(jié)中，增大臨界熱流密度系數(shù)，是指提高了臨界熱流密度和返回核態(tài)沸騰的溫度，這就使實驗段流體提早由膜態(tài)沸騰向核態(tài)沸騰轉(zhuǎn)變。本節(jié)討論臨界熱流密度系數(shù)因子分別為0.2、1.0和2.0三種情況下系統(tǒng)參數(shù)的變化(圖10)。從圖10a可看出，不同于其他兩種情況，當臨界熱流密度系數(shù)因子為2.0時，實驗段流體由過渡沸騰進入到核態(tài)沸騰區(qū)，且此時換熱系數(shù)瞬間增大，加熱棒壁面溫度瞬間下降，如圖10b、c所示。由圖10d可看到，臨界熱流密度系數(shù)的改變對實驗段冷卻劑溫度的影響不大。

5.4 氣相強制對流換熱系數(shù)的影響

由前述瞬態(tài)工況的分析可知，在膜態(tài)沸騰區(qū)時，實驗段流體氣相的換熱系數(shù)基本等于流體總的換熱系數(shù)(氣相換熱系數(shù)與液相換熱系數(shù)之和)。本節(jié)討論當氣相強制對流換熱系數(shù)(HTC)因子分別為0.2、1.0、2.0時，系統(tǒng)各熱工水力參數(shù)的變化(圖11)。從圖11可看出，以上3種條件下，實驗段流體的換熱模型轉(zhuǎn)變很相似，但冷卻劑的換熱系數(shù)卻差別很大，導致其對應的最高壁面溫度分別為708、449和438 ℃，即隨換熱系數(shù)因子的升高，最高壁面溫度降低。這是因為當換熱系數(shù)因子為0.2時，過低的換熱系數(shù)導致加熱棒在泄壓過程時出現(xiàn)了沸騰危機，這使壁面溫度在擬臨界點附近會明顯升高，一度升至大于700 ℃。而此時，實驗段內(nèi)的流體溫度也由于換熱系數(shù)太低會有一個驟降。之后隨泄壓的進行，實驗段冷卻劑的空泡份額逐漸降低，液相開始增多，氣相的換熱開始減弱，總的換熱系數(shù)開始增大。直至在2 060 s左右實驗段流體進入過渡沸騰區(qū)，液相換熱占主導地位后，這3種條件下的換熱系數(shù)基本趨于一致，壁面溫度和流體溫度也基本趨于一致。

圖10 不同臨界熱流密度系數(shù)因子下系統(tǒng)參數(shù)的變化Fig.10 System parameters of different CHF coefficient factors

圖11 不同氣相強制對流換熱系數(shù)因子下系統(tǒng)參數(shù)的變化Fig.11 System parameters on different HTC factors of forced convection to vapor

6 結(jié)論

本文使用修改后的ATHLET-SC程序?qū)WAMUP實驗裝置建模，并對回路的熱工水力特性進行分析。結(jié)果表明，被修改后的ATHLET-SC程序完全可計算跨臨界瞬態(tài)過程。計算結(jié)果顯示，在泄壓開始后的4.3 s內(nèi)，實驗段壓力降至超臨界壓力以下，實驗段內(nèi)的冷卻劑換熱模式會經(jīng)過膜態(tài)沸騰、過渡態(tài)沸騰或核態(tài)沸騰(取決于不同的計算模型)階段，最高壁面溫度發(fā)生在膜態(tài)沸騰區(qū)。敏感性分析表明：泄壓速率越大，加熱棒最高壁面溫度越高；加熱棒功率越高，最高壁面溫度越高；臨界熱流密度系數(shù)因子增大時，其換熱模型由過渡態(tài)沸騰轉(zhuǎn)變?yōu)楹藨B(tài)沸騰；氣相的換熱系數(shù)因子越小，其最高壁面溫度越高。另外需特別考慮的是，回路中換熱器的內(nèi)外壓差將達6 MPa，這在實驗過程中需予以考慮。

感謝德國GRS核安全中心提供ATHLET程序。

[1] CHENG X, LIU X J, YANG Y H. A mixed core for supercritical water-cooled reactors[J]. Nuclear Engineering and Technology, 2008, 40(2): 117-126.

[2] 程旭，劉曉晶. 超臨界水冷堆國內(nèi)外研發(fā)現(xiàn)狀與趨勢[J]. 原子能科學技術(shù)，2008，42(2):167-172.

CHENG Xu, LIU Xiaojing. Research status and prospect of supercritical water-cooled reactor[J]. Atomic Energy Science and Technology, 2008, 42(2): 167-172(in Chinese).

[3] GIF. A technology roadmap for generation Ⅳ nuclear energy systems[ED/OL]. (2002-10-22). http: ∥www.gen-4.org/Technology/roadmap. htm.

[4] SCHULENBERG T, STARFLINGER J, MARSAULT P, et al. European supercritical water cooled reactor[J]. Nuclear Engineering and Design, 2011, 241(9): 3 505-3 513.

[5] 周翀,楊燕華,程旭. 超臨界水冷堆燃料性能驗證實驗回路的冷卻劑喪失事故分析[J]. 原子能科學技術(shù)，2013，47(9):1 554-1 559.

ZHOU Chong, YANG Yanhua, CHENG Xu. Loss of coolant accident analysis of supercritical water-cooled reactor fuel qualification test loop[J]. Atomic Energy Science and Technology, 2013, 47(9): 1 544-1 559(in Chinese).

[6] ZHAO M, GU H Y, CHENG X. Experimental study on heat transfer of supercritical water flowing downward in circular tubes[J]. Annals of Nuclear Energy, 2014, 63: 339-349.

[7] FU S W, LIU X J, ZHOU C, et al. Modification and application of the system analysis code ATHLET to trans-critical simulations[J]. Annals of Nuclear Energy, 2012, 44: 40-49.

[8] ZHOU C, YANG Y H, CHENG X. Feasibility analysis of the modified ATHLET code for supercritical water cooled systems[J]. Nuclear Engineering and Design, 2012, 250: 600-612.

Thermal-hydraulic Analysis of SWAMUP Test Facility

WANG Zi-di, CAO Zhen, LIU Xiao-jing, CHENG Xu

(SchoolofNuclearScienceandEngineering,ShanghaiJiaoTongUniversity,Shanghai200240,China)

Trans-critical transients, i.e. the pressure in the reactor system undergoing a rapid decrease from the supercritical pressure to the subcritical condition, are of crucial importance for the LOCA analysis of supercritical water cooled reactor (SCWR). To obtain more knowledge about this process, the supercritical water multipurpose loop (SWAMUP) test facility in Shanghai Jiao Tong University (SJTU) will be applied to provide test data for the process. Some pre-test calculations are necessary to show the feasibility of the experiment. In this study, the trans-critical transient analysis was performed for the SWAMUP test facility with the system code ATHLET-SC which was modified in SJTU. The system behaviors including system pressure, coolant mass flow, and wall temperature of heating rod during the depressurization were presented. The effects of some important parameters such as depressurization rate, heating power of heating rod on the system characteristics were also investigated. The results indicate that the revised system code ATHLET-SC is capable of simulating thermal-hydraulic behaviors during the trans-critical transient. According to the results, the wall temperature of heating rod during the transient is kept below the design limit. However, the pressure difference between the heat exchanger inside and outside after depressurization could reach 6 MPa, which should be considered in the test.

SWAMUP test facility; trans-critical transient; ATHLET-SC code; pre-test calculation

2014-03-26;

2014-05-09

高等學校博士學科點專項科研基金資助項目(20110073120045)

汪子迪(1991—)，男，江西婺源人，碩士研究生，核能科學與工程專業(yè)

TL333

A

1000-6931(2015)07-1191-09

10.7538/yzk.2015.49.07.1191