(中國核動力研究設(shè)計院 反應(yīng)堆系統(tǒng)設(shè)計技術(shù)重點實驗室,四川 成都 610041)

AP1000核電廠的穩(wěn)壓器頂部安裝兩臺彈簧式安全閥,提供反應(yīng)堆冷卻劑系統(tǒng)的超壓保護,防止系統(tǒng)壓力超過設(shè)計值,確保反應(yīng)堆冷卻劑壓力邊界的完整性。

穩(wěn)壓器安全閥的排放能力在某些事故工況下直接影響核電廠的安全性。安全閥的動作時間、整定值壓力、回座壓力等參數(shù),目前均由試驗進行測量并整定,但安全閥的排放能力及排放特性,一般是基于系統(tǒng)運行和工程設(shè)計經(jīng)驗確定,但對安裝在電廠中的回路是否適用需進行詳細的論證和分析。同時,AP1000核電廠的穩(wěn)壓器安全閥口徑非常大,一般工業(yè)試驗臺架不能滿足此安全閥的試驗條件,而CFD分析則很好地解決了這一難題。

AP1000穩(wěn)壓器安全閥在排放的介質(zhì)為水蒸氣,在排放過程,由于壓力的變化,水蒸氣部分冷凝為液態(tài)水,使得排放的介質(zhì)從單相的水蒸氣變化復雜多變的兩相流。

國內(nèi)外關(guān)于管閥CFD三維精確仿真分析的應(yīng)用已經(jīng)非常普遍也比較成熟,但大都還局限于管閥的定常、單相流仿真分析,關(guān)于管閥非定常CFD動態(tài)效應(yīng)和多相流效應(yīng)的研究還較少。

1 基本控制方程

穩(wěn)壓器安全閥的排放介質(zhì)初始為過熱蒸汽,隨著時間的推進,系統(tǒng)壓力的不斷變化,以及喉部流速的動蕩,導致后期排放過程中,介質(zhì)由單相的蒸汽變?yōu)轱柡驼羝?、汽液兩相等。排放的介質(zhì)與系統(tǒng)初始壓力和壓力變化有關(guān),本文假設(shè)排放介質(zhì)為過熱蒸汽。

模型的基本控制方程為：

連續(xù)性方程：

動量方程：

能量方程：

采用的求解方法為有限體積法,湍流模擬采用標準的k-ε兩方程模型。

2 動網(wǎng)格及其二次開發(fā)

2.1 靜態(tài)網(wǎng)格

首先對安全閥關(guān)閉時的靜態(tài)網(wǎng)格進行劃分。設(shè)置幾何模型的進出口邊界,并利用ICEM對幾何模型采用結(jié)構(gòu)化網(wǎng)格劃分方式,對閥瓣和密封面處的網(wǎng)格進行局部加密處理,以得到排放過程中閥瓣處更加精確的動網(wǎng)格。安全閥的靜態(tài)網(wǎng)格劃分后,共得到321 418個節(jié)點和95 723個面,如圖1所示。圖2為縱截面的網(wǎng)格局部加密圖。

圖1 靜態(tài)網(wǎng)格的劃分Fig.1 Static mesh division

圖2 縱截面網(wǎng)格加密Fig.2 Mesh division on the cross section

2.2 動網(wǎng)格及其二次開發(fā)

安全閥在開啟過程中,閥瓣處的截面在不斷變化,該處的網(wǎng)格需要不斷重繪,因此,需要設(shè)置動網(wǎng)格。

本文的動網(wǎng)格技術(shù)通過以下方式實現(xiàn)：

1)當安全閥開啟到1%,2%,3%,4%,5%,6%,7%,8%,9%,10%,20%,30%,40%,50%,60%,70%,80%,90%,100%狀態(tài)時,計算自動中斷;

2)進行動網(wǎng)格的二次開發(fā),編寫外部命令,提取前一狀態(tài)的所有數(shù)據(jù)點;

3)對數(shù)據(jù)進行分析,從事先準備的網(wǎng)格中選取合適的網(wǎng)格替換;

4)求解器加載定位置處的網(wǎng)格,繼續(xù)進行計算。

