CEFR中間熱交換器一次側數值模擬

馮預恒，王一新，趙 勇，周志偉，楊紅義

(中國原子能科學研究院 中國實驗快堆工程部，北京 102413)

中國實驗快堆(CEFR)的中間熱交換器(IHX)是一、二回路的邊界，一回路主冷卻系統的冷卻劑鈉通過IHX將熱量傳遞給二回路，并將放射性鈉隔離在一回路內。CEFR共布置4臺IHX，處于一回路外熱池內。IHX體積龐大，內部布置數百根換熱管，如使用傳統三維整體建模方法，計算量大且難以完成。因此IHX整體三維數值模擬很少使用。一般采用一、二維的方法，如單管換熱模型，在確定換熱管換熱面積和冷卻劑流量的條件下，通過數次迭代計算出換熱管的換熱系數，以此確定換熱量。

傳統計算方法有兩點假設：1) 所有進入IHX的流量皆完成充分換熱，且貢獻相同；2) 換熱面積即是換熱管的總面積之和，不考慮換熱管布置等因素的影響。以上假設對回路式布置的IHX影響不大，但對于快堆一回路池式一體化布局，其IHX內為彎式換熱管，若流體不考慮矢量方向，設備忽略換熱管位置，則極大降低了IHX的計算精度。因此，在工程上采用提高計算的保守性來保證設備的安全。

本工作采用大型通用商業計算流體力學(CFD)軟件CFX對CEFR一回路IHX一次側的1/6進行三維穩態模擬。

1 設備簡介

IHX[1]主要由管束、壓力室、排放室、中心管、屏蔽部件和保護套組成。換熱管束被焊接在上、下兩塊管板上，沿直線管段在管束的軸向高度上安裝有定位帶。殼體結構及換熱管布置示于圖1、2。IHX的主要設計參數列于表1。

圖1 IHX殼體結構示意圖

圖2 換熱管布置示意圖

表1 IHX的主要設計參數

IHX內由換熱管和固定拉桿組成，其中，換熱管513根，固定拉桿27根。換熱管和固定拉桿均有補償性彎曲段，用以補償換熱管、固定拉桿和內筒體之間的熱膨脹差。彎曲段以環排方式進行管束的安裝。沿直管段在管束的軸向高度上安裝有定位帶，它由平板帶和波紋帶組成并焊在拉桿上，其作用是減少換熱管的振動。換熱管的主要參數列于表2。

表2 換熱管的主要參數

本工作模擬1/6 IHX一次側，驗證其在額定運行工況下的換熱能力： 1) IHX的出口溫度是否高于354 ℃；2) 研究IHX出口處的流體分布，驗證此處的溫度測點布置是否合理；3) 研究進、口處流體的分布，為熱池計算提供更精確的三維出口邊界條件；4) 根據IHX內的溫度分布，為系統程序提供精準的換熱管熱工流體參數；5) 以此模型的穩態計算為基礎，為模擬IHX換熱管泄漏作技術準備。

2 模型和邊界條件

2.1 模型

圖3 IHX一次側1/6模型

IHX一次側1/6模型示于圖3。IHX內的換熱管和固定拉桿共9環，最內環36根，從內向外，每環增加6根。IHX的換熱管設計非常特別，每環的換熱管長度相同，中部設計成同軸彎曲，每環的各根換熱管彎曲并設計成同環面。由內環向外環，彎曲度逐漸減小，具體結構示于圖4。

圖4 換熱管結構

IHX一次側入口開有189個φ50 mm的開孔，孔間距75 mm。出口開有231個φ50 mm的開孔，孔間距75 mm。模型的進口為32個，出口為39個，即進、出口分別比實際情況多0.5個開孔。為節省模型的計算網格數目，IHX的定位帶被省略。計算采用k-ε湍流模型，壁面對紊流的影響采用標準情況，考慮重力的影響，各流量下流道內初始流速以理論計算為準。其他輸入和模型選擇默認值。在模擬計算時，采用光滑壁面，自動生成四面體非結構化網格，同時在結構復雜和流動間隙處采用密集網格。在保證計算精度的前提下，采用網格一體化、模型分區等技術，并盡量減少網格數量。

2.2 邊界條件

1) 入口邊界條件

選擇距離IHX一次側入口200 mm為計算的入口邊界。入口流量為78.95 kg/s，入口鈉溫為515 ℃。

2) 出口邊界條件

選擇距離IHX一次側出口下方和外方200 mm處為計算的出口邊界。

3) 傳熱管條件

傳熱管的壁厚被簡化，壁面加熱系數取自IHX熱工計算說明書。壁面溫度的溫降從入口(514 ℃)到出口(310 ℃)簡化為均勻線性變化。

4) 傳熱管一次側管壁傳熱系數

計算說明書取值6.1 kW/(m2·℃)，本工作取值5 kW/(m2·℃)。

5) 物性參數

