本案例演示Fluent动网格中网格光顺的控制方法。

1 案例1：弹簧光顺

以2D、Double Precision方式启动Fluent
读取网格文件metal-sheet.msh

初始网格如下图所示。

1.1 分割计算区域

变形区域需要进行分割。

点击按钮Setting Up Domain Zones → Separate → Faces...打开如下图所示对话框，按图所示设置参数，点击按钮Separate分割边界

1.2 General设置

激活Transient计算

1.3 编译UDF

点击按钮User Defined → User Defined → Functions → Compiled...打开对话框，编译UDF

源代码如下：

 #include "udf.h"
 #include "dynamesh_tools.h"
 #define AMPL .1
 #define FREQ 10.
  
 DEFINE_GRID_MOTION(vibration1, domain, dt, time, dtime)
 {
   face_t f;
   Thread *tf = DT_THREAD ((Dynamic_Thread *)dt);
   Node *v;
   real x;
   int n;
   /* set deforming flags */
   SET_DEFORMING_THREAD_FLAG (tf->t0);
  
   begin_f_loop (f, tf)
    {
       f_node_loop (f, tf, n)
    {
       v = F_NODE (f, tf, n);
       x=NODE_X(v);
       NODE_Y(v)=+AMPL*sin(M_PI*x/2)*sin(2*M_PI*FREQ*time)+0.03;
    }
       Update_Face_Metrics (f, tf);
    }
   end_f_loop (f,tf)
   printf("n x_velocity n");
 }
  
 DEFINE_GRID_MOTION(vibration2, domain, dt, time, dtime)
 {
   face_t f;
   Thread *tf = DT_THREAD ((Dynamic_Thread *)dt);
   Node *v;
   real x;
   int n;
   /* set deforming flags */
   SET_DEFORMING_THREAD_FLAG (tf->t0);
  
   begin_f_loop (f, tf)
    {
       f_node_loop (f, tf, n)
    {
       v = F_NODE (f, tf, n);
       x=NODE_X(v);
       NODE_Y(v)=+AMPL*sin(M_PI*x/2)*sin(2*M_PI*FREQ*time)-0.03;
    }
       Update_Face_Metrics (f, tf);
    }
   end_f_loop (f,tf)
  
 }