核主泵小流量工況壓力脈動(dòng)特性

朱榮生，龍 云，付 強(qiáng)，袁壽其，王秀禮

（江蘇大學(xué) 流體機(jī)械工程技術(shù)研究中心，江蘇 鎮(zhèn)江 212013）

核反應(yīng)堆冷卻劑主循環(huán)泵簡(jiǎn)稱(chēng)核主泵，是確保核電站安全和可靠運(yùn)行的最關(guān)鍵動(dòng)力設(shè)備，屬于核I級(jí)泵，是核島內(nèi)唯一高速旋轉(zhuǎn)設(shè)備，也是一回路的壓力邊界之一。核主泵長(zhǎng)期穩(wěn)定安全地運(yùn)行對(duì)冷卻堆芯、以及防止核電站事故的發(fā)生尤為重要，因此核主泵常被喻為核電站的心臟［1］。從瞬變工況下減少熱應(yīng)力和主泵運(yùn)行工況特殊性的觀點(diǎn)出發(fā)，輕水壓水堆的主泵大都采用軸對(duì)稱(chēng)的桶型或者準(zhǔn)球型泵殼，這使得主泵的結(jié)構(gòu)明顯區(qū)別常規(guī)離心泵。

核主泵在小流量工況運(yùn)行時(shí)產(chǎn)生不穩(wěn)定流動(dòng)狀態(tài)，葉輪高速旋轉(zhuǎn)導(dǎo)致的葉輪和導(dǎo)葉間的動(dòng)靜干涉作用、邊界層分離、二次回流，使泵產(chǎn)生振動(dòng)、噪聲、壓力脈動(dòng)等現(xiàn)象，嚴(yán)重時(shí)甚至?xí)p害設(shè)備［2－4］，核主泵運(yùn)行不穩(wěn)定極易引發(fā)核事故，甚至造成核泄漏。因此針對(duì)小流量工況下核主泵內(nèi)部壓力脈動(dòng)的研究，對(duì)降低泵的振動(dòng)和提高反應(yīng)堆系統(tǒng)的穩(wěn)定性有實(shí)際意義。

國(guó)內(nèi)外學(xué)者針對(duì)泵的振動(dòng)已有一些相關(guān)研究。陳向陽(yáng)等［6］對(duì)國(guó)內(nèi)某300MW核電站主泵壓力脈動(dòng)進(jìn)行了研究，分析了葉輪和導(dǎo)葉附近的壓力脈動(dòng)特性。朱榮生等［7］通過(guò)對(duì)13種壓水室出口收縮角的研究，分析了核主泵流道內(nèi)的回流對(duì)壓力脈動(dòng)的影響。Yuan等［8］認(rèn)為葉輪流道與蝸殼在黏滯邊界層發(fā)生耦合，為了正確研究葉輪與蝸殼的耦合作用，必須同時(shí)考慮葉輪流道與蝸殼中的流動(dòng)情況。Kazem等［9］對(duì)德黑蘭反應(yīng)堆泵的啟動(dòng)瞬態(tài)進(jìn)行了模擬，通過(guò)對(duì)比分析其性能曲線(xiàn)，證明了主泵性能的優(yōu)越性。In Soo Koo等［10］對(duì)反應(yīng)堆主泵的振動(dòng)監(jiān)測(cè)、診斷系統(tǒng)進(jìn)行了分析，為查找主泵振動(dòng)的原因提供了方法。Wang等［11］試驗(yàn)研究了導(dǎo)葉在非設(shè)計(jì)點(diǎn)的不穩(wěn)定情況，得出其主要頻率的壓力脈動(dòng)出現(xiàn)在葉頻倍頻處。

本文應(yīng)用CFX軟件對(duì)AP1000核主泵進(jìn)行數(shù)值模擬分析，在偏小流量的多工況下同時(shí)監(jiān)測(cè)葉輪與導(dǎo)葉流道內(nèi)壓力脈動(dòng)的情況，以揭示壓力脈動(dòng)與泵性能之間的關(guān)系，為進(jìn)一步預(yù)測(cè)核主泵小流量非定常流動(dòng)誘發(fā)壓力脈動(dòng)提供依據(jù)。

1 數(shù)值計(jì)算方法

1．1 核主泵的基本參數(shù)

