CEFR堆頂固定屏蔽冷卻系統流動特性數值研究

馬 曉，張東輝

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

?

CEFR堆頂固定屏蔽冷卻系統流動特性數值研究

馬 曉，張東輝

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

堆頂固定屏蔽在中國實驗快堆(CEFR)中承擔著重要功能，對其進行充分冷卻極其重要。本文采用CFD方法對該設備的冷卻系統進行了三維數值研究，詳細分析了該冷卻系統的流動特性和水力學設計，并對設計中的不足提出了優化建議。研究表明，該冷卻系統基本可滿足要求，但部分環節需要優化。將調節閥盡量均勻布置可改善水平風道流場分布；入口處設置兩道通風孔可提高豎直風道內空氣流動的均勻性；調節閥開度應適當增加以進一步滿足流量分配需求。該研究可為CEFR運行安全和類似冷卻系統的設計提供參考。

中國實驗快堆；堆頂固定屏蔽；流動特性；數值研究

中國實驗快堆(CEFR)堆頂固定屏蔽位于反應堆堆坑頂部，承擔著屏蔽中子、γ輻射和熱輻射的重要功能，并且作為支承平臺為相關設備提供支承。反應堆功率運行時，堆頂固定屏蔽處在一高溫工作環境中，承受反應堆容器和各種貫穿件的散熱。為保證該設備具有足夠的結構強度，以滿足各種預期功能，其冷卻系統必須滿足熱工要求。堆頂固定屏蔽冷卻系統風道狹窄、結構復雜，難以通過理論計算和實驗測量掌握內部各位置詳細流動情況。為了詳細研究冷卻系統內部的具體流動和流量分配，全面考察堆頂固定屏蔽冷卻系統的設計能否滿足要求，本文采用計算流體力學(CFD)方法對該冷卻系統的流動特性進行計算分析。

1 堆頂固定屏蔽及其冷卻系統簡介

CEFR堆頂固定屏蔽位于反應堆堆坑頂部，由匯流箱、屏蔽箱體、屏蔽塊、環形支承裙板和密封組件組成(圖1)，主要功能是屏蔽中子和γ輻射、屏蔽熱輻射以及作為工作平臺[1]。

a——進風管；b——匯流箱；c——空氣調節閥；d——第2層水平風道；e——第3層水平風道；f——支承裙板；g——豎直風道；h——屏蔽箱體1～8分別為第2層水平風道8個區塊的編號

堆頂固定屏蔽是一大尺寸的箱式金屬結構件，最大直徑為10 480 mm，總高度為2 200 mm，其中匯流箱高度為150 mm，箱體厚度為1 255 mm，下環形裙板高度為795 mm。

堆頂固定屏蔽由多層材料組成，主要由鋼板搭建設備的主體結構，縱向有8層鋼板、4層混凝土、1層礦渣棉和3層供冷卻空氣流動的水平風道，其中水平風道高度與其徑向尺寸相比非常狹窄。在屏蔽箱體的橫截面上有沿徑向成輻射狀分布的8條筋板，將屏蔽箱體分成8個區(各區編號見圖1)，13個空氣調節閥為每個區域分配不等的流量。

為配合反應堆實際需要，在堆頂固定屏蔽上開有27個貫穿孔，一回路主循環泵、中間熱交換器等設備經由這些孔道從堆容器中伸到堆頂固定屏蔽上部，主容器支承頸也穿過堆頂固定屏蔽。冷卻系統為帶走這些部件的熱量專門設置了一系列豎直風道。另外，該設備上還存在21個盲孔。由以上分析可發現堆頂固定屏蔽是一結構復雜、空間擁擠的大型設備[2-3]。

由于堆坑內存在較多余熱，且各種貫穿件也會散發出較多熱量，因此堆頂固定屏蔽處在一高溫的工作環境，為保證該設備具有足夠的結構強度，完成各種設計功能，需對其進行足夠的通風冷卻。

堆頂固定屏蔽冷卻系統設計上既要保證整體冷卻效果，又保持一定的獨立性。為便于冷卻，該冷卻系統設計中采用了多次流量分配、水平風道與豎直風道交錯冷卻、水平風道分區塊、利用閥門調節空氣流量等多種手段，最終形成的冷卻系統非常復雜，其冷卻系統流程圖如圖2所示。

2 計算模型及方法

CFD方法是流體力學的一種有效研究手段，能夠充分利用計算機強大的計算能力解決各種復雜的流動問題，具有適用范圍廣、研究成本低、結果全面和直觀的特點[4]。由于堆頂固定屏蔽冷卻通風流道復雜、重要細節較多，難以進行理論分析求解和布置測點，因而適合采用CFD方法對其流場進行數值研究。本工作采用CFX軟件進行計算分析。

