賴喜德，程海，葉道星，杜江，陳小明

(1. 西華大學(xué)能源與動(dòng)力工程學(xué)院，四川 成都 610039； 2. 東方法馬通核泵有限責(zé)任公司，四川 德陽(yáng) 618000)

核主泵(反應(yīng)堆冷卻劑循環(huán)泵)作為核島內(nèi)控制冷卻劑循環(huán)的核一級(jí)關(guān)鍵設(shè)備，其性能直接決定了反應(yīng)堆冷卻劑系統(tǒng)運(yùn)行的穩(wěn)定性和可靠性[1].在壓水堆(PWR)核電廠的反應(yīng)堆一回路系統(tǒng)中，核主泵處于并聯(lián)運(yùn)行狀態(tài)，啟動(dòng)過(guò)程中就會(huì)產(chǎn)生反向流動(dòng)現(xiàn)象，在事故工況下，核主泵會(huì)存在多種不同的非正常運(yùn)行工況[2].為了滿足核電站系統(tǒng)設(shè)計(jì)需求，核主泵制造商必須提供主泵的運(yùn)行全特性(包括正常和非正常運(yùn)行工況)曲線，即四象限特性曲線.目前核主泵的四象限特性曲線主要是通過(guò)試驗(yàn)獲得[3-4].然而，搭建核主泵四象限特性的試驗(yàn)臺(tái)非常昂貴，要完成一個(gè)核主泵的全特性試驗(yàn)，不僅耗時(shí)，而且成本很高.近年來(lái)，采用CFD技術(shù)對(duì)泵在第一象限運(yùn)行，特別在設(shè)計(jì)工況點(diǎn)附近的數(shù)值模擬預(yù)測(cè)結(jié)果與試驗(yàn)結(jié)果之間的偏差已非常小，其數(shù)值計(jì)算方法已被證明準(zhǔn)確、可靠[5-6].GROS等[7]和COUZINET等[8]分別通過(guò)數(shù)值方法探索離心泵的四象限性能以及正轉(zhuǎn)反轉(zhuǎn)工況下的徑向力變化.楊從新等[9]對(duì)專用堆型對(duì)應(yīng)的核主泵正反轉(zhuǎn)流動(dòng)進(jìn)行了數(shù)值模擬.數(shù)值模擬方法可以預(yù)測(cè)核主泵的外特性，并掌握其內(nèi)部流場(chǎng)結(jié)構(gòu).由于核主泵多為導(dǎo)葉式混流泵，在偏離最優(yōu)工況較遠(yuǎn)的區(qū)域或非正常運(yùn)轉(zhuǎn)區(qū)域運(yùn)行時(shí)，其內(nèi)部流動(dòng)較為復(fù)雜，要準(zhǔn)確地預(yù)測(cè)其性能仍面臨很多挑戰(zhàn)[10-11].

文中采用全流道內(nèi)流場(chǎng)數(shù)值模擬方法，針對(duì)某軸封式核主泵的四象限運(yùn)行特性進(jìn)行預(yù)測(cè)，探索基于數(shù)值模擬預(yù)測(cè)四象限特性曲線的可靠性，以實(shí)現(xiàn)在核主泵產(chǎn)品開發(fā)過(guò)程中對(duì)四象限特性進(jìn)行預(yù)測(cè)與設(shè)計(jì)方案快速可靠評(píng)估，為基于全特性要求的主泵優(yōu)化設(shè)計(jì)提供一定參考.

1 主泵運(yùn)行方式與四象限特性曲線

1.1 核主泵運(yùn)行方式

