為了更準(zhǔn)確地分析軸流泵的壓力脈動特性，基于結(jié)構(gòu)化網(wǎng)格，在考慮葉頂間隙的情況下，以時均N-S方程為基本控制方程，采用標(biāo)準(zhǔn)k-ε雙方程湍流模型和SIMPLEC算法，對軸流泵模型進(jìn)行了不同的采樣頻率和采樣時間下的全流場非定常數(shù)值計算，并分析其對計算結(jié)果的影響.

　　結(jié)果表明：在進(jìn)行壓力脈動分析時，若取靜壓作為傅里葉變換的縱坐標(biāo)，頻域圖中頻率為0時的物理意義即該點平均靜壓值；采樣頻率和采樣時間對軸流泵壓力脈動時域特性和頻域特性具有重要影響，指出了因采樣頻率過小而導(dǎo)致的混疊現(xiàn)象在時域和頻域中的表現(xiàn)；并推導(dǎo)出當(dāng)保證采樣時間與葉輪周期和葉片數(shù)乘積的比值達(dá)到整數(shù)倍關(guān)系時，基頻處不會發(fā)生頻譜泄漏，并建立了普遍適用軸流泵壓力脈動數(shù)值計算采樣時間的選取公式，進(jìn)而從公式得出針對不同轉(zhuǎn)速的軸流泵，應(yīng)選擇不同的采樣時間.結(jié)果為進(jìn)一步分析軸流泵壓力脈動提供了一定參考.

　　軸流泵流量大、揚程低，在農(nóng)田灌溉、大型調(diào)水工程、市政給排水、污水處理以及電廠循環(huán)水工程等領(lǐng)域被廣泛應(yīng)用，真空技術(shù)網(wǎng)(http://www.mp99x.cn/)認(rèn)為其運行穩(wěn)定性是整個機(jī)組正常工作的必要條件之一.壓力脈動對軸流泵的運行穩(wěn)定性有重要影響，如果流道內(nèi)壓力脈動的主要頻率與系統(tǒng)固有頻率相同或相近，將引起共振現(xiàn)象，使整個機(jī)組產(chǎn)生強(qiáng)烈振動，長時間則會導(dǎo)致葉片疲勞損傷或軸承損壞，嚴(yán)重影響機(jī)組的安全運行.因此，準(zhǔn)確描述壓力脈動規(guī)律對避免系統(tǒng)共振具有重要的指導(dǎo)意義.

　　盡管國內(nèi)外學(xué)者對軸流泵的壓力脈動特性做過相關(guān)研究，但主要集中于流道內(nèi)不同位置的壓力脈動規(guī)律，無論是數(shù)值模擬還是在實際測中，采樣頻率和采樣時間的選取都較大程度依賴前人的經(jīng)驗.王福軍等對不同采樣時間下的軸流泵壓力脈動頻域圖進(jìn)行了對比研究，提出8T時得到的計算結(jié)果及計算時間是可以被工程上所接受的，但沒有就采樣頻率的取值進(jìn)行研究.姚志峰等對離心泵壓力脈動采樣頻率及采樣時間的選取進(jìn)行了分析，得到采樣頻率主導(dǎo)頻域圖的頻率范圍及采樣時間主導(dǎo)頻域圖的頻率分辨率，并推導(dǎo)出離心泵壓力脈動測試采樣頻率和采樣時間的公式，但該公式是建立在頻率分辨率給定的前提下.Kaupert等[6]采用6000Hz的采樣頻率和270T的采樣時間對離心泵進(jìn)行壓力脈動測試研究.Shi等在計算中分別對時間步長為旋轉(zhuǎn)周期的1/290和1/580進(jìn)行了對比研究，認(rèn)為時間步長選取對壓力脈動沒有明顯的影響.

　　文中以南水北調(diào)工程天津同臺試驗的軸流泵模型(模型編號為TJ-ZL-02，比轉(zhuǎn)數(shù)ns=700)為研究對象，基于結(jié)構(gòu)化網(wǎng)格，對其進(jìn)行全流場非定常數(shù)值計算，重點分析采樣頻率和采樣時間對軸流泵壓力脈動特性時序和頻域分析的影響，并對采樣時間的計算公式進(jìn)行推導(dǎo).

1、計算模型

　　計算模型包括進(jìn)口段、轉(zhuǎn)輪段、導(dǎo)葉段以及出口段，其主要幾何參數(shù)分別為葉輪直徑D2=300mm，葉片數(shù)Z=4，導(dǎo)葉葉片數(shù)Zg=7；性能參數(shù)分別為設(shè)計流量Qd=367.50L/s，設(shè)計揚程Hd=7.262m，轉(zhuǎn)速n=1450r/min.計算中為了更接近實際模型，取葉頂間隙大小與試驗泵相同，為0.5mm.考慮到網(wǎng)格質(zhì)量對計算精度有較大影響，采用正六面體網(wǎng)格對計算模型進(jìn)行劃分，通過O型網(wǎng)格控制葉片表面邊界層和其附近網(wǎng)格質(zhì)量，并對葉輪區(qū)域進(jìn)行加密處理，最終網(wǎng)格數(shù)量約為1000000，模型網(wǎng)格如圖1所示.

圖1 模型網(wǎng)格

　　計算以時均N-S方程為基本控制方程，采用標(biāo)準(zhǔn)k-ε雙方程湍流模型，進(jìn)口邊界條件設(shè)為velocityinlet，出口邊界條件設(shè)為opening.首先進(jìn)行定常計算，得到不同工況下模型泵的外特性曲線，并與試驗結(jié)果進(jìn)行對比，如圖2所示.可以看出，數(shù)值計算得到的外特性曲線與試驗曲線總體變化趨勢一致，這表明所采用的計算模型及方法是可靠的.

圖2 泵特性曲線

　　為研究采樣頻率和采樣時間對壓力脈動特性分析的影響，在葉輪前部及后部各選取1個截面，每個截面布置2排監(jiān)測點，分別為P11-P15，P21-P25.選取設(shè)計工況下，r=130mm圓周面上的監(jiān)測點P12和P22進(jìn)行分析，其位置如圖3所示.

圖3 監(jiān)測點位置

　　結(jié)論

　　通過對不同采樣頻率及采樣時間下的軸流泵模型進(jìn)行非定常數(shù)值計算與分析，得到結(jié)論如下：

　　1)在壓力脈動的頻域圖上，若取監(jiān)測點壓力作為縱坐標(biāo)，則在頻率為0處，監(jiān)測點的幅值實質(zhì)上是該測點靜壓的平均值.

　　2)從時域圖和頻域圖2個方面分析了采樣頻率對壓力脈動特性分析的影響，并指出了混疊現(xiàn)象在時域和頻域中的表現(xiàn).

　　3)推導(dǎo)了在保證采樣時間與葉輪周期和葉片數(shù)之比達(dá)到整數(shù)倍關(guān)系時，在葉頻及其倍頻處不會出現(xiàn)頻譜泄漏，即對于不同轉(zhuǎn)速的軸流泵，不能給定同一頻率分辨率.