陸 榮, 袁建平, 朱鈺雯, 付燕霞, 周幫倫

(江蘇大學(xué) 國家水泵及系統(tǒng)工程技術(shù)研究中心，江蘇 鎮(zhèn)江 212013)

余熱排出泵內(nèi)部不穩(wěn)定流動數(shù)值模擬與實驗研究

陸 榮, 袁建平, 朱鈺雯, 付燕霞, 周幫倫

(江蘇大學(xué) 國家水泵及系統(tǒng)工程技術(shù)研究中心，江蘇 鎮(zhèn)江 212013)

為了研究余熱排出泵在多工況下內(nèi)部流動特性，基于ANSYS CFX軟件，采用SST湍流模型，對模型泵進行三維非定常數(shù)值模擬，獲得了不同工況下余熱排出泵的水力性能、內(nèi)部流場結(jié)構(gòu)和壓力脈動特性，同時展開壓力脈動實驗研究，并與計算結(jié)果進行對比。研究結(jié)果表明：大流量(1.2Qd)和設(shè)計流量(1.0Qd)工況下，葉輪和導(dǎo)葉內(nèi)部流動比較穩(wěn)定，隨著流量的減小葉片進口背面附近開始形成失速旋渦，流道內(nèi)均發(fā)生不同程度的流動分離，且沿著流道向出口處發(fā)展；葉輪出口壓力脈動主頻為7fz，受導(dǎo)葉葉片數(shù)影響；導(dǎo)葉和蝸殼出口的主頻均為5fz，主要由葉頻決定；設(shè)計流量下各監(jiān)測點處壓力脈動系數(shù)幅值最小，越往小流量工況，幅值越大；說明在小流量工況下余熱排出泵內(nèi)部出現(xiàn)了不穩(wěn)定流動現(xiàn)象。

余熱排出泵；壓力脈動；流動分離；數(shù)值模擬；實驗研究

余熱排出泵是余熱排出系統(tǒng)的主要組成部分，是除核主泵之外唯一布置在核島之內(nèi)的二級泵，是關(guān)系到核島能否安全停堆的核心裝備[1]，余熱排出泵運行范圍較大，在非設(shè)計工況下運行時，泵內(nèi)部容易產(chǎn)生如回流、失速、空化等不穩(wěn)定流動現(xiàn)象[2]，這些不穩(wěn)定流動會破壞流場結(jié)構(gòu)限制泵的高效運行范圍，降低泵運行的穩(wěn)定性，嚴(yán)重時可引發(fā)機組的共振[3-5]。張人會等[6]對葉輪進口附近流動在偏離設(shè)計工況下的回流和空化特性進行了系統(tǒng)的描述, 闡述了它們的誘發(fā)機理及其對泵性能的影響。得出：葉輪進口附近的逆壓力梯度是導(dǎo)致回流形成的根本原因，減小葉片角分布系數(shù), 可使葉片最大載荷系數(shù)向葉輪出口方向移動，降低回流的誘發(fā)流量。JOHNSON等[7]運用LDV測試技術(shù)對離心泵葉輪內(nèi)流場的流動情況進行測量，研究發(fā)現(xiàn)：小流量工況下，在葉輪旋轉(zhuǎn)一個周期過程中，葉輪流道內(nèi)的失速狀態(tài)和非失速狀態(tài)交替出現(xiàn)。張睿等[8]基于非定常的 RANS方程，采用濾波器湍流模型對某一軸流泵模型在不同工況下的內(nèi)部流場進行數(shù)值計算，結(jié)果表明：在小流量工況下泵內(nèi)部發(fā)生回流和失速，并提出了改善失速工況下軸流泵水力性能的方法。

本文以余熱排出泵水力樣機為研究對象，對其進行數(shù)值模擬和實驗研究，分析在各個工況下葉輪、導(dǎo)葉內(nèi)的流動現(xiàn)象，重點掌握不穩(wěn)定流動現(xiàn)象隨流量和時間的變化規(guī)律，為提高余熱排出泵運行穩(wěn)定性和安全可靠性提供一定的參考。

1 幾何模型及計算方法

