離心泵壓力脈動(dòng)的特性分析和測(cè)試分析研究機(jī)械制造專業(yè)

1、摘要離心泵在產(chǎn)生靜態(tài)壓力分量的同時(shí)還產(chǎn)生動(dòng)態(tài)壓力分量，也就是壓力脈動(dòng)，壓力脈動(dòng)象交流信號(hào)一樣疊加在靜態(tài)壓力分量上。離心泵大多采用螺旋型蝸殼 ，蝸殼與隔舌的不對(duì)稱性造成葉輪出口流動(dòng)的環(huán)形畸變，尤其是在非設(shè)計(jì)工況點(diǎn)時(shí)，畸變更為劇烈，因此在設(shè)計(jì)主軸時(shí)，必須考慮由于壓力分布不平衡引起的脈動(dòng)特性。離心泵內(nèi)的壓力脈動(dòng)是極其復(fù)雜的，對(duì)同一臺(tái)泵，不同工況呈現(xiàn)不同性質(zhì)的壓力脈動(dòng)。泵內(nèi)部流場(chǎng)非定常的湍流特性、葉輪與進(jìn)水管的動(dòng)靜干擾、渦流及回流等都可能引起壓力脈動(dòng)，該壓力脈動(dòng)激發(fā)泵體和管道的機(jī)械振動(dòng)，影響泵的正常運(yùn)轉(zhuǎn)，降低了泵的使用壽命；如果這種激勵(lì)振動(dòng)的頻率接近泵的主頻，還會(huì)發(fā)生共振，振動(dòng)會(huì)在相鄰的環(huán)境介質(zhì)內(nèi)產(chǎn)生

2、聲壓波動(dòng)，產(chǎn)生很大的噪聲?？梢姡脙?nèi)壓力脈動(dòng)是影響水泵運(yùn)行穩(wěn)定性的重要因素，也是導(dǎo)致噪聲污染的重要原因。因而，對(duì)泵內(nèi)部壓力脈動(dòng)進(jìn)行研究，降低泵的振動(dòng)噪聲成為當(dāng)前研究的重要任務(wù)之一。本文首先應(yīng)用CFD（Computational Fluid Dynamics）技術(shù)Fluent軟件 ，在多工況下同時(shí)監(jiān)測(cè)葉輪流道與蝸殼內(nèi)壓力脈動(dòng)的情況，以揭示壓力脈動(dòng)與離心泵性能之間的關(guān)系，為進(jìn)一步預(yù)測(cè)離心泵非定常流動(dòng)誘發(fā)壓力脈動(dòng)提供依據(jù)。為了研究離心泵在不同工況下以及不同位置間壓力脈動(dòng)的變化規(guī)律及其關(guān)系，基于標(biāo)準(zhǔn) k-湍流模型和滑移網(wǎng)格模型，應(yīng)用Fluent軟件對(duì)離心泵內(nèi)部流場(chǎng)先后進(jìn)行定常與非定常三維模擬，得到大致

3、的內(nèi)部流場(chǎng)特性及監(jiān)測(cè)點(diǎn)的壓力脈動(dòng)情況。計(jì)算結(jié)果表明：在小流量工況下，由于內(nèi)部流場(chǎng)的不穩(wěn)定運(yùn)行各監(jiān)測(cè)點(diǎn)處的壓力脈動(dòng)都比較大且不均勻；在葉輪流道中，接近葉輪出口處的壓力脈動(dòng)明顯高于葉輪流道其他區(qū)域壓力脈動(dòng)。大流量工況下在蝸殼隔舌和出口處會(huì)產(chǎn)生一定回流導(dǎo)致蝸殼該處附近壓力脈動(dòng)表現(xiàn)不穩(wěn)定。蝸殼隔舌處由于該位置流域比較復(fù)雜壓力脈動(dòng)較大；蝸殼出口處的波動(dòng)情況非常平緩，這主要是因?yàn)槲仛こ隹诙蔚木彌_與穩(wěn)流作用。其次進(jìn)行測(cè)量系統(tǒng)的方案設(shè)計(jì)，先對(duì)壓力變送器和數(shù)據(jù)采集模塊的參數(shù)進(jìn)行充分的了解，完成二者的硬件選型以及電路接線圖的構(gòu)造。再利用Labview軟件編寫上位機(jī)界面，通過設(shè)置采樣頻率、采樣數(shù)據(jù)點(diǎn)、采集量程等操作

4、完成8條通道數(shù)據(jù)的采集和曲線的顯示，且保存數(shù)據(jù)。關(guān)鍵詞：離心泵，壓力脈動(dòng)，定常模擬，非定常模擬，上位機(jī)界面55AbstractCentrifugal pump produces static pressure components and produces dynamic pressure components，that is，pressure pulsation. Pressure pulsation is superimposed on static pressure components just like AC signals. Most of the spiral volute ce

5、ntrifugal pump，volute ring distortion and asymmetry of the tongue caused by impeller exit flow，especially in the non design condition，the distortion is more violent，so in the design of the spindle，the pressure must be due to the unbalanced distribution of fluctuation caused by considering. The press

6、ure pulsation in the centrifugal pump is extremely complex，and the pressure pulsation of the same pump is different in different working conditions. Pump internal flow field and unsteady turbulent characteristics，the impeller and the water inlet pipe of the rotor stator interaction，eddy current and

7、reflux may cause pressure pulsation，mechanical vibration of the pressure pulsation excitation pump and pipeline，affect the normal operation of the pump，reduce the service life of the pump； if the excitation frequency is close to the frequency of the pump. The resonance vibration will occur，will prod

8、uce pressure fluctuations in the ambient medium adjacent，have a great noise. It can be seen that the pressure pulsation in the pump is an important factor affecting the stability of the pump operation，and is also an important cause of noise pollution. Therefore，it is one of the important tasks to st

9、udy the pressure pulsation in the pump and reduce the vibration noise of the pump.This paper first applies CFD technology FLUENT software to monitor the pressure pulsation in impeller and volute simultaneously under multiple conditions，so as to reveal the relationship between pressure pulsation and

10、centrifugal pump performance，and provide a basis for further prediction of unsteady flow induced pressure pulsation in centrifugal pumps. In order to study the centrifugal pump under different working conditions and the variation of pressure fluctuation between different locations and their relation

11、ship，standard k- turbulence model and sliding mesh model based on the application of FLUENT software on the internal flow field of centrifugal pump has steady and unsteady three-dimensional simulation，get the pressure pulsation of internal flow characteristics and monitoring points of the roughly. T

12、he calculation results show that under the condition of small flow rate，the pressure pulsation at each monitoring point is large and uneven due to the unstable internal flow field，and the pressure pulsation near the impeller outlet is significantly higher than that in other regions of impeller passa

13、ge. Under the large flow condition，there will be a certain reflux in the tongue and outlet of the volute，which results in the instability of the pressure fluctuation near the worm shell. Because of the complex pressure pulsation at the location of the basin，the fluctuation of the outlet of the volut

14、e is very gentle，which is mainly due to the buffering and steady flow of the outlet section of the volute.On the other hand，on the design of the measurement system，we fully understand the parameters of the pressure transmitter and data acquisition module，and complete the hardware selection of the tw

15、o parts and the construction of the wiring diagram. Then we use LabVIEW software to write the upper machine interface. We set up sampling frequency，sampling data points and collecting ranges to complete the 8 channel data collection and curve display，and save data.Keywords：centrifugal pump, pressure