多數(shù)壓水堆核電廠反應(yīng)堆一回路系統(tǒng)(見(jiàn)圖1)中的核主泵(如美國(guó)西屋和法國(guó)法馬通的裝機(jī)900 MW的單堆都采用3臺(tái)主泵)[1]都處于并聯(lián)運(yùn)行狀態(tài)，因此，在評(píng)價(jià)核主泵的適用性和安全性時(shí)應(yīng)對(duì)整個(gè)一回路系統(tǒng)進(jìn)行分析.在壓水堆系統(tǒng)中，啟動(dòng)過(guò)程會(huì)產(chǎn)生反向流動(dòng)現(xiàn)象，在沒(méi)有任何預(yù)防反向流動(dòng)措施時(shí)，如果出現(xiàn)一個(gè)反方向的壓頭，所需的轉(zhuǎn)矩就有可能超過(guò)正常轉(zhuǎn)矩，容易導(dǎo)致電動(dòng)機(jī)損壞或失效.因此，在壓水堆一回路系統(tǒng)中的主泵上通常會(huì)設(shè)置防反轉(zhuǎn)裝置.在事故工況下，核主泵可能處于不同的狀態(tài).例如，當(dāng)防反轉(zhuǎn)裝置失效或者在大泄漏處有大量的冷卻劑流動(dòng)而斷裂失效的情況下，主泵將進(jìn)入耗能工況并進(jìn)一步變化到水輪機(jī)工況，同時(shí)伴隨旋轉(zhuǎn)方向和流動(dòng)方向的改變.與在抽水蓄能電站的水泵水輪機(jī)類似，按照運(yùn)行的轉(zhuǎn)速、流量、揚(yáng)程的正反方向，核主泵可能在8個(gè)工況區(qū)運(yùn)行，但在分析主冷卻劑系統(tǒng)設(shè)計(jì)、運(yùn)行和安全性時(shí)，主要需要考慮主泵的啟動(dòng)過(guò)程、工況過(guò)渡過(guò)程、事故斷電、惰轉(zhuǎn)特性及反向流動(dòng)特性等.

圖1 壓水堆中一回路主要設(shè)備布置

1.2 核主泵四象限特性的類比曲線

對(duì)于一般用途的泵，只需給出在正轉(zhuǎn)正流運(yùn)轉(zhuǎn)條件下的第一象限特性曲線，而核主泵必須考慮在非正常運(yùn)行條件下可能發(fā)生的包括逆流、反轉(zhuǎn)等在內(nèi)的所有工況下的揚(yáng)程與流量、轉(zhuǎn)矩與流量等之間關(guān)系，即要獲得泵的全特性.因要描述核主泵的全特性將涉及直角坐標(biāo)的4個(gè)象限，故被稱為四象限特性曲線.四象限特性曲線一般由泵出廠試驗(yàn)得到，用于掌握泵的運(yùn)行特性，通常所指的四象限特性試驗(yàn)是在泵的額定轉(zhuǎn)速下進(jìn)行的.為便于反應(yīng)堆主冷卻系統(tǒng)設(shè)計(jì)和安全校核計(jì)算，行業(yè)中要求采用四象限特性曲線的類比(同源)曲線表達(dá)方式.四象限特性試驗(yàn)結(jié)果經(jīng)過(guò)量綱一化處理后，便可以得到泵的類比曲線.

在對(duì)核主泵的四象限特性試驗(yàn)過(guò)程中通過(guò)改變流量大小和葉輪轉(zhuǎn)動(dòng)方向，測(cè)得對(duì)應(yīng)揚(yáng)程與轉(zhuǎn)矩等泵的全特性試驗(yàn)數(shù)據(jù).首先以任意工況下的泵的運(yùn)行參數(shù)(H,Q,M,n)與主泵額定工況點(diǎn)的參數(shù)(Hd,Qd,Md,nd)之比得到量綱一的揚(yáng)程h、流量q、轉(zhuǎn)矩m和轉(zhuǎn)速N，即

h=H/Hd,q=Q/Qd,m=M/Md,N=n/nd.

(1)

圖2 核主泵四象限特性類比曲線圖

2 四象限特性曲線預(yù)測(cè)方法

2.1 研究對(duì)象

以用于改進(jìn)的壓水堆(見(jiàn)圖1)的某型軸封式核主泵為研究對(duì)象，該核主泵(RCP100)的水力單元如圖3所示，其主要設(shè)計(jì)性能參數(shù)分別為流量Qd=23 790 m3/h，揚(yáng)程Hd=98 m，轉(zhuǎn)速nd=1 485 r/min，進(jìn)口壓力pin=15.16 MPa，工作溫度T=293 ℃；主要幾何參數(shù)分別為葉輪出口最大直徑D2max=800 mm，葉輪葉片數(shù)Zb=7，導(dǎo)葉數(shù)Zg=12，葉輪葉片和空間導(dǎo)葉均為空間雕塑曲面體.冷卻劑介質(zhì)物性參數(shù)分別為冷態(tài)密度1 000 kg/m3，冷態(tài)運(yùn)動(dòng)黏度8.93×10-7m2/s.熱態(tài)密度742 kg/m3，熱態(tài)運(yùn)動(dòng)黏度9.42×10-8m2/s.全流道包括直管進(jìn)水流道、混流式葉輪、空間導(dǎo)葉、類球形壓水室及出口管道部分.

