劉志超　孔繁余　王洋　謝山峰　趙立峰

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

摘要：為了減小旋渦自吸泵壓力脈動幅值及泵運(yùn)行時的噪聲，采用非等距葉片的分布方式，設(shè)計了3種葉片分布葉輪，分別對旋渦自吸泵進(jìn)行穩(wěn)態(tài)性能和壓力脈動分析。利用RNG k－ε模型求解旋渦自吸泵內(nèi)部非定常流動和性能，并獲得旋渦自吸泵內(nèi)部徑向和軸向間隙處壓力脈動數(shù)據(jù)。數(shù)值結(jié)果表明：非等距葉片分布會對旋渦自吸泵的性能產(chǎn)生一定的影響。當(dāng)調(diào)制角較小時，其對旋渦自吸泵的性能影響較??；總體上，調(diào)制角對自吸泵效率影響不大。與等距葉片相比，非等距葉片的壓力分布較小。另外，非等距葉片分布能夠有效地減小旋渦自吸泵的壓力脈動幅值和改變脈動的主頻，同時調(diào)制角的變化會影響壓力脈動特性，當(dāng)調(diào)制角為5°時，主頻幅值下降幅度最大?？刂普{(diào)制角在合適范圍內(nèi)可以獲得更好的壓力脈動特性。

關(guān)鍵詞：旋渦自吸泵；壓力脈動；非等距葉片；數(shù)值分析

Effect of unequal spacing blade distribution on pressure fluctuation of self-priming vortex pump

LIU Zhichao，KONG Fanyu，WANG Yang，XIE Shanfeng，ZHAO Lifeng

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

Abstract:In order to reduce pressure pulsation amplitude and the noise of vortex pump when the pump is running， three kinds of impeller blades are designed based on the unequal spacing blade distribution，studying the steady-state performance and pressure fluctuation of self-priming vortex pump． Internal radial and axial gap pressure data of self-priming vortex pump are obtained by computing vortex

self-priming unsteady flow and the performance of the pump with the RNG k－ε model． The numerical results show that the unequal spacing blade distribution will have a certain effect on the performance of self-priming vortex pump． When the modulation angle is small，it has little effect on the efficiency of the self suction pump． Compared with the equal spacing blades，the pressure distribution of the unequal spacing blades are much smaller． In addition，the unequal spacing blade can effectively reduce pressure pulsation amplitude of the self-priming vortex pump and change pulsating frequency distribution，and the changes of modulation angle will affect the pressure fluctuation amplitude．When the modulation angle is 5°，pressure fluctuation amplitude of main frequency declines a lot． The better pressure fluctuation characteristics can be obtained by controlling the modulation angle in the proper range．

Key words：self-priming vortex pump；pressure fluctuation；unequal spacing blade；numerical analysis

旋渦泵是一種小流量、高揚(yáng)程、比轉(zhuǎn)數(shù)很小的泵，廣泛應(yīng)用于工農(nóng)業(yè)、醫(yī)藥和航空航天領(lǐng)域。近年來由于低比轉(zhuǎn)數(shù)泵的需求日益增大，對旋渦泵的研究逐漸增多。ENGEDA等^［1］提出了一個數(shù)學(xué)模型來分析徑向葉片旋渦泵復(fù)雜的內(nèi)部流動。MEAKLMIL等^［2］利用CFD方法提出了一種開式葉輪扭曲葉片的方法來改善旋渦泵性能。QUAIL等^［3］對軸向進(jìn)口徑向出口的旋渦泵進(jìn)行了三維數(shù)值模擬并與一維理論計算進(jìn)行對比，提出了一種旋渦泵新的設(shè)計思路。CHOI等^［4］通過試驗研究了葉輪葉片角度和形狀對旋渦泵性能的影響。王洋等^［5］通過數(shù)值模擬的方法討論了葉輪葉片位置對旋渦泵性能的影響。上述研究主要集中于旋渦泵穩(wěn)態(tài)性能，而在實際應(yīng)用中，對其運(yùn)行中的振動和噪聲問題也頗為關(guān)注。

壓力脈動主要分為渦流脈動、葉頻倍頻脈動、軸頻倍頻脈動。目前國內(nèi)外關(guān)于旋渦泵的壓力脈動研究較少，但對于離心泵和混流泵的壓力脈動，國內(nèi)外學(xué)者已做了大量的研究工作，可以參考相關(guān)的研究方法．馮建軍等^［6］通過數(shù)值模擬研究了徑向?qū)~泵葉輪和導(dǎo)葉動靜干涉現(xiàn)象，研究不同參數(shù)如葉片數(shù)和徑向間隙等對壓力脈動幅值影響。黎義斌等^［7］采用RNG k －ε湍流模型和滑移網(wǎng)格技術(shù)，計算了混流泵的葉輪進(jìn)口截面、葉輪與導(dǎo)葉動靜耦合面、導(dǎo)葉出口截面的壓力脈動，進(jìn)行快速傅里葉變換，得到了不同特征截面的壓力脈動的頻率和幅值。以上相關(guān)研究表明，采用數(shù)值模擬的方法能夠比較準(zhǔn)確地預(yù)測泵內(nèi)壓力脈動特征。

