高溫高壓屏蔽電泵靜態(tài)阻力系數(shù)模擬與試驗(yàn)分析

杜洪偉，何蘭香

(佳木斯電機(jī)股份有限公司，黑龍江佳木斯154002)

0 引言

屏蔽電泵作為一種低噪聲、完全無(wú)泄漏泵，在核電和化工領(lǐng)域應(yīng)用極其廣泛。高溫高壓屏蔽電泵是屏蔽電泵的一種特殊產(chǎn)品，設(shè)計(jì)溫度350℃，設(shè)計(jì)壓力21MPa 以上。泵部分由泵體、導(dǎo)葉、葉輪組成。為了提高泵的水力效率，改善泵的性能，需要相應(yīng)的減小泵內(nèi)的水力損失。高溫高壓屏蔽電泵內(nèi)液體流動(dòng)復(fù)雜，本文以75-75-160扭曲葉片式高溫高壓屏蔽泵進(jìn)行三維流場(chǎng)阻力系數(shù)模擬與試驗(yàn)進(jìn)行對(duì)比分析。

1 基本理論

基本方程是計(jì)算機(jī)模擬計(jì)算流動(dòng)特性的基本理論和出發(fā)點(diǎn)。高溫高壓屏蔽電泵的流動(dòng)遵守質(zhì)量守恒、動(dòng)量守恒和能量守恒定律，其內(nèi)部流動(dòng)是三維、粘性、非定常的極其復(fù)雜的流動(dòng)。

1.1 流動(dòng)控制方程

假定該泵運(yùn)轉(zhuǎn)時(shí)轉(zhuǎn)速恒定，流動(dòng)控制方程為笛卡兒坐標(biāo)系下的Reynolds 時(shí)均Navier-Stockes方程，具體形式如下。

連續(xù)方程為

式中，u、μt—層流和紊流粘性系數(shù)，可以看出在Boussinesq 假設(shè)中μt被假定為各向同性;k—紊流動(dòng)能;δij—Kronecker 算子。μt和k 通過(guò)求解紊流模型方程得到。

1.2 湍流模型

采用在工程應(yīng)用中最為廣泛的兩方程模型，最基本的兩方程模型是標(biāo)準(zhǔn)的k-ε 模型進(jìn)行計(jì)算，即分別引入了k 湍動(dòng)能和湍動(dòng)能耗散率ε 的方程。其方程具體形式如下

2 靜態(tài)阻力系數(shù)模擬方法

2.1 計(jì)算區(qū)域

設(shè)計(jì)工況下泵的運(yùn)行參數(shù)見(jiàn)表1。

表1 泵的運(yùn)行參數(shù)表

計(jì)算區(qū)域包括葉輪、導(dǎo)葉和泵體的全部流道。本文采用Solidworks 三維軟件建模，通過(guò)該軟件中的布爾運(yùn)算，得到全流道實(shí)體模型如圖1 所示。

圖1 全流道三維模型圖

2.2 網(wǎng)格劃分

將水力模型以.stl 格式導(dǎo)入PumpLinx，由于計(jì)算流體域復(fù)雜，分別對(duì)葉輪、蝸殼(包括導(dǎo)葉)和入口段進(jìn)行網(wǎng)格劃分，采用專有的幾何等角自適應(yīng)二元樹(shù)(geometry Conformal Adaptative Binarytree)算法，即CAB 算法，在由封閉表面構(gòu)成的體域生成迪卡爾六面體網(wǎng)格。在靠近幾何邊界，CAB 自動(dòng)調(diào)整網(wǎng)格來(lái)適應(yīng)幾何曲面和幾何邊界線。利用最小的網(wǎng)格分辨細(xì)節(jié)特征的方法，對(duì)旋轉(zhuǎn)部件(葉輪)表面進(jìn)行了局部網(wǎng)格加密，總體網(wǎng)格劃分見(jiàn)圖2 所示。

圖2 流道網(wǎng)格劃分模型圖

2.3 計(jì)算方法

進(jìn)口邊界條件采用總壓進(jìn)口，假設(shè)在進(jìn)口截面上壓力均勻分布，通過(guò)調(diào)節(jié)進(jìn)口總壓來(lái)控制泵內(nèi)部汽蝕發(fā)生的程度;出口條件給定出口流量;壁面粗糙度設(shè)置為10μm;湍流模型為標(biāo)準(zhǔn)k-ε 模型;在葉輪、蝸殼和入口之間設(shè)置交互面，收斂精度設(shè)置為1×10-3，進(jìn)口處水的體積分?jǐn)?shù)設(shè)置為1，氣體體積分?jǐn)?shù)設(shè)置為0;計(jì)算格式:速度一、二階混合迎風(fēng)格式，其他方程采用一階迎風(fēng)格式。本文主要通過(guò)PumpLinx 軟件模擬泵運(yùn)轉(zhuǎn)時(shí)的三維流場(chǎng)情況。

2.4 靜態(tài)阻力系數(shù)計(jì)算公式

阻力系數(shù)的公式

式中，ΔP—泵進(jìn)出口的壓差，MPa;ρ—密度，kg/m3;v—出口流速，m/s;Q—流量，m3/h;D2—泵出口直徑，設(shè)計(jì)值0.075m;ξ—阻力系數(shù);Re—雷諾數(shù);ν—介質(zhì)粘度，m2/s，ν水=1×10-6m2/s。

3 模擬結(jié)果及分析

3.1 計(jì)算收斂結(jié)果

