開式自然循環回路流動不穩定性及其抑制措施研究

石雪垚，王輝，黃政，李精精，鄭云濤，黃樹亮

中國核電工程有限公司，北京 100840

利用密度差作為驅動力的自然循環回路作為一種常見的非能動安全系統，在核電廠中常被用來進行事故后的非能動熱量導出?！叭A龍一號”作為自主研發的第三代核電技術，采用了開式自然循環方案的安全殼非能動熱量導出系統(passive containment heat removal system，PCS)[1-2]，即回路出口、入口分別與開放的水箱相連。該系統無需外部電源和人為干預，具有較高的可靠性，能夠在事故后非能動地將安全殼內的熱量導出到環境，防止安全殼升溫升壓導致的破壞，提高了事故后安全殼包容放射性物質的能力，使“華龍一號”滿足了國際第三代核電的安全標準[3-6]。

針對采用開式設計的自然循環回路，Guo 等[7]基于兩相均勻流模型開發了模擬開式自然循環回路的瞬態分析程序，并對回路的循環流量、壓降、溫度等瞬態熱工參數進行了分析研究。葛魁等[8]針對“華龍一號”非能動安全殼熱量導出系統，對自然循環回路的內部熱工水力特性進行了分析研究。王輝等[9]基于“華龍一號”非能動安全殼熱量導出系統綜合性能實驗裝置的實驗結果，采用漂移流模型開發了計算分析程序，并通過實驗結果對程序進行了驗證。哈爾濱工程大學科研團隊對低壓自然循環流動特性進行了研究，研究了在閃蒸段設置內插物、注入不可凝氣體以實現消除閃蒸流動不穩定性現象[10-12]。

自然循環回路主要由密度差產生的驅動壓頭克服流動阻力來驅動流體流動。在此過程中，系統流量與若干熱工參數形成耦合，會產生流動不穩定現象。為了進一步對開式自然循環回路流動過程中的不穩定現象進行研究，本文采用RELAP5程序對開式自然循環實驗回路進行了詳細的模擬，通過實驗結果對模型進行了驗證；在該模型的基礎上進一步對回路運行過程中的流動不穩定特點進行了研究，并針對自然循環回路流動不穩定的抑制措施進行了數值分析與驗證。

1 開式自然循環實驗回路模型

為了研究開式自然循環回路的流動不穩定性并對其抑制措施進行驗證，本研究采用RELAP5程序，建立了開式自然循環實驗回路的詳細模型。

1.1 計算模型

開式自然循環實驗回路如圖1 所示，整個回路由冷管段、換熱器、熱管段、汽水分離器以及若干閥門組成。

圖1 開式自然循環實驗裝置三維示意

回路中換熱器位于實驗臺架下部的冷凝罐中[9]，冷管段入口與汽水分離器位于上部的冷卻水箱中，為了防止冷卻水箱超壓，在水箱頂部設有蒸汽排放裝置。實驗過程中，通過向冷凝罐充入水蒸氣來模擬事故后安全殼內的壓力、溫度環境。

RELAP5 程序是高度通用的熱工水力分析程序，能夠模擬復雜的氣-液流動與傳熱現象[13]。為了更加精細地復現開式自然循環回路的流動特點，本研究采用RELAP5 程序對實驗回路進行了詳細模擬。

基于RELAP5 的開式自然循環回路模型的節點劃分如圖2 所示。模型中換熱水箱、冷管段、換熱器、熱管段及汽水分離器均采用PIPE 水力學構件進行模擬，換熱水箱外部環境采用SINGLVOL模擬，并與1 個與時間相關的控制體相連來保證空間壓力的穩定。

圖2 RELAP5 計算模型節點劃分

模型中冷卻水箱劃分了20 個節點；下降段共劃分了26 個節點；換熱器劃分了10 個節點；上升段與汽水分離器采用更加精細的節點劃分，其中上升段共劃分了63 個節點，汽水分離器劃分了20 個節點。通過在回路模型各個節點設置局部阻力系數來模擬系統回路中的閥門、彎頭等。

