高溫氣冷堆進(jìn)氣事故中自然循環(huán)建立時間及影響因素分析

劉 鵬，鄭艷華，石 磊

（清華大學(xué) 核能與新能源技術(shù)研究院，北京 100084）

高溫氣冷堆事故分析中存在一種特有的事故，即進(jìn)氣事故。它是由于一回路壓力邊界發(fā)生破口造成的。雖然發(fā)生的概率極低（低于10-8／（堆·年）），但由于其可能造成嚴(yán)重的后果，因此在HTR-PM 的安全分析中必須對其進(jìn)行深入分析。針對此類事故，清華大學(xué)核能與新能源技術(shù)研究院應(yīng)核安全局的要求對HTR-10進(jìn)行了分析研究［1-2］。分析中主要考慮了擴(kuò)散和自然循環(huán)兩種現(xiàn)象，認(rèn)為空氣僅在擴(kuò)散作用下緩慢進(jìn)入堆芯和反射層流道。該種分析所使用的驗證實驗為JAERI的倒U 型管實驗［1，3］。以此分析為依據(jù)，該研究對反應(yīng)堆模型進(jìn)行了一維簡化，將所有結(jié)構(gòu)都簡化為一維管道。使用此模型計算得出的結(jié)論為：在熱氣導(dǎo)管雙端斷裂的進(jìn)氣事故中，自然循環(huán)建立的時間約為11.3h，穩(wěn)定后流量為0.04kg／s。

但這一模型在分析堆芯流動時過于簡單。由于球床堆芯內(nèi)中心區(qū)域溫度遠(yuǎn)高于外側(cè)區(qū)域，因此在無強(qiáng)迫流動時，堆芯內(nèi)會在浮力作用下形成局部的自然對流。JAERI認(rèn)為這一現(xiàn)象會加快自然循環(huán)的建立［3-4］。為了驗證這一理論，日本山梨大學(xué)搭建了新的倒U型槽實驗平臺［5］。本文針對該實驗，使用CFD 方法對其進(jìn)行模擬計算，驗證CFD 方法解決相關(guān)問題的可行性。另外，利用CFD 方法靈活的特性對影響自然循環(huán)建立時間的因素進(jìn)行更進(jìn)一步的分析，為進(jìn)氣事故的分析提供依據(jù)。

1 實驗裝置

實驗裝置包括1個矩形氣體艙室、2個豎直狹槽以及連接它們的水平通道，具體尺寸示于圖1a。為了產(chǎn)生局部自然對流，左側(cè)狹槽的兩個壁面分別設(shè)計成冷壁和熱壁。熱壁由銅板、電加熱器、絕熱層組成，通過電加熱方式調(diào)節(jié)溫度；冷壁由銅板、冷卻水管道、絕熱層組成，通過冷卻水帶走熱量，試驗中采用的冷卻水溫度不可調(diào)節(jié)。左側(cè)狹槽內(nèi)排布方式如圖1b所示。由于這些結(jié)構(gòu)的存在，左側(cè)狹槽內(nèi)實際流道寬度為20 mm。相比之下，右側(cè)狹槽未做特殊處理。

圖1 倒U 型槽實驗裝置示意圖［4-5］Fig.1 Scheme of reverse U-shaped slots experiment equipment［4-5］

實驗中，用來監(jiān)測自然循環(huán)建立情況的物理量為氣體的溫度和摩爾組分。氣體摩爾組分的測量是通過測量當(dāng)?shù)芈曀俚姆椒ㄩg接求得的，其測點位于兩個狹槽的頂端（圖1b）。氣體溫度通過位于狹槽中的K 型熱電偶測量，其分布位置列于表1。

實驗過程為：1）向裝置內(nèi)充滿高密度氣體（如Ar）；2）關(guān)閉艙室和槽之間的隔板，向槽內(nèi)充滿低密度氣體（如N2）；3）熱壁和冷壁開始工作，持續(xù)一段時間達(dá)到穩(wěn)態(tài)，此過程氣體壓強(qiáng)保持為大氣壓；4）打開隔板，實驗開始。

表1 氣體溫度測點的位置Table 1 Location of thermocouples

為得到自然循環(huán)建立時間與氣體種類的關(guān)系，實驗選取了N2／Ar、Ne／Ar、He／N2、He／Ar 4組氣體組合，其密度比分別為7／10、5／10、1.4／10、1／10。為了研究自然循環(huán)建立時間與局部自然對流強(qiáng)度的關(guān)系，左側(cè)狹槽冷熱壁面溫差選為0、10、30、50、70、100K。

2 CFD模擬驗證