16、 fluctuation, Constant simulation, Unsteady simulation, Unsteady simulation目錄目錄摘要IAbstractII第一章 緒論11.1 課題的研究背景及意義11.2 國(guó)內(nèi)外研究現(xiàn)狀21.2.1 國(guó)外研究現(xiàn)狀21.2.2 國(guó)內(nèi)研究進(jìn)展31.3 論文工作的主要內(nèi)容4第二章 泵內(nèi)流動(dòng)數(shù)值計(jì)算方法及壓力脈動(dòng)的相關(guān)理論52.1 離心泵的三維建模52.1.1 葉輪實(shí)體模型的建立52.1.2 蝸殼實(shí)體模型的建立62.1.3 進(jìn)水管實(shí)體模型的建立72.1.4 元件的裝配82.2 數(shù)值計(jì)算方法概述92.2.1 流動(dòng)控制方程102.2.2 旋轉(zhuǎn)

17、區(qū)域中流動(dòng)問題的建模122.2.3 湍流模型132.2.4 求解控制142.3 離心泵的網(wǎng)格劃分152.3.1 網(wǎng)格相關(guān)問題的概述152.3.2 網(wǎng)格劃分162.4 離心泵壓力脈動(dòng)的分類182.5 離心泵壓力脈動(dòng)的產(chǎn)生192.5.1 動(dòng)靜干擾引起的壓力脈動(dòng)192.5.2 二次流引起的壓力脈動(dòng)202.5.3 汽蝕引起的壓力脈動(dòng)212.6 壓力脈動(dòng)的分析方法212.6.1 時(shí)域法212.6.2 頻域法222.7 本章小結(jié)23第三章 離心泵內(nèi)三維數(shù)值特性分析243.1 定常數(shù)值模擬243.1.1 網(wǎng)格檢查243.1.2 設(shè)置求解器與計(jì)算模型253.1.3 設(shè)置邊界條件253.1.4 設(shè)置檢測(cè)器與計(jì)算

18、263.1.5 不同工況下數(shù)值模擬結(jié)果分析273.2 非定常數(shù)值模擬313.2.1 計(jì)算方法313.2.2 壓力脈動(dòng)特征參數(shù)的確定323.2.3 壓力脈動(dòng)監(jiān)測(cè)點(diǎn)位置的設(shè)定323.2.4 時(shí)間步長(zhǎng)及采樣頻率的確定333.2.5 不同工況下數(shù)值模擬結(jié)果分析34第四章 測(cè)量系統(tǒng)的方案設(shè)計(jì)444.1 系統(tǒng)要求444.1.1 總體要求444.1.2 系統(tǒng)功能要求454.1.3 系統(tǒng)技術(shù)指標(biāo)454.2 硬件設(shè)計(jì)464.2.1 變送器464.2.2 信號(hào)調(diào)理模塊474.2.3 數(shù)據(jù)采集卡484.2.4 電路接線圖494.3 軟件設(shè)計(jì)504.4 USB數(shù)據(jù)通信514.5 界面設(shè)計(jì)和程序打包53第五章 結(jié)論與展

19、望555.1 結(jié)論565.2 展望57第1章 緒論第1章 緒論離心泵主要用來在某一流量下產(chǎn)生一個(gè)穩(wěn)定的壓力差以達(dá)到輸送介質(zhì)的目的 ，理想的流量與揚(yáng)程曲線給人的感覺是在某一流量點(diǎn)揚(yáng)程是靜態(tài)的。事實(shí)上，泵在產(chǎn)生靜態(tài)壓力分量的同時(shí)還產(chǎn)生動(dòng)態(tài)壓力分量，也就是壓力脈動(dòng)。壓力脈動(dòng)象交流信號(hào)一樣疊加在靜態(tài)壓力分量上。離心泵內(nèi)的壓力脈動(dòng)是極其復(fù)雜的，不同的離心泵可能表現(xiàn)不同種類的壓力脈動(dòng)，即使對(duì)同一臺(tái)離心泵，當(dāng)處于不同工況時(shí)也可能呈現(xiàn)不同性質(zhì)的壓力脈動(dòng)。因此，離心泵內(nèi)壓力脈動(dòng)的定義必須包含隨機(jī)脈動(dòng)和周期脈動(dòng)兩種。設(shè)p為離心泵內(nèi)某點(diǎn)的壓力，若對(duì)于任意的時(shí)間t和充分大的周期T，有平均壓力基本不隨時(shí)間而變化。令，則

20、即為離心泵內(nèi)某點(diǎn)的壓力脈動(dòng)。這些壓力脈動(dòng)是由葉輪進(jìn)出口回流、汽蝕等非期望的流動(dòng)特性形成的。離心泵工作過程中，一定后緣厚度、一定數(shù)量的葉片會(huì)產(chǎn)生離散頻率的壓力脈動(dòng)。而且偏工況時(shí)的脫流所引起的紊流，也會(huì)產(chǎn)生寬頻帶的壓力脈動(dòng) 。寬頻帶的和離散頻率的壓力脈動(dòng)都和進(jìn)口管、葉輪、靜子（導(dǎo)葉或蝸殼）的水力設(shè)計(jì)有復(fù)雜的關(guān)系。至今沒有一個(gè)精確的理論來預(yù)測(cè)壓力脈動(dòng)的大小。1.1 課題的研究背景及意義隨著科學(xué)技術(shù)的發(fā)展，人們對(duì)泵的可靠性和噪聲等方面提出了更高的要求，而泵的運(yùn)行不穩(wěn)定和噪聲主要是由其內(nèi)部流動(dòng)情況決定的。泵內(nèi)部流場(chǎng)非定常的湍流特性、葉輪與蝸殼的動(dòng)靜干擾、渦流及回流等都可能引起壓力脈動(dòng)，該壓力脈動(dòng)激發(fā)泵體

21、和管道的機(jī)械振動(dòng)，影響泵的正常運(yùn)轉(zhuǎn)，降低了泵的使用壽命；如果這種激勵(lì)振動(dòng)的頻率接近泵的主頻，還會(huì)發(fā)生共振，振動(dòng)會(huì)在相鄰的環(huán)境介質(zhì)內(nèi)產(chǎn)生聲壓波動(dòng)，產(chǎn)生很大的噪聲1。可見，泵內(nèi)壓力脈動(dòng)是影響水泵運(yùn)行穩(wěn)定性的重要因素，也是導(dǎo)致噪聲污染的重要原因。因而，對(duì)泵內(nèi)部壓力脈動(dòng)進(jìn)行研究，降低泵的振動(dòng)噪聲成為當(dāng)前研究的重要任務(wù)之一。離心泵廣泛應(yīng)用于生活和生產(chǎn)設(shè)備中，對(duì)其內(nèi)部流動(dòng)壓力脈動(dòng)的研究越來越受到人們的重視。因此，針對(duì)離心泵壓力脈動(dòng)的研究就具有極其重要的學(xué)術(shù)價(jià)值和應(yīng)用價(jià)值，將有助于認(rèn)識(shí)壓力脈動(dòng)特性，為控制和降低壓力脈動(dòng)提供有益的參考。本文選取離心泵為研究對(duì)象，在前人對(duì)泵內(nèi)壓力脈動(dòng)的研究的基礎(chǔ)上，對(duì)其內(nèi)部流動(dòng)

