CARR輻照壓水堆小組件熱工水力分析

尹 皓，鄒 耀，劉興民

(中國原子能科學研究院 反應堆工程研究設計所，北京 102413)

?

CARR輻照壓水堆小組件熱工水力分析

尹 皓，鄒 耀，劉興民

(中國原子能科學研究院 反應堆工程研究設計所，北京 102413)

分析壓水堆4×4小組件在CARR高溫高壓回路中進行輻照考驗時的熱工水力問題。利用計算流體動力學(CFD)軟件對其進行三維數值模擬，以獲得詳細的熱工水力參數。首先，模擬簡化的燃料棒束模型，得出三維溫度與速度分布，并分析了傳熱過程。然后，模擬全尺寸小組件，與棒束模型所得的結果進行對比分析，著重研究其流動，并分析了格架的攪混特性，得出可應用于一維熱工水力程序的攪混因子。結果表明，燃料棒最高溫度可滿足安全性要求，且格架的攪混作用明顯。

壓水堆小組件；熱工水力分析；CFD軟件

中國先進研究堆(CARR)是一座高性能、多用途的反應堆，其主要用途包括核燃料和材料的堆內試驗。當前，CARR已列入國家基地建設項目，國務院已批準為CARR建造高溫高壓試驗回路和相應的功率調節裝置，以便將來進行核燃料元件的堆內試驗，研究核燃料元件在穩態或瞬態工況下的性能，為我國核電大發展開發高性能燃料元件或研發如MOX燃料、環形燃料等新型燃料元件，提高我國核電的安全性和經濟性[1]。

利用CARR開展小組件輻照考驗，獲得輻照性能數據，是實現工程應用的重要環節。此次所研究的壓水堆小組件在模擬壓水堆運行條件下輻照考驗，并在熱室內完成燃料棒的輻照后檢驗，獲得燃料元件在模擬壓水堆環境下的輻照性能數據。本文著重研究其在高溫高壓回路內考驗時的熱工水力問題，分別對比分析簡化燃料棒束模型和全尺寸小組件模型的三維溫度場和速度場。

1 小組件結構

小組件的輻照裝置結構如圖1所示，由吊架管組件、上壓緊板組件、頂部格架、中間格架、底部格架、下壓緊板組件等組成。冷卻劑進出口與考驗回路相連，回路運行時，從主回路來的冷卻劑經冷卻劑進口管進入壓力管，從壓力管和下壓緊板之間的間隙自下向上流動，在底部進入考驗燃料棒束，冷卻燃料棒，由出口管流出，進入考驗回路的堆外部分。

圖1 小組件結構Fig.1 Small assembly structure

2 簡化棒束模擬

在進行棒束模型模擬時，簡化掉上、下壓緊板組件，吊架管組件和4個格架，著重研究燃料棒內部的導熱與棒束和冷卻劑之間的對流換熱。故只建立12根燃料棒和4根導向管在壓力管內的三維模型，簡化棒束模型示于圖2。

圖2 簡化棒束模型Fig.2 Simplified rod bundle model

模型建立后，對其進行網格劃分，由于燃料棒包含很薄的鋯包殼和氦氣間隙，故其內部網格劃分相對較密，其網格劃分橫截面如圖3所示。

圖3 簡化棒束模型網格橫截面Fig.3 Mesh cross sectionof simplified rod bundle model

由于小組件放置于堆芯的一側進行輻照考驗，故其每根燃料棒的功率分布不均勻，靠近CARR堆芯的一側釋熱大，遠離的一側釋熱小。經物理計算每根燃料的功率分布不均勻因子列于表1。棒束編號示于圖4，平均線功率為31 kW/m。由表1可見，3號燃料棒與7號燃料棒最靠近CARR堆芯，故其功率分布不均勻因子最大。計算前進行邊界條件的設置(表2)。

表1 功率分布不均勻因子Table 1 Power distribution uneven factor

圖4 棒束編號Fig.4 Number of rod bundle

表2 簡化棒束模型邊界條件Table 2 Boundary condition of simplified rod bundle model

圖5 棒束溫度分布Fig.5 Temperature distribution of rod bundle

經過計算得到棒束溫度分布(圖5)。通過分析，燃料棒溫度由內而外逐漸降低。3號燃料棒與7號燃料棒最靠近CARR堆芯，故其溫度最高，其芯部最高溫度可達2 282 K，遠低于其相應燃耗下的熔點3 138 K。此外，經計算包殼處溫度的最高值為741.157 59 K，足以保證其安全性。

3 全尺寸小組件熱工水力數值模擬

在全尺寸小組件的數值模擬中主要針對最為復雜的格架部分對小組件進行模擬分析。由于格架結構非常復雜，尤其是彈簧、剛凸和攪混翼。本文簡化掉燃料棒束，在燃料棒外壁處施加熱流密度，以此來等效代替燃料棒的產熱，此簡化既符合實際的熱工條件，又能節省網格數量，能保證模擬結果的準確性。