AP1000核主泵基本參數(shù)為：額定流量Qn＝17 886 m3／h，額定揚(yáng)程 H＝111．3 m，轉(zhuǎn)速 n＝1 750 r／min，比轉(zhuǎn)速ns＝344。對(duì)葉輪、導(dǎo)葉、泵體等水力部件進(jìn)行設(shè)計(jì)。采用混流式葉輪，扭曲空間導(dǎo)葉?？紤]耐壓和運(yùn)行安全性，泵殼設(shè)計(jì)成類(lèi)似球形［12］。利用Pro／E軟件建立葉輪、導(dǎo)葉、泵殼的幾何模型，主泵三維結(jié)構(gòu)與網(wǎng)格劃分如圖1所示。

圖1 三維造型與網(wǎng)格劃分Fig．1 Three dimensional design and meshing

1．2 網(wǎng)格及計(jì)算區(qū)域劃分

根據(jù)核主泵的流動(dòng)特性將流動(dòng)區(qū)域分為4個(gè)區(qū)域：靜止區(qū)1（進(jìn)口流道），靜止區(qū)2（類(lèi)球形壓水室），靜止區(qū)3（導(dǎo)葉），旋轉(zhuǎn)區(qū)（葉輪）。

為了獲得更穩(wěn)定的流態(tài)，延長(zhǎng)一定程度的葉輪進(jìn)口段。在Workbench中對(duì)三維模型進(jìn)行網(wǎng)格劃分，在確保網(wǎng)格的計(jì)算精度和計(jì)算結(jié)果準(zhǔn)確性的基礎(chǔ)上，由于混合網(wǎng)格技術(shù)具有結(jié)構(gòu)化與非結(jié)構(gòu)化網(wǎng)格兼有的優(yōu)點(diǎn)，并且生成方便、快速，采用自動(dòng)劃分法對(duì)計(jì)算區(qū)域進(jìn)行網(wǎng)格劃分，對(duì)葉輪調(diào)整單元尺寸加密劃分，結(jié)果如下：進(jìn)口流道的網(wǎng)格數(shù)為 75 840，蝸殼網(wǎng)格數(shù)為660 446，葉輪網(wǎng)格數(shù)為1 722 717，導(dǎo)葉網(wǎng)格數(shù)883 358，網(wǎng)格無(wú)關(guān)性驗(yàn)證良好。三維網(wǎng)格如圖1（b）所示。

1．3 計(jì)算方法及邊界條件

在流量一定時(shí)，為了得到更加準(zhǔn)確的速度和壓力梯度，進(jìn)口采用壓力進(jìn)口條件，出口采用質(zhì)量出流邊界條件。壁面采用無(wú)滑移壁面邊界條件。為了更好的處理流動(dòng)邊界層，在近壁區(qū)域采用標(biāo)準(zhǔn)壁面函數(shù)。輸送介質(zhì)為清水［13］。計(jì)算過(guò)程中的亞松弛因子均采用CFX軟件的默認(rèn)值，殘差收斂精度設(shè)置為10－4。

1．4 非定常模擬設(shè)置及監(jiān)測(cè)點(diǎn)的選取

葉輪流道內(nèi)的水體為旋轉(zhuǎn)體，蝸殼內(nèi)水體為非旋轉(zhuǎn)體，以定常的收斂解作為非定常計(jì)算初始條件。非定常計(jì)算中的交界面設(shè)置為T(mén)ransient Rotor-Stator模式，該交界面對(duì)于兩部分水體間的動(dòng)靜干涉有著重要作用。葉輪旋轉(zhuǎn)4個(gè)周期，總計(jì)算時(shí)間0．137 1 s，葉輪每轉(zhuǎn)2°作為一個(gè)時(shí)間步長(zhǎng)，時(shí)間步長(zhǎng)為2×10－4s。每經(jīng)過(guò)171個(gè)時(shí)間步長(zhǎng)，葉輪旋轉(zhuǎn)一周。選取第4個(gè)周期的結(jié)果用于分析。葉輪葉片數(shù)為5，導(dǎo)葉葉片數(shù)為11，葉輪的轉(zhuǎn)動(dòng)頻率f＝29．2 Hz，葉頻為T(mén)＝146 Hz。為了監(jiān)測(cè)不同工況下核主泵內(nèi)部壓力脈動(dòng)，沿水流方向在葉輪流道內(nèi)依次選取監(jiān)測(cè)點(diǎn)Y1、Y2、Y3、Y4及在導(dǎo)葉流道內(nèi)依次選取監(jiān)測(cè)點(diǎn)G1、G2、G3、G4。在離葉輪出口1．5 mm和離導(dǎo)葉出口邊10 mm圓上每隔15°取一點(diǎn)，各取24點(diǎn)，分別記為A1～A24、B1～B24。所有的點(diǎn)均處在中截面上。各監(jiān)測(cè)點(diǎn)具體設(shè)置位置如圖2所示。