22、及其壓力脈動(dòng)進(jìn)行研究。1.2 國(guó)內(nèi)外研究現(xiàn)狀1.2.1 國(guó)外研究現(xiàn)狀自上世紀(jì)五十年代，水泵在運(yùn)行過程中所產(chǎn)生的壓力脈動(dòng)就引起了人們的注意，此時(shí)的研究將壓力脈動(dòng)看作隨機(jī)信號(hào)，如清華大學(xué)的吳仁卿、劉冀生等；80年代初，日本學(xué)者宮代裕、近藤正遭等在對(duì)壓力脈動(dòng)進(jìn)行研究時(shí)，把壓力脈動(dòng)看作是一種有規(guī)律的周期脈動(dòng)。此后，隨著測(cè)量技術(shù)的進(jìn)步以及計(jì)算流體動(dòng)力學(xué)CFD的飛速發(fā)展，壓力脈動(dòng)的研究進(jìn)入了一個(gè)新的階段。用測(cè)試的方法了解離心泵內(nèi)部壓力脈動(dòng)情況被認(rèn)為是最基本和最可靠的手段。Arndt等試驗(yàn)研究了離心泵葉輪和蝸殼間的流動(dòng)干擾，發(fā)現(xiàn)吸力面的壓力波動(dòng)大于壓力面，最大的壓力波動(dòng)出現(xiàn)在吸力面進(jìn)口邊，當(dāng)葉輪與蝸殼導(dǎo)葉間

23、的徑向間隙增大時(shí)該壓力波動(dòng)減小。Tatsuji等通過試驗(yàn)對(duì)葉輪出口速度、壓力的波動(dòng)特性進(jìn)行了研究，指出速度的波動(dòng)性深受尾跡的影響，而靜壓的波動(dòng)性受其影響要小得多，并對(duì)波動(dòng)的基波及其諧波特性進(jìn)行了描述。Wang等通過試驗(yàn)研究了導(dǎo)葉內(nèi)的壓力脈動(dòng)，得出其主頻的壓力脈動(dòng)出現(xiàn)在非設(shè)計(jì)點(diǎn)時(shí)的葉頻倍頻處。Kelder等對(duì)一低比轉(zhuǎn)速離心泵蝸殼內(nèi)的非定常流動(dòng)進(jìn)行了理論分析和試驗(yàn)研究。分析計(jì)算了在設(shè)計(jì)工況下由二次流引起的壓力脈動(dòng)，發(fā)現(xiàn)即使是在設(shè)計(jì)點(diǎn)附近，由二次流引起的壓力脈動(dòng)已經(jīng)達(dá)到了靜態(tài)揚(yáng)程的一半甚至還要大。Dong、Chu等使用帶有噪聲和壓力測(cè)量裝置的PIV（Particle Image Velocimetr

24、y）對(duì)不同隔舌間隙（間隙為葉輪半徑的7%28%）的蝸殼內(nèi)的速度場(chǎng)進(jìn)行了測(cè)量，并計(jì)算了蝸殼隔舌附近的非定常壓力場(chǎng)。研究發(fā)現(xiàn)，當(dāng)隔舌和葉輪之間的間隙約等于葉輪半徑的20%時(shí)，離心泵的振動(dòng)噪聲將會(huì)明顯地降低；在葉片間插入短葉片后，原始尾流有所減弱，但出現(xiàn)了另一射流一尾流現(xiàn)象。由于成本過高，采用試驗(yàn)手段監(jiān)測(cè)泵內(nèi)（特別是轉(zhuǎn)動(dòng)的葉片上）的壓力脈動(dòng)通常比較困難，從而迫使葉輪內(nèi)流動(dòng)的數(shù)值模擬技術(shù)獲得了更大的發(fā)展。1995年Fortes-Patella等應(yīng)用重疊網(wǎng)格和相位滯后的周期邊界條件，對(duì)泵內(nèi)葉輪和蝸殼的動(dòng)靜干擾進(jìn)行了二維非定常計(jì)算。Gonzfilez等應(yīng)用Navier-Stokes方程和滑移網(wǎng)格技術(shù)，對(duì)離

25、心泵全流道進(jìn)行了三維非定常數(shù)值模擬，討論了蝸殼隔舌附近由于動(dòng)靜干擾引起的壓力脈動(dòng)，并與試驗(yàn)值進(jìn)行了對(duì)比，給出了預(yù)測(cè)葉輪出口和蝸殼隔舌動(dòng)態(tài)相互作用的數(shù)學(xué)模型。在此前研究的基礎(chǔ)上，2006年Gonzfilez等對(duì)流場(chǎng)結(jié)構(gòu)與全局變量之間的關(guān)系作了進(jìn)一步的研究，指出蝸殼內(nèi)二次流與不同的壓力脈動(dòng)類型有關(guān)。Qin和Tsukamoto分別用奇點(diǎn)分布法和RANS（Reynolds-averaged Navier-Stokes equations）模型計(jì)算了導(dǎo)葉泵內(nèi)的葉輪-導(dǎo)葉干涉引起的非定常流，通過把兩個(gè)結(jié)果進(jìn)行對(duì)比來區(qū)別不穩(wěn)定壓力中的勢(shì)流作用，從而確定了不同來源的壓力脈動(dòng)。渦方法因?yàn)樗臒o網(wǎng)格性所以適合應(yīng)用

26、在復(fù)雜幾何體內(nèi)的流動(dòng)，可以用它來計(jì)算導(dǎo)葉泵內(nèi)復(fù)雜的非定常流的壓力脈動(dòng)。Zhu用離散渦方法對(duì)額定流量下的動(dòng)靜干擾引起的壓力脈動(dòng)進(jìn)行了計(jì)算，但是在改變工況點(diǎn)時(shí)這種方法就不適用了。Wang等在動(dòng)靜干擾問題上用邊界集成方法引出一個(gè)能改進(jìn)非穩(wěn)定壓力評(píng)估的新方法，其采用的渦方法可在變工況下應(yīng)用。隨著研究的深入，除了對(duì)所獲得的壓力脈動(dòng)現(xiàn)象進(jìn)行分析之外，人們致力于把壓力脈動(dòng)與紊流理論聯(lián)系起來深入到壓力脈動(dòng)內(nèi)在機(jī)理問題的研究。Dring等指出動(dòng)靜干擾的兩個(gè)直接機(jī)理：勢(shì)流作用和尾流作用。一方面，由于葉輪葉片和靜止的導(dǎo)葉之間的相對(duì)運(yùn)動(dòng)引起的無粘性相互作用產(chǎn)生了勢(shì)流的相互作用；另一方面，葉輪流道內(nèi)產(chǎn)生的尾跡進(jìn)入后續(xù)的

27、導(dǎo)葉流道時(shí)，由于尾跡的沖擊和傳播形成了尾跡的相互作用。此外，許多學(xué)者從離心泵的水力設(shè)計(jì)方面著手，致力于降低離心泵壓力脈動(dòng)以及振動(dòng)、噪聲的研究。大量研究表明，離心泵最強(qiáng)烈的壓力脈動(dòng)出現(xiàn)在隔舌附近。為此，許多學(xué)者試圖通過改變?nèi)~輪和隔舌參數(shù)來降低壓力脈動(dòng)水平，如優(yōu)化葉輪和隔舌的間隙、改變?nèi)~輪結(jié)構(gòu)形式（分流葉片、雙葉片葉輪）、改變隔舌形狀（階梯型隔舌）、改變?nèi)~輪葉片數(shù)等。綜上，國(guó)外對(duì)離心泵壓力脈動(dòng)的研究取得了較多的成果，有很好的借鑒價(jià)值，特別是關(guān)于壓力脈動(dòng)機(jī)理的分析以及現(xiàn)象的描述，對(duì)降低壓力脈動(dòng)的工程研究有很好的幫助。1.2.2 國(guó)內(nèi)研究現(xiàn)狀目各高校針對(duì)工程實(shí)際用軸流泵、貫流泵等內(nèi)部流動(dòng)誘導(dǎo)壓力脈動(dòng)進(jìn)