為了尋找一種減小自吸泵壓力脈動幅值和噪聲的設(shè)計方案，參考大量文獻(xiàn)發(fā)現(xiàn)，在風(fēng)機(jī)中，為了降低離散噪聲的基頻峰值，國內(nèi)外學(xué)者提出采用不等節(jié)距葉片將基頻處的噪聲峰值分散到其他頻率上，從而降低總的噪聲級^［8］。根據(jù)LIGHTHILL^［9］的氣動聲學(xué)理論，EWALD等^［10］和WRIGHT^［11］從不同角度導(dǎo)出了等節(jié)距葉片轉(zhuǎn)子的聲輻射公式。而國內(nèi)學(xué)者孫曉峰［12］又在WRIGHT的研究基礎(chǔ)上，推導(dǎo)出了不等節(jié)距葉片轉(zhuǎn)子噪聲的聲輻射公式，并通過試驗說明葉片可以在相當(dāng)大的裕度內(nèi)選擇分布方式，而不會使風(fēng)機(jī)的氣動性能惡化。馬健峰等^［13］根據(jù)文獻(xiàn)［12］葉片分組自平衡的分配方案進(jìn)行葉輪葉片的不等距設(shè)計，并進(jìn)行了數(shù)值模擬和試驗研究，證明采用不等距葉片可以降低離心通風(fēng)機(jī)基頻噪聲的峰值，但是總的噪聲級基本不變。

文中基于數(shù)值模擬，參考風(fēng)機(jī)中一種非等距葉片分布規(guī)律，設(shè)計3種非等距葉輪，研究非等距葉片對旋渦自吸泵壓力脈動的主頻和幅值的影響，從而為降低旋渦自吸泵在運(yùn)行中的脈動幅值提供一種新方法。

1　數(shù)值模型

1.1　幾何模型

文中研究的旋渦自吸泵基本參數(shù)為設(shè)計流量Q = 2 m³/h，揚(yáng)程H = 24 m，轉(zhuǎn)速n 為2 900 r/min，葉片數(shù)Z為40，葉輪直徑D為76 mm，葉輪寬度B為8 mm，兩側(cè)葉片交錯分布，徑向間隙0.10 mm，軸向間隙為0.15 mm。流體域結(jié)構(gòu)如圖1所示。

圖1　泵體和葉輪結(jié)構(gòu)

Fig.1　Self-priming vortex pump structure

1.2　葉片分布方式

文中所研究非等距葉片為葉片間周向夾角不相等。非等距葉片有許多種分布方式，如隨機(jī)分布和按一定規(guī)律分布。葉片隨機(jī)不等距分布，需要不斷試驗才能得到滿意的葉片分布，工作量很大。而按一定規(guī)律分布，只需要改變幾個參數(shù)就可以得到所需葉片分布方式。文中所研究的葉片分布方式為葉片分布角按正弦調(diào)制曲線分布，正弦調(diào)制函數(shù)為

