赫榮輝　蘇應(yīng)斌

【摘 要】壓水堆核電廠乏燃料水池靠乏燃料水池冷卻系統(tǒng)冷卻。在長(zhǎng)期運(yùn)行以后，池水中含有一定量的雜質(zhì)顆粒，水池水質(zhì)的變化對(duì)乏燃料水池冷卻系統(tǒng)產(chǎn)生一定的影響。本文利用CFD對(duì)乏燃料水池冷卻系統(tǒng)的蝶閥及相連接的管路受固體雜質(zhì)顆粒沖蝕問(wèn)題進(jìn)行了模擬分析，計(jì)算得出了在不同流速下固體顆粒對(duì)閥門及管道的沖蝕量，通過(guò)對(duì)模擬數(shù)據(jù)的處理確定了受侵蝕嚴(yán)重的部位，計(jì)算表明流速是導(dǎo)致沖蝕現(xiàn)象的主要原因。

【關(guān)鍵詞】CFD；沖蝕；流速

0 引言

壓水堆核電廠乏燃料水池靠乏燃料水池冷卻系統(tǒng)冷卻，在長(zhǎng)期運(yùn)行之后水池中會(huì)聚集一定量的固體雜質(zhì)顆粒，在乏燃料水池冷卻系統(tǒng)啟動(dòng)后，雜質(zhì)顆粒會(huì)隨著冷卻劑進(jìn)入系統(tǒng)管道中，在閥門的局部空間，由于截流作用流固混合物會(huì)對(duì)閥體結(jié)構(gòu)會(huì)產(chǎn)生一定程度的沖蝕從而影響閥門部件的壽命。在乏燃料水池冷卻系統(tǒng)設(shè)計(jì)時(shí)，系統(tǒng)流速的選擇將影響著閥門的使用壽命，因此在確定系統(tǒng)參數(shù)時(shí)有必要對(duì)潛在發(fā)生固體顆粒沖蝕的部位進(jìn)行分析。

本文利用CFD計(jì)算流體力學(xué)軟件[1]對(duì)PWR乏燃料水池冷卻系統(tǒng)帶有蝶閥的圓管內(nèi)固體顆粒沖蝕問(wèn)題進(jìn)行了數(shù)值模擬，計(jì)算了在不同流速下固體顆粒物對(duì)管道及閥門的沖蝕量，確定了受侵蝕嚴(yán)重的部位，為管路系統(tǒng)的優(yōu)化設(shè)計(jì)提供參考。

1 離散相模型

乏燃料水池冷卻系統(tǒng)部分工質(zhì)是固液混合物，從微觀角度看，由于雜質(zhì)顆粒間的間距很大，顆粒間的相互碰撞可以忽略，從而滿足離子負(fù)載流動(dòng)模式，本文采用歐拉-拉格朗日追蹤法方法的離散相模型。

2 計(jì)算模型及邊界條件

計(jì)算模型是一個(gè)帶有蝶閥的管路，管道內(nèi)徑為D=108mm，長(zhǎng)度為L(zhǎng)=1000mm，蝶閥位于管道的中部位置。由于模型的對(duì)稱特殊性，在計(jì)算中對(duì)模型進(jìn)行了簡(jiǎn)化，沿模型的一半作為計(jì)算模型，計(jì)算模型如圖1所示。

管道入口為速度入口邊界，固液兩相均以恒定的初始速度進(jìn)入流場(chǎng)，湍流強(qiáng)度為5%。

假定排放管路三個(gè)初始速度分別為w=6m/s、w=3m/s、w=1m/s。

計(jì)算中應(yīng)用以下經(jīng)驗(yàn)公式各質(zhì)點(diǎn)的運(yùn)動(dòng)速度：

U=Wmax（1-r/Rmax）1/7

式中，Wmax為進(jìn)口流速最大值，Rmax為管道內(nèi)徑，r為隨機(jī)質(zhì)點(diǎn)處對(duì)應(yīng)的幾何當(dāng)量半徑。

計(jì)算時(shí)固體顆粒的質(zhì)量流量為0.01kg/s，直徑為0.1mm，管道壁面等值粗糙度取0.1mm。

管道出口設(shè)為壓力出口，管道壁面和閥門壁面為固體壁面，顆粒與壁面的的邊界為彈性反射。

采用Simple算法對(duì)流場(chǎng)進(jìn)行數(shù)值計(jì)算，湍流模型為標(biāo)準(zhǔn)k-ε模型，動(dòng)量方程采用一階迎風(fēng)格式離散，采用標(biāo)準(zhǔn)壁面函數(shù)，壁面無(wú)滑移。