28、行了數(shù)值計(jì)算與試驗(yàn)研究。中國(guó)農(nóng)業(yè)大學(xué)的張玲采用大渦模擬對(duì)軸流泵多個(gè)工況下的壓力脈動(dòng)進(jìn)行了預(yù)測(cè)，得到了監(jiān)測(cè)點(diǎn)的時(shí)域特性和頻域特性。耿少娟等應(yīng)用商業(yè)軟件Fluent分析了無短葉片、長(zhǎng)短葉片和短短葉片三種葉輪的單級(jí)離心泵整機(jī)由于動(dòng)靜干擾引起的葉輪進(jìn)口和蝸殼出口的壓力脈動(dòng)。清華大學(xué)的徐朝暉對(duì)高速泵動(dòng)靜葉柵間采用滑移網(wǎng)格技術(shù)建立交互界面，利用RNG（Random Number Generator）湍流模型對(duì)包括誘導(dǎo)輪在內(nèi)的全流道進(jìn)行了非定常數(shù)值模擬，并對(duì)頻譜進(jìn)行了FFT分析。陳黨民等應(yīng)用商業(yè)軟件STAR-CD對(duì)三種工況下的部分流泵整機(jī)內(nèi)的非定常流場(chǎng)進(jìn)行了分析，對(duì)各葉輪流道內(nèi)壓力波動(dòng)以及揚(yáng)程的瞬變波動(dòng)進(jìn)行

29、了描述。施衛(wèi)東等通過對(duì)軸流泵壓力脈動(dòng)特性的研究，表明不同轉(zhuǎn)速下壓力脈動(dòng)的主頻發(fā)生偏移，說明壓力脈動(dòng)主頻與轉(zhuǎn)速呈非線性關(guān)系。王福軍提出軸流泵葉輪進(jìn)口前壓力脈動(dòng)的幅值最大，壓力脈動(dòng)頻率主要受葉輪轉(zhuǎn)頻影響，并提出應(yīng)避免大型機(jī)組固有模態(tài)與葉片通過頻率成整數(shù)倍。鄭源等基于不同軸流泵葉片安放角度和揚(yáng)程，探討了壓力脈動(dòng)的變化規(guī)律及特性。王洪杰等通過對(duì)渦輪壓力脈動(dòng)的研究，表明在設(shè)計(jì)工況附近，葉輪入口處空化誘發(fā)連續(xù)的高頻壓力脈動(dòng)，而動(dòng)靜干涉效應(yīng)誘發(fā)壓力脈動(dòng)以轉(zhuǎn)頻、葉頻以及兩者的倍頻為主。姚志峰等通過對(duì)種方案的雙吸離心泵葉輪進(jìn)行壓力脈動(dòng)同臺(tái)試驗(yàn)，提出相對(duì)于傳統(tǒng)雙吸葉輪，兩側(cè)葉片交錯(cuò)布置后可改善壓力脈動(dòng)特性，減輕動(dòng)

30、靜干涉的危害。潘中永通過對(duì)離心泵葉輪與鍋殼動(dòng)靜干涉的研究，提出由動(dòng)靜干涉引起的壓力脈動(dòng)的主頻與葉片通過頻率相同，最終給出了一種可以測(cè)量轉(zhuǎn)速的新方法。徐朝暉等對(duì)泵內(nèi)動(dòng)靜干涉作用進(jìn)行了數(shù)值模擬，認(rèn)為定子與動(dòng)子間各點(diǎn)的壓力脈動(dòng)頻率相同，且幅值沿周向呈現(xiàn)出正弦變化規(guī)律。1.3 論文工作的主要內(nèi)容1.使用Pro/E分別建立葉輪、蝸殼和進(jìn)水管的三維模型，建立成功后裝配成離心泵的整體結(jié)構(gòu)；再使用Gambit軟件進(jìn)行網(wǎng)格劃分，檢查劃分后的網(wǎng)格質(zhì)量是否符合要求。2.使用Fluent軟件，通過在流場(chǎng)中布置監(jiān)測(cè)點(diǎn)的方式，研究壓力脈動(dòng)與監(jiān)測(cè)點(diǎn)位置、運(yùn)行工況之間的關(guān)系，分析各監(jiān)測(cè)點(diǎn)在不同工況下壓力脈動(dòng)值的時(shí)域特性以及內(nèi)

31、部流場(chǎng)的分布情況。完成壓力變送器、數(shù)據(jù)采集模塊的器件選型和電路連接方式的測(cè)量系統(tǒng)方案設(shè)計(jì)；使用Labview軟件繪制出上位機(jī)界面，采集所在布置監(jiān)測(cè)點(diǎn)的壓力脈動(dòng)數(shù)據(jù)，顯示壓力脈動(dòng)曲線、保存實(shí)驗(yàn)數(shù)據(jù)。本章先通過定常數(shù)值模擬，得到離心泵內(nèi)部的靜壓力場(chǎng)、速度場(chǎng)，進(jìn)一步分析了隔舌處的壓力及速度分布情況。采用定常計(jì)算結(jié)果作為非定常計(jì)算的初始條件，通過對(duì)離心泵進(jìn)行三維非定常數(shù)值模擬，數(shù)值分析了不同工況下、不同位置處和不同時(shí)刻時(shí)離心泵的壓力脈動(dòng)特性。第2章 泵內(nèi)流動(dòng)數(shù)值計(jì)算方法及壓力脈動(dòng)的相關(guān)理論第2章 泵內(nèi)流動(dòng)數(shù)值計(jì)算方法及壓力脈動(dòng)的相關(guān)理論本文是以某離心泵作為研究對(duì)象，其中進(jìn)水管、葉輪和蝸殼三個(gè)主要部件作

32、為設(shè)計(jì)模型的計(jì)算區(qū)域，為了能夠更加精確地分析流動(dòng)部件內(nèi)流場(chǎng)的分布規(guī)律，可以忽略一些對(duì)流動(dòng)影響不大的部件如部件壁厚、聯(lián)接零部件等，這樣為計(jì)算模型的建模及數(shù)值模擬提供了極大的方便。本章詳細(xì)介紹了利用Pro/E軟件對(duì)計(jì)算模型進(jìn)行三維造型過程及其計(jì)算模型計(jì)算域的網(wǎng)格劃分過程。2.1 離心泵的三維建模該離心泵的轉(zhuǎn)速為1500r/min，其三個(gè)主要過流部件包括進(jìn)水管、葉輪和蝸殼。其中，葉輪是旋轉(zhuǎn)運(yùn)動(dòng)的，進(jìn)水管和蝸殼是靜止不動(dòng)的，通過葉輪與進(jìn)水管、蝸殼的相互動(dòng)靜耦合，對(duì)計(jì)算模型全流道進(jìn)行模擬計(jì)算。本文根據(jù)已知的計(jì)算模型泵的水力圖，采用功能強(qiáng)大的實(shí)體建模Pro/E軟件分別對(duì)該離心泵的三個(gè)主要過流部件進(jìn)行三維全