式中:A為調(diào)制角，A值越大，表示葉片分布越不均勻；N為調(diào)制量的循環(huán)次數(shù)；φn為均勻分布時第n個葉片相對于起始位置角度；φ’_n為非等距分布時第n個葉片相對于起始位置角度。圖2所示為非等距葉片分布示意圖。

圖2　非等距葉片分布示意圖

Fig.2　Schematic diagram of unequal spacing blade distribution

文中所研究3種非等距葉片葉輪，3種不等距葉片調(diào)制角A分別為2°，5°和8°，N均為1，等距和3種非等距葉輪葉片周向夾角如圖3所示，圖中x為葉片數(shù)，y為葉距角。

圖3　葉片周向夾角分布

Fig.3　Blade peripheral intersection angle distribution

1.3　網(wǎng)格劃分

旋渦自吸泵的計算域分為進(jìn)口段、泵體、葉輪、出口段4部分，對進(jìn)口段、泵體、葉輪以及出口段進(jìn)行非結(jié)構(gòu)網(wǎng)格劃分，并對泵體的徑向和軸向間隙處薄壁結(jié)構(gòu)進(jìn)行網(wǎng)格加密，總網(wǎng)格數(shù)為1 853 895。泵體以及葉輪網(wǎng)格如圖4所示。

圖4　旋渦自吸泵泵體和葉輪網(wǎng)格模型

Fig.4　Mesh structure of self-priming vortex pump

2　數(shù)值計算方法和監(jiān)測點分布

2.1　數(shù)值計算方法

大渦模擬是對紊流脈動的一種空間平均，即通過某種濾波函數(shù)將大尺度的渦和小尺度的渦分離開，大尺度渦直接模擬，小尺度渦用模型封閉，它可以很好地模擬旋渦自吸泵的壓力脈動，但文中主要研究非等距葉片對旋渦自吸泵壓力脈動的影響，因此在非定常計算過程中選擇了計算時間短和占用計算資源少的RNG k－ε模型。在定常計算收斂基礎(chǔ)上進(jìn)行非定常計算，取標(biāo)準(zhǔn)壁面函數(shù)。計算時，給定進(jìn)口邊界條件為壓力，出口邊界條件為質(zhì)量流量。每個非定常計算時間步長為Δt = 5.747 126 4 ×10^-5s，即每步為葉輪旋轉(zhuǎn)1°，完成一個計算周期需要360個時間步長，總計算時間為t = 0.124 137 91 s，即計算6個周期。當(dāng)每個計算時間步長內(nèi)的迭代次數(shù)達(dá)10次或控制方程變量的絕對殘差均小于10^－4時，進(jìn)入下一個時間步的迭代。

2.2　監(jiān)測點分布

在旋渦自吸泵徑向間隙、軸向間隙、流道壁面設(shè)置31個監(jiān)測點，z= 0平面處，徑向間隙和流道壁面處設(shè)置13個監(jiān)測點；z= 2 mm平面處，只在徑向間隙處設(shè)置6個監(jiān)測點；在軸向間隙處設(shè)置12個監(jiān)測點。監(jiān)測點分布如圖5所示。

圖5　旋渦自吸泵監(jiān)測點分布圖

Fig.5　Monitoring points distribution of self-priming vortex pump

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

3.1　非等距葉片對穩(wěn)態(tài)性能的影響

定常計算時，邊界條件設(shè)置為進(jìn)口壓力，出口為質(zhì)量流量，參考壓力為101.325 kPa。先對不同A值的葉輪在設(shè)計點(Q=2 m³/h)進(jìn)行定常數(shù)值模擬，得到其在設(shè)計點的揚(yáng)程和效率；按照旋渦泵相關(guān)試驗標(biāo)準(zhǔn)，保證徑向間隙0.10 mm 和軸向間隙0.15 mm，泵體其他各個參數(shù)不變，葉輪為A = 0°等距葉片，其他葉輪參數(shù)不變，進(jìn)行試驗，得到在設(shè)計點(Q = 2 m³/h) 的揚(yáng)程為24.3 m和效率0.253 4。不同A值葉輪的旋渦自吸泵在設(shè)計點的數(shù)值模擬揚(yáng)程如表1所示。