流質為鈉，流體內壓力為常壓。為真實模擬額定工況下冷池的流動及溫度的影響，各物性參數皆隨溫度而改變。液態鈉的密度ρ、黏度μ、比熱容c、導熱系數λ具體計算公式如下。

ρ=16.018 5(59.356 6-7.750 4×

10-3(1.8t+32)-0.287 2×10-6(1.8t+

32)2+0.060 3×10-9(1.8t+32)3)

μ=0.123 5ρ1/3exp(697ρ/1.8t)t≤500 ℃

c=41.86(0.389 352-1.105 99×

10-4(1.8t)+3.411 78×10-8(1.8t)2)

λ=1.729 58(54.306-1.878×10-2(1.8t+

32)+2.091 4×10-6(1.8t+32)3)

式中，t為材料的溫度。

計算區域外邊界假設為絕熱，即未考慮換熱器殼體與冷、熱池間的換熱。

3 計算結果及分析

3.1 IHX一次側溫度分布

IHX一次側溫度分布示于圖5。一回路內515 ℃的高溫鈉進入IHX，在換熱管壁外被冷卻為356 ℃。沿著低溫鈉流動方向，一次側溫度分布由高到低，在彎管段冷卻效果最好。主要冷卻工作段在下管段,根據其效率可分為5部分。

圖5 IHX換熱管溫度分布

1) 入口腔室

515 ℃的高溫鈉在進入入口腔室后被冷卻，部分沿換熱管下沉，進入IHX。剩余大部分繼續在入口腔室內向前流動，進入IHX內圈換熱管，這導致入口腔室內部溫度低于外部區域溫度。入口腔室平均溫降為6 ℃。

2) 上段直管區

上段直管區是IHX的三大工作區之一。上段直管總長0.77 m，熱鈉在換熱管壁外向下流動并冷卻，在此區域鈉溫下降約21 ℃。

3) 中部彎管區

彎管區水平高度1.074 m，進入彎管段后，熱鈉被迅速冷卻，平均溫降40 ℃。彎管區的冷卻效率較直管區平均提高35%。

4) 下段直管區

下段直管總長1.23 m，是IHX的主要工作區，一次鈉平均溫降達83 ℃，占整體溫降的50%。一次鈉流出此區域后，IHX內的平均鈉溫已與冷池的平均溫度基本相同，為360 ℃左右。

5) 出口腔室

出口腔室總長0.425 m，比入口腔室略大。在其外殼體，由上至下均勻開有39個φ50 mm的開孔。鈉流出下段直管區進入底部出口腔室后，大部分流體在管道縫隙間繼續向下流動，只有部分冷鈉從出口腔室的上排開孔流出。

出口腔室內的徑向溫度分布為外冷、內熱，與入口腔室的溫度分布相反，液態鈉在出口腔室的平均溫降為9 ℃。

3.2 IHX一次側速度分布

IHX一次側速度分布示于圖6。流體在進入IHX后，其速度分布基本規律與換熱管保持一致，但各段間、軸向和周向的速度分布各有特點。

圖6 IHX一次測流場

1) IHX一次側軸向分布

因一回路流體流經換熱管外側，傳熱管對流動的影響非常大，流體的流動方向與傳熱管基本保持一致。只是流經下彎管段時，在流動慣性的作用下，流體旋轉的角度略小于傳熱管。

2) IHX一次側周向分布

外殼體與最外排換熱管間的流速明顯高于9排傳熱管間的區域，因流體的進、出口在此區域，且流動行程相對最短，形阻較其他區域小。

3) 入口腔室

在入口腔室的外殼體開有32個φ50 mm的開孔，每個開孔的流量基本相同，但流動方向由上至下存在較大差異。通過最上層開孔進入入口腔室的流體，流動方向存在一向上40°的仰角，并在外殼體和最外3層換熱管上部之間形成一漩渦。在漩渦區鈉溫較高,降低了換熱器入口腔室內的換熱效率。通過中部開孔的流體，流動方向基本與開孔垂直。通過下3排開孔的流體，進入入口腔室后流動方向改變90°,向下流出入口腔室進入換熱器的上直管區。在入口腔室區域，60%的流體的流動方向與換熱管垂直，即管間的橫向繞流是入口腔室的主要特點。

4) 出口腔室

IHX一次側出口流場示于圖7。流體進入出口腔室后繼續沿換熱管方向流動，最終在換熱器底部流動方向改變90°，大部分流體從最下兩排開孔流出IHX，這與傳統觀點完全不同。因IHX一次側各出口的流量分配由流體的體積力決定，因此，在布置IHX一次側出口流量和溫度測點時，應考慮出口腔室的流動特點。

圖7 IHX一次側出口流場

4 結論

本工作應用CFX程序計算獲得了IHX一次側速度場、換熱管束各部分的詳細溫度分布及變化規律。計算結果表明，出口鈉溫低于354 ℃，與俄羅斯提供的計算結果較一致，IHX可滿足設計要求，在CEFR滿功率運行時，IHX可將反應堆產生的熱量傳遞給二回路，保證主熱傳輸系統的正常工作。

參考文獻：

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