IVR條件下外部冷卻流道變形的影響因素研究

田欣鷺,溫 爽，郭 超，魏 超，李鐵萍，*,王 凱

(1.環境保護部核與輻射安全中心, 北京 100082;2.西安交通大學 能源與動力工程學院, 陜西西安 710049)

IVR條件下外部冷卻流道變形的影響因素研究

田欣鷺1,溫 爽1，郭 超1，魏 超1，李鐵萍1，*,王 凱2

(1.環境保護部核與輻射安全中心, 北京 100082;2.西安交通大學 能源與動力工程學院, 陜西西安 710049)

當發生堆芯熔化事故時，壓力容器外部冷卻是保持壓力容器完整性及實現熔融物堆內滯留(In-Vessel Retention，簡稱IVR)的一項重要策略。在高溫熔融物的熱載荷和內部壓力的共同作用下，壓力容器外壁面和保溫層之間的冷卻流道可能發生變形，造成冷卻能力的降低，進而威脅到壓力容器的完整性。因此，有必要分析IVR條件下壓力容器冷卻流道變形的影響因素。結果表明，熱膨脹是造成冷卻流道變形的主要因素。在IVR策略成功的前提下，內壓和熱流密度對流道變形的影響有限。

熔融物堆內滯留；外部冷卻流道；熱膨脹；堆芯熔毀

當核電廠發生堆芯燃料棒熔化事故時，熔融物會落入壓力容器(Reactor Pressure Vessel，簡稱RPV)下封頭內,并形成高溫熔池。在熔融物熱載荷的沖擊下，RPV壁面有可能被熔穿，從而造成一回路壓力邊界完整性的破壞和放射性物質的外泄[1]。

針對這一類嚴重事故現象，越來越多的核電廠采用了熔融物堆內滯留策略來緩解事故的后果[2]，而壓力容器外部冷卻(External Reactor Vessel Cooling，簡稱ERVC)是IVR策略成功的前提之一。IVR-ERVC措施主要包括RPV內部卸壓和外部冷卻。RPV壁面的外部冷卻過程大致如下：當發生堆芯熔毀事故時，冷卻水進入到由RPV外壁面和保溫層之間的冷卻流道(如圖1所示)，對RPV壁面進行冷卻后，產生的氣液兩相流向上通過排氣孔釋放。

圖1 外部冷卻流道示意圖Fig.1 External coolant channel

目前，國內外針對IVR條件下RPV完整性的研究主要圍繞以下兩個問題開展：

第一類是RPV壁面換熱和外部冷卻能力的研究(如圖2所示)。只有當冷卻水具有足夠的冷卻能力時，才能迅速帶走衰變熱，使RPV外壁保持在較低的溫度。而這也是RPV壓力邊界保持結構完整性和包容性的必要條件。

圖2 IVR傳熱模型示意圖Fig.2 Heat transfer model of IVR

第二類是RPV壁面在高溫下的結構力學研究。在力學載荷和熱載荷的共同作用下，RPV壁面很容易產生高溫蠕變變形，而蠕變現象正是造成高溫承壓容器失效的主要因素[3]。避免RPV過度變形和斷裂失效，是RPV保持結構完整性和包容性的充分條件。近年來，我國和一些核電發達國家圍繞這部分內容開展了一系列研究[4]。這些研究包括針對不同牌號RPV用鋼(16MND5[5]，SA533B1，SA508[6]等)開展高溫蠕變性能的試驗研究[7]；有些是利用縮比模型[8]和相似準則[9]對假想的嚴重性超過設計基準事故并造成堆芯明顯惡化的事故工況進行模擬試驗，并利用試驗結果對原型進行研究；還有一些則采用有限元數值模擬結合試驗數據的方式對RPV的高溫蠕變變形進行分析[8]。一些文獻對IVR條件下RPV完整性的研究現狀進行了分析和總結[10]。

如上文所述，外部冷卻流道由RPV外壁面和保溫層構成，而保溫層的變形可以忽略，因此冷卻流道的變形被認為只與RPV外壁面相關。雖然冷卻流道的設計尺寸在在設計上可以滿足冷卻的要求，但由于事故工況下冷卻流道的尺寸可能變窄，因此有必要對流道變形后的冷卻能力進行評估。本文在上述各項研究的基礎上，對IVR條件下RPV外壁面的變形進行數值模擬計算，并分析了影響RPV壁面變形的因素，為找出傳熱能力惡化最嚴重的冷卻流道提供必要的參數。