1.1 幾何模型

本次選用的余熱排出泵為臥式、單級單吸、后開門式的離心泵，其結(jié)構(gòu)示意圖如圖1所示。其主要設(shè)計參數(shù)：流量Qd=312 m3/h，揚程Hd=38 m，轉(zhuǎn)速n=1 490 r/min，比轉(zhuǎn)速ns=105。

通過三維造型軟件Pro/E對余熱排出泵水力樣機進行三維造型，整個計算區(qū)域包括進口段、葉輪、導(dǎo)葉、蝸殼及出口段。為了使流體在進入葉輪前能得到充分發(fā)展并較準(zhǔn)確地模擬流體的流動狀況，對進、出口段進行了適當(dāng)延長，如圖2所示。

圖1 余熱排出泵結(jié)構(gòu)圖Fig.1 Structure sketch drawing of residual heat removal pump

圖2 余熱排出泵三維模型Fig.2 3D model of residual heat removal pump

1.2 網(wǎng)格劃分及無關(guān)性檢驗

網(wǎng)格劃分是重要的數(shù)值計算前處理步驟，網(wǎng)格質(zhì)量和數(shù)量直接影響到計算結(jié)果準(zhǔn)確性。本文運用ANSYS ICEM軟件對計算域進行網(wǎng)格劃分，選用適合形貌復(fù)雜、質(zhì)量較高的六面體結(jié)構(gòu)化網(wǎng)格[9]，并進行網(wǎng)格無關(guān)性分析，其對應(yīng)的設(shè)計工況下的計算結(jié)果如表1、圖3所示。

表1 網(wǎng)格數(shù)目對比Tab.1 Comparison of meshes

圖3 網(wǎng)格無關(guān)性分析Fig.3 Check of grid independence

由圖3可知，網(wǎng)格數(shù)大于390萬后揚程趨于穩(wěn)定，且與設(shè)計值的誤差小于5%，為了節(jié)約計算時間及資源選用第四套的網(wǎng)格數(shù)進行計算。各過流部件的網(wǎng)格信息如表2所示，計算域的網(wǎng)格劃分如圖4所示。

表2 過流部件的網(wǎng)格基本信息Tab.2 Geometrical parameters of flow passage components

圖4 計算域網(wǎng)格Fig.4 Mesh of computational domain

1.3 湍流模型及邊界條件

SST湍流模型考慮了湍流剪切力的運輸，能夠適應(yīng)逆壓梯度變化的流動現(xiàn)象，能精確地預(yù)測流動的開始和逆壓力梯度下流體的分離[10-12]。所以文中選用SST湍流模型和連續(xù)方程使動量方程封閉，采用SIMPLEC算法求解三維不可壓縮雷諾平均N-S方程邊界件分別設(shè)置為進口為全壓進口，出口為質(zhì)量流量出口；壁面為無滑移、光滑、絕熱壁面；葉輪流道區(qū)域設(shè)置為旋轉(zhuǎn)坐標(biāo)系，進口段、導(dǎo)葉以及蝸殼區(qū)域設(shè)置為靜止坐標(biāo)系，計算收斂精度設(shè)置為1.00×10-5，網(wǎng)格節(jié)點之間的匹配方式設(shè)置為GGI模式[13]。

2 外特性預(yù)測及實驗驗證

為獲得余熱排出泵水力樣機的性能，對其進行外特性試驗，進而對比試驗值和計算值之間的誤差，驗證所采用數(shù)值計算方法的正確性和準(zhǔn)確性。本次試驗在江蘇大學(xué)流體機械研究中心實驗室開式試驗臺上進行性，試驗方法參照GB/T 3216—2005《回轉(zhuǎn)動力泵水力性能驗收試驗1級和2級》的相關(guān)規(guī)定。試驗裝置簡圖如圖5所示，試驗滿足2級精度要求。通過試驗獲得全流量范圍內(nèi)泵揚程H和效率η的真實數(shù)據(jù)，將試驗結(jié)果與計算結(jié)果進行對比，如圖6所示。