33、流道建模，然后裝配成整體的全流道實(shí)體模型。液體在泵內(nèi)運(yùn)動(dòng)情況可通過泵的水力特性外在表現(xiàn)出來，而內(nèi)部流體的流動(dòng)情況是預(yù)測(cè)泵性能的基礎(chǔ)。2.1.1 葉輪實(shí)體模型的建立離心泵的葉輪是完成液體傳輸?shù)闹饕考?。本文的葉輪結(jié)構(gòu)對(duì)稱，其具體建模過程如下：（1）在圓柱坐標(biāo)系中，使用偏移坐標(biāo)系基準(zhǔn)點(diǎn)命令，分別輸入所給葉片軸面圖中的每條軸面截線坐標(biāo)點(diǎn)；使用插入基準(zhǔn)曲線命令，將同一軸面截線上的坐標(biāo)點(diǎn)連接起來；再使用邊界混合工具將各個(gè)曲線連接起來，便生成了葉片的各個(gè)面；再使用曲面合并命令，將葉片六個(gè)面兩兩合并，最終構(gòu)成一個(gè)封閉的曲面；選中形成的閉合曲面，利用實(shí)體化命令將閉合曲面變?yōu)閷?shí)體；參照選擇葉輪葉片的工作面和背面

34、，驅(qū)動(dòng)曲面選擇葉輪葉片的內(nèi)側(cè)邊界混合曲面，使用完全倒圓角工具對(duì)葉片模型加以修飾。（2）建組從邊界混合操作開始到最后倒圓角操作結(jié)束，右鍵組選擇陣列操作，選定旋轉(zhuǎn)軸，第一方向的陣列成員數(shù)輸入6，再生選項(xiàng)選一般即通過假設(shè)所有成員形狀各異且可能彼此相交來計(jì)算成員幾何，使用陣列工具創(chuàng)造出另外5個(gè)葉輪葉片。（3）使用草繪工具繪制出葉輪的軸面投影圖，旋轉(zhuǎn)軸選擇內(nèi)部CL，角度選擇在第一個(gè)方向從草繪平面以指定的角度作為實(shí)體旋轉(zhuǎn)，使用旋轉(zhuǎn)工具生成出不完整的葉輪模型。（4）使用插入-共享數(shù)據(jù)-合并/繼承命令，將參照類型設(shè)置為外部，選取葉片實(shí)體作為參照模型，定義兩個(gè)對(duì)齊約束類型。選擇約束類型為對(duì)齊：分別選定葉輪模型和

35、葉片模型中的基準(zhǔn)軸；選擇約束類型為對(duì)齊：選定葉輪模型和葉片模型中的基準(zhǔn)平面，偏移選擇重合即將原件置于和組件參照重合的位置。移除材料、選擇從屬關(guān)系，指定約束時(shí)在組件窗口中顯示元件，實(shí)現(xiàn)外部切除。生成葉輪的實(shí)體模型保存為108-16-yl.prt，模型如圖2.1所示。圖2.1 離心泵的葉輪模型Fig.2.1 Impeller of centrifugal pump2.1.2 蝸殼實(shí)體模型的建立蝸殼是離心泵中的又一重要過流部件，其具體建模過程如下：（1） 利用草繪工具繪制出離心泵蝸殼的截面圖。為方便后續(xù)做掃描混合的時(shí)候使用，使用基準(zhǔn)點(diǎn)工具創(chuàng)建出蝸殼八個(gè)斷面與基圓的交點(diǎn)。按照給定的蝸殼的水力模型圖，再

36、利用草繪工具分別做出八個(gè)斷面的形狀。對(duì)于其他斷面，參照基準(zhǔn)坐標(biāo)系CS0，類型選擇圓柱坐標(biāo)系，使用偏移坐標(biāo)系工具，輸入蝸殼水力模型圖上的斷面數(shù)據(jù)，選擇創(chuàng)建出的基準(zhǔn)點(diǎn)，使用插入基準(zhǔn)曲線工具，繪制出其他斷面的形狀。（2） 使用偏移、鏡像和草繪工具繪制出引導(dǎo)線，在第一方向上選擇相鄰兩個(gè)斷面作為鏈，在第二方向上選擇引導(dǎo)線作為鏈，使用邊界混合工具繪制出相鄰兩個(gè)斷面的過渡曲面。在蝸殼的隔舌處，選中繪制處的樣條曲線，利用拉伸工具，從中心向兩側(cè)拉升做一拉伸曲面，深度要超過渦室進(jìn)口寬度，此步驟是方便后續(xù)的補(bǔ)面之用。參照選擇修剪的面組為鏡像步驟的一側(cè)曲面、修剪對(duì)象選擇拉伸出的曲面，保留修剪曲面，使用修剪工具進(jìn)行修剪

37、；使用相同的方法對(duì)鏡像步驟的另一側(cè)曲面進(jìn)行修剪，最后使用邊界混合工具生成蝸殼隔舌處的過渡曲面。（3） 利用曲面合并命令，兩兩合并逐步將所有的曲面合并為一個(gè)閉合的曲面。選中合并后的閉合曲面，利用實(shí)體化命令將閉合曲面變?yōu)閷?shí)體。蝸殼的實(shí)體模型成功建立保存為108-16-wk.prt，模型如圖2.2所示。圖2.2 離心泵的蝸殼模型Fig.2.2 Volute of centrifugal pump2.1.3 進(jìn)水管實(shí)體模型的建立進(jìn)水管模型的建立過程比較簡(jiǎn)單，使用草繪工具繪制出葉輪進(jìn)口面曲線，從草繪平面以指定的深度值拉伸為實(shí)體即可獲得進(jìn)水管的實(shí)體模型，保存為108-16-jsg.prt，如圖2.3所示。

38、圖2.3 離心泵的進(jìn)水管模型Fig.2.3 Inlet pipe of centrifugal pump2.1.4 元件的裝配（1）使用裝配工具將元件108-16-wk.prt添加到組件中，使用用戶定義約束集，約束類型選擇在缺省位置裝配元件；運(yùn)動(dòng)類型選擇平移定義方向，在視圖平面中相對(duì)來調(diào)整或選取參照；指定約束時(shí)在組件窗口中顯示元件。（2）使用裝配工具將元件108-16-yl.prt添加到組件中，定義兩個(gè)對(duì)齊約束類型。選擇約束類型為對(duì)齊：選定葉輪模型中的基準(zhǔn)軸A_1和蝸殼模型中的基準(zhǔn)軸A_1；選擇約束類型為對(duì)齊：選定葉輪模型中的FRONT基準(zhǔn)平面和ASM_FRONT基準(zhǔn)平面，偏移選擇重合即將原件

39、置于和組件參照重合的位置。指定約束時(shí)在組件窗口中顯示元件。（3）使用裝配工具將元件108-16-jsg.prt添加到組件中，仍定義兩個(gè)對(duì)齊約束類型。選擇約束類型為對(duì)齊：選定進(jìn)口延伸段模型中的基準(zhǔn)軸A_1和葉輪模型中的基準(zhǔn)軸A_1；選擇約束類型為對(duì)齊：選定進(jìn)口延伸段模型中的FRONT基準(zhǔn)平面和葉輪模型中的進(jìn)口圓面，偏移選擇重合即將原件置于和組件參照重合的位置。指定約束時(shí)在組件窗口中顯示元件。裝配完成后的離心泵實(shí)體模型組件如圖2.4所示。至此，整個(gè)離心泵的三維建模工作完成，模型保存為108-16-zhuangpei.asm。為方便后續(xù)導(dǎo)入Gambit軟件里進(jìn)行網(wǎng)格劃分，將模型108-16-zhua