計(jì)算結(jié)果中流速、壓力、湍流動(dòng)能、湍流耗散率和蒸汽質(zhì)量分?jǐn)?shù)均具有良好的收斂性，并且收斂精度小于10-4。

3.2 內(nèi)部流動(dòng)的特性圖

決定內(nèi)部流動(dòng)特性的主要因素是壓力和速度，因此本文僅以靜態(tài)阻力系數(shù)模擬過(guò)程中額定流量點(diǎn)40m3/h 為例，模擬得出相應(yīng)的總壓、速度流線分布(詳見(jiàn)圖3 和圖4)。

圖3 總壓流線圖

圖4 速度流線圖

由圖3 可見(jiàn)，水力模型在葉輪出口處與導(dǎo)葉進(jìn)口處壓力變化較大，并且導(dǎo)葉出口處與泵體內(nèi)腔銜接處變化突出;由圖4 可見(jiàn)葉輪出口處與導(dǎo)葉進(jìn)口處速度驟變，其他流體區(qū)域速度變化相對(duì)均勻。以上兩圖可直觀看出，葉輪出口與導(dǎo)葉入口處流場(chǎng)變化較大，效率損失嚴(yán)重，在后續(xù)的改進(jìn)設(shè)計(jì)中應(yīng)重點(diǎn)考慮葉輪與導(dǎo)葉的匹配設(shè)計(jì)，從而提高水力模型效率。

4 阻力系數(shù)計(jì)算結(jié)果

根據(jù)上述模擬得出不同流量點(diǎn)的出、入口壓差(即揚(yáng)程)，計(jì)算靜態(tài)阻力系數(shù)，參數(shù)見(jiàn)表2。通過(guò)計(jì)算結(jié)果繪制流量、雷諾數(shù)與阻力系數(shù)曲線見(jiàn)圖5、圖6。

表2 阻力系數(shù)表(計(jì)算值)

圖5 流量-阻力系數(shù)曲線圖(計(jì)算值)

圖6 雷諾數(shù)-阻力系數(shù)曲線圖(計(jì)算值)

5 阻力系數(shù)試驗(yàn)結(jié)果

對(duì)75-75-160 模型泵進(jìn)行阻力特性試驗(yàn)，在電機(jī)不供電的冷態(tài)停車狀態(tài)下，電泵過(guò)流流量為正向(從進(jìn)口到出口)，測(cè)定水力部件過(guò)流流量與泵進(jìn)出口壓差(水頭損失)的關(guān)系，最終計(jì)算電泵在不同流量下的阻力系數(shù)。通過(guò)試驗(yàn)得到的結(jié)果見(jiàn)表3，并繪制流量、雷諾數(shù)與阻力系數(shù)曲線見(jiàn)圖7、圖8。

表3 阻力系數(shù)表(試驗(yàn)值)

圖7 流量-阻力系數(shù)曲線圖(試驗(yàn)值)

圖8 雷諾數(shù)-阻力系數(shù)曲線圖(試驗(yàn)值)

6 結(jié)語(yǔ)

本文所建立的水力模型是完全按照?qǐng)D紙?jiān)O(shè)計(jì)的，但在實(shí)際生產(chǎn)過(guò)程中水力部件的加工(如鑄造、鍛壓、車削等)對(duì)水力性能都會(huì)有一定的影響。因此，模擬的結(jié)果是理想狀態(tài)下水力模型內(nèi)部流場(chǎng)反映的靜態(tài)阻力和阻力系數(shù)的大小及變化趨勢(shì)。本次計(jì)算的理論計(jì)算值和試驗(yàn)數(shù)據(jù)相差不大，變化趨勢(shì)相似，模擬的結(jié)果是合理的。采用PumpLinx 軟件，基于Reynolds 時(shí)均N-S 方程和RNG k-ε 湍流模型，進(jìn)行三維流場(chǎng)數(shù)值模擬，分析流體在離心泵內(nèi)的速度和壓力分布，揭示其內(nèi)部的流動(dòng)規(guī)律，為產(chǎn)品的設(shè)計(jì)或改進(jìn)優(yōu)化設(shè)計(jì)提供重要依據(jù)。改善離心泵的水力性能，提高效率是離心泵優(yōu)化設(shè)計(jì)的重要研究?jī)?nèi)容。通過(guò)對(duì)比模擬結(jié)果與試驗(yàn)結(jié)果，流場(chǎng)分析結(jié)果可以有效預(yù)測(cè)性能變化趨勢(shì)，利用三維流場(chǎng)分析可以替代那些不可能實(shí)現(xiàn)的或者實(shí)現(xiàn)起來(lái)相當(dāng)復(fù)雜或費(fèi)時(shí)費(fèi)力的試驗(yàn)，不僅節(jié)約試驗(yàn)費(fèi)用，節(jié)約時(shí)間，而且可以定性地分析內(nèi)部流場(chǎng)的變化規(guī)律。

［1］ 關(guān)醒凡.現(xiàn)代泵設(shè)計(jì)手冊(cè)［M］.北京:宇航出版社，1995.9.(2).

［2］ 查森.離心式和軸流式水泵［M］.北京:中國(guó)工業(yè)出版社，1961.(12).

［3］ 施衛(wèi)東.流體力學(xué)教程［M］.四川:西南交通大學(xué)出版社，1996.(1).

［4］ Kurokawa J，Jing J.·Performance Prediction of CentrifugalImpeller and Scale Effects［A］.Proc.IAHR-Beijing 89Symposium，1989.(5).

［5］ 施衛(wèi)東.流體機(jī)械［M］.四川:西南交通大學(xué)出版社，1996.(1).