圖3 核主泵水力單元流道三維模型

2.2 三維流動(dòng)數(shù)值模擬方法

2.2.1 湍流模型及求解方法

在實(shí)際運(yùn)行中，核主泵內(nèi)部為三維非定常流動(dòng)，采用有限體積法進(jìn)行求解，其控制方程[12-13]的通用形式為

(2)

式中：ρ為密度；Г為廣義擴(kuò)散系數(shù)；S為廣義源項(xiàng);φ為描述流場(chǎng)的通用變量，可以表示p，ui，T等物理量.φ不同，對(duì)應(yīng)于連續(xù)方程、動(dòng)量守恒方程、能量方程等流體狀態(tài)方程.

對(duì)式(2)進(jìn)行雷諾時(shí)均化，并采用SSTk-ω湍流模型封閉方程組[14].SSTk-ω模型的雙方程為

(3)

(4)

式中：Pk為湍流脈動(dòng)動(dòng)能k的生成項(xiàng)；Pω為湍流脈動(dòng)頻率ω的生成項(xiàng)；Γk，Γω分別為k和ω的有效擴(kuò)散系數(shù)；Yk，Yω分別為k和ω的耗散項(xiàng)；Dω為正交擴(kuò)散項(xiàng).

SSTk-ω模型是一個(gè)自適應(yīng)湍流模型，它在完全湍流區(qū)運(yùn)用k-ε模型進(jìn)行求解，在近壁區(qū)運(yùn)用k-ω模型求解.對(duì)水力單元過(guò)流部件流道幾何離散化，采用ANSYS-CFX?求解.流體在壁面處的流動(dòng)模式采用無(wú)滑移邊界條件，在CFX?中選用基于k-ω方程的自動(dòng)壁面處理模型[9].該模型根據(jù)網(wǎng)格密度的大小自動(dòng)判斷和選擇處理邊界層內(nèi)的流動(dòng)，不需要嚴(yán)格控制y+<20，其y+值只提供近壁分辨率信息.

2.2.2 全流道離散

對(duì)如圖4所示的全部流體域網(wǎng)格用ANSYS/ICEM CFD?進(jìn)行離散，其中進(jìn)口段、葉輪和空間導(dǎo)葉流道內(nèi)的網(wǎng)格為六面體結(jié)構(gòu)化網(wǎng)格，殼體和出口延伸段內(nèi)網(wǎng)格為混合網(wǎng)格.由于整個(gè)流體域中有空間三維扭曲的葉輪和空間導(dǎo)葉的葉片，網(wǎng)格很難都滿足y+<20，考慮到文中研究的重點(diǎn)是外特性預(yù)測(cè)，而不是邊界層內(nèi)的流動(dòng)規(guī)律，在CFX?中一般要求y+≤90即可滿足非邊界層內(nèi)流動(dòng)的計(jì)算要求[9].對(duì)進(jìn)口段、葉輪、空間導(dǎo)葉流道內(nèi)的近壁區(qū)網(wǎng)格進(jìn)行特別的加密處理，使過(guò)流表面絕大部分區(qū)域的y+值控制在12～20，局部在21～65，前后蓋板與葉片的過(guò)渡局部區(qū)域y+<12，導(dǎo)葉區(qū)域絕大部分區(qū)域的y+值控制在12～20，局部在20～30，能滿足基于k-ω方程的自動(dòng)壁面處理模型的要求.

圖4 核主泵的全流道流體域模型

2.2.3 邊界條件設(shè)置

全流道流場(chǎng)數(shù)值模擬為所有過(guò)流部件耦合計(jì)算，葉輪與空間導(dǎo)葉之間的交界面采用STAGE模型級(jí)間連接.對(duì)流項(xiàng)和湍流項(xiàng)數(shù)值精度均為2階.計(jì)算時(shí)按設(shè)定工況以質(zhì)量流量作為進(jìn)口邊界條件.