構造幾何模型時要著重進行簡化。由于機械制造與加工的要求，格架上的結構十分復雜，其中包括多處的窄縫、定位凸起、焊點等，這些結構會給網格劃分帶來很多困難。根據經驗可知這些細小的結構對整個組件的熱工水力影響非常小，為此將這些結構簡化掉。將彈簧的倒角與圓角全部拉直，并填充彈簧兩邊的小空隙，彈簧與條帶之間的線接觸簡化為面接觸，如圖6所示。

圖6 簡化前、后的倒角Fig.6 Chamfer before and after simplification

格架處網格劃分最大的困難在于彈簧、剛突與燃料棒及導向管間的接觸空隙趨近于零，不能生成計算網格。目前定位格架模型簡化主要有3種方法：1) 將彈簧、剛突與燃料棒間的線接觸改為面接觸；2) 在彈簧、剛突與燃料棒間留一小空隙；3) 完全忽略彈簧與剛突的影響。

根據實際結構，同時盡量減小模型簡化對計算的影響，在彈簧、剛突與燃料棒間采用第2種方法，留出一微小縫隙，而在彈簧、剛凸與導向管間采取第1種方法，將線接觸改為面接觸。這是由于燃料棒有發熱率，改為間隙是比較合理的，而導向管不發熱，故改為面接觸較為簡單。具體簡化情況如下。

改變燃料棒與彈簧、剛凸的接觸方式，使燃料棒與彈簧、剛凸之間產生0.3 mm左右的縫隙。在導向管的四周添加4個非常小的凸臺，將線接觸改為面接觸。間隙處簡化前、后示于圖7。

圖7 間隙處簡化前、后示意圖Fig.7 Schemes before and after gap simplification

就此完成了整個組件模型的簡化，然后建立壓力管幾何模型并流出進出口段。用長方體流道對小組件的實體部分求差，得到的流體區域即為計算域。將流體域沿流動方向劃分為4部分，對每部分進行單獨的網格劃分。使用非結構化四面體網格加拉伸網格的形式，格架部分由于結構復雜使用非結構化網格，流道部分使用拉伸出的三棱柱網格。在無滑移壁面劃分有附面層網格。具體步驟示于圖8。

圖8 網格劃分流程Fig.8 Flow chart of meshing

劃分網格時要調整參數，使間隙處網格密，其他區域網格稀疏。同時，附面層處要生成至少3層網格，燃料棒與格架間隙處網格示于圖9。經統計最終的網格數量達1 599萬。

由于建立的幾何結構模型是一體的，故4段流體域沿z軸的坐標是確定的。所以將4段流體域導入CFD中，4段流體域是銜接起來的。只需將4段流體域間的交界面設為interface，并將4段流體域分別設置物性參數即可。邊界條件列于表3。

圖9 燃料棒與格架間隙處網格Fig.9 Mesh between fuel rod and grid gap

表3 CFD邊界條件Table 3 Boundary condition in CFD

設置好參數后，在CFD中進行計算并進行結果處理。z=0.9 m與中間格架處的溫度分布云圖示于圖10。

速度矢量的沿程分布示于圖11。由圖11可看出，入口處的速度是沿流動方向均勻分布的，當遇到下壓緊板時速度方向發生改變，也有突增，局部甚至可達12 m/s。此后，速度的變化較均勻，在格架處有少許增加，隨后趨于穩定。

圖10 z=0.9 m和中間格架處的溫度分布云圖Fig.10 Temperature contours at z=0.9 m and middle grid

圖11 速度矢量沿程分布Fig.11 Velocity vector distribution along fuel rod

雖然局部區域的速度可達12 m/s，但經計算，其沿x方向與y方向分量的最大值僅為2.5 m/s。此外，速度的最大值出現在下壓緊板處，即小組件的一端，此處是有固定即有力來約束的，故并不會給組件帶來很大的橫向力，不會發生流致振動。

大量文獻報道了格架攪混翼可在下游的一定距離內產生并保持一定強度的渦流與交叉流[2-5]。形式不同的評價因子被定義來研究攪混翼所引起的流動特性，雖然形式上略有差別，但所描述的流動本質是相同的。常采用的是渦流攪混因子和交叉流攪混因子。

渦流攪混因子S可用來衡量渦流對流體的攪混作用，為無量綱數，定義如下：

其中：R為子通道中心到燃料棒表面的距離，m；r為通道截面上流域內各點到中心的距離，m；vt為橫截面上引起渦流的橫向速度，m/s；ua為橫截面上各點的軸向速度，m/s。

交叉流攪混因子F用來衡量通道間交叉攪混的強度，為無量綱數，定義如下：