1 計算模型和材料參數

1.1 計算模型和假設

本文利用通用有限元計算軟件ANSYS對某堆型的RPV的壁面進行了全尺寸建模、網格剖分和模擬計算。由于計算過程涉及到蠕變等非線性特征，為了降低計算量，根據RPV的結構和載荷特征，在建模時通過必要的假設將其簡化為二維軸對稱模型。為了消除邊界條件對計算的影響，在模型中加入了RPV下筒體部分。網格以四邊形網格為主，只有在幾何形態變化較大的部分采用少量三角形網格。

本文對分析過程進行了解耦處理，分兩步進行計算：首先計算出高溫熔融物和冷卻水共同作用下RPV壁面的穩態溫度場；之后將溫度場導入到力學模型中，計算在高溫和內壓(包括熔融物的載荷)共同作用下RPV壁面的變形。本文的計算基于以下假設：

(1)RPV下封頭為半球狀，堆內最終形成兩層穩定結構熔融池。在堆芯熔融物完全坍塌到下封頭之前，堆腔水完全淹沒RPV下封頭，并且堆腔注水系統的補水量大于堆腔中水的蒸發量。外壁冷卻水始終處于沸騰危機之下。

(2)器壁無內熱源，熔化時消耗的相變潛熱暫不考慮，器壁材料的熔點為1600K，溫度超過熔點的那部分器壁直接按熔化處理，不考慮器壁熔化消耗的相變潛熱，不考慮器壁以及堆芯熔融物的輻射傳熱。

1.2 材料參數

本文研究的RPV下封頭所使用的材料牌號是SA-508GR.3CL.2。該材料常規的熱學和機械性能參數可以在ASME和RCC-M等規范中查到。但該材料的高溫性能和蠕變性能只能參考相關文獻。

SA-508GR.3CL.2的高溫熱學和機械性能參考了文獻中法國牌號16MND5鋼的數據[11]。事實上，雖然這兩種材料在不同國家被賦予了不同的牌號，但二者化學成分幾乎完全相同。

而蠕變性能的確定卻相對復雜。標準和文獻中均沒有SA-508的蠕變參數。因此只能參考其他RPV用鋼材料的蠕變數據。目前，文獻中有蠕變試驗數據可查的，且與SA-508化學成分最為相似的材料是SA-533B1。該材料有兩組高溫蠕變性能的試驗數據，蠕變數據1[12]采用的是Modified Time Hardening(primary)模型，蠕變數據2[13]采用的是Time Hardening模型。

為驗證數據的準確性，本文選取Rempe[14]試驗對這兩組數據進行對比。作為驗證性的試驗，Rempe試驗采用了1m×0.2m的矩形試驗件，溫度T=1150K，其中一端施加固支約束，另一端施加26.5 MPa的壓力。選用不同的蠕變參數進行計算的結果如圖3所示。

圖3 蠕變數據對比Fig.3 Creep data comparison

如圖3所示，蠕變數據1[12]與試驗結果最為接近，最大相差不超過5%。蠕變數據2[13]與試驗結果相差大約4倍左右。本文選擇與試驗數據接近的蠕變數據1作為最終的輸入數據。但在需要保證計算結果保守性時，選擇蠕變數據2更加合理。

2 RPV壁面溫度場分析

溫度場的分析基于對有限元模型設置合理的邊界條件和熱載荷。這里施加的邊界條件和載荷主要包括溫度和熱流密度：

(1)對于壓力容器外壁面，設置400 K的溫度約束。IVR措施成功的前提之一是RPV外部被充分冷卻，此時400K的溫度約束相對于冷卻水的飽和溫度和過熱度來說，是比較保守的。