1-入口閥；2-進口壓力傳感器；3-水力樣機；4-彈性聯(lián)軸器；5-電機；6-出口壓力傳感器；7-渦輪流量計；8-出口閥圖5 試驗臺裝置簡圖Fig.5 Schematic of the pump test system

圖6 水力性能對比Fig.6 Comparison of hydraulic performance

從圖6中可以看出：數(shù)值計算得到的性能曲線和試驗測得的性能曲線的變化趨勢一致；小流量工況下，計算效率值與試驗效率值吻合較好；隨著流量的增大，計算揚程值與試驗揚程值吻合較好；揚程的最大誤差為4.453%，效率的最大誤差為4.735%，揚程誤差和效率誤差均小于5%，因此，本次所采用的數(shù)值計算方法能夠較準(zhǔn)確地預(yù)測余熱排出泵的外特性，具有較高的可信度。

2.1 葉輪內(nèi)部流場分析

圖7為不同流量工況下葉輪中截面速度流線圖?？梢钥闯?，不同流量工況下，葉片背面的流速均大于工作面。大流量(1.2Qd)和設(shè)計流量工況下，葉輪各流道內(nèi)的流動較規(guī)律，流體的相對速度基本沿葉片型線方向，只有部分流道內(nèi)有輕微的脫流現(xiàn)象。隨著流量的減小，葉片背面流動分離加劇，當(dāng)流量減小至0.6Qd時流道內(nèi)開始形成失速旋渦，并在葉輪出口產(chǎn)生與葉輪旋轉(zhuǎn)方向相反的軸向漩渦，隨著流量的減小，旋渦數(shù)量逐漸增加，尺寸不斷增大，以至堵塞部分流道。當(dāng)流量減小到0.4Qd時，多個流道同時出現(xiàn)了失速旋渦和軸向旋渦，且嚴(yán)重阻塞流道，在葉輪出口處具有明顯的射流-尾跡現(xiàn)象。

2.2 導(dǎo)葉內(nèi)部流場分析

圖8中，從大流量和設(shè)計流量工況下導(dǎo)葉流道內(nèi)速度分布規(guī)律基本一致，導(dǎo)葉進口無流動分離和回流現(xiàn)象產(chǎn)生，導(dǎo)葉部分流道的出口處存在分離渦；隨著流量的減小，導(dǎo)葉流道進口處的部分流體開始偏離主流方向，導(dǎo)葉各流道出口的流動分離不斷加劇。當(dāng)流量降低至0.6Qd時，導(dǎo)葉各流道進口均產(chǎn)生回流，且出口的分離旋渦尺寸擴大，當(dāng)流量減小至0.4Qd，在流道的進口處產(chǎn)生了失速旋渦，且出口處的流動分離得到進一步的發(fā)展，使流線變得十分紊亂，嚴(yán)重阻塞了流道。

圖7 不同流量工況下葉輪中截面速度流線圖Fig.7 The streamline of velocity on the middle section of the impeller at different flow rates

圖8 不同流量工況下導(dǎo)葉中截面速度流線圖Fig.8 The streamline of velocity on the middle section of the diffuser of at different flow rates

通過對多個工況下葉輪及導(dǎo)葉內(nèi)流動情況的分析可知，大流量下及設(shè)計流量下其內(nèi)部流動較規(guī)律，隨著流量的減小，葉輪和導(dǎo)葉流道內(nèi)相繼出現(xiàn)失速和流動分離等不穩(wěn)定流動現(xiàn)象，這些不穩(wěn)定流動現(xiàn)象容易引起泵性能下降、噪聲、振動等危害[14]，所以因盡量避免余熱排出泵在小流量工況下運行。

3 壓力脈動數(shù)值模擬與實驗研究

為研究余熱排出泵內(nèi)的壓力脈動特性，在泵內(nèi)中截面上設(shè)置了8個監(jiān)測點并進行非定常計算，為獲得較為穩(wěn)定的結(jié)果設(shè)置計算總時間為葉輪旋轉(zhuǎn)10個周期，并取最后一個周期內(nèi)的數(shù)據(jù)作壓力脈動分析，監(jiān)測點分布情況如圖9(a)所示。重點對葉輪出口P3，導(dǎo)葉出口P6，及蝸殼出口P8分析并對P8處設(shè)計工況下的壓力脈動進行實驗驗證。試驗以清水為介質(zhì)，由虛擬儀器采集系統(tǒng)獲得壓力脈動試驗數(shù)據(jù)，傳感器安裝位置如圖9(b)所示。