表1　不同A值時葉輪模擬揚(yáng)程和效率

Tab.1　Simulation head and efficiency of impeller with different A values

由表1可知，在設(shè)計點，A = 0°的定常數(shù)值模擬揚(yáng)程與A = 0°的葉輪試驗揚(yáng)程相差不大，在相對誤差( 5％)范圍內(nèi)，模擬效率和試驗相差也在相對誤差范圍內(nèi)，數(shù)值模擬的結(jié)果有一定可信度。由表1可以看出，在A值較小時，非等距葉片對旋渦自吸泵的影響較??；在A值偏大時，對旋渦自吸泵的揚(yáng)程影響較大；調(diào)制角A對旋渦自吸泵的效率影響不大。因此，存在一個合理的A值，對旋渦自吸泵揚(yáng)程影響比較小，在一定范圍內(nèi)選擇A值不會明顯降低旋渦自吸泵揚(yáng)程；總體來說，調(diào)制角A對旋渦自吸泵的效率影響不大。

不同A值對旋渦自吸泵流道壓力影響如圖6所示。由圖6可知，泵體內(nèi)壓力最低的位置發(fā)生變化，在A = 0°時，泵體內(nèi)壓力最低的位置在靠近流道進(jìn)口處，而不等距葉片分布時，壓力最低處靠近流道進(jìn)口和整個進(jìn)口流道。A從0°到5°，進(jìn)口流道的壓力逐漸降低，而A = 8°時，進(jìn)口流道壓力比A = 5°時有所上升，但仍比A = 0°時要??；無論A = 2°，5°，8°，其他位置壓力都比A = 0°時要低，當(dāng)A從2°到8°變化時，其他位置壓力變化不大。非等距葉片對旋渦自吸泵的進(jìn)口流道壓力有影響，而且不同A值對泵體進(jìn)口流道處的壓力影響不同，而在其他位置，等距葉片分布要比非等距分布的壓力要大，非等距分布時，調(diào)制角A的變化對壓力影響不大。

圖6　不同A時流道壓力分布

Fig.6　Pressure distribution of flow passage with different A values

3.2　非等距葉片對非定常流動壓力脈動特性的影響

文中主要研究在設(shè)計點(Q = 2 m³/h)旋渦自吸泵的壓力脈動。以定常計算的結(jié)果作為初始條件，對不同A值的旋渦自吸泵進(jìn)行非定常計算，得到31個監(jiān)測點的壓力脈動數(shù)據(jù)，進(jìn)行傅里葉變換得到壓力脈動的頻域數(shù)據(jù)。

3.2.1　等距葉片壓力脈動特性

圖7為等距葉片壓力脈動頻域特性。由圖7可知，徑向間隙z= 0和z = 2 mm處的壓力脈動主要集中在葉片通過頻率以及倍頻，幅值最大的位置在靠近出口流道的監(jiān)測點P14處，監(jiān)測點從P14到P19和從P1到P6，脈動幅值均逐漸減小。軸向間隙處的壓力脈動也主要集中于葉片通過頻率以及倍頻處，幅值最大的位置也是在靠近出口流道點P20處，從點P20到P25脈動幅值逐漸減小，從點P31到P30脈動幅值增大，從點P30到P26，脈動幅值逐漸減小。在壁面流道z = 0處，壓力脈動也是主要集中在葉片通過頻率以及倍頻，但幅值遠(yuǎn)小于徑向和軸向間隙處的脈動幅值。以上脈動主要集中在葉片通過頻率及其倍頻處，主要是由于葉片等距分布，每個葉輪水體通過隔舌的時間相同，產(chǎn)生波動在葉片通過頻率及其倍頻處疊加，使這些頻率處脈動更強(qiáng)烈。在軸向間隙和徑向間隙處，壓力脈動幅值最大的位置都在靠近流道出口處。在徑向間隙z = 0處，葉頻處的脈動幅值比較小，主要是兩側(cè)葉片交錯分布，削弱了葉輪與泵體隔舌之間的動靜干涉。在z =2.0 mm 和z = 3.9 mm處，葉輪兩側(cè)交錯葉片分布對動靜干涉作用影響不大。

