（機(jī)械電子工程專(zhuān)業(yè)論文）離心泵內(nèi)部流場(chǎng)的數(shù)值模擬與實(shí)驗(yàn)分析.pdf

分類(lèi)號(hào) 幽 u d c 工學(xué)碩士學(xué)位論文 離心泵內(nèi)部流場(chǎng)的數(shù)值模擬與實(shí)驗(yàn)分析 作者姓名 薛冰軍 指導(dǎo)教師 魏效玲教授 申請(qǐng)學(xué)位級(jí)別 工學(xué)碩士 學(xué)科專(zhuān)業(yè) 機(jī)械設(shè)計(jì)及理論 所在單位 機(jī)電工程學(xué)院 授予學(xué)位單位 河北工程大學(xué) ad i s s e r t a t i o ns u b m i t t e dt o h e b e iu n i v e r s i t yo fe n g i n e e r i n g f o rt h ea c a d e m i cd e g r e eo fm a s t e ro f e n g i n e e r i n g t h en u m e r i c a ls i m u l a t i o na n d e x p e r i m e n t a l a n a l y s i so fi n n e rf l o w f i e l do fc e n t r i f u g a l p u m p c a n d i d a t e s u p e r v i s o r a c a d e m i c d e g r e ea p p l i e df o r s p e c i a l t y c o l l e g e d e p a r t m e n t x u eb i n g i u n p r o f w e ix i a o l i n g m a s t e ro f e n g i n e e r i n g m e c h a n i c a ld e s i g na n dt h e o r y c o l l e g eo fm e c h a n i c a l e l e c t r i c a h e b e iu n i v e r s i t yo f e n g i n e e r i n g 一1 1一 j u n e 2 0 1 2 獨(dú)創(chuàng)性聲明 本人鄭重聲明 所呈交的學(xué)位論文 是本人在導(dǎo)師的指導(dǎo)下 獨(dú)立進(jìn)行研 究工作所取得的成果 除文中已經(jīng)注明引用的內(nèi)容外 本論文不含任何其他個(gè)人 或集體已經(jīng)發(fā)表或撰寫(xiě)過(guò)的研究成果 也不包含為獲得河北工程大堂或其他教 育機(jī)構(gòu)的學(xué)位或證書(shū)而使用過(guò)的材料 對(duì)本文的研究做出重要貢獻(xiàn)的個(gè)人和集體 均已在論文中作了明確的說(shuō)明并表示了謝意 本人完全意識(shí)到本聲明的法律結(jié)果 由本人承擔(dān) 學(xué)位論文作者簽名 謬劾夠竄 簽字同期 年 多月丁同 學(xué)位論文版權(quán)使用授權(quán)書(shū) 本學(xué)位論文作者完全了解塑蘭堡王猩盤(pán)鱟有關(guān)保留 使用學(xué)位論文的規(guī) 定 特授權(quán)塑皇堡三猩盤(pán)堂可以將學(xué)位論文的全部或部分內(nèi)容編入有關(guān)數(shù)據(jù)庫(kù) 進(jìn)行檢索 并采用影印 縮印或掃描等復(fù)制手段保存 匯編以供查閱和借閱 同 意學(xué)校向國(guó)家有關(guān)部門(mén)或機(jī)構(gòu)送交論文的復(fù)印件和電子文檔 保密的學(xué)位論文在解密后適用本授權(quán)說(shuō)明 學(xué)位論文作者簽名 櫪倜旆 簽字日期 力 年鈿曠日 導(dǎo)師簽名 粼途 簽字同期 如 1 年6 月r 同 摘要 摘要 離心泵是一種應(yīng)用較廣泛的能耗大效率低的設(shè)備 近年來(lái) 隨著計(jì)算機(jī)技術(shù) 的快速發(fā)展 利用f l u e n t 實(shí)現(xiàn)離心泵的內(nèi)部流場(chǎng)數(shù)值模擬為了解和掌握離心泵 復(fù)雜的內(nèi)部流動(dòng)開(kāi)辟了新的途徑 彌補(bǔ)了傳統(tǒng)設(shè)計(jì)方法的不足 對(duì)離心泵的優(yōu)化 設(shè)計(jì) 水力性能的改善 提高其效率具有非常重要的意義 本文首先采用p r o e 軟件建立離心泵的三維實(shí)體模型 然后采用計(jì)算流體動(dòng) 力學(xué)軟件f l u e n t 對(duì)r y 型離心泵進(jìn)行內(nèi)部流場(chǎng)的數(shù)值模擬 得到不同工況下離 心泵內(nèi)部流場(chǎng)不同時(shí)刻的速度分布和壓力分布 分析了離心泵內(nèi)部的流動(dòng)規(guī)律 并進(jìn)行性能預(yù)測(cè) 在此基礎(chǔ)上對(duì)r y 型離心泵進(jìn)行優(yōu)化改進(jìn) 并對(duì)優(yōu)化后的離心 泵進(jìn)行不同工況下內(nèi)部流場(chǎng)的數(shù)值模擬 通過(guò)計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果的比較驗(yàn)證該 方法計(jì)算離心泵內(nèi)部流動(dòng)的可行性和正確性 本文主要作了以下幾方面的研究工作 1 介紹了離心泵內(nèi)部流場(chǎng)數(shù)值模擬的國(guó)內(nèi)外發(fā)展研究現(xiàn)狀及商用c f d 軟 件的發(fā)展 從控制方程 湍流模型問(wèn)題的提出及解決方法 湍流模型的后處理等 多方面介紹了湍流流動(dòng)數(shù)值模擬理論 2 介紹了計(jì)算流體動(dòng)力學(xué)的基本理論知識(shí) 闡述了流體與流動(dòng)的基本特性 給出了流體動(dòng)力學(xué)的基本控制方程 3 n 用p r o e 軟件對(duì)r y 型離心泵過(guò)流部件進(jìn)行實(shí)體造型 建立離心泵的三 維實(shí)體模型 4 介紹了f l u e n t 環(huán)境下離心泵的數(shù)值計(jì)算過(guò)程 應(yīng)用滑移網(wǎng)格法對(duì)離心泵 內(nèi)部湍流流動(dòng)進(jìn)行非定常數(shù)植模擬 得到不同工況下各個(gè)時(shí)刻的速度分布和壓力 分布 通過(guò)計(jì)算結(jié)果比較分析 揭示了離心泵內(nèi)部的流動(dòng)規(guī)律 5 根據(jù)數(shù)值模擬結(jié)果對(duì)現(xiàn)有r y 型離心泵進(jìn)行優(yōu)化設(shè)計(jì) 并以此為基礎(chǔ)進(jìn)行 離心泵的性能預(yù)測(cè) 將預(yù)測(cè)結(jié)果與已有的實(shí)驗(yàn)測(cè)試結(jié)果進(jìn)行分析比較 驗(yàn)證了所 采用方法的可行性和 f 確性 關(guān)鍵詞 離心泵 三維建模 計(jì)算流體動(dòng)力學(xué) 流場(chǎng)分析 數(shù)值模擬 滑移網(wǎng)格 a b s t r a c t a b s t r a c t c e n t r i f u g a lp u m p i sw i d e l ya p p l i e de q u i p m e n tw h i c hn e e d sm o r ee n e r g yb u th a s l o we f f i c i e n c y i nr e c e n ty e a r s a l o n gw i t ht h er a p i dd e v e l o p m e n to fc o m p u t e r t e c h n o l o g y w ec a na c h i e v ei n t e m a lf l o wf i e l dv a l u es i m u l a t i o no fc e n t r i f u g a lp u m pb y f l u e n t i to p e n san e wa v e n u em a s t e r i n gi n t e m a lf l o wf i e l do fc e n t r i f u g a lp u m pa n d m a k e su pf o r t h el a c ko fd e s i g n i th a sv e r yi m p o r t a n ts i g n i f i c a n c eo nt h eo p t i m i z a t i o n d e s i g no ft h ec e n t r i f u g a lp u m p h y d r a u l i cp e r f o r m a n c ei m p r o v e m e n t i m p r o v i n g e f f i c i e n c y t h i sp a p e rf i r s t l ys e tu pa3 de n t i t ym o d e lb a s e do nt h ep r o e t h e nt h en u m e r i c a l s i m u l a t i o no ft h ec e n t r i f u g a lp u m pf o ri n t e r n a lf l o wf i e l db yu s i n gt h es o f t w a r eo f f l u e n tc a ng e tt h ec e n t r i f u g a lp u m pi n t e m a lf l o wf i e l do fv e l o c i t yd i s t r i b u t i o na n d s t r e s sd i s t r i b u t i o no fd i f f e r e n tw o r k i n gc o n d i t i o n sa n dt i m e t h r o u g ht h i s w ec a n a n a l y z et h ei n t e m a ll a wo ft h ec e n t r i f u g a lp u m pa n dc a l c u l a t et h ep e r f o r m a n c e b a s e d o ni t w eo p t i m i z er yc e n t r i f u g a lp u m pa n dp e r f o r mt h en u m e r i c a ls i m u l a t i o no ft h e f l o wf i e l du n d e rd i f f e r e n tw o r k i n gc o n d i t i o n s t h r o u g ht h ec o m p a r i s o nb e t w e e nt h e c a l c u l a t i o nr e s u l t sa n de x p e r i m e n t a lr e s u l t s v e r i f y i n gt h ec o r r e c t n e s sa n df e a s i b i l i t yo f t h ew a yt oc a l c u l a t ei n t e m a lf l o wo fc e n t r i f u g a lp u m p t h i sp a p e rm a i n l yd o e st h ef o l l o w i n gj o b s 1 i n t r o d u c i n gc u r r e n tr e s e a r c hd e v e l o p m e n ta th o m ea n d a b r o a da b o u tc e n t r i f u g a l p u m p i ta l s oi n t r o d u c e st h ed e v e l o p m e n to ft h es o f t w a r eo nc f d a n dt h et u r b u l e n tf l o w n u m e r i c a ls i m u l a t i o nt h e o r yb a s e do nt h ec o n t r o le q u a t i o n s p u t t i n gf o r w a r dt u r b u l e n c e m o d e lo ft h ep r o b l e ma n dt h es o l v i n gm e t h o d s t u r b u l e n tm o d e lp o s t p r o c e s s i n ga n ds o o n 2 i n t r o d u c i n gt h eb a s i ct h e o r yk n o w l e d g eo fc o m p u t a t i o n a l f l u i dd y n a m i c s e x p o u n d i n gb a s i cc h a r a c t e r i s t i c so ft h ef l o wo ff l u i d a n dg i v i n gt h eb a s i cc o n t r o l e q u a t i o n so ft h ef l o wo ff l u i d 3 e s t a b l i s h i n gt h et h r e e d i m e n s i o n a le n t i t ym o d e lo fc e n t r i f u g a lp u m pu s i n g p r o e 4 i n t r o d u c i n gn u m e r i c a lc a l c u l a t i o no fc e n t r i f u g a lp u m pb a s e do nf l u e n t w e c a ng e tv e l o c i t yd i s t r i b u t i o na n ds t r e s sd i s t r i b u t i o no fd i f f e r e n tw o r k i n gc o n d i t i o n sa n d t i m eb yp e r f o r m i n gt h en u m e r i c a ls i m u l a t i o no ft h ef l o wf i e l du s i n gs l i d i n gg r i d ii a b s t r a c t t e c h n o l o g y t h r o u g ht h ef l o wf i e l dc a l c u l a t i o nr e s u l t s w ec a ng e tf l o wl a w 5 b a s e do nt h er e s u l to fn u m e r i c a ls i m u l a t i o n w eo p t i m i z er yc e n t r i f u g a lp u m p a n dt h e nt e s ti t s p e r f o r m a n c e a n a l y s i s i n gf o r e c a s t e dr e s u l t sa n dt e s t r e s u l t st o d e m o n s t r a t et h ec o r r e c t n e s sa n df e a s i b i l i t yo ft h em e t h o d k e y w o r d s c e n t r i f u g a lp u m p 3 dm o d e l i n g c o m p u t a t i o n a lf l u i dd y n a m i c s f o wf i e l d a n a l y s i s n u m e r i c a ls i m u l a t i o n s l i d i n gg r i d 目錄 目錄 摘要 i a b s t r a c t i i 第1 章緒論 1 1 1 本課題的研究背景及研究意義 1 1 2 國(guó)內(nèi)外相關(guān)技術(shù)的發(fā)展現(xiàn)狀 2 1 2 1 國(guó)外相關(guān)技術(shù)發(fā)展現(xiàn)狀 2 1 2 2 國(guó)內(nèi)相關(guān)技術(shù)發(fā)展現(xiàn)狀 2 1 3 離心泵內(nèi)部流場(chǎng)數(shù)值模擬的研究現(xiàn)狀及發(fā)展趨勢(shì) 3 1 4 計(jì)算流體動(dòng)力學(xué) c f d 技術(shù)的發(fā)展 6 1 5 本課題研究?jī)?nèi)容 6 1 6 本章小結(jié) 7 第2 章離心泵內(nèi)部流場(chǎng)數(shù)值模擬理論基礎(chǔ) 8 2 1 計(jì)算流體動(dòng)力學(xué)理論基礎(chǔ) 8 2 1 1 計(jì)算流體動(dòng)力學(xué)基本思想 8 2 1 2 流體與流動(dòng)的基本特性 9 2 1 3 流體動(dòng)力學(xué)的基本控制方程 1 0 2 2 湍流流動(dòng)的數(shù)值模擬方法 1 3 2 3 常見(jiàn)湍流模型及應(yīng)用分析 14 2 4 離心泵內(nèi)部流場(chǎng)數(shù)值模擬方法 2 0 2 5 本章小結(jié) 2 2 第3 章離心泵的三維實(shí)體建模 2 3 3 1 p r o e 造型軟件介紹 2 3 3 2 葉輪三維造型設(shè)計(jì) 2 4 3 3 蝸殼三維造型設(shè)計(jì) 2 7 3 4 本章小結(jié) 3 0 第4 章離心泵內(nèi)部流場(chǎng)的數(shù)值模擬 3 l 4 1f l u e n t 軟件介紹 31 4 2 幾何模型的簡(jiǎn)化生成與網(wǎng)格的劃分 3 2 4 2 1 幾何模型簡(jiǎn)化與生成 3 3 4 2 2 網(wǎng)格的劃分 3 3 4 3 旋轉(zhuǎn)葉輪和靜止蝸殼間的耦合模型 3 4 4 4 邊界條件 3 5 4 4 1 進(jìn)口邊界條件 3 6 4 4 2 出口邊界條件 3 7 4 4 3 固體壁面邊界條件 3 7 4 5 離心泵內(nèi)部流場(chǎng)的數(shù)值模擬計(jì)算 3 8 4 5 1 設(shè)置計(jì)算模型 3 8 4 5 2 設(shè)置邊界條件 3 8 4 5 3 設(shè)置求解控制參數(shù) 3 8 4 5 4 迭代計(jì)算 3 9 4 6 迭代結(jié)果與分析 3 9 目錄 4 7 本章小結(jié) 一4 2 第5 章離心泵數(shù)值模擬結(jié)果與實(shí)驗(yàn)值比較 4 3 5 1 離心泵性能實(shí)驗(yàn) 4 3 5 1 1 基本原理 4 3 5 1 2 實(shí)驗(yàn)數(shù)據(jù) 4 4 5 1 3 繪制離心泵工作性能曲線(xiàn) 4 6 5 2 優(yōu)化后離心泵數(shù)值模擬 4 7 5 2 1 幾何模型的建立 4 7 5 2 2 數(shù)值模擬結(jié)果分析 4 8 5 3 優(yōu)化前后數(shù)值模擬結(jié)果比較 4 9 5 3 1 同一時(shí)刻的速度分布比較 4 9 5 3 2 同一時(shí)刻的壓力分布比較 5 0 5 3 3 同一時(shí)刻的湍動(dòng)能分布比較 5 1 5 4 性能預(yù)測(cè)與實(shí)驗(yàn)值的比較 5 2 5 4 1f l u e n t 數(shù)值報(bào)告介紹 5 2 5 4 2 揚(yáng)程的預(yù)測(cè) 5 3 5 4 3 效率的預(yù)測(cè) 5 3 5 4 4 預(yù)測(cè)值與實(shí)驗(yàn)值的比較 一 5 4 5 5 本章小結(jié) 5 6 結(jié)論 5 7 致謝 5 9 參考文獻(xiàn) 6 0 作者簡(jiǎn)介 6 3 攻讀碩士學(xué)位期間發(fā)表的學(xué)術(shù)論文 6 3 第1 章緒論 第1 章緒論 泵與風(fēng)機(jī)屬于通用機(jī)械的范疇 在石油 化工 造船 水利 輕工及電力等國(guó)民 經(jīng)濟(jì)的各個(gè)領(lǐng)域中有著廣泛的應(yīng)用 同時(shí) 泵與風(fēng)機(jī)又是消耗能量最多的通用機(jī)械之 一 據(jù)統(tǒng)計(jì) 我國(guó)水泵與風(fēng)機(jī)裝機(jī)總功率已達(dá)1 5 9 億千瓦 其中風(fēng)機(jī)0 4 9 億千瓦 水 泵1 1 億千瓦 年總耗電量達(dá)到3 2 0 0 億千瓦時(shí) 大約消耗當(dāng)年全國(guó)電力消耗總量的 1 3 占工業(yè)用電量的4 0 我國(guó)泵與風(fēng)機(jī)的生產(chǎn)廠家眾多 產(chǎn)品質(zhì)量良莠不齊 實(shí)際 運(yùn)行效率更低 僅為4 0 左右 據(jù)有關(guān)部門(mén)計(jì)算 如果能夠提高泵與風(fēng)機(jī)的設(shè)計(jì)效率 和運(yùn)行效率 僅在電廠的風(fēng)機(jī)和水泵節(jié)約電量就達(dá)到2 5 6 9 億千瓦時(shí) 年 因此 改進(jìn) 和提高泵與風(fēng)機(jī)的設(shè)計(jì)和運(yùn)行效率能夠節(jié)省大量的能源 對(duì)社會(huì)經(jīng)濟(jì)的可持續(xù)發(fā)展具 有十分重要的意義u 離心泵結(jié)構(gòu)的復(fù)雜性決定了其內(nèi)部流場(chǎng)的分布相當(dāng)復(fù)雜 這也給離心泵內(nèi)部流動(dòng) 情況的分析研究帶來(lái)了一定困難 近年來(lái) 國(guó)內(nèi)外對(duì)這方面的研究取得了一定成果 最早的是采用實(shí)驗(yàn)方法 測(cè)得流場(chǎng)的分布情況 從而判斷離心泵結(jié)構(gòu)是否合理 近幾 年來(lái) l d v p d v p l v p i v 等技術(shù)的發(fā)展進(jìn)一步豐富了離心泵流場(chǎng)的測(cè)量 但這 些測(cè)試手段需要花費(fèi)大量地人力物力并且需要熟練的實(shí)驗(yàn)技術(shù) 離心泵內(nèi)部流場(chǎng)數(shù)值 模擬的研究已經(jīng)變得相當(dāng)活躍 實(shí)現(xiàn)離心泵的內(nèi)部流場(chǎng)數(shù)值模擬對(duì)于離心泵的優(yōu)化設(shè) 計(jì) 改善離心泵的水力性能以達(dá)到增效節(jié)能的目的具有重要的現(xiàn)實(shí)意義 而且 基于 c f d 數(shù)值模擬的現(xiàn)代設(shè)計(jì)方法還可以彌補(bǔ)傳統(tǒng)設(shè)計(jì)方法的不足從而可以設(shè)計(jì)出高效 節(jié)能的離心泵 這些都將對(duì)泵與風(fēng)機(jī)的設(shè)計(jì)研發(fā)帶來(lái)積極的影響 1 1 本課題的研究背景及研究意義 在石油工業(yè)中 離心泵是油田原油短距離集中輸送 石油管道長(zhǎng)距離輸送和石油 化工流程中的重要流體輸送設(shè)備 具有結(jié)構(gòu)簡(jiǎn)單 操作方便 運(yùn)轉(zhuǎn)平穩(wěn) 便于維護(hù)等 優(yōu)點(diǎn) 因此 它在石油礦場(chǎng)上 儲(chǔ)運(yùn)過(guò)程中和煉油廠內(nèi)應(yīng)用廣泛 石油 石油化工行 業(yè)中所使用的離心泵大多為蝸殼式離心泵 其中中低比轉(zhuǎn)速泵占很大的比例 另外 離心泵還具有壓力高 功率大 耗電高等參數(shù)特點(diǎn) 其中 離心泵單機(jī)消耗最大功率 已接近1 0 0 0 k w 啦 據(jù)統(tǒng)計(jì) 目前 油田泵的耗電量約占油田總耗電電量的1 2 其中 離心泵的消耗占有很大比重 通過(guò)對(duì)大慶油田輸油泵的調(diào)查發(fā)現(xiàn) 對(duì)任意選定的十臺(tái)離 心泵進(jìn)行調(diào)查 分析 計(jì)算 結(jié)果表明平均每臺(tái)離心泵浪費(fèi)功率8 1 2 k w 油田僅此一 項(xiàng) 按保守的估算每年消耗于出口節(jié)流閥門(mén)的功率就達(dá)4 5 0 0 萬(wàn)度電 隨著現(xiàn)代石油 工業(yè)的快速發(fā)展 離心泵的重要性越來(lái)越明顯 但是目前市場(chǎng) lj 的離心泵存在 個(gè)普 河北工程大學(xué)碩士學(xué)位論文 遍的閥題就是泵的效率低下 因此 提高離心泵運(yùn)行效率 降低離心泵電能消耗對(duì)提 i 高企業(yè)經(jīng)濟(jì)效益 國(guó)民經(jīng)濟(jì)的可持續(xù)發(fā)展都有十分重要的現(xiàn)實(shí)意義 蔓曩 j 1 2 國(guó)內(nèi)外相關(guān)技術(shù)的發(fā)展現(xiàn)狀 1 2 1 國(guó)外相關(guān)技術(shù)發(fā)展現(xiàn)狀 早在19 2 0 年左右 就有學(xué)者關(guān)注于離心泵內(nèi)部流場(chǎng)研究 f i s c h e r 和t h o m a 用傳 統(tǒng)的流場(chǎng)顯示技術(shù)證明了離心泵葉輪內(nèi)部的真實(shí)流動(dòng)與基于理想流體的流動(dòng)分析 有明顯差別n 3 1 9 3 3 年 英國(guó)科學(xué)家t h o r n 應(yīng)用手搖計(jì)算機(jī)完成了對(duì)一個(gè)外掠圓柱流 動(dòng)的數(shù)值計(jì)算嘲 1 9 5 7 年 a c o s t a 和b o w e r m a n 研究表明 在相當(dāng)大的范圍內(nèi) 葉片的 吸力面?zhèn)仁潜萌~輪內(nèi)部流動(dòng)損失相對(duì)集中的區(qū)域嘲 f o w l e r 于2 0 世紀(jì)6 0 年代發(fā)現(xiàn) 了流道出口附近壓力面?zhèn)攘魉倜黠@大于吸力面?zhèn)?與勢(shì)流預(yù)測(cè)結(jié)論正好相反h 3 j o h n s o n 和m o o r e 實(shí)驗(yàn)研究了離心泵葉輪流量對(duì)尾跡的影響嘲 a m d t 等實(shí)驗(yàn)研究了 離心泵葉輪與蝸殼間的流動(dòng)干擾 發(fā)現(xiàn)壓力面的壓力波動(dòng)小于吸力面 最大的壓力 波動(dòng)出現(xiàn)在吸力面進(jìn)口邊嘲 s t o f f e l 等人先后研究了蝸殼對(duì)葉輪內(nèi)流的影響 實(shí)驗(yàn) 結(jié)果表明 蝸殼與葉輪間的徑向間隙大到一定值后 蝸殼對(duì)葉輪內(nèi)流的影響不再顯 著 但葉輪出流仍不對(duì)稱(chēng)n0 m 觚i s hs i n h a 等利用p i v 技術(shù)對(duì)帶有葉輪擴(kuò)壓器的離 心泵內(nèi)部的轉(zhuǎn)子與靜子的相互干涉 失速現(xiàn)象進(jìn)行了分析 計(jì)算和實(shí)驗(yàn)研究u 1 1 1 2 2 國(guó)內(nèi)相關(guān)技術(shù)發(fā)展現(xiàn)狀 國(guó)內(nèi)科研人員也離心泵內(nèi)部流動(dòng)測(cè)量進(jìn)行了分析研究 吳仲華教授于1 9 5 2 年提 出了基于兩類(lèi)相對(duì)流面的通用理論 在工業(yè)界被廣泛應(yīng)用 引 這大大影響了葉輪機(jī)械 內(nèi)部流動(dòng)的數(shù)值模擬研究 薛敦松等用l d v 測(cè)量葉輪內(nèi)流動(dòng)實(shí)驗(yàn)表明較大的包角比 較小的出口角對(duì)內(nèi)流均壓化有利 3 1 李森虎等利用閃光測(cè)速儀研究了離心泵和誘導(dǎo) 輪內(nèi)部汽蝕空泡的產(chǎn)生和發(fā)展n 李文廣等利用l d v 對(duì)離心泵內(nèi)部流動(dòng)測(cè)量實(shí)驗(yàn) 實(shí)驗(yàn)表明離心泵葉輪出口處流動(dòng)是主流一尾流模型n 置1 6 1 f l u e n t 軟件廣泛應(yīng)用以后 第1 章緒論 國(guó)內(nèi)很多學(xué)者對(duì)水力機(jī)械的內(nèi)部流場(chǎng)進(jìn)行了分析研究 國(guó)內(nèi)的戴江 吳玉林等計(jì)算過(guò) 泵葉輪內(nèi)的二維或三維粘性流動(dòng) 7 18 江蘇大學(xué)袁壽其計(jì)算過(guò)低比轉(zhuǎn)速離心泵葉 輪內(nèi)的流動(dòng)n9 李海峰 吳玉林運(yùn)用壓力修正法對(duì)兩種離心泵的設(shè)計(jì)方案葉輪內(nèi)紊 流進(jìn)行了數(shù)值模擬 通過(guò)把計(jì)算結(jié)果和實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比 可以為葉輪方案的比 較和選擇提供有價(jià)值的信息圓1 國(guó)內(nèi)外的很多專(zhuān)家學(xué)者對(duì)離心泵的內(nèi)部流場(chǎng)進(jìn)行了大量分析研究 大多對(duì)離心泵 內(nèi)部流場(chǎng)進(jìn)行了定常計(jì)算 從而研究流體的壓力 速度分布以及流動(dòng)中的各種現(xiàn)象 而采用非定常計(jì)算的文獻(xiàn)相對(duì)較少 本文對(duì)r y 型離心泵進(jìn)行非定常分析計(jì)算 可以 得到離心泵的瞬時(shí)流場(chǎng)分布 1 3 離心泵內(nèi)部流場(chǎng)數(shù)值模擬的研究現(xiàn)狀及發(fā)展趨勢(shì) 計(jì)算流體動(dòng)力學(xué)隨著計(jì)算機(jī)技術(shù)的不斷更新而飛速發(fā)展 同時(shí)由于水力機(jī)械葉輪 內(nèi)部流場(chǎng)復(fù)雜多變不易用實(shí)驗(yàn)手段測(cè)量 從而使得泵葉輪內(nèi)部流場(chǎng)的數(shù)值模擬技術(shù)得 到了更大的發(fā)展 特別是近些年來(lái) 葉輪機(jī)械的數(shù)值模擬已由原來(lái)的無(wú)粘性發(fā)展到考 慮流體粘性 原來(lái)的二維流場(chǎng)分析到三維流場(chǎng)數(shù)值模擬 c f d 作為主要流體分析方法 能推出多種優(yōu)化的物理模型 如定常流動(dòng)與非定常流動(dòng) 不可壓縮流動(dòng)與可壓縮流動(dòng) 紊流 層流 化學(xué)反應(yīng) 傳熱等 應(yīng)用c f d 軟件可以在計(jì)算機(jī)上真實(shí)反映離心泵內(nèi)流 體復(fù)雜的二維及三維流動(dòng)情況 從而能分析泵內(nèi)速度及壓力等的分布 為泵的優(yōu)化設(shè) 計(jì) 改進(jìn)提供理論依據(jù) 隨著國(guó)民經(jīng)濟(jì)的發(fā)展 工業(yè)生產(chǎn)對(duì)葉輪機(jī)械的設(shè)計(jì)和研發(fā)提出了更加嚴(yán)格的要 求 只有研究了離心泵內(nèi)部流場(chǎng)分布情況和能量損失產(chǎn)生的原因 才能對(duì)葉輪機(jī)械的 結(jié)構(gòu)進(jìn)行優(yōu)化設(shè)計(jì) 從而提高葉輪機(jī)械性能 因此提高葉輪機(jī)械設(shè)計(jì)水平的關(guān)鍵是研 究葉輪機(jī)械的內(nèi)部流場(chǎng) 目前理論分析 實(shí)驗(yàn)研究和數(shù)值模擬是用于研究分析葉輪機(jī) 械內(nèi)部流場(chǎng)的主要方法 數(shù)值模擬在葉輪機(jī)械的發(fā)展中起著越來(lái)越重要的作用 其顯 著特點(diǎn)是高效率 小投資和周期短等 葉輪機(jī)械內(nèi)部流場(chǎng)的數(shù)值模擬已逐漸發(fā)展為計(jì) 算流體動(dòng)力學(xué)重要的應(yīng)用領(lǐng)域 葉輪機(jī)械內(nèi)部流場(chǎng)的數(shù)值模擬發(fā)展時(shí)期如下晗 1 無(wú)粘流數(shù)值模擬時(shí)期 無(wú)粘性數(shù)值模擬主要集中在上世紀(jì)的5 0 年代到8 0 年代 該時(shí)期數(shù)值模擬計(jì)算以 e u l e r 方程 流函數(shù)或勢(shì)函數(shù)作為控制方程進(jìn)行求解 其中 吳仲華教授于1 9 5 2 年提 出了s l s 2 兩類(lèi)相對(duì)流面理論 這大大影響了葉輪機(jī)械內(nèi)部流動(dòng)的數(shù)值模擬研究陛2 此后 在普遍采用這種理論來(lái)模擬計(jì)算葉輪內(nèi)部流動(dòng)的同時(shí) 產(chǎn)生了 些新的數(shù)值模 擬方法 比如準(zhǔn) f 交面法 流線(xiàn)曲率法等 河北工程大學(xué)碩士學(xué)位論文 1 流線(xiàn)曲率法 其基本思想是 在泵的流道內(nèi)假定近似的流線(xiàn) 借助正交線(xiàn)上的 常微分速度梯度方程求解流動(dòng)量新的迭代值 再依據(jù)流量等值反插方法調(diào)整各流線(xiàn)的 位置 迭代收斂直至達(dá)到給定的精度 陳勝利 吳達(dá)人分析得出了離心泵葉輪內(nèi)平均 相對(duì)流面的流線(xiàn)分布位置及速度分布 最后借助快速近似法得出了葉片表面的速度 壓力分布 最終實(shí)現(xiàn)了離心泵葉輪內(nèi)部流場(chǎng)的數(shù)值模擬瞳3 2 準(zhǔn)正交面法 其基本思想是 在泵葉輪流道中從各個(gè)準(zhǔn)正交面上分別沿s 1 s 2 流面與準(zhǔn)正交面的交線(xiàn)求解速度梯度方程的積分 得到新的相對(duì)速度w 值 通過(guò)準(zhǔn)正 交面流量相等原則修正相對(duì)速度w 值 積分之前在流道內(nèi)選定初始流線(xiàn)節(jié)點(diǎn)和相對(duì)速 度w 的分布 接著反復(fù)迭代直至得到滿(mǎn)足流量條件的相對(duì)速度w 分布 然后反插等分 流量點(diǎn) 求得流線(xiàn)節(jié)點(diǎn)新坐標(biāo) 再求新的相對(duì)速度w 分布和新的流線(xiàn)節(jié)點(diǎn) 反復(fù)迭代 至相對(duì)速度w 收斂 2 有粘流數(shù)值模擬時(shí)期 上世紀(jì)8 0 年代到9 0 年代 葉輪機(jī)械內(nèi)部流場(chǎng)的數(shù)值模擬開(kāi)始綜合考慮內(nèi)部流場(chǎng) 的粘性 回流及旋渦對(duì)內(nèi)部流場(chǎng)的影響 更復(fù)雜的數(shù)值模擬計(jì)算方法隨著計(jì)算機(jī)技術(shù) 的快速發(fā)展也開(kāi)始出現(xiàn) 主要方法有射流一尾流模型 渦量一流函數(shù)法 勢(shì)流一邊界層的 迭代解法等瞰1 這些方法在近期國(guó)內(nèi)外仍有著比較廣泛的應(yīng)用 特別是勢(shì)流一邊界層迭 代求解法 許多學(xué)者對(duì)于該方法在葉輪內(nèi)部流場(chǎng)計(jì)算中的應(yīng)用也作了不少研究分析 2 5 2 7 o 3 完全粘性流數(shù)值模擬時(shí)期 1 9 9 0 年開(kāi)始 很多學(xué)者對(duì)葉輪內(nèi)部流場(chǎng)計(jì)算方法的研究開(kāi)始考慮湍流運(yùn)動(dòng)的影 響 由此離心泵葉輪內(nèi)部流場(chǎng)數(shù)值模擬開(kāi)始進(jìn)入三維粘性數(shù)值模擬時(shí)期 該時(shí)期葉輪 內(nèi)部流場(chǎng)數(shù)值模擬的主流是直接求解雷諾時(shí)均化的n a v i e r s t o k e s 方程組 然后結(jié)合湍 流模型計(jì)算葉輪內(nèi)的三維粘性流動(dòng) 該時(shí)期數(shù)值模擬方法主要有 1 壓力修正法 壓力修正法最早是由c h o r i n 于1 9 6 7 年和1 9 6 8 年提出的瞳8 馴 基本思想是 對(duì)給定 的壓力場(chǎng)按順序求解速度代數(shù)方程 得到的速度場(chǎng)不一定滿(mǎn)足連續(xù)性方程 這就需要 對(duì)給定壓力場(chǎng)進(jìn)行修正 為此將動(dòng)量方程的離散形式所規(guī)定的壓力與速度的關(guān)系代入 連續(xù)性方程的離散形式 得到壓力修 f 方程 根據(jù)壓力修 f 方程得出壓力修 f 值 進(jìn) 而去改進(jìn)速度 最終得到這一迭代層次上滿(mǎn)足連續(xù)性方程的解 然后用新的速度值去 修正動(dòng)量離散方程的系數(shù) 開(kāi)始下一迭代層次運(yùn)算 反復(fù)直至收斂 該方法自其問(wèn)世 以來(lái)在葉輪機(jī)械的內(nèi)部流場(chǎng)計(jì)算中得到了極廣泛的應(yīng)用 江蘇大學(xué)袁壽其借助壓力修 f 法分析研究了過(guò)低比轉(zhuǎn)速離心泵葉輪內(nèi)的流動(dòng) 李海峰 吳玉林運(yùn)用壓力修 f 法 對(duì)兩種離心泵日十輪內(nèi)紊流進(jìn)行了數(shù)值模擬 根據(jù)訓(xùn) 算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)的對(duì)比 為葉輪 方案的設(shè)計(jì)和比較選擇提供了有價(jià)值的理論信息瞳 2 第1 章緒論 2 時(shí)間相關(guān)法 時(shí)間相關(guān)法分為顯式和隱式兩種 顯式法優(yōu)點(diǎn)是計(jì)算簡(jiǎn)單 便于矢量化 邊界條 件也容易實(shí)現(xiàn) 但條件穩(wěn)定限制 時(shí)間步長(zhǎng)限制 限制了它的應(yīng)用范圍 隱式法允許采 用的時(shí)間步長(zhǎng)較大 收斂速度也比較快 對(duì)穩(wěn)定及非穩(wěn)定問(wèn)題的求解都很有效 3 擬可壓縮方法 葉輪機(jī)械的內(nèi)部流場(chǎng)介質(zhì)一般為不可壓縮介質(zhì)或介質(zhì)可壓縮性很小 甚至可以忽 略 不可壓縮流動(dòng)的n a v i e r s t o k e s 方程為橢圓型方程 該方程計(jì)算速度比較慢 不能 使用推進(jìn)法 為此可以通過(guò)數(shù)學(xué)處理把橢圓型方程轉(zhuǎn)變成雙曲線(xiàn)型方程 雙曲線(xiàn)型方 程格式多 計(jì)算精確 效率高 而且可以使用時(shí)間推進(jìn)法求解 而擬可壓縮方法采用 高度發(fā)達(dá)的時(shí)間推進(jìn)法求解 求解速度較壓力修正法快 可以求解分離流 大流動(dòng)梯 度等高復(fù)雜流動(dòng)問(wèn)題 但是 擬可壓縮方法的最大缺點(diǎn)是其穩(wěn)定性和收斂性對(duì)參數(shù)設(shè) 置太敏感 需要比較高的存儲(chǔ)能力 在近3 0 年內(nèi) p h o e n i c s c f d f l u e n t 等大型商用計(jì)算軟件相繼問(wèn)世 已 經(jīng)廣泛應(yīng)用于包括航空航天領(lǐng)域 離心泵在內(nèi)的諸多領(lǐng)域 李文廣應(yīng)用p h o e n i c s 軟件對(duì)離心泵葉輪內(nèi)部流場(chǎng)進(jìn)行了三維紊流計(jì)算 對(duì)比分析了計(jì)算結(jié)果與l d v 的測(cè) 量結(jié)果跚1 清華大學(xué)唐輝 何楓使用f l u e n t 軟件對(duì)離心泵內(nèi)部流場(chǎng)進(jìn)行了數(shù)值模擬 比較了三種計(jì)算模型的數(shù)值模擬結(jié)果 3 l 譚宗柒等使用c f d 采用多重坐標(biāo)系法 m u l t i p l er e f e r e n c ef r a m e m r f 模擬了離心泵內(nèi)流體的二維流動(dòng)情況口2 1 呂培文 以煤漿離心泵內(nèi)部流場(chǎng)為研究對(duì)象 采用m r f 方法對(duì)離心泵內(nèi)部流場(chǎng)三維流動(dòng)進(jìn)行了 定常數(shù)值計(jì)算 得到了內(nèi)部流場(chǎng)的流動(dòng)速度 壓力分布圖并計(jì)算預(yù)測(cè)了離心泵揚(yáng)程口3 張崢采用隱式修正s m p l e c 算法 使用m r f 方法對(duì)i b 型離心泵的內(nèi)部流場(chǎng)進(jìn)行了 三維數(shù)值模擬 并根據(jù)分析結(jié)果分析總結(jié)了離心泵內(nèi)的流體流動(dòng)規(guī)律 揭示了離心泵 內(nèi)流動(dòng)存在的漩渦 二次流等現(xiàn)象 圳 借助c f d 流體數(shù)值分析軟件 工作人員可以有 更多的精力考慮所計(jì)算的內(nèi)部流場(chǎng)問(wèn)題的物理本質(zhì) 邊界條件的設(shè)置以及計(jì)算結(jié)果是 否合理等方面問(wèn)題 離心泵葉輪內(nèi)部流場(chǎng)的研究發(fā)展趨勢(shì)如下 1 進(jìn)一步研究具體應(yīng)用情況下的數(shù)學(xué)模型 在標(biāo)準(zhǔn)湍流模型的基礎(chǔ)上 著重考 慮離心泵葉輪的曲率和旋轉(zhuǎn)對(duì)流動(dòng)的影響 進(jìn)一步尋求適合于離心泵葉輪內(nèi)流計(jì)算的 數(shù)學(xué)模型 迄今還不存在通用的模型 因此不同的情況需要不同的湍流模型 2 1 優(yōu)化設(shè)計(jì) 研究離心泵內(nèi)部流場(chǎng)的最終目的是為了提高離心泵的使用效率 節(jié) 省能源 3 計(jì)算網(wǎng)格劃分的進(jìn)一步研究 在數(shù)值模擬計(jì)算過(guò)程中 網(wǎng)格生成占用計(jì)算任務(wù) 的大部分時(shí)問(wèn) 而c f d 計(jì)算對(duì)計(jì)算網(wǎng)格的疏密程度要求很高 因此進(jìn)一步研究討 算網(wǎng) 格自動(dòng)或半自動(dòng)的生成方法 可以進(jìn)一步提高復(fù)雜邊界的網(wǎng)格生成質(zhì)鼉 河北工程大學(xué)碩士學(xué)位論文 4 更有效的算法研究 目前算法研究的重點(diǎn)之一尋求更有效地算法 1 4 計(jì)算流體動(dòng)力學(xué) c f d 技術(shù)的發(fā)展 計(jì)算流體動(dòng)力學(xué) c o m p u t a t i o n a lf l u i dd y n a m i c s 簡(jiǎn)稱(chēng)c f d 是近代流體動(dòng)力學(xué) 數(shù)值計(jì)算技術(shù)和計(jì)算機(jī)科學(xué)技術(shù)相結(jié)合的產(chǎn)物 主要作用是借助計(jì)算機(jī)對(duì)含有流體流 動(dòng) 熱傳遞等物理現(xiàn)象的機(jī)械系統(tǒng)數(shù)值進(jìn)行數(shù)值模擬 流動(dòng)虛擬仿真并最終在計(jì)算機(jī) 上實(shí)現(xiàn)流場(chǎng)可視化 1 9 4 6 年 美國(guó)賓夕法尼亞大學(xué)研制成功了世界上第一臺(tái)電子計(jì)算機(jī)e n i a c 與 此同時(shí) 計(jì)算機(jī)之父 v o n n e u m a n n 預(yù)言 解決流體非線(xiàn)性問(wèn)題的常用方 法一解析方法將被數(shù)值方法取代 在此后的十多年中 眾多的專(zhuān)家學(xué)者投入了大量的 精力在多方面展開(kāi)了對(duì)數(shù)值方法和理論的研究 涌現(xiàn)出了大量的算法 為計(jì)算流體動(dòng) 力學(xué)的形成打下了堅(jiān)實(shí)的理論基礎(chǔ) 1 9 6 0 年以后 大容量 高速度 多功能計(jì)算機(jī)的研制及其推廣應(yīng)用 尤其是計(jì)算 技術(shù)的發(fā)展 為計(jì)算流體動(dòng)力學(xué)的形成創(chuàng)造了有利條件 在此之后 相繼產(chǎn)生了不少 流體動(dòng)力學(xué)的數(shù)值方法 主要有有限元法 邊界元法 裝配法 算子分裂法 分?jǐn)?shù)步 法以及激波捕捉法等 與此同時(shí) 數(shù)學(xué)模型研究 離散化方法的理論分析研究和差分 方法的定性分析理論也取得了重大進(jìn)展 從而為計(jì)算流體動(dòng)力學(xué)數(shù)值方法的設(shè)計(jì) 選 擇和應(yīng)用提供了理論依據(jù) 對(duì)計(jì)算流體動(dòng)力學(xué)的發(fā)展起到了積極的作用 如今 計(jì)算流體動(dòng)力學(xué) c f d 己進(jìn)入計(jì)算機(jī)數(shù)值模擬 數(shù)值實(shí)驗(yàn)和計(jì)算機(jī)設(shè)計(jì) 時(shí)代 其大致可分為有限差分法 f i n i t ed i f f e r e n c em e t h o d f d m 有限元法 f i n i t e e l e m e n tm e t h o d f e m 和有限體積法 f i n i t ev o l u m em e t h o d f v m 三個(gè)分支 其 中有限體積法被目前的大部分商用c f d 軟件所采用 1 9 8 0 年 s v p a t a n k e r 在其專(zhuān) 著 n u m e r i c a lh e a tt r a n s f e ra n df l u i df l o w 中對(duì)f v m 作了全面闡述 此后該方法得 到了廣泛應(yīng)用1 35 目前針對(duì)f v m 的分析研究和應(yīng)用拓展也在進(jìn)一步深化 如p c h o w 等學(xué)者提出的擴(kuò)展有限體積法適用于任意多邊形非結(jié)構(gòu)網(wǎng)格 目前 常用的商用c f d 軟件有p h o e n i c s s t a r c d f i d a p c f x f l u e n t 其中功能最強(qiáng)大 適用性最 廣 國(guó)外內(nèi)使用最為廣泛是f l u e n t 1 5 本課題研究?jī)?nèi)容 本課題研究?jī)?nèi)容主要主要有以下幾方面 第1 章緒論 1 介紹離心泵內(nèi)部流場(chǎng)數(shù)值模擬的研究現(xiàn)狀及計(jì)算流體動(dòng)力學(xué)技術(shù)的發(fā)展 從 控制方程 湍流流動(dòng)的數(shù)值模擬方法 湍流模型后處理等方面介紹了湍流流動(dòng)數(shù)值模 擬理論 2 介紹了計(jì)算流體動(dòng)力學(xué)的基本理論知識(shí) 闡述了流體與流動(dòng)的基本特性 給出了流體動(dòng)力學(xué)的基本控制方程 3 運(yùn)用p r o e 造型軟件對(duì)r y 型離心泵過(guò)流部件進(jìn)行實(shí)體造型 建立離心泵的 三維實(shí)體模型 4 介紹了f l u e n t 對(duì)離心泵內(nèi)部流場(chǎng)的數(shù)值計(jì)算過(guò)程 應(yīng)用滑移網(wǎng)格模型對(duì)離 心泵內(nèi)部湍流流動(dòng)進(jìn)行了非定常數(shù)植模擬 得到不同工況下各個(gè)時(shí)刻的速度分布 壓 力分布及流場(chǎng)分布情況 通過(guò)對(duì)數(shù)值計(jì)算結(jié)果的分析比較 揭示了離心泵內(nèi)部的流動(dòng) 規(guī)律 5 根據(jù)數(shù)值模擬結(jié)果對(duì)原型號(hào)泵進(jìn)行優(yōu)化設(shè)計(jì)改進(jìn) 對(duì)優(yōu)化后結(jié)構(gòu)進(jìn)行數(shù)值模 擬 并以此為基礎(chǔ)進(jìn)行離心泵的性能預(yù)測(cè) 將數(shù)值計(jì)算的預(yù)測(cè)結(jié)果與已有的實(shí)驗(yàn)測(cè)試 結(jié)果進(jìn)行分析比較 驗(yàn)證了所采用方法的可行性和正確性 為離心泵的優(yōu)化設(shè)計(jì)研究 提供參考 1 6 本章小結(jié) 本章介紹了本課題的研究背景和研究?jī)?nèi)容的必要性及研究意義 闡述了離心泵內(nèi) 部流動(dòng)數(shù)值模擬的國(guó)內(nèi)外發(fā)展現(xiàn)狀 給出了本課題研究所涉及到的內(nèi)容 河北工程大學(xué)碩士學(xué)位論文 第2 章離心泵內(nèi)部流場(chǎng)數(shù)值模擬理論基礎(chǔ) 隨著科技的進(jìn)步和經(jīng)濟(jì)的發(fā)展 各行各業(yè)對(duì)高性能離心泵的需求越來(lái)越多 但是 傳統(tǒng)設(shè)計(jì)方法并不能滿(mǎn)足這個(gè)需求 這就要求設(shè)計(jì)者采用現(xiàn)代設(shè)計(jì)方法和理論 詳盡 地掌握離心泵性能和內(nèi)部流場(chǎng)的分布狀況 而計(jì)算流體動(dòng)力學(xué) c o m p u t a t i o n a lf l u i d d y n a m i c s c f d 將數(shù)值計(jì)算方法和數(shù)據(jù)可視化技術(shù)有機(jī)結(jié)合起來(lái) 對(duì)流動(dòng) 換熱等 物理現(xiàn)象進(jìn)行模擬分析 是解決流動(dòng)和換熱問(wèn)題的全新技術(shù)手段 隨著計(jì)算機(jī)技術(shù)的 快速發(fā)展 目前計(jì)算流體動(dòng)力學(xué)已經(jīng)廣泛應(yīng)用到現(xiàn)代科學(xué)研究和工程應(yīng)用當(dāng)中 2 1 計(jì)算流體動(dòng)力學(xué)理論基礎(chǔ)瞄一胡 2 1 1 計(jì)算流體動(dòng)力學(xué)基本思想 c f d 的基本思想是 用一系列離散點(diǎn)的變量值的集合來(lái)替代原來(lái)在時(shí)間域和空間 域上連續(xù)的物理量場(chǎng) 通過(guò)一定的原則和方式建立能夠反映這些離散點(diǎn)上場(chǎng)變量關(guān)系 的代數(shù)方程組 最后借助于求解代數(shù)方程組獲得場(chǎng)變量的近似值 c f d 對(duì)流動(dòng)過(guò)程進(jìn) 行的數(shù)值模擬受流動(dòng)基本方程控制 最終可以得到復(fù)雜流場(chǎng)內(nèi)各個(gè)位置上的基本物理 量 壓力 速度 溫度等 的分布 非定常計(jì)算情況下還可以得到這些物理量隨時(shí)間 變化的分布情況 另外 根據(jù)c f d 數(shù)值模擬結(jié)果結(jié)合c a d 軟件可以對(duì)現(xiàn)有結(jié)構(gòu)進(jìn)行 優(yōu)化設(shè)計(jì) 流動(dòng)控制方程一般是非線(xiàn)性的 自變量多 求解域的幾何形狀和邊界條件都比較 復(fù)雜 想得到解析解很困難 而利用c f d 進(jìn)行數(shù)值模擬不但可以得到滿(mǎn)足需要的數(shù)值 解 還可以在計(jì)算機(jī)上形象地再現(xiàn)流動(dòng)情景 另外 可以在c f d 數(shù)值模擬運(yùn)算時(shí)選擇 不同的流動(dòng)參數(shù)進(jìn)行各種數(shù)值實(shí)驗(yàn) 這樣方案對(duì)比不會(huì)受到物理模型和實(shí)驗(yàn)?zāi)Ｐ偷南?制 方便簡(jiǎn)單 靈活性好 得到的資料完整詳細(xì) 與實(shí)驗(yàn)室進(jìn)行實(shí)驗(yàn)相比 省錢(qián)省時(shí) 節(jié)約成本 當(dāng)然c f d 也存在一定的局限性 首先 它得到的計(jì)算結(jié)果是有限離散點(diǎn)上的數(shù)值 解 由于采用的求解方法離散近似 因此結(jié)果必定存在一定的計(jì)算誤差 其次 它需 要借助理論分析或者模型實(shí)驗(yàn)提供某些流動(dòng)參數(shù) 最后還需要對(duì)計(jì)算得出的數(shù)值解進(jìn) 行驗(yàn)證 再者 操作者的經(jīng)驗(yàn)與操作技巧決定了資料的收集 整理與使用 另外 c f d 進(jìn)行數(shù)值計(jì)算對(duì)計(jì)算機(jī)軟硬件配置的要求比較高 而且由于數(shù)值計(jì)算方法的原因還可 能導(dǎo)致討1 算誤差大 造成結(jié)果不真實(shí) 第2 章離心泵內(nèi)部流場(chǎng)數(shù)值模擬理論基礎(chǔ) 2 1 2 流體與流動(dòng)的基本特性 c f d 的主要研究對(duì)象是流體 流體性質(zhì)和流動(dòng)狀態(tài)決定著c f d 中計(jì)算模型的選 擇設(shè)置以及計(jì)算方法的選擇 也決定數(shù)值計(jì)算中各個(gè)物理量的分布情況 1 理想流體與黏性流體 黏性應(yīng)力是指流體在運(yùn)動(dòng)時(shí)對(duì)相鄰兩層流體間相對(duì)運(yùn)動(dòng)的抵抗力 黏性是指流體 抵抗變形的性質(zhì) 即流體具有的這種抵抗兩層流體間相對(duì)滑動(dòng)的能力 流體性質(zhì)決定了流行黏性的大小 并且流體黏性大小與流體溫度有關(guān) 隨著流體 溫度的變化而顯著變化 無(wú)黏流體 也叫做理性流體是指流體無(wú)黏性 即在流體的黏 性較小 運(yùn)動(dòng)的相對(duì)速度也不大情況下 流體的黏性應(yīng)力相比其他類(lèi)型力 比如慣性 力 可以忽略不計(jì) 這時(shí)可以近似認(rèn)為流體無(wú)黏性 相反 黏性較大的流體 則稱(chēng)為 黏性流體 2 牛頓流體與非牛頓流體 流體劃分為牛頓流體和非牛頓流體的依據(jù)是即流體是否滿(mǎn)足牛頓內(nèi)摩擦定律 符合 牛頓內(nèi)摩擦定律 切應(yīng)力與速度梯度成正比的流體稱(chēng)為牛頓流體 相反 流體內(nèi)部切 應(yīng)力與速度梯度不是線(xiàn)性分布 不符合牛頓內(nèi)摩擦定律的流體稱(chēng)為非牛頓流體 3 可壓流體與不可壓流體 流體劃分為可壓流體與不可壓流體的依據(jù)是流體壓縮性的大小 流體密度隨壓強(qiáng) 變化較大且不可視為常數(shù)的流體稱(chēng)為可壓流體 相反 流體密度隨壓強(qiáng)變化較小甚至 可以看作常數(shù)的流體稱(chēng)為不可壓流體 一般情況下流體運(yùn)動(dòng)都可視為不可壓流體進(jìn)行分析計(jì)算 水擊現(xiàn)象除外 低速 氣體在v 5 0 m s 情況下也可視為不可壓流體 這是因?yàn)榈退贇怏w的壓強(qiáng)與溫度在流動(dòng) 過(guò)程中變化很小 它的密度變化也比較小 可以近似看作常數(shù) 但是高速氣流則需要 作為可壓流體來(lái)處理 4 定常流動(dòng)與非定常流動(dòng) 流體流動(dòng)的物理量 如速度 壓力等 是否隨時(shí)間變化是流動(dòng)分為定常流動(dòng)和非 定常流動(dòng)的劃分依據(jù) 定常流動(dòng)情況下流體中各個(gè)物理量不隨時(shí)間變化而改變 因此 定常流動(dòng)也叫做恒定流動(dòng)或穩(wěn)態(tài)流動(dòng) 非定常流動(dòng)情況下流體中的物理量隨時(shí)間變化 而變化 因此 非定常流動(dòng)也叫做非恒定流動(dòng)或瞬態(tài)流動(dòng) 流體機(jī)械在正常運(yùn)轉(zhuǎn)情況 下可以看作定常流動(dòng) 但是在啟停機(jī)時(shí) 流體流動(dòng)一 般是非定常流動(dòng) 5 層流與湍流 流體流動(dòng)時(shí)在自然界中的存在形式主要是層流和湍流 一般情況下 層流只是特 殊情況 而湍流則是普遍存在 河北工程大學(xué)碩士學(xué)位論文 2 1 3 流體動(dòng)力學(xué)的基本控制方程 流體流動(dòng)所遵循的基本守恒定律包括 質(zhì)量守恒定律 動(dòng)量守恒定律 能量守恒 定律 這些守恒定律用控制方程來(lái)描述 1 質(zhì)量守叵方程 任何流動(dòng)問(wèn)題都滿(mǎn)足質(zhì)量守恒定律 基本含義是指單位時(shí)間內(nèi)流體微元體中質(zhì)量 的增加等于同一時(shí)間間隔內(nèi)流入該微元體的凈質(zhì)量 控制方程如下 i o p 譬 挈 挈 0 2 1 8 t 瓠跏娩 定義矢量符號(hào)斫v 口 冬 要 冬 上式 2 1 可改寫(xiě)成 等 d v 0 2 2 在上述式子中 p 密度 r 時(shí)間 u 速度矢量 方程 2 1 是瞬態(tài)三維可壓流體的連續(xù)性方程 如果流體屬于不可壓縮流體 則 密度p 為常數(shù) 式 2 1 變?yōu)?罷 譬 娑 0 0 2 3 i l 優(yōu) o y c 7 z 如果流體密度不隨時(shí)間變化而改變 式 2 1 變?yōu)?掣 挈 挈 0 2 4 a 一 7 2 動(dòng)量守恒方程 動(dòng)量守恒定律實(shí)際上就是牛頓第二定律 基本含義是指外界作用在微元體上的各 種力之和等于該微元體中流體動(dòng)量對(duì)時(shí)間的變化率 因此 可以寫(xiě)出x y 和z 方向 上的動(dòng)量守恒方程 即 警砌 訓(xùn) 一罷 孥 等 誓 c 泣5 a 等脅 訓(xùn) 一考 誓 孥 誓 e 沼5 b 百o p u d i v 小一詈 魯 等 誓 c 2 5 c 上式中 p 誡體中微元體上的壓力 p a r w r 彳 因分子黏性作用產(chǎn)生的作用在微元體表面上的黏性應(yīng) 力f 的分量 f f f 微元體上3 個(gè)方向的體積力 若體積力只受重力 r z 軸豎直向上 則有e 0 e 0 t p g 式 2 5 a 2 5 b 2 5 c 的動(dòng)量守恒方程對(duì)任何類(lèi)型的流體 包括非牛頓流 體 都成立 對(duì)于牛頓流體 有 鏟2 喀 肋v l0 x 礦2 孝 樅1 丁 2 瓦o w 勵(lì)v f 掣可 叫i 萬(wàn) 瓦o v j r 拋 i 鏟鏟 老 蕓 t y z t z y 叫i 瓦o v 萬(wàn)o w j f1 式中 動(dòng)力黏度 五 第二黏度 一般情況下取兄 一2 3 將式 2 6 帶入 2 5 可得 等砌 甜 瑚v 咿咖 一罷城 等砌v 擾 瑚v o t g r a 咖 一考 s 等 咖 疣v 礦日咖 一罷o z s c 2 6 2 7 a 2 7 b 2 7 c 式中 礦酬 掣 掣 業(yè)8 z 符號(hào)趴s 那 代表動(dòng)量方程的廣義源項(xiàng) s f s s f sv s 二 c 5 二 其 p s s 和s 二可用如下公式表示 河北工程大學(xué)碩士學(xué)位論文 j 吳f 熹1 熹f 票 曇f 豢1 曇 砌v 甜 2 8 a 虬2 瓦l 瓦 j 萬(wàn)l 瓦 l 瓦l 瓦j 瓦 砌例 曠甜針甜撲 豺?qū)Ｃ{ s y2 瓦l 萬(wàn)j 萬(wàn)l 萬(wàn)j 瓦 萬(wàn)j 萬(wàn)呦 j j 曇 老 曇 老 丟卜老 曇 勵(lì)v 甜 2 瓦 瓦j 萬(wàn)i 瓦j 瓦p 瓦j 瓦愀例 2 8 b 2 8 c 通常情況下 s s y 和s 為小量 特別對(duì)于黏性為常數(shù)的不可壓縮流體 有 s s s u 7 5 程t 2 4 7 明展升彤瓦如卜 煎型 皇幽 a p u v 墊 o t8 x a y 8 z 丟 罷 茜 孝 魯 老 一塞 甌 型 塹塵 墮型 a c o o t舐 砂8z 曇 喀 專(zhuān) 磅 魯 老 一考 s y 塹 煎型 煎型 趔 o t8 x a y 8 z 曇 蕓 號(hào)卜茅 魯 罷 一考 疋瓦l 瓦 j 萬(wàn)l 萬(wàn)j 瓦p 瓦j 一舌檔z 式 2 7 及 2 9 為動(dòng)量方程 也稱(chēng)作n a v i e r s t o k e s n s 方程 2 9 a 2 9 b 2 9 c 3 能量守恒方程 具