圖9 監(jiān)測點和壓力傳感器位置Fig.9 Location of monitoring points and pressure transmitter

為了直觀地反映壓力脈動的幅度，以便進行不同對象、不同位置間壓力脈動大小的比較，現(xiàn)引入壓力脈動系數(shù)這個參數(shù)，壓力脈動系數(shù)公式[15]為

(1)

為了分析葉輪內(nèi)的壓力脈動特性，選取葉輪旋轉(zhuǎn)第10個周期內(nèi)的壓力脈動數(shù)據(jù)進行頻譜分析，運用快速傅里葉變換(FFT)計算得到不同流量下各監(jiān)測點處處壓力脈動的頻域圖。葉輪的轉(zhuǎn)速為n=1 490 r/min，因此，定義葉輪的轉(zhuǎn)動頻率為fz=1 490/60 Hz=24.833 Hz，葉頻fb=124.165 Hz=5fz。

3.1 葉輪、導(dǎo)葉出口壓力脈動

圖10(a)表明，在4種不同流量下葉輪出口處壓力脈動的主要頻率均為7fz，等于導(dǎo)葉葉頻，這表明葉輪出口處監(jiān)測點的壓力脈動主要由轉(zhuǎn)頻決定，且受導(dǎo)葉的動靜干涉作用影響較大；隨著流量的減小，葉輪出口壓力脈動中的高頻成分逐漸減少，其脈動強度逐漸降低，而低頻脈動成分逐漸增多，其脈動強度也逐漸增強。

圖10(b)顯示，當(dāng)流量降低至0.4Qd時，壓力脈動的主頻為2fz，其余三種不同流量工況下的主頻均為葉頻fb，這說明導(dǎo)葉出口處的壓力脈動主要由葉頻決定；大流量工況下，導(dǎo)葉出口處壓力脈動的次主頻為35fz，設(shè)計流量和0.6Qd工況下壓力脈動的次主頻為15fz；當(dāng)流量降至0.4Qd時，次主頻則為葉頻fb。因此，導(dǎo)葉出口處壓力脈動的主頻隨流量的降低逐漸向低頻區(qū)靠近；同時隨著流量的減小，高頻壓力脈動成分逐漸減少，低頻壓力脈動強度逐漸增強，越偏離設(shè)計流量，壓力脈動強度越增強。

圖10 不同流量下葉輪、導(dǎo)葉出口壓力脈動頻域圖Fig.10 Frequency domain comparison at the outlet of impeller,diffuser under multi-conditions

在設(shè)計流量下葉輪與導(dǎo)葉出口點的壓力脈動系數(shù)幅值均達到最小，越往小流量點偏移，其幅值越大，說明余熱排出泵在設(shè)計流量下運行較安全，而在小流量工況下運行時，內(nèi)部流動十分不穩(wěn)定，由此引發(fā)的高幅壓力脈動易對機組造成損壞。

3.2 蝸殼出口壓力脈動

為了驗證數(shù)值模擬的準(zhǔn)確性，對設(shè)計流量下蝸殼出口P8的壓力脈動進行實驗對比，如圖11所示。實驗測得蝸殼出口管主頻頻率為5fz，由于采用頻率的原因模擬結(jié)果與實驗結(jié)果在高頻段稍有差別，但在低頻段兩者均采集得到了葉輪葉頻和轉(zhuǎn)頻處的頻率，結(jié)果基本吻合，說明蝸殼內(nèi)的壓力脈動主要由葉頻決定。因此，在余熱排出泵結(jié)構(gòu)設(shè)計時應(yīng)注意分析水泵主要結(jié)構(gòu)的固有模態(tài)是否與轉(zhuǎn)頻、葉頻成整數(shù)倍，并避免共振的發(fā)生。

