Fluent动网格｜光顺

详细信息

作者优秀: 优秀教师/意见领袖/博士学历/特邀专家
平台推荐: 内容稀缺

8月前浏览10552

本文摘要（由AI生成）：

本文介绍了Fluent动网格中网格光顺的控制方法，通过两个案例演示了弹簧光顺和扩散光顺的实现过程。弹簧光顺采用Smoothing方法，使用UDF定义网格运动，通过设置变形区域、动网格设置、动网格预览等步骤实现。扩散光顺采用Diffusion方法，同样使用UDF定义网格运动，通过设置General参数、读取Profile文件、设置动网格、指定区域运动类型、网格预览等步骤实现。两个案例均采用2D、Double Precision方式启动Fluent，并读取相应的网格文件。

本案例演示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)
  
 }