40、ngpei.asm保存副本為STEP類型文件即108-16-zhuangpei_asm.stp。圖2.4 離心泵的整體三維模型Fig.2.4 Three dimensional model of centrifugal pump2.2 數(shù)值計(jì)算方法概述近年來，隨著計(jì)算機(jī)技術(shù)以及數(shù)值計(jì)算理論的更新?lián)Q代，針對(duì)流場(chǎng)進(jìn)數(shù)值行模擬和預(yù)測(cè)的CFD技術(shù)應(yīng)運(yùn)而生。利用CFD技術(shù)，可以在短時(shí)間內(nèi)預(yù)測(cè)流場(chǎng)性能，并通過改變流場(chǎng)中參數(shù)設(shè)置反復(fù)計(jì)算，使泵產(chǎn)品達(dá)到最佳設(shè)計(jì)效果，同時(shí)獲得一些實(shí)驗(yàn)中無法獲得的流場(chǎng)分布情況。經(jīng)過多年的發(fā)展，CFD技術(shù)完全可以準(zhǔn)確地反映流體機(jī)械內(nèi)部的流場(chǎng)分布情況，因此在新產(chǎn)品開發(fā)方面應(yīng)用廣泛，在

41、有效減少設(shè)計(jì)的工作量并提高設(shè)計(jì)的可靠性的同時(shí)節(jié)省了實(shí)驗(yàn)資源。在基于Fluent軟件平臺(tái)，對(duì)離心泵的內(nèi)部流場(chǎng)進(jìn)行數(shù)值計(jì)算的過程中，需要分別對(duì)求解器中參數(shù)進(jìn)行設(shè)置。每一個(gè)參數(shù)都有它的意義和適用范圍，并且求解的模型不同，各參數(shù)不同的匹配方式得到的計(jì)算結(jié)果不同。因此，為了正確、高效運(yùn)用Fluent軟件解決具體問題，對(duì)CFD技術(shù)的基礎(chǔ)知識(shí)做充分的了解是十分必要的。下面在針離心泵內(nèi)部流動(dòng)數(shù)值模擬時(shí)的應(yīng)用過程，對(duì)CFD技術(shù)做簡(jiǎn)單介紹。2.2.1 流動(dòng)控制方程通常情況下，泵可以認(rèn)為是一種工作介質(zhì)為絕熱的不可壓縮流體的工作機(jī)，液體在泵內(nèi)的流動(dòng)狀態(tài)主要表現(xiàn)為兩個(gè)過程。第一旋轉(zhuǎn)的葉輪葉片對(duì)液體做功，第二液體自身動(dòng)能

42、與勢(shì)能之間相互轉(zhuǎn)化。因此，在運(yùn)用CFD技術(shù)對(duì)離心泵流動(dòng)問題數(shù)值模擬的過程中，并不需要考慮其能量方程，而只需要考慮連續(xù)性方程和動(dòng)量方程。具體控制方程如下：（1） 慣性坐標(biāo)系下的控制方程（a） 連續(xù)性方程連續(xù)性方程是基于質(zhì)量守恒定律，也就是說同一流體質(zhì)點(diǎn)的質(zhì)量在運(yùn)動(dòng)過程中始終不變。不可壓縮流體連續(xù)性方程為： （2.1）其中為哈密頓算子，。直角坐標(biāo)系中連續(xù)性方程的表達(dá)式為： （2.2）其中，u、v、w分別為速度在x、y、z3個(gè)坐標(biāo)軸方向上的分量。（b） 動(dòng)量方程（Navier-Stokes方程）動(dòng)量方程即納威-斯托克斯方程，即體積V中流體質(zhì)點(diǎn)動(dòng)量的變化率等于作用在該體積上的質(zhì)量力與表面力以及粘性力三

43、者之和。不可壓縮流體的動(dòng)量方程表達(dá)式為： （2.3）式中，t為時(shí)間，為作用在單位質(zhì)量流體上的質(zhì)量力，為流體的密度，p為壓力，v為運(yùn)動(dòng)粘度，為拉普拉斯算子。直角坐標(biāo)系中動(dòng)量方程的表達(dá)式為： （2.4a） （2.4b） （2.4c）其中，、為分別為質(zhì)量力在x、y、z三個(gè)方向上的分量。（2）旋轉(zhuǎn)坐標(biāo)系中的控制方程慣性坐標(biāo)系下，由于葉輪內(nèi)流體隨著葉輪葉片一起旋轉(zhuǎn)，認(rèn)為流動(dòng)處于周期性地非定常狀態(tài)?？梢越⒁粋€(gè)與葉輪以同樣角速度旋轉(zhuǎn)的坐標(biāo)系，此時(shí)葉輪內(nèi)部的流動(dòng)可以認(rèn)為是定常狀態(tài)，從而使得求解過程得以簡(jiǎn)化。相對(duì)速度和絕對(duì)速度之間的關(guān)系如下： （2.5）其中，為遷移速度，為葉輪的旋轉(zhuǎn)角速度，為質(zhì)點(diǎn)矢徑。絕對(duì)加

44、速度與相對(duì)加速度及牽連加速度和科氏加速度之間的表達(dá)式為： （2.6）其中，。（a） 連續(xù)性方程 （2.7）其中，下標(biāo)r表示旋轉(zhuǎn)坐標(biāo)系和中的參數(shù)。在相對(duì)運(yùn)動(dòng)中直角坐標(biāo)系下的表達(dá)式為： （2.8）其中，、分別為在x、y、z方向上的分量。（b） 動(dòng)量方程當(dāng)轉(zhuǎn)速恒定，系統(tǒng)角加速度慣性力項(xiàng)為零，所以動(dòng)量方程為： （2.9）其中： （2.10）表示作用在單位質(zhì)量流體上的科氏力，表示作用在單位質(zhì)量流體上的離心力。動(dòng)量方程在相對(duì)運(yùn)動(dòng)直角坐標(biāo)系中的表達(dá)式為： （2.11a） （2.11b） （2.11c）其中，、分別為在x、y、z三個(gè)坐標(biāo)軸上的分量： （2.12a） （2.12b） （2.12c）其中：、分別為

45、在x、y、z三個(gè)坐標(biāo)軸上的分量，、分別為在x、y、z三個(gè)坐標(biāo)軸上的分量。2.2.2 旋轉(zhuǎn)區(qū)域中流動(dòng)問題的建模（1）多參考系模型多參考系模型MRF（Multi-Reference Frame）為不同旋轉(zhuǎn)或者移動(dòng)速度的每一個(gè)細(xì)小單元體的穩(wěn)態(tài)近似。MRF模型是通過創(chuàng)建多個(gè)參考系來實(shí)現(xiàn)動(dòng)靜耦合，其動(dòng)靜交界面上的網(wǎng)格沒有移動(dòng)，也就是對(duì)于每個(gè)參考系的交界面的速度（絕對(duì)速度）必須相同。該方法比較適合于均勻混合時(shí)（邊界上流動(dòng)區(qū)域幾乎相一致）。顯然多參考系模型方法是簡(jiǎn)化了的，它可以在很多時(shí)均流動(dòng)的應(yīng)用場(chǎng)合里提供合理的模型。一般來說，當(dāng)轉(zhuǎn)子和定子之間干涉作用相對(duì)較弱、并無大范圍瞬態(tài)影響時(shí)，MRF模型完全可以滿足模