2.2.4 網(wǎng)格無(wú)關(guān)性驗(yàn)證

為了驗(yàn)證網(wǎng)格精度對(duì)數(shù)值計(jì)算結(jié)果的影響，采用基于流場(chǎng)數(shù)值模擬預(yù)測(cè)出的外特性與該核主泵在冷態(tài)試驗(yàn)曲線[9]相比較，進(jìn)行網(wǎng)格無(wú)關(guān)性驗(yàn)證，如圖5所示.

圖5 水力單元預(yù)測(cè)結(jié)果與試驗(yàn)外特性對(duì)比

由圖5可以看出：當(dāng)全流道網(wǎng)格單元總數(shù)達(dá)到280萬(wàn)后，在設(shè)計(jì)工況下核主泵揚(yáng)程的相對(duì)偏差為±0.60%，效率的相對(duì)偏差為±0.85%；在0.2Qd～1.2Qd流量范圍內(nèi)，揚(yáng)程的最大相對(duì)偏差為1.38%，效率的最大相對(duì)偏差為1.41%，計(jì)算結(jié)果滿足工程實(shí)際的需要.最終確定數(shù)值計(jì)算的網(wǎng)格單元數(shù)總數(shù)約為285萬(wàn).

2.3 外特性計(jì)算與四象限特性預(yù)測(cè)方法

2.3.1 外特性計(jì)算方法

在熱態(tài)工況下對(duì)不同流量的核主泵內(nèi)流場(chǎng)進(jìn)行三維非定常流動(dòng)數(shù)值模擬，得到全流場(chǎng)的速度矢量分布和壓力分布.提取相關(guān)數(shù)據(jù)[8]，核主泵外特性計(jì)算公式為

(5)

(6)

式中：pin，pout分別為泵進(jìn)口、出口質(zhì)量加權(quán)平均總壓；M為壓力或黏性力對(duì)葉輪回轉(zhuǎn)軸的力矩；ω′為葉輪的角速度；ρ為流體密度；g為重力加速度.

2.3.2 四象限特性的預(yù)測(cè)計(jì)算工況確定

參考圖2所示的核主泵類比曲線運(yùn)行范圍，根據(jù)四象限特性曲線繪制數(shù)據(jù)要求，類似于試驗(yàn)過(guò)程的參數(shù)調(diào)節(jié)，按定轉(zhuǎn)速對(duì)核主泵水力單元進(jìn)行全流道定常流場(chǎng)數(shù)值模擬，即保持轉(zhuǎn)速大小不變(按區(qū)域分正轉(zhuǎn)和反轉(zhuǎn))，改變流量(從負(fù)流量到正流量).數(shù)值計(jì)算流量從±0.1Qd～±1.6Qd內(nèi)，每隔0.1Qd取1個(gè)計(jì)算工況，共計(jì)64個(gè)工況點(diǎn)進(jìn)行全流道數(shù)值計(jì)算，提取相關(guān)流場(chǎng)數(shù)據(jù)并計(jì)算出核主泵外特性.

2.3.3 外特性預(yù)測(cè)方法的驗(yàn)證

為了驗(yàn)證基于數(shù)值模擬預(yù)測(cè)方法的可靠性，首先對(duì)核主泵四象限特性曲線中的正轉(zhuǎn)正流工況進(jìn)行計(jì)算，采用上述方法用CFX?軟件對(duì)RCP100型混流式核主泵進(jìn)行不同流量泵工況下的數(shù)值模擬，計(jì)算包括0.1Qd～1.6Qd在內(nèi)共16個(gè)工況點(diǎn).對(duì)每個(gè)工況下的揚(yáng)程和轉(zhuǎn)矩性能曲線進(jìn)行計(jì)算，并與試驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比，如圖6所示.

圖6 核主泵正轉(zhuǎn)正流工況數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果對(duì)比

由圖6可以看出，數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果一致性較好，數(shù)值計(jì)算結(jié)果與試驗(yàn)測(cè)出揚(yáng)程和轉(zhuǎn)矩的偏差大部分點(diǎn)小于±3%，個(gè)別點(diǎn)在±5%以內(nèi)，這表明文中所采用的數(shù)值計(jì)算方法是可靠的.

