竇立波, 朱玉蘭, 陸紀(jì)根, 居杰, 王友宏

（江蘇亞梅泵業(yè)集團(tuán)有限公司，江蘇 靖江214500）

0 引 言

自吸式旋渦泵是一種廣泛應(yīng)用于石化、航空業(yè)、農(nóng)業(yè)和市政工程等領(lǐng)域的通用機(jī)械產(chǎn)品，具有揚(yáng)程高、流量小、結(jié)構(gòu)簡(jiǎn)單、質(zhì)量輕、體積小和能自吸等特點(diǎn)，其水力性能特征介乎于普通離心泵和容積泵之間，具有良好的穩(wěn)定性和可靠性[1]。自吸式旋渦泵主要由葉輪、泵體、氣水分離室、電磁閥和進(jìn)、出口管等組成，其中葉輪、泵體和儲(chǔ)水箱是影響其性能的核心關(guān)鍵水力部件，其葉輪分為開式葉輪和閉式葉輪2類，其泵體流道可分為開式流道、閉式流道和半開式流道3種。目前常用的自吸式旋渦泵結(jié)構(gòu)主要采用閉式葉輪和開式流道的結(jié)構(gòu)型式，其具有較好的自吸性能和效率，適用于工業(yè)、農(nóng)業(yè)等領(lǐng)域。在以往的研究中，主要采用經(jīng)驗(yàn)設(shè)計(jì)和試驗(yàn)方法進(jìn)行自吸式旋渦泵產(chǎn)品的設(shè)計(jì)，對(duì)于旋渦泵內(nèi)流場(chǎng)分布和泵體流道內(nèi)部增壓規(guī)律缺少量化研究，近幾年來，國(guó)內(nèi)外學(xué)者利用CFD(Computational Fluid Dynamics)數(shù)值仿真技術(shù)對(duì)旋渦泵進(jìn)行了研究，取得了一定的進(jìn)展[2]。

沙毅、唐濤等[3-4]采用經(jīng)驗(yàn)設(shè)計(jì)和試驗(yàn)方法開展了自吸式旋渦泵的效率分析、水力設(shè)計(jì)與試驗(yàn)研究，分析了閉式旋渦泵的水力設(shè)計(jì)模型和方法，并通過模型試驗(yàn)進(jìn)行了驗(yàn)證，對(duì)于葉輪和泵體模型的設(shè)計(jì)優(yōu)化提供了理論和試驗(yàn)基礎(chǔ)。朱祖超、王樂勤等[5-6]開展了小流量、高揚(yáng)程旋渦泵的設(shè)計(jì)理論與試驗(yàn)方法研究，建立了小流量旋渦泵的水力模型設(shè)計(jì)理論，并進(jìn)行了試驗(yàn)驗(yàn)證研究，為后續(xù)的特種旋渦泵產(chǎn)品開發(fā)提供了參考。施衛(wèi)東等[7]采用CFD數(shù)值模擬方法研究了不同流道截面形狀對(duì)于其性能影響的規(guī)律，并對(duì)于旋渦泵內(nèi)流場(chǎng)的分布特征進(jìn)行了比較研究。劉志超等[8]采用CFD非定常數(shù)值模擬方法研究了旋渦泵葉片結(jié)構(gòu)和分布對(duì)于其內(nèi)部壓力脈動(dòng)和能量分布的影響，分析了葉片分布對(duì)于旋渦泵噪聲和振動(dòng)特性的影響規(guī)律，并總結(jié)了旋渦泵的優(yōu)化設(shè)計(jì)方法。本文以某工業(yè)用自吸式旋渦泵產(chǎn)品為目標(biāo)，開展數(shù)值模擬和試驗(yàn)驗(yàn)證研究，分析了自吸式旋渦泵的水力性能和內(nèi)部流場(chǎng)分布規(guī)律，為后續(xù)的水力模型改進(jìn)設(shè)計(jì)提供有效參考。

1 設(shè)計(jì)參數(shù)和三維模型

以某公司的工業(yè)用自吸式旋渦泵為研究對(duì)象，其主要設(shè)計(jì)參數(shù)為：設(shè)計(jì)工況點(diǎn)的流量Q0= 8 m3/h，揚(yáng)程H0=24 m，轉(zhuǎn)速n=1450 r/min。其葉輪采用閉式葉輪，葉輪葉片數(shù)為80，泵體流道采用開式流道結(jié)構(gòu)，泵體流道型式為環(huán)形流道。采用三維造型軟件建造獲得的自吸式旋渦泵內(nèi)部流場(chǎng)三維模型（如圖1），其主要包括進(jìn)口管、出口管、氣水分離室、葉輪和泵體流道等部分，為了增強(qiáng)模型的穩(wěn)定性，三維模型的進(jìn)口管和出口管進(jìn)行了一定的延伸[9]。