46、擬要求。通常在旋轉(zhuǎn)機(jī)械數(shù)值模擬當(dāng)中，可以先采用MRF模型進(jìn)行流場(chǎng)計(jì)算，然后將計(jì)算結(jié)果作為瞬態(tài)滑動(dòng)網(wǎng)格模型的初始條件再繼續(xù)精確計(jì)算。（2）滑移網(wǎng)格技術(shù)滑移網(wǎng)格技術(shù)（Moving Mesh）運(yùn)用兩個(gè)或更多的單元區(qū)域。在每個(gè)單元區(qū)域中至少存在一個(gè)邊界的交界面，該交界面和另一單元區(qū)域相鄰，彼此相鄰的兩個(gè)單元的交界面形成“網(wǎng)格交界面”，因此這兩個(gè)單元彼此相對(duì)于滑動(dòng)。需要注意的是，網(wǎng)格交界面必須定義，因?yàn)閮蛇吘辛黧w通過，“網(wǎng)格交界面”必須定義于動(dòng)靜交界流體區(qū)域內(nèi)，不能位于轉(zhuǎn)子或定子邊緣的固體邊界部分。在計(jì)算過程中，相鄰單元區(qū)域的網(wǎng)格沿著“網(wǎng)格交界面”相互滑動(dòng)（旋轉(zhuǎn)或平移），此時(shí)交界面兩側(cè)的網(wǎng)格結(jié)點(diǎn)可以不

47、相互重合，但是兩側(cè)的通量一定要保證相同。對(duì)于非定常周期性的計(jì)算，其最終解并不取決于起始階段解的時(shí)間步長(zhǎng)。由于流動(dòng)起始階段的精確解對(duì)計(jì)算結(jié)果影響不大，因此可以采用“較大”的時(shí)間步長(zhǎng)開始計(jì)算，這樣能夠使流動(dòng)問題的解更快地出現(xiàn)周期性變化規(guī)律。然而，當(dāng)解出現(xiàn)周期性時(shí)，應(yīng)當(dāng)減少時(shí)間步長(zhǎng)從而獲得精確解。另外，如果計(jì)算過程中采用二階時(shí)間精度求解，那么改變時(shí)間步長(zhǎng)有可能影響當(dāng)前解的精度。若求解過程中采用“較大”時(shí)間步長(zhǎng)開始計(jì)算，隨后改變時(shí)間步長(zhǎng)的幅度不應(yīng)超過20%，否則將不能確保改變時(shí)間步長(zhǎng)后，在最后幾個(gè)周期內(nèi)的解滿足周期性狀態(tài)。2.2.3 湍流模型湍流是一種非常復(fù)雜的不穩(wěn)定三維流動(dòng)，在湍流中流場(chǎng)的各個(gè)物理量

48、均隨時(shí)間和空間隨機(jī)改變。為了研究湍流的運(yùn)動(dòng)規(guī)律，把流場(chǎng)中任一點(diǎn)的瞬時(shí)物理量近似看作是平均值和脈動(dòng)值之和，對(duì)Navier-Stokes方程作時(shí)間平均得到雷諾方程，從而產(chǎn)生新的未知量雷諾應(yīng)力項(xiàng)，為了使湍流的平均雷諾方程封閉，需要建立不同的湍流模型。因此在數(shù)值模擬中，建立正確的湍流模型對(duì)計(jì)算結(jié)果的準(zhǔn)確性顯得尤為重要。多年以來，學(xué)者對(duì)湍流相關(guān)理論進(jìn)行了深入研究，在實(shí)踐運(yùn)用中效果顯著、收益良好。由于湍流流動(dòng)是一種高度無序的復(fù)雜流動(dòng)，對(duì)其運(yùn)動(dòng)方程只能通過數(shù)值模擬的方法獲得其近似解。泵內(nèi)部的絕大多數(shù)液體均處于湍流充分發(fā)展區(qū)，為描述泵體內(nèi)湍流的瞬時(shí)運(yùn)動(dòng)，可以采用連續(xù)性方程和動(dòng)量方程。目前，湍流模擬的方法可分為

49、直接數(shù)值模擬（Direct Numerical Simulation，簡(jiǎn)稱DNS）和非直接數(shù)值模擬兩大類?；诮坪秃?jiǎn)化方法的不同，可以將非直接數(shù)值模擬方法分為大渦模擬（Large Eddy Simulation，簡(jiǎn)稱LES），分離渦模擬（Dettached Eddy Simulation，簡(jiǎn)稱DES）、雷諾平均法（RANS）和統(tǒng)計(jì)平均法等。湍流模型如圖2.5所示。雷諾時(shí)均法在關(guān)于泵內(nèi)部流場(chǎng)分析的文章中使用最為廣泛。在滿足工程實(shí)際應(yīng)用的同時(shí)避免了直接數(shù)值模擬法計(jì)算量大的問題，本文隨后選用RANS兩方程模型中的Standard k-模型進(jìn)行計(jì)算，在此重點(diǎn)介紹。Standard k-模型是由Lau

50、nder和Spalding最早提出來的，后來被工程界廣泛采納。該湍流模型基于渦黏性假設(shè)，其湍動(dòng)能和耗散率的控制方程為： （2.13） （2.14）式中：和為用戶自定義的源項(xiàng)：為湍流粘性系數(shù)，模型常數(shù)=0.09，=1.44，=1.0，=1.3。標(biāo)準(zhǔn)k-模型具有穩(wěn)定、簡(jiǎn)單和經(jīng)濟(jì)等特點(diǎn)，該湍流模型目前使用最為廣泛，它能夠在較大的工程范圍內(nèi)滿足一定的精度。圖2.5 湍流數(shù)值模擬計(jì)算方法及湍流模型Fig.2.5 Numerical simulation method and turbulence model for turbulent flow2.2.4 求解控制對(duì)于在求解域內(nèi)所建立的偏微分方程，理論上

51、是有真解的。但是，由于所處理問題自身的復(fù)雜性，很難獲得方程的解析解。因?yàn)椋谇蠼庥騼?nèi)存在著有限的因變量，為此將這些因變量值采用數(shù)值模擬方法變?yōu)槲粗縼硖幚?，并建立代?shù)方程組來表示這些未知量，然后求解所建立的代數(shù)方程，得到這些節(jié)點(diǎn)值。同樣的，為求得計(jì)算域內(nèi)其他節(jié)點(diǎn)位置上的值，可采用插值法來求解。這樣，這種用數(shù)值離散來求解變量獲得所求問題的精確解的方法稱為離散近似。三維控制方程離散后，離散方程中還包含有各種未知量，而這些變量并不能夠直接求解出來，因此，還必須對(duì)它們進(jìn)行特殊處理，對(duì)離散方程做出特定的調(diào)整。本章選用SIMPLE（Semi-Implicit Method for Pressure-Lin

52、ked Equations）算法求解離散方程，對(duì)所有變量整場(chǎng)聯(lián)立求解，同時(shí)求解連續(xù)方程和動(dòng)量方程組。SIMPLE算法為求解壓力耦合方程組的半隱式方法，該方法基于交錯(cuò)網(wǎng)格，屬于壓力修正法的一種，主要用于求解不可壓縮流場(chǎng)的數(shù)值方法，用途最為廣泛。該方法采用“猜測(cè)-修正”的過程，在交錯(cuò)網(wǎng)格的基礎(chǔ)上計(jì)算壓力場(chǎng)，從而求解Navier-Stokes方程。SIMPLE算法的基本思想可以簡(jiǎn)單的歸納為：對(duì)于給定的壓力場(chǎng)，求解離散形式的動(dòng)量方程，以此得出速度場(chǎng)。由于壓力場(chǎng)是人為假定的，不夠精確，由此得到的速度場(chǎng)一般不滿足相應(yīng)的連續(xù)方程，因此，必須對(duì)給定的壓力場(chǎng)加以修正。2.3 離心泵的網(wǎng)格劃分網(wǎng)格劃分?jǐn)?shù)值離散的基