選取商用軟件FLUENT 14.5對該實驗進(jìn)行模擬計算。由于實驗裝置具有對稱特性，故對其進(jìn)行1／2建模以便簡化計算。幾何模型以及網(wǎng)格劃分如圖2 所示，計算區(qū)域包含9 216個六面體網(wǎng)格，時間步長在0.02～1s內(nèi)變化。混合氣體選取多組分模型，混合物的密度、比熱、黏度等物理量采用非可壓縮理想氣體模型或氣體混合模型。由于氣體種類僅兩種，且暫不考慮多孔介質(zhì)問題，因此擴(kuò)散系數(shù)等于兩氣體的二元擴(kuò)散系數(shù)，計算方法為分子動理論。對于左側(cè)狹槽內(nèi)自然對流的分析，為了更準(zhǔn)確地得到黏性底層內(nèi)切應(yīng)力，文中采用低Re 的k-ω 模型，為了滿足模型要求，左側(cè)狹槽內(nèi)網(wǎng)格的y+約為1.2［6］。

圖2 實驗裝置的幾何模型及網(wǎng)格劃分Fig.2 Geometrical model and mesh of experiment equipment

對實驗的驗證，選取文獻(xiàn)［5］中給出數(shù)據(jù)最多的N2／Ar組分，冷熱壁面溫差為50 K。實驗中采用溫度變化和濃度變化兩種方式判斷自然循環(huán)的建立。二者的計算結(jié)果和模擬結(jié)果如圖3、4所示。圖3 為狹槽頂部測點所得Ar濃度的變化，伴隨擴(kuò)散過程，兩側(cè)的N2含量逐漸增加，左側(cè)由于溫度較高，擴(kuò)散過程明顯更快。在擴(kuò)散到一定階段（約140 min），形成穩(wěn)定的自然循環(huán)。從圖3可看出，計算結(jié)果雖與實驗值仍有一定的差距，但其整體趨勢以及自然循環(huán)建立的時間基本一致。圖4 為兩個槽內(nèi)12個熱電偶（表1）所測的溫度變化曲線。由于文獻(xiàn)［5］中并未測量冷卻水的溫度，所以計算所使用的冷壁面溫度只能假設(shè)為常溫（約300K），因此模擬結(jié)果與實驗結(jié)果存在一定差別也在意料之中，但本文所關(guān)注的自然循環(huán)建立時間與實驗結(jié)果是基本相符的。綜上所述，用FLUENT進(jìn)行的數(shù)值計算可較好地重現(xiàn)這一實驗的全過程，這也表明，由于實驗?zāi)M反應(yīng)堆復(fù)雜的結(jié)構(gòu)難度較大，用CFD 模擬的方式來進(jìn)一步研究分析進(jìn)氣事故是合理可行的。

圖3 氬氣摩爾組分的變化Fig.3 Mole fraction change of Ar

3 敏感因素分析

通過與實驗結(jié)果對比，驗證了所建立數(shù)值模型的合理性。在此基礎(chǔ)上，利用該模型可進(jìn)一步研究自然循環(huán)現(xiàn)象的影響因素。文獻(xiàn)［1-2］中基于一維假設(shè)已對這一現(xiàn)象做過一定的分析，提出了諸如溫差、氣體密度比、擴(kuò)散系數(shù)等影響自然循環(huán)建立的因素。但隨著局部自然對流現(xiàn)象的引入，一些新的問題需要進(jìn)行考慮。

圖4 氣體溫度變化Fig.4 Temperature change of gas

3.1 對流強(qiáng)度

文獻(xiàn)［5］研究了局部自然對流強(qiáng)度對自然循環(huán)建立的影響，得到的結(jié)論是：增大左側(cè)狹槽兩壁的溫差，可增強(qiáng)局部自然對流，有利于自然循環(huán)的建立（表2）。但增大溫差的過程中，由于冷壁（冷卻水）的溫度不變，左側(cè)的平均溫度也是升高的，而溫度也是影響自然循環(huán)建立的一個重要因素，因此僅通過這組實驗得出以上結(jié)論是不充分的。

表2 自然循環(huán)建立時間與溫差的關(guān)系［5］Table 2 Onset time of natural circulation with temperature difference［5］

在此基礎(chǔ)上，用CFD 模擬補(bǔ)充一組對比工況：左側(cè)狹槽內(nèi)平均溫度與表2 中ΔT=50K時相同，且兩壁面溫度相等，槽內(nèi)無自然對流。對左側(cè)狹槽內(nèi)氣體溫度做體積平均，所得結(jié)果作為補(bǔ)充實驗的壁面溫度，約為325K。兩組實驗中右側(cè)狹槽頂部Ar組分變化如圖5 所示。在左側(cè)狹槽內(nèi)溫度均勻時，自然循環(huán)需要230min建立；而相同平均溫度下，存在局部自然對流會使左側(cè)狹槽內(nèi)氣體更快地均勻混合，并通過上方的連接通道向右側(cè)狹槽擴(kuò)散，使兩側(cè)Ar含量均增加。因此，局部自然對流的存在會導(dǎo)致自然循環(huán)過程的加速建立。