計(jì)算時(shí)，假定了最大的跟蹤距離是管道的總長(zhǎng)度L，為了觀察到速度慢的粒子運(yùn)動(dòng)軌跡，最大跟蹤時(shí)間設(shè)為10L/υ。

網(wǎng)格采用四面體網(wǎng)格，網(wǎng)格數(shù)量為50萬(wàn)。

3 數(shù)值仿真與結(jié)果分析

3.1 固體相運(yùn)動(dòng)

從流線圖（圖3）中可以看出，流體進(jìn)入管道內(nèi)在閥門前的流動(dòng)很平穩(wěn)，固體顆粒做直線運(yùn)動(dòng)，經(jīng)過(guò)閥門后在慣性力的作用下向閥板外側(cè)的壁面聚集，這一過(guò)程中固體顆粒以不同的角度與壁面發(fā)生碰撞，角度為0°～45°。固體顆粒與壁面碰撞后運(yùn)動(dòng)方向發(fā)生劇烈變化，從粒子軌跡動(dòng)圖（圖4）中可以看出固體顆粒經(jīng)過(guò)閥門時(shí)與管壁以及閥門發(fā)生碰撞，碰撞最激烈的部位是閥門出口的閥板邊緣及與閥板相垂直的管壁上。閥門出口的閥板邊緣碰撞激烈的原因是由于閥門的截流作用引起的，閥板相垂直的管壁上碰撞激烈的原因是由于鏡面反射及固體顆粒流經(jīng)所受慣性力的作用。

3.2 流動(dòng)沖蝕問(wèn)題分析

流速為6m/s時(shí)蝶閥及管道受沖蝕計(jì)算結(jié)果見圖5-圖6，從圖中可以看出，在閥門入口前未發(fā)生沖蝕現(xiàn)象，在閥門出口的管壁上以及閥板上均發(fā)生了較為明顯的沖蝕，沖蝕率密度最大值為25kg/m2·s，從云圖中可以看出在蝶閥與管壁相近處均發(fā)生了較嚴(yán)重的沖蝕，沖蝕的最大點(diǎn)發(fā)生在與閥板相垂直的管壁上。由3.1節(jié)的分析可知，碰撞發(fā)生最激烈的部位是閥門出口的閥板邊緣與閥板相垂直的管壁上，由沖蝕率密度云圖上可以看出正是這些部位的沖蝕現(xiàn)象最嚴(yán)重，僅個(gè)別的閥板部位及閥門出口管道壁面未發(fā)生沖蝕現(xiàn)象。

為了驗(yàn)證流速對(duì)沖蝕的影響，將固液流速調(diào)整為3m/s和1m/s，重新進(jìn)行了計(jì)算，計(jì)算結(jié)果見圖7、圖8、圖9、圖10。

從圖中可以看出流速減慢后，沖蝕現(xiàn)象明顯降低，但沖蝕嚴(yán)重的部位依舊是閥門出口的閥板邊緣及與閥板相垂直的管壁上。流速為3m/s時(shí)，沖蝕率密度最大值為25kg/m2·s，流速為1m/s時(shí)，沖蝕率密度最大值僅為10kg/m2·s。當(dāng)速度為3m/s時(shí)，閥板及部分出口管道發(fā)生了沖蝕現(xiàn)象，當(dāng)速度降為1.0m/s時(shí)，沖蝕現(xiàn)象幾乎很少發(fā)生，但閥板邊緣的沖蝕現(xiàn)象仍舊存在。

經(jīng)過(guò)本文的分析可知，流速是導(dǎo)致沖蝕現(xiàn)象的最大原因之一，在帶有固液兩相流動(dòng)的系統(tǒng)設(shè)計(jì)時(shí)，一定需要避免高流速，建議工藝系統(tǒng)流速應(yīng)控制在1m/s以內(nèi)。

4 結(jié)論

（1）本文分析的帶有閥門的管路，在管內(nèi)為固液兩相流動(dòng)時(shí)，固體粒子與邊界碰撞會(huì)發(fā)生碰撞，碰撞最激烈的部位是閥門出口的閥板邊緣與閥板相垂直的管壁上。

（2）當(dāng)流速為6m/s時(shí)，閥門出口的管道及閥板上沖蝕現(xiàn)象很明顯，當(dāng)速度降低時(shí)，沖蝕現(xiàn)象會(huì)隨之降低。建議系統(tǒng)流速應(yīng)控制在1m/s以內(nèi)。

【參考文獻(xiàn)】

[1]ANSYS 12.1Help[Z].

[2]泥沙顆粒在沖擊式水輪機(jī)斗葉內(nèi)壁面的沖蝕磨損研究[J].曹勇，等.水力發(fā)電 2014（7）.

[責(zé)任編輯：楊玉潔]