圖11 蝸殼出口壓力脈動頻域圖Fig.11 Frequency domain of the monitoring point at the outlet of the volute

4 結(jié) 論

本文通過對余熱排出泵水力樣機的內(nèi)部流動特性進行數(shù)值模擬和實驗研究得到以下結(jié)論：

(1)在大流量和設(shè)計流量工況下，葉輪和導(dǎo)葉內(nèi)部流動比較穩(wěn)定，小流量工況下，葉輪和導(dǎo)葉流道內(nèi)出現(xiàn)了失速旋渦和流動分離等不穩(wěn)定流動現(xiàn)象，隨著流量的減小，其不穩(wěn)定程度加劇，嚴(yán)重影響泵的性能。

(2)各流量工況下，葉輪出口壓力脈動主要由轉(zhuǎn)頻決定，為7fz，受導(dǎo)葉葉片數(shù)影響；蝸殼和導(dǎo)葉出口壓力脈動主頻為5fz，主要由葉頻決定。

(3)葉輪和導(dǎo)葉出口壓力脈動系數(shù)幅值隨流量的減小急劇增加，應(yīng)避免泵在小流量工況下運行。

(4)壓力脈動的計算結(jié)果與實驗結(jié)果在基本吻合，證明本文的數(shù)值計算具有一定的可靠性。

[1] 洪鋒，袁建平，張金鳳，等．余熱排出泵小破口失水事故空化特性數(shù)值分析[J]．哈爾濱工程大學(xué)學(xué)報，2015，36(3)：297-301． HONG Feng，YUAN Jianping，ZHANG Jinfeng，et al．Numerical analysis of cavitating flow characteristics in residual heat removal pumps during the SBLOCA[J]．Journal of Harbin Engineering University，2015，36(3)：297-301．

[2] DUCOIN A J，ASTOLFI A，DENISET F，et al．Computational and experimental investigation of flow over a transient pitching hydrofoil[J]．European Journal of Mechanics B，2009，28(6)：728-743．

[3] 袁建平，朱鈺雯，周邦倫，等．余熱排出泵內(nèi)部壓力脈動特性分析[J]．排灌機械工程學(xué)報，2015，33(6)：475-480． YUAN Jianping，ZHU Yuwen，ZHOU Banglun，et al． Analysis on pressure fluctuation in residual heat removal pumps[J]．Journal of Drainage and Irrigation Machinery Engineering，2015，33(6)：475-480．

[4] 朱榮生，蘇保穩(wěn)，楊愛玲，等．離心泵壓力脈動特性分析[J]．農(nóng)業(yè)機械學(xué)報，2010，41(11)：43-47． ZHU Rongsheng，SU Baowen，YANG Ailing，et al．Numerical investigation of unsteady pressure fluctuation in centrifugal pump[J]．Journal of Agricultural Machinery，2010,41(11):43-47．

[5] KHALIFA A E，Al-QUTUB A M，BEN-MANSOUR R．Study of pressure fluctuations and induced vibration at blade-passing frequencies of a double volute pump[J]．Arabian Journal for Science and Engineering，2011，36(7)：1333-1345．

[6] 張人會，張學(xué)靜，楊軍虎．非設(shè)計工況下葉輪進口附近的流動及其控制[J]．甘肅工業(yè)大學(xué)學(xué)報，2003，29(4)：64-66． ZHANG Renhui，ZHANG Xuejing，YANG Junhu．Flow field around intake of impeller under off-design conditions and its control[J]．Journal of Gansu University of Technology，2003，29(4)：64-66．

[7] JOHNSON D A，PEDERSEN N，JACOBSEN C B．Measurements of rotating stall inside a centrifugal pump impeller[C]// The Fifth International Symposium on Pumping Machinery. Houston，TX，2005:1281-1288．

[8] 張睿，陳紅勛．改善失速工況下軸流泵水力性能的研究[J]．水力發(fā)電學(xué)報，2014，33(3)：292-298． ZHANG Rui，CHEN Hongxun．Study on the improvement of hydrodynamic performance of hydrodynamic performance of axial-flow pump at stall condition[J]．Journal of Hydroelectric Enginnering,2014，33(3)：292-298．