2 計(jì)算結(jié)果及分析

由于流體的黏性作用以及旋轉(zhuǎn)部件與靜止部件的動(dòng)靜干擾，使得離心泵內(nèi)的流場(chǎng)呈現(xiàn)非定常的流動(dòng)特征。這種流動(dòng)特征引起流場(chǎng)的壓力脈動(dòng)，這將引起噪聲，同時(shí)在葉片上產(chǎn)生一個(gè)交變作用力，使葉片發(fā)生振動(dòng)。

2．1 圓周方向壓力脈動(dòng)時(shí)域圖

為了分析葉輪與導(dǎo)葉動(dòng)靜干涉以及導(dǎo)葉與蝸殼對(duì)流體造成的不穩(wěn)定流動(dòng)，在葉輪和導(dǎo)葉出口邊分別取點(diǎn)A1～A24、B1～B24，監(jiān)測(cè)流體靜壓變化情況，監(jiān)測(cè)結(jié)果如圖3所示。

圖2 泵內(nèi)各監(jiān)測(cè)點(diǎn)示意Fig．2 Indicator points in the pump

圖3 靜壓分布情況Fig．3 Static pressure distributions

從圖3（a）可以看出，葉輪與導(dǎo)葉間隙存在較大的壓力脈動(dòng)，葉輪與導(dǎo)葉的相互干擾以及非設(shè)計(jì)工況點(diǎn)的二次回流對(duì)流體造成不穩(wěn)定流動(dòng)。流體壓力波動(dòng)沿圓周方向分布并不均勻，基本成正弦變化規(guī)律，各工況壓力脈動(dòng)的最大值出現(xiàn)在 45°、120°、210°、300°附近，對(duì)應(yīng)于圖2所示葉片處監(jiān)測(cè)點(diǎn)位置。設(shè)計(jì)工況點(diǎn)的壓力脈動(dòng)較為均勻，非設(shè)計(jì)工況點(diǎn)壓力變化較大，其中0．2 Qn工況的壓力脈動(dòng)最大，最大值為2．655 MPa，且波動(dòng)幅度最大，從 0．2Qn到 1．0Qn，隨著流量變大，靜壓及其波動(dòng)幅度逐漸降低，沿周向每隔約30°出現(xiàn)一次峰值，造成的原因是導(dǎo)葉葉片數(shù)影響葉輪出口壓力變化。各均值隨流量變大而降低，但0．6Qn的周向壓力均值大于0．4Qn。

從圖3（b）可以看出，導(dǎo)葉出口附近也存在壓力脈動(dòng)，但峰值和幅度比葉輪出口邊小，沿周向的壓力脈動(dòng)也略為平緩，這是因?yàn)閷?dǎo)葉起到使流動(dòng)穩(wěn)定的作用。和葉輪出口附近壓力脈動(dòng)一樣，流體壓力波動(dòng)沿圓周方向分布并不均勻，基本成正弦變化規(guī)律，每隔30°出現(xiàn)一次峰值。在240°～300°之間，各流量下的壓力脈動(dòng)存在明顯的一直上升突降的過(guò)程，畸變發(fā)生在285°，90°～240°之間的壓力脈動(dòng)較為平緩，0°即對(duì)應(yīng)于出口位置的壓力脈動(dòng)最小，這主要是因?yàn)榍蛐螇核覍?duì)稱(chēng)結(jié)構(gòu)影響了流體在壓水室流動(dòng)方向末端的出流。在圖中較為明顯地看出隨著流量變大，導(dǎo)葉出口邊壓力脈動(dòng)呈現(xiàn)梯度變化，各工況點(diǎn)的壓力脈動(dòng)變化規(guī)律近乎一致。0．2Qn的壓力脈動(dòng)最大，依次分別為 0．3Qn、0．6Qn、0．4Qn、1．0Qn，與葉輪出口壓力脈動(dòng)均值變化規(guī)律一致。