2 數(shù)值模型

2.1 模型網(wǎng)格劃分

根據(jù)圖1可知，自吸式旋渦泵的內(nèi)流場(chǎng)三維模型形狀較為復(fù)雜，為了減少數(shù)值模擬的計(jì)算工作量和準(zhǔn)確捕捉其內(nèi)部流動(dòng)特征，采用混合網(wǎng)格結(jié)構(gòu)對(duì)該內(nèi)流場(chǎng)三維模型進(jìn)行網(wǎng)格劃分。利用專業(yè)網(wǎng)格軟件ICEM進(jìn)行該模型的網(wǎng)格劃分，其進(jìn)口管、出口管和葉輪部分采用結(jié)構(gòu)化六面體網(wǎng)格結(jié)構(gòu)，氣水分離室和泵體流道等部分采用適應(yīng)性較強(qiáng)的非結(jié)構(gòu)化四面體網(wǎng)格結(jié)構(gòu)，同時(shí)在葉輪和泵體流道交界區(qū)域進(jìn)行了網(wǎng)格局部加密處理。網(wǎng)格無關(guān)性分析后，確定該網(wǎng)格模型節(jié)點(diǎn)總數(shù)約為121萬(wàn)，此時(shí)的計(jì)算誤差相對(duì)很小，最終確定模型參數(shù)如表1所示[10]，網(wǎng)格模型如圖2所示。

圖1 自吸式旋渦泵內(nèi)流場(chǎng)三維模型

表1 網(wǎng)格模型參數(shù)

2.2 模型算法和邊界

按照以往的研究文獻(xiàn)和模型分析方法，采用RNG k-ε湍流模型對(duì)自吸式旋渦泵三維流場(chǎng)模型進(jìn)行求解，葉輪旋轉(zhuǎn)區(qū)域采用多重參考系（MRF）方法進(jìn)行求解，進(jìn)口管的進(jìn)口邊界采用無旋的速度進(jìn)口條件，出口管的出口邊界條件采用壓力出口條件，模型所有動(dòng)、靜壁面采用無滑移條件[11]。對(duì)流項(xiàng)和擴(kuò)散項(xiàng)分別采用二階迎風(fēng)差分和中心差分格式，收斂殘差值為10-4，最終確定數(shù)值模型算法和邊界條件如表2所示。

為準(zhǔn)確了解流體介質(zhì)在泵體流道內(nèi)的靜壓變化過程，沿著泵體流道內(nèi)的流體流動(dòng)方向，分別選擇如圖3所示的A、B、C、D、E共5個(gè)過流斷面作為監(jiān)測(cè)截面，對(duì)泵體流道內(nèi)部的流體靜壓增加過程進(jìn)行監(jiān)測(cè)和描述。

圖2 自吸式旋渦泵網(wǎng)格劃分

3 結(jié)果分析

3.1 水力性能驗(yàn)證

按照上述自吸式旋渦泵模型的技術(shù)參數(shù)，通過其數(shù)值模型計(jì)算的工況點(diǎn)包括0.1Q0、0.5Q0、1.0Q0（設(shè)計(jì)工況點(diǎn)）和1.5Q0共4個(gè)工況點(diǎn)，同時(shí)按照J(rèn)B/T 7743-2011《旋渦泵》和GB 3216-2005《回轉(zhuǎn)動(dòng)力泵水力性能驗(yàn)收試驗(yàn)》標(biāo)準(zhǔn)要求開展了自吸式旋渦泵的水力性能試驗(yàn)，并對(duì)其水力性能試驗(yàn)結(jié)果進(jìn)行了分析。通過對(duì)數(shù)值模擬和試驗(yàn)結(jié)果進(jìn)行對(duì)比分析，獲得自吸式旋渦泵的數(shù)值模擬和水力性能試驗(yàn)的揚(yáng)程曲線對(duì)比，如圖4所示。

表2 模型算法和邊界條件

圖3 泵體流道監(jiān)測(cè)截面

圖4 揚(yáng)程曲線對(duì)比