53、礎(chǔ)，是成功實(shí)現(xiàn)數(shù)值模擬的前提，并且計(jì)算網(wǎng)格質(zhì)量的好壞和所選湍流模型是否正確對(duì)于數(shù)值模擬結(jié)果的可靠性和精度至關(guān)重要。2.3.1 網(wǎng)格相關(guān)問題的概述在對(duì)指定問題進(jìn)行CFD計(jì)算之前，首先要將計(jì)算域進(jìn)行離散化，即對(duì)空間上連續(xù)的計(jì)算區(qū)域劃分為若干個(gè)子區(qū)域，并確定每個(gè)區(qū)域中的節(jié)點(diǎn)，從而生成網(wǎng)格?？刂品匠痰碾x散是以網(wǎng)格為基礎(chǔ)的，離散化的物理量都儲(chǔ)存在網(wǎng)格節(jié)點(diǎn)上，網(wǎng)格在離散過程中起著關(guān)鍵的作用。網(wǎng)格劃分的形式及密度等將直接影響數(shù)值模擬的計(jì)算結(jié)果。因此，當(dāng)劃分網(wǎng)格時(shí)，網(wǎng)格的好壞會(huì)直接影響到計(jì)算問題的計(jì)算精度和收斂速度，同時(shí)關(guān)系到流場(chǎng)計(jì)算是否成功，其是否合理的布置對(duì)數(shù)值解的精度也有重要的影響。一般對(duì)網(wǎng)格劃分的要求

54、有如下幾點(diǎn)：（1）劃分的網(wǎng)格應(yīng)相適應(yīng)于求解區(qū)域所設(shè)置的邊界；（2）在高梯度及特定分析的區(qū)域劃分網(wǎng)格時(shí)需加密；（3）在邊界上劃分的網(wǎng)格線應(yīng)盡量正交于邊界，以便準(zhǔn)確方便的進(jìn)行邊界上的微分；（4）網(wǎng)格劃分時(shí)其間距在梯度變化時(shí)應(yīng)過渡較緩；（5）劃分網(wǎng)格線時(shí)應(yīng)盡可能與流體流動(dòng)方向重合。在數(shù)值模擬中，網(wǎng)格劃分是一個(gè)重要的前處理過程，當(dāng)計(jì)算域比較復(fù)雜時(shí)，生成網(wǎng)格所需時(shí)間也相對(duì)較長(zhǎng)。Gambit劃分網(wǎng)格軟件提供了多種網(wǎng)格單元，根據(jù)用戶的要求自動(dòng)完成劃分網(wǎng)格的工作，并且可以生成結(jié)構(gòu)網(wǎng)格、非結(jié)構(gòu)網(wǎng)格和混合網(wǎng)格等多種類型的網(wǎng)格，其能較的好匹配于Fluent數(shù)值模擬軟件，所以本文釆用Gambit軟件對(duì)離心泵進(jìn)行網(wǎng)格劃

55、分。目前，比較常用的網(wǎng)格有結(jié)構(gòu)化網(wǎng)格和非結(jié)構(gòu)化網(wǎng)格。（1）結(jié)構(gòu)化網(wǎng)格劃分的網(wǎng)格節(jié)點(diǎn)排列有序，當(dāng)給定了一個(gè)節(jié)點(diǎn)的編號(hào)時(shí)，立即可以得到其相鄰節(jié)點(diǎn)的編號(hào)，這種網(wǎng)格稱之為結(jié)構(gòu)化網(wǎng)格。它是一種傳統(tǒng)的網(wǎng)格形式，釆用結(jié)構(gòu)化網(wǎng)格有很多優(yōu)點(diǎn)：它可以很容易的實(shí)現(xiàn)區(qū)域的邊界擬合，適用于流體和表面應(yīng)力集中等方面的計(jì)算；網(wǎng)格生成的速度快；網(wǎng)格生成的質(zhì)量好；數(shù)據(jù)結(jié)構(gòu)簡(jiǎn)單；對(duì)曲面或空間的擬合大多數(shù)采用參數(shù)化或樣條插值的方法得到，區(qū)域光滑，與實(shí)際的模型更容易接近。（2）非結(jié)構(gòu)化網(wǎng)格非結(jié)構(gòu)化網(wǎng)格是指其網(wǎng)格的節(jié)點(diǎn)以一種不規(guī)則的方式布置在流場(chǎng)中。這種網(wǎng)格雖然生產(chǎn)過程比較復(fù)雜，但卻有著極大的適應(yīng)性，尤其對(duì)具有復(fù)雜邊界的流場(chǎng)計(jì)算問題特

56、別有效。其具體優(yōu)點(diǎn)是允許自由的劃分網(wǎng)格大小，劃分網(wǎng)格時(shí)比較靈活，實(shí)用性非常廣；它的缺點(diǎn)是劃分后的網(wǎng)格排列沒有規(guī)則性，需要尋找出各個(gè)單元之間的相關(guān)性，進(jìn)行數(shù)值模擬時(shí)有很大的計(jì)算量，計(jì)算所需的時(shí)間長(zhǎng)，占用計(jì)算機(jī)的內(nèi)存多。在本文計(jì)算中，由于離心泵的模型比較復(fù)雜，其中存在著很多不規(guī)則區(qū)域。因此，釆用非結(jié)構(gòu)化網(wǎng)格劃分模型網(wǎng)格，既能減少對(duì)網(wǎng)格劃分的工作量，又能很好的得到計(jì)算解。2.3.2 網(wǎng)格劃分本文釆用Gambit軟件對(duì)離心泵流道進(jìn)行網(wǎng)格的劃分。下面具體介紹網(wǎng)格劃分的步驟：（1） 導(dǎo)入建立好的離心泵實(shí)體模型108-16-zhuangpei_asm.stp并進(jìn)行檢查。（2）將三個(gè)區(qū)域分別進(jìn)行網(wǎng)格劃分，給出

57、相應(yīng)的網(wǎng)格尺寸和網(wǎng)格類型，整體采用四面體非結(jié)構(gòu)網(wǎng)格劃分網(wǎng)格。比如先對(duì)葉輪模型進(jìn)行劃分網(wǎng)格，選擇指定顯示屬性（Specify Display Attributes）工具里的體（Volumes）模塊將蝸殼和進(jìn)水管模型指定為不可見，單獨(dú)顯示葉輪模型。（3）使用合并命令，合并模型中的碎面。葉輪和蝸殼模型中的曲線和曲面數(shù)比較繁雜，先使用合并面（Merge Faces）命令來簡(jiǎn)化模型。另外對(duì)于某些狹長(zhǎng)的曲面，如果直接在上面網(wǎng)格劃分會(huì)導(dǎo)致極其尖銳的三角形網(wǎng)格。所以最好先進(jìn)行合并面命令來消除這類影響，當(dāng)然如此布置的面網(wǎng)格也相應(yīng)地降低模型的精度，某些細(xì)節(jié)地方會(huì)被抹平。（4）合并面命令完成后開始體網(wǎng)格的劃分。Mesh Volumes面板中的Elements命令提供了3種體網(wǎng)格劃分類型：Hex：六面體網(wǎng)格形式；Hex/Wedge：主要以六面體網(wǎng)格形式，在適當(dāng)位置包含楔形網(wǎng)格；Te