2．2 葉輪流道內(nèi)壓力脈動(dòng)時(shí)域與頻域分析

圖4 點(diǎn)Y1，Y2，Y3和Y4的時(shí)域和頻域特性Fig．4 Time and frequency domain characteristics of fluctuation at point Y1，Y2，Y3，Y4

圖4 為葉輪流道監(jiān)測(cè)點(diǎn)Y1～Y4壓力脈動(dòng)的時(shí)域特性和通過(guò)快速傅里葉變化（FFT）得到的對(duì)應(yīng)頻域特性。此處約定波動(dòng)幅度 CA＝（pmax－pmin）／pmax（pmax為對(duì)應(yīng)點(diǎn)的最大壓力；pmin為對(duì)應(yīng)點(diǎn)的最小壓力）。頻譜圖中，橫坐標(biāo)為頻率值，縱坐標(biāo)為各個(gè)頻率值對(duì)應(yīng)的壓力脈動(dòng)能量幅值。由圖4可以看出，在葉輪流道內(nèi)壓力波動(dòng)比較明顯，主要是由于葉輪和導(dǎo)葉的動(dòng)靜干涉以及小流量時(shí)葉輪內(nèi)的渦流造成的。葉輪內(nèi)各監(jiān)測(cè)點(diǎn)處在設(shè)計(jì)工況點(diǎn)的壓力脈動(dòng)呈現(xiàn)明顯的規(guī)律性變化，在一個(gè)周期內(nèi)出現(xiàn)11個(gè)波峰和11個(gè)波谷，這是由葉輪和導(dǎo)葉間動(dòng)靜干涉作用引起的。而葉輪內(nèi)各監(jiān)測(cè)點(diǎn)在0．2Qn～0．6Qn工況點(diǎn)的壓力脈動(dòng)雖也出現(xiàn)約11個(gè)波峰和11個(gè)波谷，但其規(guī)律性不顯著，同時(shí)其波動(dòng)幅度比設(shè)計(jì)工況點(diǎn)更為劇烈，這是因?yàn)樵谛×髁繀^(qū)域葉輪內(nèi)出現(xiàn)的二次回流及葉輪和導(dǎo)葉動(dòng)靜干涉成為流體振動(dòng)的動(dòng)力源，從而產(chǎn)生壓力脈動(dòng)。

在監(jiān)測(cè)點(diǎn) Y1點(diǎn)，0．2Qn、0．3Qn、1．0Qn的平均壓力很低，因此 pmax很小，從而導(dǎo)致 CA異常大，0．4Qn和0．6Qn的壓力脈動(dòng)均值明顯大于 0．2Qn、0．3Qn、1．0Qn，造成的原因可能是在0．2Qn和0．3Qn工況葉輪進(jìn)口產(chǎn)生回流。監(jiān)測(cè)點(diǎn)Y2和Y3壓力脈動(dòng)均值變化趨勢(shì)與Y1點(diǎn)相似，說(shuō)明在0．2Qn和 0．3Qn工況葉輪的中段渦結(jié)構(gòu)仍然存在，但此時(shí)回流結(jié)構(gòu)對(duì)Y3點(diǎn)的影響略有減弱。監(jiān)測(cè)點(diǎn)Y1～Y3在0．4Qn和0．6Qn工況的壓力脈動(dòng)遠(yuǎn)大于其它工況點(diǎn)，原因是隨著流量的減小，尤其是在0．2Qn和0．3Qn工況，葉輪入口處的回流區(qū)域逐漸加劇，導(dǎo)致葉輪入口堵塞效應(yīng)增強(qiáng)，流動(dòng)損失增加。而在監(jiān)測(cè)點(diǎn) Y4，從 0．2Qn～1．0Qn壓力脈動(dòng)均值呈遞減趨勢(shì)，這是因?yàn)閅4點(diǎn)處于葉輪流道出口中間，液流流動(dòng)通暢，受葉輪內(nèi)回流影響小，從而導(dǎo)致葉輪對(duì)不同工況下流體做功產(chǎn)生明顯的壓力梯度。

