二維離心泵流動模擬上機指導書(fluent)(共15頁)

1、 二維離心泵的數(shù)值(shz)模擬(mn)與性能預測一、實驗(shyn)目的熟悉和掌握CFD數(shù)值模擬的基本方法，能夠獨立進行簡單二維水力模型的CFD數(shù)值模擬。二、研究對象研究如圖所示的二維離心泵，該泵由旋轉(zhuǎn)的葉輪和靜止的蝸殼兩部分構(gòu)成。流體從葉輪中央的圓形進口沿徑向均勻進入葉輪，經(jīng)過旋轉(zhuǎn)的葉片作用后，得到能量，從蝸殼出口排出。已知葉輪的葉片數(shù)為6，葉輪進、出口直徑分別為120mm和220mm，葉片進口安放角（葉片與圓周方向夾角）和出口安放角分別為和，葉片厚度為3mm。蝸殼隔舌角 （葉輪中心至蝸殼螺旋線起點的連線與水平夾角）為 ，出口段擴散角為。 圖表 SEQ 圖表 * ARABIC 1 二維離心

2、泵示意圖 圖表 SEQ 圖表 * ARABIC 2 UG NX所繪二維離心泵三、計算步驟1、 利用 Gambit 對計算區(qū)域離散化和指定邊界條件類型 步驟 1：導入幾何模型 生成幾何模型的方式有許多種，如Autocad，Pro/E，UG NX等，Gambit也自帶簡單的繪圖功能，在這里UG NX繪圖。如下圖，我們將會給出繪圖文件。在UG NX繪完圖后，需將結(jié)果(ji gu)導出以便Gambit使用，這里導出為Parasolid，生成后綴名為x_t的文件。主要(zhyo)內(nèi)容：在Gambit中選擇File/Import/Parasolid命令，選取(xunq)先前生成的文件11.x_t，則二維離

3、心泵模型被裝入到Gambit。結(jié)果如圖：圖表 SEQ 圖表 * ARABIC 3 導入到Gambit的二維離心泵幾何模型步驟 2：網(wǎng)格劃分 為了對幾何區(qū)域劃分網(wǎng)格，單擊Operation / Mesh / Face / Mesh Face按鈕，彈出如圖所示的Mesh Faces對話框。在Faces列表框中選取蝸殼區(qū)域，在Elements列表框中選擇Tri(三角形單元)，在Type列表框中選拌Pave(非結(jié)構(gòu)網(wǎng)格)，選中Scheme命令組中的Apply復選框，然后，從Spacing區(qū)域的列表框中選擇Interval Size(指定網(wǎng)格間隔)，在文本框中輸入10，選中該區(qū)域中的Apply復選框，最

4、后選中Options區(qū)域中的Mesh復選框，單擊Apply按鈕，則生成蝸殼內(nèi)流體區(qū)域的網(wǎng)格。圖表(tbio) SEQ 圖表(tbio) * ARABIC 4 Mesh Faces對話框重復上述命令，從Faces列表框中選取葉輪區(qū)域，在Spacing文木框(m kun)中輸入“8，則可生成轉(zhuǎn)輪內(nèi)流體區(qū)域網(wǎng)格。離心泵的網(wǎng)格結(jié)果如圖所示。圖表 SEQ 圖表 * ARABIC 5 二維離心泵網(wǎng)格步驟(bzhu) 3：邊界條件類型(lixng)的指定 （1）指定邊界(binji)類型單擊 Operation / Zones / Specify Boundary Types按鈕，彈出 Specify Bo

