Fluent動網格【6】：部件變形案例

本案例描述使用動網格過程中處理邊界變形的問題。

案例描述

本案例幾何為一個拋物線旋轉成型的幾何體。如圖所示。

其中上壁面剛體運動引起拋物面變形。剛體運動方程為：

\[v=\left\{ \begin{array}{c} \begin{matrix} -0.3t& ,t<=1\\ \end{matrix}\\ \begin{matrix} 0.3t-0.6& ,1<t<=3\\ \end{matrix}\\ \begin{matrix} -0.3t+1.2& ,t<=4\\ \end{matrix}\\ \begin{matrix} 0& ,t>4\\ \end{matrix}\\ \end{array} \right. \]

如圖所示。

拋物面方程為：

\[x = -0.4z^2+0.5z+0.3 \]

UDF

UDF可寫成：

#include "udf.h"
#include#include "udf.h"
#include "dynamesh_tools.h"
 
/*40層,每層高0.02m，時間步長可以為0.05s*/
DEFINE_CG_MOTION(velocity,dt,vel,omega,time,dtime)
{
    if(time <= 1)
    {
        vel[2] = -0.3*time;
    }
    else if(time > 1 && time <= 3 )
    {
        vel[2] = 0.3*time-0.6;
    }
    else if(time <= 4)
    {
        vel[2] = -0.3* time+1.2;
    }
    else
    {
        vel[2]= 0;
    }
 
}
 
 
DEFINE_GEOM(parabola,domain,dt,position)
{
    real radius;
     real norm;
     real x,y,z;
 
     x = position[0];
     y = position[1];
     z = position[2];
 
     radius = -0.4 * z * z + 0.5 * z + 0.3;
     norm = sqrt(x * x + y * y );
 
     position[0] = position[0] * radius/norm;
     position[1] = position[1] * radius/norm;
}

模型創建

模型創建在ICEM CFD中完成，通過導入點坐標的方式創建拋物線，之后繞軸旋轉形成拋物面。

划分六面體網格。

動網格設置

采用layering方式進行網格處理。設置拋物面變形及頂部面剛性運動。

視頻操作

---

更多CFD資料可微信掃描下方二維碼關注微信公眾號。