葉輪流道內(nèi) Y1、Y2、Y3、Y4監(jiān)測(cè)點(diǎn)在 0．4Qn的主頻為 138 Hz（0．94 T），其他工況的主頻為 331．5 Hz（2．27 T），壓力脈動(dòng)主要集中在低頻區(qū)，各監(jiān)測(cè)點(diǎn)的主頻與理論計(jì)算所得的葉輪頻率及倍頻之間存在一定的偏差，主要是因?yàn)樾×髁繒r(shí)葉輪內(nèi)回流引起流體的運(yùn)動(dòng)。壓力脈動(dòng)主要集中在低頻區(qū)，幾乎監(jiān)測(cè)不到高頻成分的存在，這是因?yàn)樵诒脙?nèi)部無(wú)空化發(fā)生的情況下，壓力脈動(dòng)的低頻信號(hào)主要是由葉輪和導(dǎo)葉間的動(dòng)靜干涉、高速旋轉(zhuǎn)葉片表面的流動(dòng)分離以及泵內(nèi)部產(chǎn)生的旋渦引起的。在設(shè)計(jì)工況點(diǎn)主頻處的脈動(dòng)幅值遠(yuǎn)高于諧波處的脈動(dòng)幅值，因此可以認(rèn)為設(shè)計(jì)點(diǎn)葉頻在由壓力脈動(dòng)誘發(fā)的泵振動(dòng)中起主導(dǎo)作用。而在小流量區(qū)域，主次頻處的脈動(dòng)幅值大小差別并不明顯，甚至在低頻區(qū)域存在多個(gè)次高峰值，這表明在小流量區(qū)域誘導(dǎo)泵振動(dòng)的原因不僅僅是葉輪導(dǎo)葉的動(dòng)靜干涉作用，還包括葉輪流道內(nèi)的脫流、湍流、渦流，且兩者所占比重相當(dāng)。

2．3 導(dǎo)葉流道內(nèi)壓力脈動(dòng)時(shí)域與頻域分析

圖5為導(dǎo)葉流道監(jiān)測(cè)點(diǎn)G1、G3、G4壓力脈動(dòng)的時(shí)域特性和對(duì)應(yīng)頻域特性。從圖中可以看出，在導(dǎo)葉流道內(nèi)壓力波動(dòng)比較明顯，小流量的壓力波動(dòng)更為劇烈，主要是由于葉輪和導(dǎo)葉的動(dòng)靜干涉以及球形壓水室對(duì)導(dǎo)葉出流的影響。導(dǎo)葉內(nèi)各監(jiān)測(cè)點(diǎn)處在設(shè)計(jì)工況點(diǎn)的壓力脈動(dòng)呈現(xiàn)明顯的規(guī)律性變化，在一個(gè)周期內(nèi)出現(xiàn)5個(gè)大波峰和5個(gè)大波谷，每個(gè)大波峰又包含1個(gè)小波峰和一個(gè)小波谷，小波峰峰值遠(yuǎn)小于大波峰，這說(shuō)明宏觀的葉輪導(dǎo)葉動(dòng)靜干涉對(duì)壓力脈動(dòng)產(chǎn)生起主要作用。而導(dǎo)葉內(nèi)G1和G3監(jiān)測(cè)點(diǎn)在0．2Qn～0．6Qn工況點(diǎn)的壓力脈動(dòng)雖也出現(xiàn)約5個(gè)波峰和5個(gè)波谷，但其規(guī)律性減弱，同時(shí)其波動(dòng)比設(shè)計(jì)工況點(diǎn)更為劇烈，而監(jiān)測(cè)點(diǎn)G4在0．2Qn～0．6Qn工況并未表現(xiàn)出明顯的5個(gè)波峰和5個(gè)波谷，且規(guī)律性不顯著，這說(shuō)明小流量區(qū)域的壓力脈動(dòng)除了動(dòng)靜干涉影響外，導(dǎo)葉內(nèi)的二次回流也起到明顯作用，尤其是導(dǎo)葉出口附近G4點(diǎn)的二次回流。

圖5 點(diǎn)G1，G2，G3，G4的時(shí)域和頻域特性Fig．5 Time and frequency domain characteristics of fluctuation at point G1，G2，G3，G4

