1.7 UDF使用流程_ANSYS Fluent 二次开发指南-QQ阅读女生现言网

上QQ阅读APP看本书，新人免费读10天

设备和账号都新为新人

1.7　UDF使用流程

这里以一个简单的初始化案例来描述UDF的源代码编写、编译及加载过程。通过此案例可以熟悉UDF的整个使用流程。

1.7.1　Fluent中的Patch

Fluent中提供了全局初始化以及局部Patch功能。对于整体区域的全局初始化可以采用starndard及hybrid方法进行初始化，指定各种物理量的初始分布。而对于计算域中的局部区域初始化，则可以通过Patch功能来实现。

在使用Patch方法时，首先需要对要进行Patch的区域进行标记。鼠标右键选择树性节点Solution > Cell Registers，点击弹出菜单项New→Region...可弹出区域定义对话框（图1-4）。

图1-4　创建区域

可以在弹出的对话框中设置几何条件来标记（Mark）区域，如图1-5所示。

图1-5　区域定义对话框

此对话框中可定义的形状类型在2D模型中只有三种：Quad、Circle或Cylinder。在3D模型中对应的是Hex、Sphere以及Cylinder。对于更复杂的形状则无能为力，此时可以借助UDF来解决问题。

1.7.2　案例描述

如图1-6所示的矩形区域为计算区域，其初始温度为300K，计算模型尺寸如图所示。

图1-6　几何模型

图中椭圆部分为要进行初始化处理的区域，其初始温度为500K。生成计算网格如图1-7所示。

注　意

在创建几何模型时，确保几何左下角为坐标原点，否则需要更改UDF程序中的椭圆中心坐标。

图1-7　计算网格

1.7.3　编写UDF源文件

对于这种椭圆形区域的初始化，只能采用UDF来实现。利用DEFINE_INIT宏来实现这种区域的标记工作。

本案例中椭圆方程为：

因此可在文本编辑器中编写UDF代码如下：

#include ”udf.h”
DEFINE_INIT(domainInit,d)
{
    cell_t c;
    Thread *t;
    real xc[ND_ND];
   &ensp;real x;
&ensp;   real y;
    thread_loop_c(t,d)
    {
        begin_c_loop_all(c,t)
        {
            C_CENTROID(xc,c,t);
            x = xc[0];
            y = xc[1];
            
            if(pow((x-0.05),2)/(0.03*0.03)+
                pow(y-0.03,2)/pow(0.015,2)<1)
            {
                C_T(c,t) = 500;
     &ensp;      }
            else
            {
                C_T(c,t)=300;
            }
        }
        end_c_loop_all(c,t)
    }
}