文档介绍：第 10 章应用举例 边界条件 源项 物理属性 反应速率( Reacting Rates ) 用户定义标量( User_Defined Scalars ) 边界条件这部分包含了边界条件 UDFs 的两个应用。两个在 FLUENT 中都是作为解释式 UDFs 被执行的。 涡轮叶片的抛物线速度入口分布要考虑的涡轮叶片显示在 Figure 中。非结构化网格用于模拟叶片周围的流场。区域从底部周期性边界延伸到顶部周期性边界,左边是速度入口,右边是压力出口。 Figure : The Grid for the Turbine Vane Example 常数 x 速度应用于入口的流场与抛物线 x 速度应用于入口的流场作了比较。当采用分段线性分布的型线的应用是有效的对边界型线选择,多项式的详细说明只能通过用户定义函数来完成。常数速度应用于流场入口的结果显示在 Figure 和 Figure 中。当流动移动到涡轮叶片周围时初始常速度场被扭曲。 Figure : Velocity Magnitude Contours for a Constant Inlet x Velocity Figure : Velocity Vectors for a Constant Inlet x Velocity 现在入口 x 速度将用以下型线描述: 这里变量 y 在人口中心是 ,在顶部和底部其值分别延伸到 0745 .0?。这样 x 速度在入口中心为 20m/sec ,在边缘为 0。 UDF 用于传入入口上的这个抛物线分布。 C 源代码( )显示如下。函数使用了 Section 中描述的 Fluent 提供的求解器函数。/***********************************************************************/ /* */ /* UDF for specifying steady-state velocity profile boundary condition */ /***********************************************************************/ #include "" DEFINE_PROFILE(inlet_x_velocity, thread, position) { real x[ND_ND]; /* this will hold the position vector */ real y; face_t f; begin_f_loop(f, thread) { F_CENTROID(x,f,thread); y= x[1]; F_PROFILE(f, thread, position) =20. - y*y/(.0745*.0745)*20.; } end_f_loop(f, thread) } 函数,被命名为 inlet_x_velocity, 使用了 DEFINE_PROFILE 定义并且有两个自变量: thread 和 position 。 Thread 是一个指向面的 thread 的指针, position 是一个整数,它是每个循环( loop )内为变量设置的数值标签。函数通过声名变量 f 作为 face_t 的数据类型。一维数组 x 和变量 y 被定义为 real 数据类型。循环宏用于循环区域内每个面来创建型线,或数据数组。每个循环内, F_CENTROIDS 为带指标 f 的面输出面质心(数组 x)的值,指标 f在被 thread 指向的线( thread )上。存储在 x[1] 中的 y坐标用于为变量 y赋值,然后用于计算 x速度。这个值接着被分配到 F_PROFILE ,使用整数 position (在你从 Velocity Inlet 面板中选择的 UD F 作为 x速度的边界条件的基础上通过求解器传递给它)来设置内存中 x速度面值。为了在 FLUENT 中使用这个解释式 UDF ,你必须首先编译它。 Define User-Defined Functions Interpreted... 在 Interpreted UDFs 面板中,在 Source File Name 区域命名你的函数。如果必要,在 mand Name 区域输入你的 C 预处理程序类型和在 Stack Size 下输入堆栈大小。打开 Display Assembly Listing 在你的控制台窗口中能看到汇编列表,当函数编