速度随时间变化的UDF(1)

管道入口瞬态速度分布在本例中，入口速度满足关系式：，管道长1m，半径0.2m，管道内流体为空气，密度1 kg/m3，粘性系数2105 kg/m-s。入口速度分布与时间有关，随时间正弦变化。C源程序如下：/*/* unsteady.c */* UDF for specifying a transient velocity profile boundary condition */*/#include udf.hDEFINE_PROFILE(unsteady_velocity, thread, position)face_t f;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，变量flow-time存贮流场当前历经时间，函数RP_GET_REAL得到当前时间，由于本例为非稳态问题，需要选择非稳态解法器，如下激活非稳态解法器：DefineModelsSolver本例使用Interpreted型UDFs，在面板Interpreted UDFs中编译连接UDF，在面板Velocity Inlet中XVelocity下拉列表选择函数udf unsteady_velocity。进行FLUENT默认初始化之后，就可以在Iterate面板进行迭代计算。SolveIterate本例取时间步长为0.0314s，总计算次数为60，则计算总时间为0.031460s。每时间步长最多迭代20次，每次迭代前都要更新边界条件并且输出结果。计算60步之后，我们就可以检查压力出口的速度大小。如果要在计算过程中，查看相关信息，需要事先设定。打开Surface Monitors面板：SolveMonitorsSurface把Surface Monitors设为1，monitor1就可以使用了。我们可以输入新文件名代替monitor1，然后选择是Plot，Print，还是Write。Every有两个选项，分别为Iteration，Flow Time或Flow Time，用来定义监测流场相关变量的时间间隔。点击Define之后，出现Define Surface Monitor 面板：在相应下拉列表中选择Velocity和Velocity Magnitude，选择所要输出参数的面pressure outlet5。监测输出参数类型为平均值（Average），x轴取流动时间。设置好之后，就会在每一步时间步长都输出速度大小，以供计算时检测。我们也可以通过File XY Plot面板，手工绘制文件monitor-1.out的图形。PlotFile在Files里选择文件，如果需要添加文件点击按钮add，选择文件之后点击Plot按钮，可以得到文件输出的图形。图6.1.6输出不同时刻本例压力出口的速度大小。图6.1.6