FLUENT嵌入式滑移網(wǎng)格實(shí)例

[本例內(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所示。本計算域?yàn)橐粋€封閉系統(tǒng),其中包含一個十字旋轉(zhuǎn)區(qū)域(直徑0.1m),其圓心偏離計算域中心區(qū)域在x及y方向上均為0.1m。旋轉(zhuǎn)區(qū)域fluid-rotor順時針旋轉(zhuǎn)速度為2rad/s,旋轉(zhuǎn)區(qū)域fluid-circle圓心位于計算域中心,其沿逆時針旋轉(zhuǎn)速度為1rad/s。工作流體為液態(tài)水,其密度1000kg/m3,粘度0.001kg/m-s。

2、前期UDF準(zhǔn)備

本例需要UDF計算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的初始圓心位置與計算域圓心間的角度。Radr為fluid_rotor中心到計算域中心的距離。由于本例中fluid_rotor為旋轉(zhuǎn),因此該半徑為定值。

(3)thetap1為在當(dāng)前時間步旋轉(zhuǎn)的角度。*omega=omegar,此條語句的作用為將速度矢量賦值給omega以傳入求解器。

(4)ox1及oy1為計算的當(dāng)前fluid_rotor圓心位置。并將其值賦值給origin數(shù)組。

3、fluent模型設(shè)置

導(dǎo)入msh文件,選擇瞬態(tài)計算,利用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選項,注意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)選項。

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)行計算。時間步數(shù)取1800,即旋轉(zhuǎn)10圈。如圖6所示。

圖6 計算設(shè)置

7、計算結(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ì)量損失。