搜索FLUENT官方教程之嵌入式滑移網(wǎng)格實(shí)例
2016-08-11 by:CAE仿真在線 來源:互聯(lián)網(wǎng)
[本例內(nèi)容來源于fluent官方教程]
現(xiàn)實(shí)生活中經(jīng)常存在這樣的模型:物體在運(yùn)動的同時,還包含著旋轉(zhuǎn)。換句話說,即其在公轉(zhuǎn)的同時還存在自轉(zhuǎn)現(xiàn)象。例如運(yùn)動的割草機(jī),其整體在運(yùn)動,同時其割草部件在告訴旋轉(zhuǎn)。同樣的例子還有很多,比如吸塵器、旋轉(zhuǎn)飛行的足球等等。ANSYS FLUENT13.0之后的版本提供了嵌入式滑移網(wǎng)格模型,允許用戶對這一現(xiàn)象進(jìn)行模擬,從而可以替代動網(wǎng)格,因此在某種程度上保證了網(wǎng)格質(zhì)量,也不會喪失太多的精度。
圖1 幾何描述
1、幾何描述
本例幾何描述如圖1所示。本計(jì)算域?yàn)橐粋€封閉系統(tǒng),其中包含一個十字旋轉(zhuǎn)區(qū)域(直徑0.1m),其圓心偏離計(jì)算域中心區(qū)域在x及y方向上均為0.1m。旋轉(zhuǎn)區(qū)域fluid-rotor順時針旋轉(zhuǎn)速度為2rad/s,旋轉(zhuǎn)區(qū)域fluid-circle圓心位于計(jì)算域中心,其沿逆時針旋轉(zhuǎn)速度為1rad/s。工作流體為液態(tài)水,其密度1000kg/m3,粘度0.001kg/m-s。
2、前期UDF準(zhǔn)備
本例需要UDF計(jì)算fluid-rotor區(qū)域圓心坐標(biāo)。我們這里利用DEFINE_ZONE_MOTION宏。
程序代碼如下:
#include "udf.h"
#define PI 3.141592654
DEFINE_ZONE_MOTION(rotor, omega, axis, origin, velocity, time, dtime)
{
real theta0, thetap1, omegac, omegar, radr, ox1, oy1;
omegar = -2.0;
omegac = 1.0;
theta0 = PI/4.;
radr = 0.141435;
thetap1 = omegac*(time+dtime);
*omega = omegar;
ox1 = radr*cos(theta0+thetap1);
oy1 = radr*sin(theta0+thetap1);
origin[0] = ox1;
origin[1] = oy1;
origin[2] = 0.0;
}
程序解釋:
(1)omegar變量定義為fluid_rotor的旋轉(zhuǎn)角速度,omegac定義為fluid_circle區(qū)域的旋轉(zhuǎn)角速度。由右手定則確定旋轉(zhuǎn)方向。
(2)theta0為fluid_rotor的初始圓心位置與計(jì)算域圓心間的角度。Radr為fluid_rotor中心到計(jì)算域中心的距離。由于本例中fluid_rotor為旋轉(zhuǎn),因此該半徑為定值。
(3)thetap1為在當(dāng)前時間步旋轉(zhuǎn)的角度。*omega=omegar,此條語句的作用為將速度矢量賦值給omega以傳入求解器。
(4)ox1及oy1為計(jì)算的當(dāng)前fluid_rotor圓心位置。并將其值賦值給origin數(shù)組。
3、fluent模型設(shè)置
導(dǎo)入msh文件,選擇瞬態(tài)計(jì)算,利用starnd k-e湍流模型,采用標(biāo)準(zhǔn)壁面函數(shù)。
選擇材料為water-liquid,修改密度1000kg/m3,粘度0.001kg/m-s。
設(shè)置三區(qū)域介質(zhì)為water-liquid。
4、加載UDF
本例udf可以采用解釋型,也可以采用編譯型。
【define】> 【user-defined…】>【functions】>【interpreted udfs】,彈出設(shè)置對話框,如圖2所示。選擇創(chuàng)建的udf源文件。點(diǎn)擊interpret。
圖2 解釋udf
4、區(qū)域設(shè)置
Fluid-outer設(shè)置介質(zhì)為water-liquid,其他保持默認(rèn)。該區(qū)域?yàn)殪o止。
Fluid-circle設(shè)置如圖3所示。
圖3 circle區(qū)域設(shè)置
勾選mesh motion選項(xiàng),注意relative to cell zone為absolute,設(shè)置旋轉(zhuǎn)速度speed為1rad/s。注意此處選擇absolute,表示fluid-circle的絕對旋轉(zhuǎn)速度為1rad/s。
Fluid-rotor區(qū)域設(shè)置如圖4所示。
圖4 rotor區(qū)域設(shè)置
與circle區(qū)域設(shè)置類似,不過選擇relative to cell zone為fluid-circle,表示該區(qū)域是相對于fluid-circle區(qū)域運(yùn)動的。同時設(shè)置zone motion function為rotor。
5、邊界條件設(shè)置
本例為封閉系統(tǒng),為簡便起見,壁面采用默認(rèn)設(shè)置,即無滑移壁面邊界。
本例需要建立兩組interface,即fluid-rotor與fluid-circle之間,fluid-circle與fluid-outer之間。如圖5所示。
圖5 建立interface對
6、其他設(shè)置
Solution methods:壓力耦合采用coupled方法。其他采用默認(rèn)選項(xiàng)。
Solution controls:采用默認(rèn)設(shè)置。
初始化:采用默認(rèn)設(shè)置初始化即可。
動畫設(shè)置:可以錄制速度及壓力等物理量云圖動畫。
時間步長:Δt = (0.0349066 rad) /(2 rad/sec) = 0.0174533 sec,我們假定每時間步旋轉(zhuǎn)2°,即0.0249066rad,以最大速度2rad/s進(jìn)行計(jì)算。時間步數(shù)取1800,即旋轉(zhuǎn)10圈。如圖6所示。
圖6 計(jì)算設(shè)置
7、計(jì)算結(jié)果
略
8、總結(jié)
(1)這類公轉(zhuǎn)+自轉(zhuǎn)的問題,既可以使用動網(wǎng)格實(shí)現(xiàn),也可以采用fluent中的嵌入式滑移網(wǎng)格實(shí)現(xiàn)。
(2)此例中兩區(qū)域均為旋轉(zhuǎn)運(yùn)動,對于滑移網(wǎng)格特點(diǎn),平移運(yùn)動也可以實(shí)現(xiàn)。
(3)采用嵌入式滑移網(wǎng)格,避免了網(wǎng)格重構(gòu)的網(wǎng)格質(zhì)量損失。
相關(guān)標(biāo)簽搜索：FLUENT官方教程之嵌入式滑移網(wǎng)格實(shí)例 fluent滑移網(wǎng)格 fluent官方教程 Fluent培訓(xùn) Fluent流體培訓(xùn) Fluent軟件培訓(xùn) fluent技術(shù)教程 fluent在線視頻教程 fluent資料下載 fluent分析理論 fluent化學(xué)反應(yīng) fluent軟件下載 UDF編程代做
