Fluent udf 第10章 应用举例

*/
/***********************************************************************/
#include "udf.h"
DEFINE_PROFILE(inlet_x_velocity, thread, position)
{
real x[ND_ND];
增加 Surface Monitors 的指标到 1。这将使你定义 monitor-1（如果需要，你可重 命名，在 Name 下的文本输入框内）的参数。选择 Plot，所以所选的数量将作为 计算的收益被画出来。选择 Print 在控制台窗口来看所选数的变化值。选择 Write， 因而信息将被写入文件，它将被给一个名字为 monitor-1.out。（如果你改变了 monitor 的名字，这个名字将用作输出文件的前缀）。 在 Every 下面，你可选择 Iteration, Time Step, or Flow Time。为了监视每个时间 步的结果，你应选择 Time Step 选项。通过点击 Define...你可在 Define Surface Monitor 面板中指定被监视的数量。
这个图精细地显示了速度以期望的振幅值 5m/s 围绕平衡值 20m/s 上下振荡。 10.2 源项 这部分包含了源项 UDF 的应用。在 FLUENT 中它是作为解释式 UDF 被执行的。 10.2.1 给管道流动添加动量源项（Adding a Momentum Source to a Duct Flow） 当使用 UDF 模拟源项时，理解函数调用的内容是很重要的。当你增加了一个源 项，它在单元上执行全局的循环时，FLUENT 将调用你的函数。你的函数应只计 算源项并返回它到求解器。
*/
/********************************************************************
**/
#include "udf.h"
DEFINE_PROFILE(unsteady_velocity, thread, position)
{ face_t f;
Figure 10.1.4: Velocity Magnitude Contours for a Parabolic Inlet x Velocity
第 10 章 应用举例
Figure 10.1.5: Velocity Vectors for a Parabolic Inlet x Velocity
/* 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.; }
/***********************************************************************/
/* vprofile.c
*/
/* UDF for specifying steady-state velocity profile boundary condition
第 10 章 应用举例
函数通过声名变量 f 作为 face_t 的数据类型。一维数组 x 和变量 y 被定义为 real 数据类 型。循环宏用于循环区域内每个面来创建型线，或数据数组。每个循环内，F_CENTROIDS 为带指标 f 的面输出面质心（数组 x）的值，指标 f 在被 thread 指向的线（thread）上。 存储在 x[1]中的 y 坐标用于为变量 y 赋值，然后用于计算 x 速度。这个值接着被分配到 F_PROFILE，使用整数 position（在你从 Velocity Inlet 面板中选择的 UDF 作为 x 速度的 边界条件的基础上通过求解器传递给它）来设置内存中 x 速度面值。 为了在 FLUENT 中使用这个解释式 UDF，你必须首先编译它。 Define User-Defined Functions Interpreted...
Figure 10.1.2: Velocity Magnitude Contours for a Constant Inlet x Velocity
Figure 10.1.3: Velocity Vectors for a Constant Inlet x Velocity
现在入口 x 速度将用以下型线描述：
第 10 章 应用举例
在 X-Velocity 下拉列表下，选择 udf inlet_x_velocity，这个名字在上面已经命名给函数。 这个函数将被使用，而不是出现在 X-Velocity 区域的 0 值（在这个例子中）。点 OK 接受 这个新边界条件并关闭面板。 在求解运行到收敛之后，获得了一个修正的速度场如 Figure 10.1.4 and 10.1.5 所示。速度 场在入口中心显示了最大值，到边缘上降为 0。
在 Interpreted UDFs 面板中，在 Source File Name 区域命名你的函数。如果必要，在 CPP Command Name 区域输入你的 C 预处理程序类型和在 Stack Size 下输入堆栈大小。打开 Display Assembly Listing 在你的控制台窗口中能看到汇编列表，当函数编译时。点击 Compile 然后 Close 面板。 为了选择这个用户定义函数作为所选区域的速度边界条件，打开 Velocity Inlet 面板。
你可通过在 Files 下面的文本框中输入它的名字来读取输出文件和选择这个文件 并点击 Add....，然后点击 Plot，你可获得 Figure10.1.6 的图形显示。
第 10 章 应用举例
Figure 10.1.6: Average Velocity Magnitude at the Pressure Outlet
第 10 章 应用举例
在这个例子中，在 Report Of 下的下拉列表中选择 Velocity... and Velocity Magnitude。报告的区域是 pressure-outlet-5,，在 Surfaces 列表中选择它。 Area-Weighted Average 在 Report Type 的下拉列表中选择，Flow Time 在 X Axis 的下拉列表中选择。 一旦完成了你要求的迭代次数，监视器应出现在所选绘图窗口中。作为选择，你 可通过打开 File XY Plot 面板读入文件。 Plot File...
第 10 章 应用举例
第 10 章 应用举例
本章包含了 FLUENT 中 UDFs 的应用例子。 10.1 边界条件 10.2 源项 10.3 物理属性 10.4 反应速率（Reacting Rates） 10.5 用户定义标量（User_Defined Scalars） 10.1 边界条件 这部分包含了边界条件 UDFs 的两个应用。两个在 FLUENT 中都是作为解释式 UDFs 被执行的。 10.1.1 涡轮叶片的抛物线速度入口分布 要考虑的涡轮叶片显示在 Figure 10.1.1 中。非结构化网格用于模拟叶片周围的流 场。区域从底部周期性边界延伸到顶部周期性边界，左边是速度入口，右边是压 里变量 y 在人口中心是 0.0，在顶部和底部其值分别延伸到 ± 0.0745 。这样 x 速度在
入口中心为 20m/sec，在边缘为 0。 UDF 用于传入入口上的这个抛物线分布。C 源代码（vprofile.c）显示如下。函数使用 了 Section 5.3 中 描 述 的 Fluent 提 供 的 求 解 器 函 数 。
第 10 章 应用举例
begin_f_loop(f, thread) { real t = RP_Get_Real("flow-time"); F_PROFILE(f, thread, position) = 20. + 5.0*sin(10.*t); }
end_f_loop(f, thread) } 函数，被命名为 unsteady_velocity,使用 DEFINE_PROFILE macro 定义。The utility RP_Get_Real("flow-time") 用 于 访 问 real 流 动 时 间 ， 它 被 分 配 给 变 量 t 。 （RP_Get_Real 的详细内容见 Section 6.9）。 在你编译这个 UDF 之前，你必须在 Solver 面板中指定为非稳态流动计算。 Define Models Solver... 接下来你将打开 Interpreted UDFs 面板。在 Source File Name 下面的文本框内 输入函数的名字，如果必要，还有你的 C 预处理程序的名字。打开 Display Assembly Listing。点击 Compile 然后关闭面板。 Define User-Defined Functions Interpreted...
10.1.2 管内流动的瞬态速度入口分布 这个例子中，瞬时的周期性的速度边界条件使用 UDF 应用于管子的入口。速度形式为：
（10.1.1） 管子长 1m，半径为 0.2m。假设充满密度为 1kg/m3 和粘度为 2 ×10−5 kg/m-s 的空气。空
气速度关于平衡值 v0 = 20m / s 的波动振幅为 5m/s，频率为 10rad/s。
Figure 10.1.1: The Grid for the Turbine Vane Example
常数 x 速度应用于入口的流场与抛物线 x 速度应用于入口的流场作了比较。当采用分段
第 10 章 应用举例
线性分布的型线的应用是有效的对边界型线选择，多项式的详细说明只能通过用户定 义函数来完成。 常数速度应用于流场入口的结果显示在 Figure 10.1.2 和 Figure 10.1.3 中。当流动移动到 涡轮叶片周围时初始常速度场被扭曲。