5、undary Types對話框在對話框中，進行如下操作：將Type欄選為Velocity_inlet，在Entity欄選取Edges， 并選中葉輪進口圓作進口邊界。將Type欄選為Outflow，在 Entity欄選取Edges，并選中渦殼出口邊作為出口邊界。將Type欄選為Wall，在Entity欄選取Edges，并選中組成6個葉片的所有邊作為壁面邊界1.將Type欄選為Wall，在Entity欄選取Edges，并選中葉輪出口作為壁面邊界2(該壁面將在FLUENT的邊界條件的指定中改為 interior 類型）。將Type欄選為Wall，在Entity欄選取Edges，并選中除出口邊外的組成

6、渦殼的所有邊，作為壁面邊界3。 圖表(tbio) SEQ 圖表(tbio) * ARABIC 6 Specify Boundary Types對話框 圖表(tbio) SEQ 圖表 * ARABIC 7 Specify Continuum Types對話框（2）設(shè)定區(qū)域類型單擊 Operation/Zones/Specify Continuum Types 按鈕。彈出 Specify Continuum Types對話框，該對話框用于指定區(qū)域類。在Type欄選擇FLUID，在Entity欄選擇 Faces，并選中代表葉輪流體的面作為流動區(qū)域1。在確認對話框域上方的Add 復選框選中的前提下，單

7、擊Apply按鈕。在 Type欄選擇FLUID，在Entity欄選擇 Faces，并選中代表渦殼流體的面作為流動區(qū)域2。最后，選擇File / Export / Mesh命令，給定文件名（如11.msh），并選中Export 2D Mesh復選框，可將上述網(wǎng)格存盤。2. 利用(lyng)Fluent 求解(qi ji)器進行求解 步驟(bzhu) 1：網(wǎng)格文件的讀入、檢查及顯示等 選擇FileReadCase將網(wǎng)格導入，選擇DisplayGrid，顯示網(wǎng)格，如圖：圖表 SEQ 圖表 * ARABIC 8 Fluent中的網(wǎng)格顯示選擇FileCheck對網(wǎng)格文件進行檢查，這里要注意最小的網(wǎng)格體積

8、（minimum volume）值一定要大于0。 網(wǎng)格按比例縮放。在GAMBIT中，生成網(wǎng)格使用的單位是mm，在FLUENT中默認單位是m，需要縮放。為此，選擇Grid / Scale命令，在彈出的Scale Grid對話框中，在Grid Was Created In 下拉菜單中，選取mm，然后單擊Scale按鈕。選擇角速度單位(dnwi)。這主要是為了后續(xù)定義轉(zhuǎn)動參考系作準備。為此，選擇Define/units命令，在打開的Set Units對話框的Quantities下拉列表中選擇(xunz)angular-velocity，在Units下拉列表(li bio)中選擇rpm。步驟 2：選擇

9、計算模型 按下述過程選擇求解器的格式與計算模型，并設(shè)置運行環(huán)境: (1)選擇求解器。選擇Define/Models/Solver命令，彈出圖所示的Solver對話框。在Salver選項組中選擇Segregated（分離式求解器），在Space選項組中選擇2D（二維問題），在Time選項組中選擇Steady（穩(wěn)態(tài)流動），其他用默認值。 (2)設(shè)置運行(ynxng)環(huán)境。選擇Define/Operating Conditions命令(mng lng)，彈出圖7.G所示的Operating Conditions對話框。保持FLUENT默認的參考(cnko)壓力值（一個標準大氣壓）和零參考壓力的位置，

10、選中Gravity復選框，表示計及秉力，并在Y一欄輸入“-9.81 。 (4)設(shè)置(shzh)湍流模型。選擇Define/Models/Viscous命令，彈出Viscous Model對話框，如圖所示。保留圖所示的默認設(shè)置，即表示選擇(xunz)標準k-模型(mxng)，模型的系數(shù)均用默認值。步驟 3：定義物理性質(zhì)本例中流動介質(zhì)為水，F(xiàn)LUENT材料數(shù)據(jù)庫中已包含water liquid這一介質(zhì)，因此，直接復制即可。為此，選擇Define/Materials命令，彈出Materials對話框，單擊Database按鈕，打開 FLUENT材料數(shù)據(jù)庫，在新對話框的Fluid Materials下

11、拉列表中選擇wafer-1iquid，單擊Copy按鈕。單擊Close按鈕關(guān)閉Materials對話樞。步驟 4：設(shè)置邊界條件 拉下述過程設(shè)置邊界條件：（1）將兩個流體區(qū)域的交界面從wall改為interior類型。 FLUENT將彈出一個提示框，點擊Yes按鈕，F(xiàn)LUENT會自動將wall面和其shadow面合并為interior類型。（2）內(nèi)部區(qū)域定義轉(zhuǎn)動參考系。為此，選擇Define/ Boundary Conditions命令，打開Boundary Conditions對話框，然后按下述過程操作：在Material Name下拉列表中，選擇Water-liquid，保持Rotation

12、-Axis origin的默認值(0,0),該點為旋轉(zhuǎn)區(qū)域圓形邊界的曲率中心。在Motion Type下拉列表中，選擇Moving Reference Frame，下拉滾動條，將Rotational Velocity設(shè)為-1470rprn。單擊OK按鈕關(guān)閉對話框。 (3)按類似(li s)辦法，為另一個(y )Fluid區(qū)域（外部(wib)的）選擇流動介質(zhì)為Water-liquid，其他參數(shù)使用默認值。注意該區(qū)域不是轉(zhuǎn)動參考系。 (4)設(shè)置進口邊界。將進日Inlet設(shè)為velocity-inlet(速度進口邊界條件)，在Velocity Inlet 對話框的Velocity Magnitude

13、文本框中輸入“2.2” m/s（取自已知條件)。 (5)設(shè)置出口邊界。將出口Outlet設(shè)為Outflow(出流邊界條件)，參數(shù)使用默認值。 (6)定義代表葉片的壁血區(qū)域相對于轉(zhuǎn)動流體區(qū)域的轉(zhuǎn)動速度（注意內(nèi)部流體區(qū)域被指定為轉(zhuǎn)動參考系）。為此，按下述過程進行：在Boundary Conditions對話框中單擊wall_1然后單擊Set按鈕，啟動Wall對話框。在Momenturn選項中，選中Moving Wall。展開后的對話框?qū)@示運動參數(shù)。在Motion選項組中，選擇Relative to Adjacent Cell done和Rotational o將速度(sd)(相對值)Speed設(shè)