由于RELAP5 程序模擬冷凝罐中含不可凝氣體的冷凝換熱時具有較大的局限性，且本文重點關注回路內部自然循環流動的不穩定性，本模型在換熱器傳熱管外表面采用恒定熱流密度的邊界條件來模擬冷卻罐對回路換熱器的加熱。

1.2 模型驗證

表1給出了開式自然循環實驗回路2 個設計工況的實驗結果與RELAP5 模型計算結果。圖3 給出了系統運行在工況2 下，當系統出現流動不穩定現象時實驗得到的自然循環流量與模型計算結果比較。

表1 不同工況開式自然循環回路實驗與計算結果比較

圖3 工況2 自然循環流量

從表1 和圖3 的比較結果可以看出，本研究基于RELAP5 程序建立的開式自然循環回路分析模型，在2 個設計工況下能夠對回路的流動進行準確模擬，并且在工況2 中流量波動的形態與實驗結果符合較好，因此可以將該模型用于系統回路流動不穩定性研究。

2 開式自然循環流動不穩定性特點

開式自然回路在運行過程中，當換熱水箱溫度達到或接近飽和溫度時，冷卻水流經換熱器后，水溫會上升并超過換熱水箱內的飽和溫度；在上升管段中，隨著壓力的降低，在某一位置會出現閃蒸現象，在閃蒸點位置以上隨著高度的升高，壓力進一步降低，由于閃蒸造成的含汽率逐漸增大，密度逐漸下降，進一步驅動回路的循環流動。在這個過程中，當系統流量受到擾動而增大時，一方面造成循環阻力增加；但另一方面由于流量增大后，會導致換熱器出口水溫下降，熱管段水溫隨之下降，熱管段出現閃蒸的位置上升，造成回路自然循環驅動壓頭的下降，最終會減小系統流量。

從上述分析過程中可以看出，開式自然循環回路本身流動是否穩定取決于系統回路阻力隨流量的變化率和驅動壓頭隨流量的變化率。當系統回路阻力隨流量的變化率小于驅動壓頭隨流量的變化率時，系統流動處于不穩定狀態，反之則處于流動穩定狀態[14]。

目前，國內針對自然循環過程中流動不穩定性的實驗研究及理論分析主要針對閉式循環[14-15]，而本文研究的開式自然循環實驗回路在運行過程中不但存在冷卻水箱—下降段—換熱器—上升段—汽水分離器—冷卻水箱這一自然循環路徑，在換熱水箱內部，還存在冷卻水箱與汽水分離器內部之間的自然循環路徑，因此流動情況更為復雜。

綜上所述，開式自然循環回路流動受到多個因素的共同影響，由于實驗回路汽水分離器的存在，使得自然循環特性更加復雜，難以用理論推導的方法對其流動不穩定性特性進行研究。因此本文采用經過驗證的回路分析模型，對不同流量和功率下的回路流動狀態進行數值模擬，獲得實驗回路流動不穩定性的特點。

圖4給出了在不同功率狀態下的自然循環流量，圖5 給出了不同功率狀態下汽水分離器上部（圖2 中A 點）、下部（圖2 中B 點）及距離汽水分離器1 m 位置處（圖2 中C 點）的含汽率。

圖4 回路流量隨功率變化

圖5 汽水分離器及上升段含汽率隨回路功率變化

從圖4 和圖5 可以看出，當回路運行狀態功率在0～0.2 MW 時，在熱管段和汽水分離器內均未發生閃蒸，回路內流體呈單相狀態，并且流動平穩；當運行狀態在0.2～0.6 MW 時，在汽水分離器內發生閃蒸現象，導致汽水分離器出口含汽率出現周期性波動（見圖5 中A 點曲線），同時，回路內出現了小幅度的流量波動（見圖4）；當運行狀態在0.6～1.2 MW 時，汽水分離器底部和上升段上部開始出現閃蒸現象，含汽率呈周期性波動（見圖5 中B、C 點），造成回路內流量出現較大幅度的波動（見圖4）；當運行狀態在1.2 MW 功率以上時，出現閃蒸的位置逐漸穩定，上升段B、C 點含汽率和回路流量不再出現波動現象。