實際運行中，堆頂固定屏蔽下部熱、上部冷，各貫穿件周圍溫度較高，其不均勻的溫度分布會使冷卻系統各處空氣的溫度、物性、所受浮升力等變得不均勻，進而對水平風道空氣流態及系統壓降等方面產生影響。對其進行溫度場和流場耦合計算可詳細考慮溫度對流場的影響。但堆頂固定屏蔽尺寸巨大、結構極其復雜，進行溫度場和流場的整體計算時流固高度熱耦合、熱源很復雜，且需要極多的網格和計算資源。因此，本文暫不考慮溫度不均勻性對流場分布的影響和局部可能存在的自然對流。

圖2 堆頂固定屏蔽冷卻系統流程Fig.2 Process of cooling system for fixed shielding platform

本文建立了與實際結構尺寸1∶1的計算模型，如圖3所示，包括匯流箱、第2層水平風道、第3層水平風道、10類豎直風道等，對于一些影響流動和流量分配的重要細節區域(如13個空氣調節閥和豎直風道入口處的通風窗、大量通風孔等結構)也按照實際情況建模和分析，充分考慮這些結構對流動造成的影響。

在阻力環節建模方面，通風窗共有兩類，分別是中間熱交換器風道入口處均勻分布的4個100 mm×100 mm×10 mm的通風窗和獨立熱交換器風道入口處均勻分布的3個85 mm×95 mm×20 mm的通風窗，流體經過這些通風窗后進入一環形腔，然后流經通風孔進入對應的豎直風道中。各豎直風道入口處均有1圈沿圓周均勻分布、數目不等的通風孔，其直徑分為30 mm和50 mm兩種，熱交換器風道入口通風孔的大小可通過進風調節器來調節，建模時采用設計開度對應的尺寸。調節閥建模時，對于閥門入口窗、內外圓筒、螺桿、螺帽、肋板等主要結構均按照實際尺寸建模，對于螺紋等少量幾何細節由于對流動影響非常小，建模時進行了省略。調節閥的幾何模型如圖4所示。

a——進風管；b——匯流箱；c、h、i、j、l——豎直風道；d——調節閥；e——第2層水平風道；f——第3層水平風道；g、k——通風窗和進風調節器

圖4 調節閥幾何模型Fig.4 Geometry model of regulating valve

本文采用四面體非結構網格劃分方法，流動模型采用k-ε模型，壁面附近湍流處理采用壁面函數方法，并考慮重力對流動的影響。計算針對額定工況進行，此時冷卻系統空氣流量為46 692 m3/h，空氣平均溫度為45 ℃，空氣物性選取標準大氣壓、45 ℃時的空氣物性數據[5]。邊界條件選用速度入口和壓力出口，入口速度為額定工況對應的入口流速25.27 m/s，出口壓力為0 Pa，壁面采用無滑移邊界條件。

針對網格數量對計算結果的影響進行了敏感性分析，以進出口之間的壓降作為敏感性分析所考察的變量，結果列于表1。由表1可知，當網格數從527萬增加到703萬時，進出口壓降只變化了0.6%，這說明計算結果基本達到了網格無關解，因此最終的網格數取703萬。

3 計算結果和分析

3.1 冷卻系統整體流動特性

圖5為堆頂固定屏蔽各層水平風道的流場。由圖5可見：通風窗和通風孔等節流處、幾何結構較狹窄位置、送風量較集中位置空氣流速明顯較大，其余區域流速較低，流動較為平穩；匯流箱流速整體高于第2層水平風道流速和第3層水平風道流速，流動最為劇烈，第3層水平風道整體流速最慢。

表1 網格敏感性分析Table 1 Mesh sensitivity analysis

a——匯流箱；b——第2層水平風道；c——第3層水平風道

匯流箱體積較小，孔道卻有很多，這使得匯流箱流道非常狹窄和擁擠，流體在流經狹窄流道、通風窗、通風孔和調節閥時流速顯著增大。考察匯流箱外圍近乎均勻分布的6個進風管道，可以發現靠近主泵的兩個進風管由于處在匯流箱非常狹窄的角落，導致其周圍流體流速明顯較其余4個進風管附近流速大，這增大了匯流箱內流體的流動不均勻性，這兩個進風管的位置設置存在問題，需進一步優化。匯流箱內最高流速可達60 m/s以上，出現在幾何結構較狹窄位置和換熱器通風孔區域，距進風管較近的位置由于空間狹窄，導致該位置流速達30 m/s以上；第2層水平風道流通空間相對大一些，最大流速不到39 m/s，出現在主泵風道附近區域，豎直風道入口通風孔處流速較高，普遍在10 m/s以上，主泵風道和換熱器風道周圍區域出現了高流速區；第3層水平風道整體流速最慢，最高流速出現在該風道外圍通風孔附近區域，流速為5～10 m/s，其余區域流速普遍小于3 m/s，且該風道出現了多處漩渦結構。