(2)對于壓力容器內壁面，給定了沿壁面角度變化的熱流密度。這里假設所發生的事故為全場斷電(Station Blackout，簡稱SBO疊加非能動預熱導出余熱導出(Passive Residual Heat Removal，簡稱PRHR完全失效事故。

溫度場的計算結果顯示下封頭部分壁面的溫度超出了材料的熔點(1600K)。因此要對這一部分超出器壁熔點的區域進行“殺單元”處理。即當壁面溫度超過熔點時，該部分器壁直接按熔化處理，不考慮器壁熔化消耗的相變潛熱，也不考慮器壁以及堆芯熔融物的輻射傳熱。

本文通過ANSYS計算嚴重事故下的RPV下封頭壁面的厚度，并與文獻[15]中使用編程方法得到的結果進行對比，來驗證“殺單元”方法的準確性。圖4是根據文獻中的參數計算出的下封頭剩余厚度分布。在0-37℃，溫度并未超過熔點，而內壁面從37℃開始熔化，在76℃熔化程度最嚴重，壁面不足50mm。計算值與文獻[15]中的厚度符合較好，因此采用ANSYS 的“殺單元”方法是可行的。

圖4 下封頭壁厚的變化Fig.4 Variation of bottom head thickness

利用上述方法，對超過1600K的單元進行“殺單元”操作——如果一個網格4個頂點的平均溫度超過1600K，即將此單元“剔除”。“殺單元”之后剩余壁面的溫度場如圖5所示。可以看出，由于熱流密度較高，RPV內壁面大部分區域都超過了材料熔點，造成了壁面大面積的熔毀。最嚴重的部位，壁面厚度減少超過了80%。因此，在進行力學分析時不能忽略熔壁現象的影響。本文采用的方法是根據溫度場的計算結果，對熔壁后的RPV重新進行建模后，再進行力學分析。

圖5 壓力容器壁面溫度場Fig.5 Temperature fields of RPV

3 RPV壁面變形分析

在進行力學場分析時，除了需要設置對稱邊界條件和約束外，還需要施加以下幾類載荷：

(1)在RPV外壁面施加冷卻水壓力(即冷卻水對RPV的浮力)。

(2)在RPV內壁面施加0.1MPa的內壓，以及熔融物的重力。RPV內部卸壓是IVR措施成功的前提，本文假設一回路已經充分卸壓，并留有0.1MPa的裕度。

RPV壁面的位移場如圖6所示，其外表面的位移大小即是該點外部冷卻流道的變化量。

圖6 壓力容器壁面位移場Fig.6 Displacement fields of RPV

計算結果顯示，水平方向的位移在RPV下封頭與下筒體的交接處達到最大值，而垂直方向的位移則在45度左右的位置最大。因此這兩個位置上外部冷卻流道的縮減最為明顯，在進行換熱能力的評估時應予以重點關注。

4 外部冷卻流道變形分析

如上文所述，當發生堆芯熔化事故時，冷卻劑進入RPV外壁和保溫層之間的冷卻流道，對RPV壁面進行冷卻。在熱膨脹和蠕變共同作用下。RPV外壁面會發生形變，從而造成冷卻劑流道形狀的變化。在外壁保溫層幾何尺寸已知的情況下，流道的最終形狀可以通過RPV外壁面各節點的位移計算確定。外部冷卻流道隨壁面角度的變形量見表1(壓力容器底部為0度，下封頭與下筒體交接處為90度)。

表1 外部冷卻流道的變化量Table 1 Variation of external coolant channel

上述結果表明：在RPV下封頭形成熔池后，外部冷卻流道的寬度有所降低，縮減值隨角度不同在13mm至18mm之間變化，并分別在45度和90度左右出現了峰值。

雖然這個變化量相對于RPV的整體尺寸來說幾乎可以忽略，也不會造成結構的失效，但對于寬度只有90mm至150mm左右的外部冷卻流道來說，縮減量達到了約原始尺寸的1/10至1/5。這一縮減量是否會造成RPV壁面局部的傳熱惡化，還需要通過進一步計算RPV壁面的臨界熱流密度曲線來確定。

5 外部冷卻流道變形影響因素分析

5.1 熱膨脹和蠕變對冷卻流道變形的影響

造成RPV外壁面變形，即外部冷卻流道變形的因素，主要包括高溫下RPV壁面的熱膨脹和蠕變變形。

當IVR措施成功時，表明卸壓已經成功。此時因為載荷有限，所以蠕變對RPV外壁面變形場的影響幾乎可以忽略。

相反，在高溫度梯度下產生的熱應力，對RPV的力學場尤其是冷卻流道尺寸的影響卻遠遠大于其他因素的作用。如果在與上文完全相同的情況下重新進行計算，但只考慮熱膨脹而忽略蠕變現象時，冷卻流道在0度位置縮減量變為13.11 mm，與之前考慮蠕變現象的結果相比減小了0.07 mm。隨著角度的增加，二者計算結果的偏差不斷增加，在30度、60度和90度分別達到0.28 mm、0.39 mm和1.2 mm。但相對于總體縮減量來說，仍然很小。由此可見，在IVR條件下，造成RPV變形以及冷卻流道變窄的主要因素是熱膨脹，高溫蠕變的影響遠小于熱膨脹。

5.2 熱流密度和內壓對冷卻流道變形的影響

為了研究影響冷卻流道變形的因素，本節對內壓和熱流密度進行調整后，重新進行了計算，并將計算的結果與之前的結果進行對比，定性的研究了上述兩類因素對結果的影響大小。

本節涉及的計算有以下兩種：

(1)低內壓情況下的流道變形的計算：將之前施加的0.1MPa內壓改為0，其余邊界條件與之前一致。IVR措施要求在發生堆芯熔化事故時，RPV內部需要降壓，因此內壓可以認為與外界相同。相對于之前的分析，這一處理雖然喪失了保守性，但更接近于實際情況。

(2)低熱流密度情況下的流道變形的計算：將施加在RPV內壁面的熱流密度減小10%，其余邊界條件與之前一致。這一改變基于以下兩個原因，一方面，之前施加的熱流密度曲線為最不利工況時的包絡曲線，因此在進行對比分析時減小了熱流密度；另一方面，如果施加的熱流密度過小，又不能真實反應熔堆事故的實際情況。10%減幅的選擇是平衡以上兩種考慮后所確定的。

這兩類計算的結果與之前計算結果的對比見表2。

表2 不同工況外部冷卻流道變化量Table 2 Variation of external coolant channel under different conditions

其中工況1為原始工況，工況2降低了內壓，工況3降低了熱流密度。

可以看到，當內壓減小時，冷卻流道的縮減值有所降低，但并不十分明顯。這說明冷卻流道的變形隨壓力的增加而增加，但在內壓不大的情況下，壓力容器內壓力對冷卻流道變形的影響非常有限。

相對內壓而言，熱流密度對冷卻流道變形的影響較大，這主要由于以下兩個原因：

(1)RPV外壁面，即外部冷卻流道的變形受熱膨脹影響較大，而熱流密度的降低會使RPV壁面的溫度有所降低。

(2)溫度場的變化還會進一步影響到RPV壁面熔化的范圍，從而改變RPV的幾何形狀。

但通過表2可以看出，盡管熱流密度對冷卻流道的影響比內壓大，但仍然十分有限。這主要是由于：當熱流密度非常大時，即使降低了一些，但大部分壁面仍然出現了熔化現象，此時熱流密度的變化只改變了RPV壁面的剩余厚度，對溫度場的分布影響不大，RPV大部分內壁面仍處于熔點溫度。由于熱膨脹造成的流道變形主要與RPV整體尺寸和溫度場分布有關，而與局部壁面厚度關系不大，因此熱流密度降低10%，對冷卻流道的變形影響不會太大。

需要注意的是，上述結論成立的前提是施加在RPV內壁面上的熱流密度很大，足以使大部分壁面熔化。如果施加的熱流密度較小，以至于RPV幾乎不發生熔壁，那么溫度場的大小和分布便會隨熱流密度的減小而發生顯著變化，從而使熱膨脹造成的冷卻流道變形大大減少。由于本文研究的內容是熔堆情況下冷卻流道的變形，因此這種情況不在本文討論范圍以內。

綜上所述，在熔堆現象發生后，如果RPV外部成功實現冷卻，且內部成功卸壓，那么剩余的內壓對冷卻流道的影響幾乎可以忽略；而當施加在內壁的熱流密度較高，足以使RPV大部分壁面發生熔壁現象時，那么熱流密度的改變也不會對RPV變形的結果產生本質影響。

6 結論

本文通過建立RPV壁面的全尺寸模型，并施加合理的邊界條件和載荷，對IVR條件下外部冷卻流道變形及其影響因素進行了分析后，得出以下幾點結論：

(1)RPV下封頭的變形會造成外部冷卻流道尺寸的降低，縮減值在13mm至18mm之間，大約占原始流道寬度的1/5至1/10。外部冷卻流道變窄后的換熱能力需要進行重新評估。

(2)造成冷卻流道變形的最主要因素是高溫下RPV壁面的熱膨脹，蠕變變形雖然是造成RPV失效的主要因素，但對冷卻流道變形的影響不大。

(3)在發生嚴重事故熔堆時，如果內部卸壓成功，且外部冷卻能保證RPV壁面溫度保持不變，那么RPV內部壓力和熔融物熱流密度的改變對外部流道變形的影響不大。

[1] 武志瑋, 寧冬, 姚偉達. 嚴重事故下反應堆壓力容器材料高溫蠕變研究進展[J]. 核安全, 2011(2):20-24.

[2]Kym?l?inen O, Tuomisto H, Theofanous T G. In-vessel retention of corium at the Loviisaplant[J].Nuclear Engineering & Design, 1997, 169(1):109-130.

[3]李琳, 臧希年. 壓水堆核電廠嚴重事故下堆芯熔融物的冷卻研究[J]. 核安全, 2007(4):39-44.

[4]THEOFANOUS TG,LIUC,ADDITON S,et al.In-vessel coolability and retention of a core melt[J]. Nuclear Engineering and Design,1997, 169(1-3): 1-48..

[5]Devos J, Sainte C C, Poette C, et al. CEA programme to model the failure of the lower head in severe accidents[J]. Nuclear Engineering & Design, 1999, 191(1):3-15.

[6]Thinnes G L, Korth G E, Chavez S A, et al. High-temperature creep and tensile data for pressure vessel steels SA533B1 and SA508-CL2[J]. Nuclear Engineering & Design, 1994, 148(2-3):343-350.

[7]Sehgal BR, Nourgaliev R R, Dinh T N. Characterization of heat transfer processes in a melt pool convection and vessel-creep experiment[J]. Nuclear Engineering & Design, 2002, 211(2-3):173-187.

[8]H.-G. Willschütz, E. Altstadt, B.R. Sehgal, et al. Coupled thermal structural analysis of LWR vessel creep failure experiments[J]. Nuclear Engineering& Design, 2001, 208(3):265-282.

[9]Pilch M M, Ludwigsen J S, Chu T Y, et al. Creep failure of a reactor pressure vessel lower head under severe accident conditions[J]. Office of Scientific & Technical Information Technical Reports, 1998.

[10]姚彥貴,寧冬,武志瑋,等.假想堆芯熔化嚴重事故下反應堆壓力容器完整性的研究進展與建議[C]//第十七屆全國反應堆結構力學會議論文集,2012-10-15,中國上海:出版社不詳,2012:258-266.

[11]H.-G. Willschütz, E. Altstadt. Generation of a High Temperature Material Data Base and its Application to Creep Tests with French or German RPV-steel[J]. 2002.

[12]Villanueva W, Tran C T, Kudinov P. Coupled thermo-mechanical creep analysis for boiling water reactor pressure vessel lower head[J]. Nuclear Engineering & Design, 2012, 249(249):146-153.

[13]Rashid Y R. Creep considerations for the lower head[J]. Nuclear Engineering& Design, 1997, 169(1):101-108.

[14]Rempe J L, Chavez S A, Thinnes G L. Light water reactor lower head failure analysis[J]. 1993.

[15]郭濤, 包士毅, 高增梁. 嚴重事故IVR下反應堆壓力容器穩態溫度場計算[J]. 壓力容器, 2012, 29(6):34-39.

StudyonInfluencingFactorsofExternalCoolingChannelDeformationunderIVR

TIAN Xinlu1, WEN Shuang1,GUO Chao1,WEI Chao1,LI Tieping1,*，WANG Kai2

(1. Nuclear and Radiation Safety Center, MEP, Beijing 100082, China;2.School of Energy and Power Engineering, Xi’an Jiaotong Univerisity, Xi’an 710049, China)

The external reactor vessel cooling is an important strategy for maintaining the integrity of the pressure vessel and for holding the melt in the case of core melt. The cooling channel between the pressure vessel outer wall and insulation layer may be deformed under the combined action of the thermal load and the internal pressure. It maybe leads to cooling capacity reducing and then threatens the integrity of the pressure vessel. Therefore, it is necessary to analyze the influencing factors of the cooling channel deformation of the pressure vessel under IVR. The results show that the thermal dilatation is the main factor causing the deformation of the cooling channel. Under the premise of successful implementationof IVR strategy, the influence of internal pressure and heat flux on cooling channel deformation is limited.

In-Vessel Retention；external coolant channel；thermal dilatation；core melt

TH49

：A

：1672- 5360(2017)02- 0005-06

2017- 04- 14

2017- 06- 03

國家科技重大專項課題，項目編號：2015ZX06002007-003

田欣鷺(1989—)，女，遼寧沈陽人，助理工程師，本科，船舶工程專業，現主要從事核電廠設備力學研究工作