從圖4可知，自吸式旋渦泵的揚(yáng)程隨著流量的增大而近似呈線性降低，數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果趨勢(shì)一致，其中數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果在1.0Q0流量（設(shè)計(jì)工況點(diǎn)）附近時(shí)的相對(duì)誤差最小，隨著流量的增大和減小，其相對(duì)誤差逐漸增大。經(jīng)過計(jì)算可知，在設(shè)計(jì)工況點(diǎn)，數(shù)值模擬和試驗(yàn)結(jié)果的揚(yáng)程相對(duì)誤差約為3%，顯示了該數(shù)值模型計(jì)算的準(zhǔn)確性。

3.2 泵體流道內(nèi)部靜壓分析

隨著旋渦泵葉輪葉片的高速旋轉(zhuǎn)，泵體流道內(nèi)會(huì)產(chǎn)生周期性的沖擊、渦旋和剪切流動(dòng)等現(xiàn)象，流體在旋渦泵泵體流道內(nèi)不斷增壓，旋渦泵內(nèi)部流體流動(dòng)屬于非常復(fù)雜的三維空間湍流狀態(tài)，通過常規(guī)分析方法難以研究其內(nèi)部流動(dòng)規(guī)律。為了分析自吸式旋渦泵泵體流道內(nèi)的流體靜壓增壓過程和靜壓分布規(guī)律，根據(jù)上述數(shù)值模型計(jì)算結(jié)果，沿著流體流動(dòng)方向，對(duì)圖3所示的泵體流道監(jiān)測(cè)截面A、B、C、D、E進(jìn)行了靜壓分析，各監(jiān)測(cè)截面在不同工況點(diǎn)下的靜壓曲線如圖5所示（為便于分析，設(shè)定監(jiān)測(cè)截面A的靜壓值為零）。

圖5 監(jiān)測(cè)截面（過流斷面）靜壓曲線

從圖5可知，不同工況點(diǎn)條件下，泵體流道內(nèi)部沿著流動(dòng)方向的流體靜壓變化曲線基本一致，其總體近似呈線性增長(zhǎng)，從泵體流道的進(jìn)口監(jiān)測(cè)截面A開始，流體靜壓逐漸增大至泵體流道的出口監(jiān)測(cè)截面E。從圖5的靜壓曲線變化特征可知，在泵體流道的進(jìn)口和出口附近，即流體從監(jiān)測(cè)截面A→B和監(jiān)測(cè)截面D→E的流動(dòng)過程中，其靜壓增長(zhǎng)相對(duì)較為平緩，而流體從監(jiān)測(cè)截面B→D的流動(dòng)過程中，其靜壓增長(zhǎng)相對(duì)較快。

根據(jù)上述結(jié)果可知，流體在泵體流道內(nèi)部的靜壓分布規(guī)律基本一致，其總體近似呈線性增長(zhǎng)。通過詳細(xì)分析發(fā)現(xiàn)，相比較而言，流體在泵體流道的進(jìn)、出口附近增長(zhǎng)速率相對(duì)較低，其在泵體流道中間段（監(jiān)測(cè)截面B→D）增長(zhǎng)速率相對(duì)較高，該現(xiàn)象反映了泵體流道內(nèi)部流動(dòng)狀態(tài)的變化趨勢(shì)。

4 結(jié) 論

1）研究結(jié)果顯示，基于RNG k-ε湍流模型的自吸式旋渦泵數(shù)值模型可以較為準(zhǔn)確地預(yù)測(cè)旋渦泵的水力性能曲線和內(nèi)部流場(chǎng)特征，在設(shè)計(jì)工況點(diǎn)時(shí)，數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果在設(shè)計(jì)工況點(diǎn)附近最為接近，其揚(yáng)程相對(duì)誤差約為3%。

2）沿著流體流動(dòng)方向，自吸式旋渦泵泵體流道內(nèi)的靜壓近似呈線性增長(zhǎng)，不同工況下的泵體流道內(nèi)靜壓增長(zhǎng)曲線基本一致。

3）在泵體流道進(jìn)、出口流動(dòng)條件的影響下，流體在泵體流道中間段（監(jiān)測(cè)截面B-D）流動(dòng)過程中的靜壓增長(zhǎng)相對(duì)較快，表明在該流動(dòng)區(qū)域內(nèi)的葉輪與流體之間的動(dòng)能交換更為高效顯著，該結(jié)果可以為后續(xù)的泵體改進(jìn)提供有效參考。