其中：s為燃料棒表面間距，m；vcross為穿過棒間橫截面的交叉流速分量，m/s；ubulk為截面軸向流速，m/s。

采用中間兩個格架之間的全部流體域, 總長度約為240 mm，并取4×4小組件最中心處的子通道來研究攪混因子的變化。通過三維結果輸出的數據來計算，所得攪混因子沿程變化示于圖12。

由圖12可見，以當量直徑Dh為單位，渦流攪混因子剛開始時沿流動方向逐漸增大，增至8Dh左右達到峰值，然后逐漸衰減。交叉流攪混因子沿流動方向先略有減小后逐漸增大，在6Dh左右時達到最大值，隨后逐漸減小，到12Dh時下降趨勢變緩。由此可見，格架的攪混作用要先經歷一充分發展的過程再經歷一逐漸衰減的過程，總體趨勢屬于先增大至峰值后逐漸減小。

圖12 攪混因子沿程變化Fig.12 Mixing factor change along fuel rod

許多一維軟件，如COBRA，均需要格架攪混性能參數作為輸入文件，不同的軟件攪混因子的定義不同，但總體而言均是以上述兩個攪混因子作為計算依據。

4 結論

1) 燃料棒內部的最高溫度正比于燃料棒的功率分布不均勻因子，從而可驗證靠近CARR堆芯的一側燃料棒芯部溫度較大，且其芯部最高溫度與其他燃料棒芯部最高溫度之間的溫差最高可達500 K。

2) 格架與上、下壓緊板對冷卻劑溫度的影響不是太大，但對流動速度影響較大。尤其是下壓緊板，由于是流體流經的第1個組件，其對流動的影響不僅表現在改變速度方向上，還使局部最高速度達12 m/s。

3) 通過對攪混因子的計算表明格架的攪混翼對冷卻劑的攪混作用十分明顯，尤其是在攪混翼下游6Dh～8Dh之間達到最大，其后逐漸減小。

4) 本次計算的結果可為今后針對該壓水堆小組件進行的輻照考驗試驗提供很好的理論依據，燃料棒溫度分布的大致趨勢、格架對速度的影響以及格架攪混因子的變化趨勢對其他結構的燃料組件也有一定的借鑒性。此外，采用CFD軟件進行三維熱工水力模擬普遍適用于各類型的燃料棒束與燃料組件。

[1] 刁均輝. CARR堆內試驗燃料棒束熱工水力性能分析[D]. 北京：中國原子能科學研究院，2008.

[2] CUI X Z, KIM K Y. Three-dimensional analysis of turbulent heat transfer and flow through mixing vane in a subchannel of nuclear reactor[J]. Journal of Nuclear Science and Technology, 2003, 40(10): 719-724.

[3] LEE C M, CHOI Y D. Comparison of thermo-hydraulic performances of large scale vortex flow (LSVF) and small scale vortex flow (SSVF) mixing vanes in 17×17 nuclear rod bundle[J]. Nuclear Engineering and Design, 2007, 237(24): 2 322-2 331.

[4] IN W K, CHUN T H, OH D S, et al. CFD analysis of turbulent flows in rod bundles for nuclear fuel spacer design[C]∥Transactions of the 15th International Conference on Structural Mechanics in Reactor Technology. Seoul: SMiRT, 1999: 365-372.

[5] IN W K, CHUN T H, SHIN C H, et al. Numerical computation of heat transfer enhancement of a PWR rod bundle with mixing vane spacers[J]. Nuclear Technology, 2008, 161(1): 69-79.

Thermal-hydraulic Analysis of PWR Small Assembly for Irradiation Test of CARR

YIN Hao, ZOU Yao, LIU Xing-min

(ChinaInstituteofAtomicEnergy,P.O.Box275-33,Beijing102413,China)

The thermal-hydraulic behaviors of the PWR 4×4 small assembly tested in the high temperature and high pressure loop of China Advanced Research Reactor were analyzed. The CFD method was used to carry out 3D simulation of the model, thus detailed thermal-hydraulic parameters were obtained. Firstly, the simplified model was simulated to give the 3D temperature and velocity distributions and analyze the heat transfer process. Then the whole scale small assembly model was simulated and the simulation results were compared with those of simplified rod bundle model. Its flow behavior was studied and flow mixing characteristics of the grids were analyzed, and the mixing factor of the grid was calculated and can be used for further thermal-hydraulic study. It is shown that the highest temperature of the fuel rod meets the design limit and the mixing effect of the grid is obvious.

PWR small assembly; thermal-hydraulic analysis; CFD software

2014-03-20;

2014-04-16

尹 皓(1989—)，女，河北邯鄲人，碩士研究生，從事反應堆熱工水力研究

TL334

A

1000-6931(2015)06-1069-06

10.7538/yzk.2015.49.06.1069