14、為0 rad/s，即葉片壁而相對(xingdu)于fluid轉(zhuǎn)動區(qū)域的速度為完成后，轉(zhuǎn)動軸的原點(Rotation-Axis origin)將位于(wiy)(0，0)，這樣，葉片將隨周圍流體起同速轉(zhuǎn)動。步驟 5：求解設(shè)置 過程如下:選擇求解控制參數(shù)。選擇Solve / Controls/Solution命令，彈出Solution Controls對話框。在Discretization選項組中，對速度與速度與壓力禍合方式，選擇SIMPLEC，對Momentum(動量). Turbulence Kinetic Energy(湍動能)和Turbulence Dissipation Rate(耗散率)

15、，均選擇Second order Upwind（二階迎風各式），其余用默認值。（2）啟動繪制(huzh)殘差功能。選擇Solve/Monitors/Residual命令(mng lng)，彈出residual Monitors對話框。在options選項組中，選中Plot復選框，其余(qy)用默認值。（3）用進口的流動初始條件初始化整個流場的解。選擇Solve/Initialize/Initialize命令，彈出Solution Initialization對話框。在Compute From下拉中選擇inlet，在Reference Frame組合框中，選擇Absolute。(4)保存設(shè)置文件。選擇File/Write/base命令，將當前設(shè)置保存到案例文件。(5)啟動迭代計算。選擇Solve/ Iterate命令，彈出Iterate對話框，設(shè)置迭代次數(shù)為400，其余用默認值。單擊Iterate，迭代開始。步驟 6：保存結(jié)果 選擇FileWriteCase & Data，保存所有的設(shè)置和所有的數(shù)據(jù)。七、結(jié)果處理 (1)顯示壓力等值線。通過Display: Contours命令，打開Contours對話框如下：顯示(xinsh)速度矢量圖。通過(tnggu)Display:Vectors 命令(mng lng)，顯示湍流