[9] 謝龍漢，趙新宇，張炯明．ANSYS CFX流體分析及仿真[M]．北京：電子工業(yè)出版社，2012．

[10] CORPORATION H P．Numerical simulation on the effect of turbulence models on impingement cooling of double chamber model[J]．Mathematical Problems in Engineering，2013，70(3)：1005-1021．

[11] 張德勝，吳蘇青，施衛(wèi)東，等．不同湍流模型在軸流泵葉頂泄漏渦模擬中的應(yīng)用與驗證[J]．農(nóng)業(yè)工程學(xué)報，2013，29(13)：46-53． ZHANG Desheng，WU Suqing，SHI Weidong，et al．Application and experiment of different turbulence models for simulating tip leakage vortex in axial flow pump[J]．Transactions of the Chinese Society of Agricultural Engineering，2013，29(13)：46-53．

[12] 施衛(wèi)東，徐磊，王川，等．蝸殼式離心泵內(nèi)部非定常數(shù)值計算與分析[J]．農(nóng)業(yè)機械學(xué)報，2014，45(3)：49-53． SHI Weidong，XU Lei，WANG Chuan，et al．Numerical calculation and analysis of inner unsteady flow for volute centrifugal pump[J]．Transactions of the Chinese Society for Agricultural Machinery，2013，29(13)：46-53．

[13] 劉厚林，劉明明，白羽，等．基于GCI的離心泵網(wǎng)格的收斂性[J]．江蘇大學(xué)學(xué)報(自然科學(xué)版)，2014，35(3)：279-283． LIU Houlin，LIU Mingming，BAI Yu，et al．Grid convergence based on GCI for centrifugal pump[J]．Journal of Jiangsu University(Natural Science Edition)，2014，35(3)：279-283．

[14] LUCIUS A， BRENNER G．Numerical simulation and evaluation of velocity fluctuations during rotating stall of a centrifugal pump[J]．Journal of Fluids Engineering，2011，133(8)：081101-081108．

[15] SPENCE R，AMARAL-TEIXEIRA J．Investigation into pressure pulsations in a centrifugal pump using numerical methods supported by industrial tests [J]．Comput Fluids，2008，37(6)：690-704．

Numerical and experimental study on flow instabilities in residual heat removal pumps

LU Rong, YUAN Jianping, ZHU Yuwen, FU Yanxia, ZHOU Banglun

(National Research Center of Pumps，Jiangsu University，Zhenjiang 212013，China)

In order to investigate the flow characteristics of residual heat removal pumps under multi conditions，three-dimensional unsteady numerical simulations were conducted by using the SST turbulence model in ANSYS CFX．The performance, inner flow fields, and pressure fluctuation of the residual heat removal pumps under multi conditions was obtained and compared with experimental results．The results show that：the inner flow keep regular under 1.2Qdand 1.0Qd．However，with the decrease of the flow rate，the stall vortex generates near the suction side at the impeller and diffuser inlet associated with the onset of varying degrees of flow separation in the flow channels，which is intensified toward the outlet．The main frequency of pressure fluctuation at the outlet of the impeller affected by the number of vanes is 7fz．The main frequencies at the outlet of both the diffuser and volute is 5fz，which are the blade passing frequency．The amplitudes of fluctuating pressure coefficient in all the monitoring points reach its minimum value under design condition，which increase gradually with the decrease of flow rate．It is concluded that the flow instabilities inside the pump occurs under small flow rate condition．

residual heat removal pump；pressure fluctuation；flow separation；numerical simulation；experiment

國家科技支撐計劃項目(2011BAF14B04)；國家自然科學(xué)基金資助項目(51349004)；江蘇高校自然科學(xué)研究項目(09KJB570001)

2015-07-14 修改稿收到日期：2015-12-01

陸榮 男，碩士生，1990年生

袁建平 男，研究員，博士生導(dǎo)師，1970年生

TH212；TH213.3

A

10.13465/j.cnki.jvs.2016.24.006