3 核主泵的四象限特性曲線預(yù)測(cè)

采用在正轉(zhuǎn)正流工況流場(chǎng)數(shù)值模擬及外特性預(yù)測(cè)方法對(duì)核主泵其他3個(gè)象限工況(包括正轉(zhuǎn)逆流、反轉(zhuǎn)正流和反轉(zhuǎn)逆流工況)進(jìn)行數(shù)值計(jì)算，計(jì)算包括±0.1Qd～±1.6Qd在內(nèi)共48個(gè)工況點(diǎn).

按照上述的換算方法，首先將核主泵的類比曲線換算成常規(guī)四象限特性曲線，再將數(shù)值模擬預(yù)測(cè)計(jì)算的性能曲線換算成常規(guī)四象限特性曲線，并將數(shù)值計(jì)算與試驗(yàn)得到的常規(guī)四象限特性曲線進(jìn)行對(duì)比，如圖7所示.

圖7 數(shù)值預(yù)測(cè)結(jié)果與試驗(yàn)數(shù)據(jù)曲線對(duì)比

由圖7可以看出：核主泵的正轉(zhuǎn)逆流工況及反轉(zhuǎn)逆流工況預(yù)測(cè)結(jié)果與試驗(yàn)結(jié)果一致性好，各曲線變化趨勢(shì)相同且重合度高，大部分工況點(diǎn)的揚(yáng)程誤差和轉(zhuǎn)矩偏差在±3%以內(nèi)，個(gè)別點(diǎn)在±5%以內(nèi)；核主泵在反轉(zhuǎn)正流工況(第四象限)數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果的各曲線變化趨勢(shì)相同，但預(yù)測(cè)與試驗(yàn)曲線間的重合度不高，大部分工況點(diǎn)的揚(yáng)程和轉(zhuǎn)矩偏差在±10%以內(nèi)，個(gè)別點(diǎn)在±20%內(nèi).雖然此工況區(qū)采用數(shù)值模擬預(yù)測(cè)計(jì)算結(jié)果與試驗(yàn)相差較大，但性能曲線的變化趨勢(shì)是相同的，這說(shuō)明采用基于數(shù)值模擬的預(yù)測(cè)方法能夠替代大部分流量范圍的試驗(yàn)來(lái)獲得可靠的特性曲線，但對(duì)核主泵在反轉(zhuǎn)正流工況的準(zhǔn)確數(shù)值模擬方法還需要進(jìn)一步研究.

4 結(jié) 論

在定轉(zhuǎn)速工況下，對(duì)某軸封式核主泵四象限運(yùn)行特性進(jìn)行了全流道內(nèi)流場(chǎng)數(shù)值模擬并預(yù)測(cè)了其四象限特性曲線，將計(jì)算結(jié)果與核主泵的試驗(yàn)結(jié)果進(jìn)行對(duì)比，得到結(jié)論如下：

1) 在正轉(zhuǎn)逆流工況及反轉(zhuǎn)逆流工況，基于數(shù)值模擬的預(yù)測(cè)結(jié)果與試驗(yàn)結(jié)果一致性較好，大部分工況點(diǎn)的揚(yáng)程誤差和轉(zhuǎn)矩偏差在±3%以內(nèi)，個(gè)別點(diǎn)在±5%以內(nèi).在反轉(zhuǎn)正流工況預(yù)測(cè)結(jié)果與試驗(yàn)結(jié)果的各曲線變化趨勢(shì)相同，但兩者偏差較大，大部分工況點(diǎn)的揚(yáng)程和轉(zhuǎn)矩曲線偏差在±10%以內(nèi)，小部分點(diǎn)在±20%以內(nèi).

2) 在核主泵研發(fā)過(guò)程中，對(duì)核主泵性能的最終評(píng)估，大部分工況可采用數(shù)值模擬的方法來(lái)預(yù)測(cè)外特性數(shù)據(jù)，對(duì)一些極端工況可依靠試驗(yàn)，從而可節(jié)約試驗(yàn)成本，縮短研發(fā)周期.

3) 對(duì)一些極端工況(超小和超大流量)以及核主泵反轉(zhuǎn)正流工況的數(shù)值模擬方法還需要進(jìn)一步研究，以提高預(yù)測(cè)結(jié)果的準(zhǔn)確性.