圖7　等距葉片壓力脈動頻域特性

Fig.7　Frequency spectrum characteristics of equal spacing blade

3.2.2　非等距葉片對壓力脈動的影響

由于旋渦自吸泵的徑向和軸向間隙很小，徑向和軸向間隙處由于動靜干涉引起的壓力脈動最劇烈。而且，由上文研究可知，徑向和軸向間隙處P14和P20的壓力幅值最大，下文主要研究監(jiān)測點P14和P20處脈動。

圖8為葉片非等距分布時，徑向和軸向間隙P14和P20處壓力脈動頻域圖。

圖8　不等距葉片壓力脈動頻域特性

Fig.8　Frequency spectrum characteristics of unequal spacing blade

由圖8可知，葉片非等距分布時，脈動主頻和幅值發(fā)生變化，脈動幅值不再集中于葉片通過頻率及其倍頻處，而是分散到其他頻率處，A值越大，幅值分散越明顯。以上現(xiàn)象主要是由于葉片非等距分布，每個葉輪水體通過隔舌的時間不再相同，從而使每個葉輪產(chǎn)生的脈動不再疊加，而是有一定的離散。在監(jiān)測點P14和P20處，當(dāng)調(diào)制角A = 2°時，脈動主頻還是葉片通過頻率(1 929.33 Hz)；在調(diào)制角A = 5°和8°時，脈動主頻不再是葉片通過頻率，分別為1 833.40 Hz和1 736.91 Hz；在監(jiān)測點P14和P20處，3種非等距葉片，都使壓力脈動幅值有所降低，而A = 5°時的主頻幅值下降最多，點P14和P20的壓力脈動幅值分別從99.3 kPa和84.5 kPa下降至54.3 kPa和47.6 kPa。由此可知，當(dāng)調(diào)制角A很小時，非等距葉片僅能使主頻的幅值降低，而主頻不變；當(dāng)調(diào)制角A比較大時，能使脈動的主頻和幅值下降。但并不是調(diào)制角A越大，主頻的幅值越小，要使徑向和軸向間隙處的脈動幅值最小，調(diào)制角A不宜過大。

由圖8可知，隨著調(diào)制角A的不斷增大，葉片通過頻率處及其倍頻處的脈動幅值不斷減小。因此，非等距分布葉片能使通過頻率以及倍頻處的脈動幅值得到調(diào)制，調(diào)制角A越大，對葉片通過頻率以及倍頻處的調(diào)制越明顯。但隨著調(diào)制角A的增大，其他頻率處的脈動幅值上升很快，特別是低于葉片通過頻率處的脈動幅值上升很快。為了控制總體的脈動幅值，應(yīng)該把調(diào)制角A控制在一定范圍。

4　結(jié)論

1)非等距葉片分布對旋渦自吸泵穩(wěn)態(tài)性能有一定影響，當(dāng)調(diào)制角A較小時，對旋渦自吸泵的性能影響不大，當(dāng)調(diào)制角A過大時，對性能影響比較大。在不影響旋渦自吸泵性能的情況下，可以在一定范圍內(nèi)選擇葉片的分布方式，A在0°～5°對泵的性能影響不大。

2)非等距葉片分布能使旋渦自吸泵的脈動主頻和幅值發(fā)生變化，脈動幅值不再集中于葉片通過頻率及其倍頻處，而是分散到其他頻率處。對葉片通過頻率以及倍頻處脈動有調(diào)制作用，A值越大，調(diào)制越明顯。A值在一定范圍內(nèi)，能使脈動幅值減小到最小，A = 5°時，脈動幅值最小。

參考文獻(xiàn)