小流量工況下導(dǎo)葉流道內(nèi)Y1、Y2、Y3、Y4點(diǎn)的壓力脈動(dòng)隨著流量的增加，壓力脈動(dòng)均值不斷減小，但G1和G2點(diǎn)在0．4Qn略有畸變，造成的原因是0．2Qn和0．3Qn工況流道內(nèi)的回流產(chǎn)生較大的能量損失。而除了G4點(diǎn)外，其它各點(diǎn)波動(dòng)幅度也基本隨著流量增加呈下降趨勢(shì)，說(shuō)明隨著流量的增加，葉輪和導(dǎo)葉內(nèi)的回流區(qū)域減小，流體的通順能力逐漸增強(qiáng)。在頻譜圖中，導(dǎo)葉流道內(nèi)各監(jiān)測(cè)點(diǎn)的壓力脈動(dòng)主頻為138 Hz（0．94T），這表明小流量區(qū)域的壓力脈動(dòng)雖也受導(dǎo)葉內(nèi)的回流影響，但葉輪和導(dǎo)葉的動(dòng)靜干涉占主要作用。

2．4 葉輪出口邊周向壓力脈動(dòng)

圖6為葉輪出口周向監(jiān)測(cè)點(diǎn) A1、A7、A13和 A19的壓力脈動(dòng)的時(shí)域特性與對(duì)應(yīng)頻域特性。由圖6可以看出，葉輪出口周向壓力波動(dòng)比較明顯，各監(jiān)測(cè)點(diǎn)處在設(shè)計(jì)工況點(diǎn)的壓力脈動(dòng)呈現(xiàn)規(guī)律性變化，A1點(diǎn)在一個(gè)周期內(nèi)出現(xiàn)11個(gè)波峰和11個(gè)波谷，A7、A13和A19點(diǎn)則在一個(gè)周期內(nèi)出現(xiàn)5個(gè)大波峰和5個(gè)大波谷，每個(gè)波峰又包含1個(gè)小波峰和1個(gè)大波峰，這主要是由葉輪和導(dǎo)葉間動(dòng)靜干涉作用以及葉輪出口射流尾跡引起的。監(jiān)測(cè)點(diǎn)A1在0．2Qn～0．6Qn工況點(diǎn)的壓力脈動(dòng)也出現(xiàn)約11個(gè)波峰和11個(gè)波谷，其周期性規(guī)律減弱。監(jiān)測(cè)點(diǎn) A7、A13和 A19在 0．2Qn～0．6Qn工況點(diǎn)的壓力脈動(dòng)基本呈現(xiàn)5個(gè)波峰和5個(gè)波谷，其規(guī)律性不明顯，波動(dòng)更為劇烈。

圖6 點(diǎn)A1，A7，A13和A19的時(shí)域和頻域特性Fig．6 Time and frequency domain characteristics of fluctuation at point A1，A7，A13，A19

各監(jiān)測(cè)點(diǎn)的壓力脈動(dòng)均值隨著流量的增加呈現(xiàn)下降趨勢(shì)，但A1和A7點(diǎn)在0．4Qn略有畸變，波動(dòng)幅度隨著流量的增加呈減小趨勢(shì)，但A1和A19點(diǎn)在1．0Qn略增，造成的原因是小流量區(qū)域流道內(nèi)的回流產(chǎn)生較大的能量損失，葉輪出口周向監(jiān)測(cè)點(diǎn)A1、A7、A13和A19的壓力脈動(dòng)峰值信號(hào)主要產(chǎn)生在葉頻及倍葉頻處，壓力脈動(dòng)主要集中在低頻區(qū)。但在A7點(diǎn)小流量壓力脈動(dòng)規(guī)律性差且主頻混亂，這表明球形壓水室導(dǎo)致葉輪出口液流的不均勻。

2．5 導(dǎo)葉出口邊周向壓力脈動(dòng)