形成上述規律的主要原因是回路內的自然循環由溫差導致的密度差來驅動，并且開式自然循環回路具有一定的自適應特性。隨著排熱功率的上升，回路冷、熱管段的溫差逐漸增大，分別超過了A、B、C 點的飽和壓力，在局部發生閃蒸。而在B 點剛發生閃蒸時，含氣率的變化率對流量的影響更加明顯，導致在0.6～1.2 MW 時流動不穩定性較為明顯。

3 開式自然循環回路流動不穩定抑制

自然循環管道內的流動不穩定性可能導致系統管道震動、水錘等現象，雖然實驗過程中未發現可能影響回路安全運行的現象，但為了確?；芈烽L期運行的穩定性，仍需要對流動不穩定的抑制措施進行研究。

從前文針對流動不穩定性特點的數值分析可以看出，在一定功率水平下，上升管道內冷卻劑閃蒸點位置、含汽率的變化造成了驅動壓頭和流動阻力的周期變化，最終造成流動不穩定性現象。

本研究通過增加上升管段含汽率的方式以提高流動的穩定性，在設計上通過向回路上升管道注入不可凝氣體的方式實現，該措施在工程上易于實施，不會對回路整體功能造成影響。

為了驗證該措施對抑制流動不穩定的作用，分析中分別假設在汽水分離器底部B 點、上升管道距離冷卻水箱底部1 m 位置處C 點和上升管道距離冷卻水箱底部2 m 位置處D 點（見圖2）以恒定速率注入氮氣，對設計工況2 的流動狀態進行計算。

圖6分別給出了在不同位置以不同速率注入氮氣時在工況2 下的流動情況。表2 給出了不同注入位置和注入流量下回路自然循環流量。

表2 不同氮氣注入情況下的系統循環流量kg/s

圖6 不同注入點注入氮氣時工況2 自然循環流量

從計算結果可以看出，當在汽水分離器底部注入氮氣時（見圖6(a)），需要至少5 g/s 的注入速率才能對系統的流量波動產生良好的抑制效果。隨著注入位置的降低，所需要的氮氣注入速率逐漸下降。在距離冷卻水箱2 m 位置處的D 點，只需要1 g/s 的注入速率即可滿足對流量波動的抑制需求（見圖6(c)）。因此，通過在較低位置注入不可凝氣體，能夠提高系統回路的驅動壓頭以及驅動壓頭的穩定性，從而抑制系統流量的波動。

此外，注入不可凝氣體在抑制流量波動的同時還提高了系統自然循環的速率，有利于提高回路的冷卻能力，為核電站提供更加充分的安全裕量。

4 結論

本文對開式自然循環回路流動不穩定現象進行了進行了詳細的數值模擬。

1）通過對不同工作狀態下的回路流動不穩定現象進行數值計算，對回路流動不穩定性的特點及其形成原因進行了分析。

2）對采用局部注入不可凝氣體來抑制回路流動不穩定現象的有效性進行了論證。結果表明，在回路上升段距離冷卻水箱底部一定范圍內，以恒定流量注入不可凝氣體，能夠對自然循環回路運行過程中的流動不穩定現象有較好的抑制作用，在距離冷卻水箱底部2 m 位置處的上升段，僅需1 g/s 的不可凝氣體注入速率便能夠產生較好的抑制效果。同時，該措施還可以提高自然循環回路的循環流量，從而提高了系統的帶熱能力。

本研究結果可以為同樣采用開式自然循環設計的“華龍一號”PCS 系統優化和性能提升提供依據。