［1］ENGEDA A，RAHEEL M M． Systematic design approach for radial blade regenerative turbomachines［J］． Journal of propulsion and power， 2005， 21( 5) : 884 － 892．

［2］MEAKHAIL T，PARK S O． An improved theory for regenerative pump performance［J］．Journal power and energy， 2005， 219( 3) : 213－222．

［3］QUAIL F J，SCANLON T，BAUMGARTNER A． Design study of a regenerative pump using one-dimensional and three dimensional numerical techniques［J］． European journal of mechanics B/fluids， 2012， 31: 181－187．

［4］CHOI W C，YOO I S，PARK M R，et al． Experimental study on the effect of blade angle on regenerative pump performance［J］．Journal of power and energy，2013， 227( 5) : 585－592．

［5］王洋，傅劍輝，蔣其松．閉式葉輪葉片位置對旋禍泵性能的影響［J］． 農(nóng)機(jī)化研究， 2010， 32( 9) : 150－154．

WANG Yang，F(xiàn)U Jianhui， JIANG Qisong． The influence of the vortex pump performance to the position of a closed impeller blade［J］． Journal of agricultural mechanization

research， 2010， 32( 9) : 150－154．( in Chinese)

［6］FENG J J，BENRA F，DOHMEN H J． Numerical investigation on pressure fluctuations for different configurations of vaned diffuser pumps［J］． International journal of rotating machinery，2007， 2007: 1－10．

［7］黎義斌，李仁年，王秀勇，等．低比轉(zhuǎn)數(shù)混流泵壓力脈動特性的數(shù)值模擬［J］．排灌機(jī)械工程學(xué)報， 2013， 31( 3) : 205－209．

LI Yibin，LI Rennian，WANG Xiuyong， et al． Numerical analysis of pressure fluctuation in low specific speed mixed-flow pump［J］． Journal of drainage and irrigation

machinery engineering，2013，31 ( 3 ) : 205－209．( in Chinese)

［8］智乃剛，蕭濱詩． 風(fēng)機(jī)噪聲控制技術(shù)［M］． 北京: 機(jī)械工業(yè)出版社，1985．

［9］LIGHTHILL M J． On sound generated aerodynamically.I.General theory［J］．Proceedings of the royal society A， 1952，211( 1107) : 564 － 587．

［10］EWALD D，PAVLOVIC A，BOLLINGER J G． Noise reduction by applying modulation principles［J］． Journal of the acoustical society of America， 1971， 49: 1385．

［11］WRIGHT S E．Sound radiation from a lifting rotor generated by asymmetric disk loading［J］． Journal of sound＆ vibration，1969，9( 2) : 223－226．

［12］孫曉峰．不等節(jié)距葉片風(fēng)機(jī)氣動聲學(xué)特性的研究［J］． 北京航空航天大學(xué)學(xué)報， 1986( 4) : 137 － 145．

SUN Xiaofeng． The aeroacoustic nature of unequally spaced fan［J］． Journal of Beijing Institute of Aeronautics and Astronautics， 1986( 4) : 137－145．( in Chinese)

［13］馬健峰，袁民建，劉秋洪，等． 不等距葉片離心風(fēng)機(jī)氣動噪聲的數(shù)值與實驗研究［J］． 噪聲與振動控制，2008( 3) : 100－103．

MA Jianfeng，YUAN Minjian，LIU Qiuhong，et al． Experimental and numerical study of aeroacoustics of a centrifugal fan with nonisometric blade impeller［J］． Noise and vibration control，2008 ( 3 ) : 100－103．( in Chinese)

基金項目: “十二五”國家科技支撐計劃項目( 2011BAF14B01)；江蘇高校優(yōu)勢學(xué)科建設(shè)工程項目

作者簡介：

劉志超( 1988—) ，男，河南盧氏人，碩士研究生( 2008bigsea@163.com) ，主要從事流體機(jī)械及工程研究

孔繁余( 1956—) ，男，江蘇揚(yáng)州人，研究員，博士生導(dǎo)師( kongfy291s@ sohu.com) ，主要從事石化、軍工等特種泵的設(shè)計、試驗研究