圖7為小流量工況監(jiān)測(cè)點(diǎn)B1、B7、B13和B19的壓力脈動(dòng)的時(shí)域特性與對(duì)應(yīng)頻域特性。由圖7可以看出，導(dǎo)葉出口周向壓力波動(dòng)比較明顯，各監(jiān)測(cè)點(diǎn)處在設(shè)計(jì)工況點(diǎn)的壓力脈動(dòng)呈現(xiàn)規(guī)律性變化，在一個(gè)周期內(nèi)出現(xiàn)5個(gè)大波峰和5個(gè)大波谷、5個(gè)小波峰和5個(gè)小波谷，這表明在設(shè)計(jì)點(diǎn)葉輪和導(dǎo)葉間動(dòng)靜干涉作用占主導(dǎo)作用。各監(jiān)測(cè)點(diǎn)在0．2Qn壓力脈動(dòng)出現(xiàn)5個(gè)波峰和5個(gè)波谷，其周期性規(guī)律比設(shè)計(jì)點(diǎn)減弱；在0．4Qn工況各監(jiān)測(cè)點(diǎn)的壓力脈動(dòng)基本呈現(xiàn)2個(gè)大波峰和2個(gè)大波谷，近似成兩個(gè)拱形，每個(gè)拱形又包含若干波峰和波谷；在0．6Qn工況各監(jiān)測(cè)點(diǎn)的拱形大波峰被抑制，呈現(xiàn)多個(gè)無(wú)規(guī)律的波形，這是因?yàn)槭軐?dǎo)葉末端回流影響和球形壓水室束縛，在0．4Qn工況導(dǎo)葉出口環(huán)向產(chǎn)生不均勻的類(lèi)似喘息的脈動(dòng)，從而導(dǎo)致壓力脈動(dòng)呈2個(gè)大拱形，而隨著流量的增加，導(dǎo)葉出口回流區(qū)域減小，流道變得通順，這種喘振被抑制。

小流量各工況下導(dǎo)葉出口周向壓力脈動(dòng)主要集中在低頻區(qū)，峰值信號(hào)主要產(chǎn)生在葉頻處，0．3Qn、0．4Qn、0．6Qn的壓力脈動(dòng)特性和設(shè)計(jì)點(diǎn)的區(qū)別很大，說(shuō)明小流量區(qū)域?qū)~出口周向的流動(dòng)極不穩(wěn)定，造成的原因主要是受導(dǎo)葉出口回流和球形壓水室束縛。從B7點(diǎn)到B19點(diǎn)，壓力脈動(dòng)均值不斷增加，這是因?yàn)槲仛な占黧w造成的。導(dǎo)葉出口周向的壓力脈動(dòng)幅度明顯小于葉輪周向出口，說(shuō)明導(dǎo)葉起到了很好的導(dǎo)流作用，對(duì)降低泵的振動(dòng)起到很好的作用。

圖7 小流量各工況下點(diǎn)B1，B7，B13和B19的時(shí)域和頻域特性Fig．7 Time and frequency domain characteristics of fluctuation at point B1，B7，B13，B19 in small flux

3 流線(xiàn)分布

核主泵內(nèi)部產(chǎn)生的壓力脈動(dòng)，主要原因?yàn)槿~輪與蝸殼存在的動(dòng)靜干涉，球形壓水室非對(duì)稱(chēng)結(jié)構(gòu)，及流道內(nèi)的回流，因此，分析不同流量下流體的流動(dòng)狀態(tài)對(duì)壓力脈動(dòng)分析有重要意義。圖8為小流量工況下中截面的流線(xiàn)分布，從圖中可以看出設(shè)計(jì)工況流線(xiàn)非常光滑，而在0．2Qn～0．6Qn工況存在明顯的回流，主要存在于葉輪和導(dǎo)葉進(jìn)出口處，且偏離設(shè)計(jì)工況越大，流線(xiàn)扭曲越嚴(yán)重，回流面積增大，回流造成流道堵塞，速度方向產(chǎn)生突變，這種突變必然在流道內(nèi)產(chǎn)生一組方向相反的正負(fù)壓力波，從而導(dǎo)致小流量工況壓力脈動(dòng)波動(dòng)劇烈且周期性差。小流量工況下球形壓水室出口出現(xiàn)明顯的漩渦，這表明對(duì)稱(chēng)式的類(lèi)球形壓水室不利于液體出流。

圖8 小流量區(qū)域流線(xiàn)分布Fig．8 Streamline distribution at small flow rate conditions

4 結(jié) 論

（1）核主泵內(nèi)壓力脈動(dòng)明顯，葉頻在由壓力脈動(dòng)誘發(fā)的振動(dòng)中起主導(dǎo)作用，主要表現(xiàn)為葉輪和導(dǎo)葉間的動(dòng)靜干涉。

（2）葉輪導(dǎo)葉流道內(nèi)的回流造成了小流量工況葉輪和導(dǎo)葉流道及其周向的不穩(wěn)定壓力脈動(dòng)，回流主要存在于葉輪和導(dǎo)葉進(jìn)出口位置，因此該區(qū)域的壓力波動(dòng)劇烈且周期性差。