圖6 冷卻系統壓力分布Fig.6 Pressure distribution of cooling system

冷卻系統壓力分布如圖6所示，匯流箱內壓力最高，流體流經通風孔后壓力明顯降低。入口處平均壓力為2 828 Pa，出口壓力為0 Pa，系統壓降為2 828 Pa。

該冷卻系統在設計上應用了分層冷卻的方式，通過兩次流量分配，由水平風道為各豎直風道分配流量。3層水平風道、兩次流量分配的設計既保證了對整個固定屏蔽的整體冷卻，又增加了空氣流量調節的靈活性。

該冷卻系統將大尺寸風道化整為零，分別進行控制，通過將第2層水平風道分為8個區塊并設置數目不等的空氣調節閥的方式，提高了對大區塊流量分配進行調節的效率和精細程度。另外，分成8個區塊可將冷卻空氣限制在各區塊內部流動，減少大范圍橫向流動的出現，使流場更為優化，減弱對設備穩定性的影響。

堆頂固定屏蔽冷卻系統通過兩次流量分配、分層、分塊、獨立的調節閥和數量極多的通風孔實現了對該設備內部結構(尤其是混凝土)和各貫穿件進行多層次、全方位的有效冷卻。

3.2 調節閥位置分布對流場的影響

圖7為堆頂固定屏蔽第2層水平風道第2區和第6區的流場。第2層水平風道第2區和第6區幾何結構和尺寸相同，唯一的差別是3個空氣調節閥的布置方式不同：第2區的3個閥門集中布置在主泵風道左上方，而第6區調節閥分布相對均勻一些，兩個位于主泵風道左下角，另一個位于主泵風道右下角。從圖7可看到，空氣調節閥位置對水平風道流場形態有顯著影響：較集中的調節閥布置方式造成流場出現大范圍高速區，這不利于設備的穩定性(圖7a)；較均勻的調節閥布置方式有利于形成更優的流場，可避免出現大范圍高速區，有利于設備穩定性(圖7b)。通過以上分析，建議盡量均勻布置調節閥，將流動對設備穩定性的影響減至最小。

3.3 通風孔對流場的影響

圖8為主泵風道和中間熱交換器風道的局部流場。從圖8可清楚看出，主泵風道豎直段空氣流動發生了偏離豎直向下方向的情況，而在中間熱交換器風道，流體在風道豎直段保持豎直向下流動，沒有發生偏離。兩類風道入口通風孔附近的流動情況示于圖9。由于中間熱交換器風道內設置了兩道通風孔，導致空氣在連續兩次流量分配后變得更加均勻，空氣在豎直段均向下均勻流動，沒有偏離。而主泵風道只設一道通風孔，不足以使空氣經過一次流量分配后變得均勻，因而空氣流動發生偏離。與之相似的其他豎直風道也出現了類似現象。

圖7 第2層水平風道第2區(a)和第6區(b)的流場Fig.7 Flow fields of part 2 (a) and part 6 (b) in the second horizontal channel

圖8 主泵風道(a)和中間熱交換器風道(b)的流場Fig.8 Flow fields of primary pump channel (a) and intermediate heat exchanger channel (b)

圖9 主泵風道入口(a)和中間熱交換器風道入口(b)的局部流場Fig.9 Local flow fields of inlet region for primary pump channel (a) and inlet region for intermediate heat exchanger channel (b)

圖10為中間熱交換器風道流場截面圖，從圖10可清楚看到通風窗和通風孔對流體的節流作用，流體經過這些結構時流速迅速增大，而進入豎直風道后流動變得均勻緩慢。另外，通風窗外側區域存在明顯漩渦流動，在從通風孔進入的高速流體上下區域出現了兩個明顯的漩渦結構，這與流體向低壓區高速流動夾帶周圍流體轉動的預期相吻合。流入中間熱交換器風道的流體就是這樣經過兩道節流結構的作用而變得更加均勻。

圖10 中間熱交換器風道流場截面圖Fig.10 Flow field on section plane of intermediate heat exchanger channel

通過以上分析可發現，采用兩道節流孔能增強豎直風道空氣流動的均勻性。由于豎直風道內空氣均勻流動有利于對設備進行冷卻和減小對設備穩定性的影響，建議今后類似冷卻系統盡量采用兩道節流裝置，以提高豎直風道空氣流動均勻性。