3 計算工況、邊界條件

模擬全廠失電工況(SBO)下,穩(wěn)壓器安全閥開啟這一最惡劣的設(shè)計擴展工況。

穩(wěn)壓器安全閥入口的壓力隨時間的變化示意圖如圖3所示,通過編寫自定義程序給出。

圖3 入口壓力隨時間變化的示意圖Fig.3 Inlet pressure variation versus time

設(shè)定安全閥排放的初始入口溫度為343 ℃,初始流速為0 m/s,初始含汽率為1,安全閥在入口壓力達到17.23 MPa時開啟。

4 計算結(jié)果分析

4.1 閥門內(nèi)部壓力分布和速度分布

圖4 0.002 s時刻閥門內(nèi)部壓力場合速度場Fig.4 Pressure and speed feed at 0.002 s

圖5 0.026 s時刻閥門內(nèi)部壓力場合速度場Fig.5 Pressure and speed feed at 0.026 s

圖6 0.056 s時刻閥門內(nèi)部壓力場合速度場Fig.6 Pressure and speed feed at 0.056 s

圖7 0.096 s時刻閥門內(nèi)部壓力場合速度場Fig.7 Pressure and speed feed at 0.096 s

圖4、圖5、圖6和圖7給出了0.002 s、0.026 s、0.056 s和0.096 s的閥門內(nèi)部壓力分布和速度分布?？梢钥闯?當安全閥入口壓力達到開啟整定壓力。0.002 s時,閥門開啟角度較小,從速度分布來看,存在明顯的射流現(xiàn)象,即向四周急速的噴射,射流角度較小,隨著閥門開度不斷增大,射流角度逐漸增大,對閥座的沖擊力也逐漸增大;閥門的壓力在閥門開啟后迅速降低,而入口壓力隨著閥門排量的限制和蓄壓的原因,先減小,后增大,在0.096 s時,由于給定的入口壓力增大,導致閥門入口壓力增大。

4.2 閥位變化仿真結(jié)果

閥門隨時間的開啟高度的變化如圖8所示,由于沒有機械的阻尼機制,所以閥門動能的消耗只通過湍動能的耗散變成熱能。故其在跳起后呈現(xiàn)了出了圖8中所示的,先跳起一個比較大的高度,然后又回彈的運動過程。

圖8 閥門位移隨時間的變化Fig.8 Disposition change versus time

從圖8可以看出,閥門從開始開啟到完全開啟,由于閥門入口壓力大,閥門呈線性開啟,隨后,由于閥門入口壓力下降,閥門的逐漸關(guān)閉,并在之后出現(xiàn)了閥門的開啟和關(guān)閉的交替變化,說明AP1000彈簧式安全閥的入口壓力較小時,會出現(xiàn)頻跳現(xiàn)象,容易造成閥瓣的損壞。

4.3 閥瓣受力變化仿真結(jié)果

閥門受力主要為氣體對閥門的氣動力和彈簧給予閥門的力。圖9所示為這兩個力隨時間的變化。

圖9 閥門受力隨時間的變化Fig.9 Forces on the valve versus time

閥門開啟之初,氣動力和彈簧力線性升高至最大值,閥門線性全開。隨后,由于閥門的頻跳現(xiàn)象,導致氣動力和彈簧力出現(xiàn)了一定程度的動蕩。

5 結(jié) 論

本文在AP1000穩(wěn)壓器安全閥內(nèi)部的靜態(tài)CFX網(wǎng)格的基礎(chǔ)上,進行了動網(wǎng)格的二次開發(fā),通過編寫外部命令,提取不同閥門開度時的靜態(tài)網(wǎng)格,實現(xiàn)了對AP1000穩(wěn)壓器安全閥動態(tài)排放時閥門各個截面處的網(wǎng)格的模擬,通過CFD計算防震,得到了AP1000穩(wěn)壓器安全閥在排放過程中的流動特性,一定程度上復現(xiàn)了安全閥開啟過程中壓力、流速、開啟高度、閥瓣受力的動態(tài)變化特性,可以為安全閥的性能評估以及安全閥及其下游管道的載荷計算提供一定的參考輸入。