（3）核主泵的振動(dòng)，不利于核電站的安全穩(wěn)定運(yùn)行，通過(guò)對(duì)小流量工況的壓力脈動(dòng)分析，對(duì)預(yù)測(cè)核主泵在極端工況下的動(dòng)態(tài)特性和推進(jìn)核主泵國(guó)產(chǎn)化具有十分重要的意義。

［1］蔡龍，張麗平．淺談壓水堆核電站主泵［J］．水泵技術(shù)，2007（4）：1－9．CAI Long，ZHANG Li-ping．Discussion of PWR main pump［J］．Pump Technology，2007（4）：1－9．

［2］黃經(jīng)國(guó)．壓水堆核電廠冷卻劑主循環(huán)泵技術(shù)歷程和發(fā)展［J］．水泵技術(shù)，2009（4）：1－8．HUANG Jing-guo．PWR nuclear power plant primary coolant circulation pump technology process and development［J］．Pump Technology，2009（4）：1－8．

［3］Runkel J，Stegemann D，Vortriede A．Operating experience with an on-line vibration control system for PWR main coolant pumps［J］．Nuclear Engineering Design，1998，183：157－167．

［4］Jose G，Carlos S．Unsteady flow structure and global variables in a centrifugal pump［J］．Journal of Fluids Engineering，2006，128（9）：937－946．

［5］袁壽其，薛菲，袁建平，等．離心泵壓力脈動(dòng)對(duì)流動(dòng)噪聲影響的試驗(yàn)研究［J］．排灌機(jī)械，2009，27（5）：287－290．YUAN Shou-qi， XUE Fei， YUAN Jian-ping， et al．Experimental study on impact of pressure fluctuation on flow-noise centrifugal pump［J］． Drainage and Irrigation Machinery，2009，27（5）：287－290．

［6］陳向陽(yáng)，袁丹青，楊敏官，等．300MW級(jí)核電站主泵壓力脈動(dòng)研究 ［J］．核動(dòng)力工程，2010，31（3）：78－82．CHEN Xiang-yang，YUAN Dan-qing，YANG Min-guan，et al．Study on pressure fluctuation in main pump of 300 MW nuclear power plants［J］．Nuclear Power Engineering，2010，31（3）：78－82．

［7］朱榮生，李小龍，袁壽其，等．反應(yīng)堆主泵壓水室出口收縮角對(duì)水力性能的影響［J］．核動(dòng)力工程，2012，33（2）：97－102．ZHU Rong-sheng，LI Xiao-long， YUAN Shou-qi， et al．Effect of pumping chamber outlet contraction angle on hydraulic performance of main nuclear reactor pump［J］．Nuclear Power Engineering，2012，33（2）：97－102．

［8］ Yuan Shou-qi，Ni Yong-yan．Unsteady turbulent simulation and pressure fluctuation analysis for centrifugal pumps［J］．Chinese Journal of Mechanical Engineering，2009，22（1）：1－6．

［9］Kazem F，Anis Bousbia-salah，F(xiàn)ransceso D A．A model for the analysis of pump start-up transients in Tehran research reactor［J］．Progress in Nuclear Energy，2007，49（7）：499－510．

［10］ In Soo Koo，Whan Woo Kim．The development of reactor coolant pump vibration monitoring and a diagnostic system in the nuclear power plant［J］．ISA Transactions，2000，39（3）：309－316．

［11］Wang H，Tsukamoto H．Experimental and numerical study of unsteady flow in a diffuser pump at off-design conditions［J］．Journal of Fluids Engineering，2003，125（9）：767－778．

［12］王春林，彭娜，康燦，等．模型泵內(nèi)部流場(chǎng)的數(shù)值模擬［J］．核動(dòng)力工程，2009，30（4）：81－85．WANG Chun-lin，PENG Na，KANG Can，et al．Numerical simulation of interior flow field of nuclear model pump［J］．Nuclear Power Engineering，2009，30（4）：81－85．

［13］張野，王曉放，介紅恩．壓水堆冷卻劑中硼酸濃度對(duì)核主泵性能影響研究［J］．核動(dòng)力工程，2011，32（4）：95－98．ZHANG Ye，WANGXiao-fang，JIE Hong-en．Effect of boric acid concentration on reactor coolant pump performance in PWRs［J］．Nuclear Power Engineering，2011，32（4）：95－98．