豎直風道空氣流動的均勻性還與第2層水平風道流動均勻性相關，增強流動均勻性還可從優化6個進風管和13個空氣調節閥的位置分布來實現，使空氣在水平風道中分布更加均勻合理，避免某些進風管道或閥門位置不合理導致水平風道出現流動極不均勻的情況。

3.4 流量分配情況

表2、3分別列出第2層水平風道各區塊和6個換熱器風道的流量分配情況。經計算發現，流入全部區塊流量較所需流量低13.8%(從而換熱器風道以外的其他豎直風道流量較所需流量低13.8%)，流入6個換熱器風道的流量較所需流量高20.6%。這說明按照設計文件中空氣調節閥的開度[6]會導致第2層水平風道流量顯著偏低、換熱器風道流量顯著偏高，這對于依靠從第2層水平風道取風來冷卻的很多貫穿件是十分不利的。

表2 各區塊的流量分配Table 2 Distribution of flow rate for each part

在實際運行中，為保證各風道流量盡量與所需流量一致，應適當增大調節閥的開度，以增加流入第2層水平風道各區塊及對應豎直風道的冷卻空氣流量。

表3 6個換熱器風道的流量分配Table 3 Distribution of flow rate for six heat exchanger channels

4 結論

本文利用CFD方法對額定工況下CEFR堆頂固定屏蔽冷卻系統的流場進行了三維數值模擬研究，主要結論如下。

1) 3層水平風道內空氣平均流速自上而下逐漸降低，第3層水平風道內空氣流動尤其緩慢，并存在很多局部漩渦；

2) 豎直風道流場呈現出不同的特點，部分豎直風道空氣均勻豎直向下流動，另一些豎直風道內空氣流動出現了偏離和不均勻的情況，在豎直風道入口處應盡量采用兩道通風孔的設計，提高流動均勻性，增強冷卻效果；

3) 調節閥位置分布對流場有明顯影響，不合理的位置分布會造成大范圍的高流速區，對設備穩定性有不利影響，應當使調節閥盡量均勻布置，改善流場分布；

4) 按照設計文件中的閥門開度，第2層水平風道流量低于設計流量，換熱器豎直風道流量高于設計流量，實際運行中為了更好地達到冷卻效果，應合理增大調節閥開度，提高第2層水平風道各區塊的冷卻空氣流量。

堆頂固定屏蔽在CEFR中承擔重要功能，需要對其進行足夠的冷卻。本工作對其冷卻系統的流動特性進行了詳細分析，可對今后類似冷卻系統的設計提供一定參考，對于快堆運行安全也具有重要參考意義。

[1] 中國原子能科學研究院中國實驗快堆工程部. 中國實驗快堆最終安全分析報告[R]. 北京：中國原子能科學研究院，2008.

[2] 中國原子能科學研究院中國實驗快堆工程部. 堆頂固定屏蔽技術規格書[R]. 北京：中國原子能科學研究院，2003.

[3] 中國原子能科學研究院中國實驗快堆工程部. 堆頂固定屏蔽施工圖冊[R]. 北京：中國原子能科學研究院，2006.

[4] 王福軍. 計算流體動力學分析[M]. 北京：清華大學出版社，2004.

[5] 楊世銘，陶文銓. 傳熱學[M]. 北京：高等教育出版社，2006.

[6] 雒曉衛. 堆頂固定屏蔽通風系統流量調節閥水力設計[R]. 北京：清華大學核能與新能源技術研究院，2001.

Numerical Research on Flow Characteristic of Cooling Syste for CEFR Fixed Shielding Platform

MA Xiao, ZHANG Dong-hui

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

The fixed shielding platform plays important function in China Experimental Fast Reactor (CEFR), and it’s extremely important to sufficiently cool it. In this paper, CFD method was used to perform three-dimensional numerical research on the cooling system for fixed shielding platform. Flow characteristic and hydraulic design of the cooling system were analyzed detailedly. The optimization suggestions were put forward on the defects in the design. The results show that the cooling system can basically meet the requirements, but some parts need to be optimized. Flow field distribution in horizontal channels will be improved with more uniformly distributed regulating valves. The uniformity of airflow in vertical channels can be improved with two sets of throttles at inlets. The openness of regulating valves should be properly increased to further meet the requirements of flow distribution. This work can provide reference for operation safety of CEFR and design of similar cooling systems in the future.

China Experimental Fast Reactor; fixed shielding platform; flow characteristic; numerical research

2014-03-20;

2014-11-25

馬 曉(1988—)，男，山西長治人，碩士研究生，從事反應堆熱工水力研究

TL333

A

1000-6931(2015)07-1220-07

10.7538/yzk.2015.49.07.1220