讓網(wǎng)格動起來 | (十三)FLUENT區(qū)域運動【轉發(fā)】

2017-06-01  by:CAE仿真在線  來源:互聯(lián)網(wǎng)



本次使用的是滑移網(wǎng)格,動網(wǎng)格實現(xiàn)以后再討論。要實現(xiàn)的運動如下圖所示。

讓網(wǎng)格動起來 | (十三)FLUENT區(qū)域運動【轉發(fā)】fluent結果圖片1

杯子中裝滿水,現(xiàn)在以速度1rad/s延續(xù)1s鐘使杯子傾斜1rad,觀察5s鐘內(nèi)水的變化情況。

本例可以用滑移網(wǎng)格或動網(wǎng)格實現(xiàn),但是使用滑移網(wǎng)格能夠保持較好的網(wǎng)格質量。本例使用滑移網(wǎng)格。

FLUENT中的滑移網(wǎng)格可以最大限度的代替動網(wǎng)格,尤其對于一些分界面確定的計算模型。滑移網(wǎng)格可以用于瞬態(tài)模擬中(其它如MRF,SRF,MP則很少用于瞬態(tài)模擬中)。

本例涉及到的內(nèi)容包括:

(1)分界面幾何模型的建立。涉及到多幾何體的創(chuàng)建。在workbench的DM模塊中很容易解決此類問題。但是如果要在ICEM CFD中創(chuàng)建網(wǎng)格的話,則需要進行一些特別的處理。主要是各部分模型網(wǎng)格的組裝問題。

(2)區(qū)域運動的指定。在本例中主要是指定運動區(qū)域的旋轉速度。需要注意的是旋轉中心與旋轉方向的設定。

(3)多相流的使用。本例中使用的是VOF模型。

1、幾何模型

本例的幾何模型如下圖所示。尺寸是隨便定的。

讓網(wǎng)格動起來 | (十三)FLUENT區(qū)域運動【轉發(fā)】fluent結果圖片2

如圖所示,計算模型分為三個區(qū)域,1、2、3分別對應名稱為braket、zone與cup,如前所述,區(qū)域zone與cup為運動區(qū)域,運動方式用UDF進行定義。

2、邊界類型

一共有兩對interface,分別位于zone區(qū)域與cup區(qū)域,zone區(qū)域與braket區(qū)域。如下圖所示。braket區(qū)域的左側、右側及下方邊界均為wall類型,上方邊界為pressure_outlet類型,內(nèi)部邊界為interface類型。cup區(qū)域左、右、下側邊界為wall類型,上方邊界為interface類型。在mesh interface中設置interface對,一共是兩對。

讓網(wǎng)格動起來 | (十三)FLUENT區(qū)域運動【轉發(fā)】fluent結果圖片3

最終組合后的網(wǎng)格如下圖所示。

讓網(wǎng)格動起來 | (十三)FLUENT區(qū)域運動【轉發(fā)】fluent仿真分析圖片4

3、求解參數(shù)設置

導入模型后,在scale面板中進行必要的scale操作以使模型尺度滿足要求。在general面板中設置使用瞬態(tài)模擬,并設置重力加速度。選擇RNG K-Epsilon湍流模型。設置多相流模型為VOF模型,如右下圖所示。

讓網(wǎng)格動起來 | (十三)FLUENT區(qū)域運動【轉發(fā)】fluent仿真分析圖片5 讓網(wǎng)格動起來 | (十三)FLUENT區(qū)域運動【轉發(fā)】fluent仿真分析圖片6

4、添加材料

在材料數(shù)據(jù)庫中添加水至當前材料中。設置空氣為主相,水為第二相。VOF中主相與第二相的設置并沒太多的限制,但是一般將需進行patch的相設置為第二相,或者將含量較少的相設置為第二相。當然沒有明確的規(guī)定。設置水-空氣的表面張力系數(shù)為0.07。

5、區(qū)域運動UDF

這里使用宏DEFINE_TRANSIENT_PROFILE宏,設置當時間小于1s時,旋轉速度1rad/s,其它時間保持靜止。編譯此UDF文件。代碼如下:

#include "udf.h"
DEFINE_TRANSIENT_PROFILE(rotation_rate_ramp,time)
{
real rotation_rate = 0;
if(time <= 1)
{
rotation_rate = 1.0;
}else
{
rotation_rate = 0.0;
}
return rotation_rate;
}

6、區(qū)域設置

主要設置區(qū)域cup與zone的運動。FLUENT中區(qū)域默認為靜止。以區(qū)域cup為例,設置運動如下圖所示:

讓網(wǎng)格動起來 | (十三)FLUENT區(qū)域運動【轉發(fā)】fluent仿真分析圖片7

zone區(qū)域設置與此完全相同。

7、邊界條件

需要特別說明的是,屬于某區(qū)域的wall需要設置相對速度,與區(qū)域的相對速度為0。主要是針對區(qū)域cup與zone。

設置braket區(qū)域中的pressure_outlet邊界靜壓為0,并設定湍流條件。

8、求解器控制參數(shù)

采用默認。

9、初始化

選擇all zone進行初始化。patch區(qū)域cup的第二相(水)體積分數(shù)為1。如下圖所示。

讓網(wǎng)格動起來 | (十三)FLUENT區(qū)域運動【轉發(fā)】fluent流體分析圖片8

10、設置迭代進行計算

設置時間步長0.1s,時間步數(shù)50步進行迭代。如下圖所示。

讓網(wǎng)格動起來 | (十三)FLUENT區(qū)域運動【轉發(fā)】fluent流體分析圖片9

11、計算結果

如本篇開始所示。

——流沙


開放分享:優(yōu)質有限元技術文章,助你自學成才

相關標簽搜索:讓網(wǎng)格動起來 | (十三)FLUENT區(qū)域運動【轉發(fā)】 Fluent培訓 Fluent流體培訓 Fluent軟件培訓 fluent技術教程 fluent在線視頻教程 fluent資料下載 fluent分析理論 fluent化學反應 fluent軟件下載 UDF編程代做 Fluent、CFX流體分析 HFSS電磁分析 

編輯
在線報名:
  • 客服在線請直接聯(lián)系我們的客服,您也可以通過下面的方式進行在線報名,我們會及時給您回復電話,謝謝!
驗證碼

全國服務熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓QQ咨詢:點擊咨詢 點擊咨詢
項目QQ咨詢:點擊咨詢
email:kf@1cae.com