圖5 右側(cè)狹槽頂部Ar摩爾組分變化Fig.5 Mole fraction change of Ar in right slot

3.2 阻力

球床堆芯、上反射層、下反射層都是會產(chǎn)生較大流動阻力的區(qū)域，這些區(qū)域的阻力是否會對自然循環(huán)的建立造成影響也是研究的內(nèi)容之一。對阻力的研究包括兩個方面：一是局部對流區(qū)域外的阻力；二是局部對流區(qū)域內(nèi)部阻力。

針對第1種情況，本文在左側(cè)狹槽的頂部、底部以及右側(cè)狹槽的底部設(shè)置局部阻力，這3種情況分別記為Case 1、Case 2、Case 3。阻力系數(shù)的設(shè)置原則為：使穩(wěn)定循環(huán)后的局部壓降約為沿程壓降的20%。計算結(jié)果列于表3，局部阻力的存在以及位置均不會顯著影響自然循環(huán)建立的速度。因為在局部對流區(qū)域之外，傳質(zhì)的主要方式是擴(kuò)散，而阻力不影響擴(kuò)散作用。

表3 不同局部阻力情況下自然循環(huán)建立時間Table 3 Onset time of natural circulation in different cases

自然對流存在于球床結(jié)構(gòu)中，數(shù)值模擬中通常使用多孔介質(zhì)模型對球床進(jìn)行簡化。用多孔介質(zhì)模擬左側(cè)狹槽，并設(shè)置不同的各向同性的黏性阻力系數(shù)，分別為104、105、106。在不同阻力系數(shù)下，H=0.3 m 處的水平截面的速度分布如圖6所示?？煽闯?，隨黏性阻力系數(shù)的增大，對流的速度減小，強(qiáng)度變?nèi)?。?.1節(jié)的結(jié)論可知，對流強(qiáng)度會影響自然循環(huán)建立的時間，因此對流區(qū)域內(nèi)阻力會影響自然循環(huán)建立時間。

圖6 不同黏度阻力下自然對流速度分布Fig.6 Velocity of natural convection with different viscosities

4 結(jié)論

本文對日本山梨大學(xué)的U 型槽實驗進(jìn)行了數(shù)值驗證，數(shù)值計算結(jié)果與實驗值符合較好。這不僅驗證了FLUENT 計算模型的合理性，也方便于針對高溫氣冷堆進(jìn)氣事故做一些初步分析。主要結(jié)論如下：

1）通過實驗驗證，F(xiàn)LUENT 程序可很好地模擬由隔板打開到自然循環(huán)建立這一過程。

2）局部自然對流會影響整體自然循環(huán)建立的時間。對流強(qiáng)度越強(qiáng)，自然循環(huán)建立所需時間越短。

3）在對流存在的區(qū)域之外，局部阻力不會對自然循環(huán)建立的時間造成影響。但在自然對流存在的區(qū)域內(nèi)，由于阻力會降低自然對流的強(qiáng)度，因此會推延自然循環(huán)的建立。

以上研究均是建立在簡單模型的基礎(chǔ)上的，對于高溫氣冷堆堆芯的這種復(fù)雜結(jié)構(gòu)還需要更進(jìn)一步研究。

［1］ 劉寶亭.10 MW 高溫氣冷試驗堆進(jìn)氣事故研究［D］.北京：清華大學(xué)，1998.

［2］ YANHUA Z，F(xiàn)UBING C，LEI S.Analysis of diffusion process and influence factors in the air ingress accident of the HTR-PM［J］.Nuclear Engineering and Design，2014，271：397-403.

［3］ TAKEDA T.Research and development on prevention of air ingress during the primary-pipe rupture accident in the HTTR［J］.Nuclear Engineering and Design，2004，233（1-3）：197-209.

［4］ TAKEDA T.Air ingress phenomena in a depressurization accident of the very-high temperature reactor［J］.Nuclear Engineering and Design，2010，240（10）：2 443-2 450.

［5］ TAKEDA T，MIZUNO H，F(xiàn)UNATANI S.Study on mixing process of two component gases in a vertical fluid layer［J］.Nuclear Engineering and Design，2013，271：424-430.

［6］ 陶文銓.數(shù)值傳熱學(xué)［M］.西安：西安交通大學(xué)出版社，2004：385-391.