新手压力入口问题编程

花花家园

压力在0.01s之前是随时间变化，0.01s之后就一直维持在0.01s之时的压力值。这个需要用到udf吗？下面的这个udf对吗？错的话还请高手指点！
#include"udf.h“
DEFINE_PROFILE(unsteady_pressure,thread,position)
{
cell_t_c;
begin_c_loop(c,thread)
{
realt=RP_Get_Real("flow-time");
t<=0.01;
C_PROFILE(c,thread,position)=(-60000*t*t+120*t)*1000000;
t>=0.01;
C_PROFILE(c,thread,position)=600000;
}
end_f_loop(c,thread)
}

gbd20022002

#include"udf.h“
DEFINE_PROFILE(unsteady_pressure,thread,position)
{
face_t f;
real t=RP_Get_Real("flow-time");
begin_f_loop(f,thread)
{
t<=0.01;
F_PROFILE(f,thread,position)=(-60000.*t*t+120.*t)*1000000.;
t>=0.01;
F_PROFILE(f,thread,position)=600000.;
}
end_f_loop(f,thread)
}

[ 本帖最后由 gbd20022002 于 2009-9-4 11:06 编辑 ]

花花家园

那你刚说的整型实型是怎么回事？程序中需要说明下吗！

gbd20022002

我重新修改了下，你看看吧，整实混合运算不影响大局，只是细节问题，系统会转换的。
复制计算一下

花花家园

还有一个，我那个是个单元的，也就是三维的，如果换做你的，都用face的不就不合适了嘛！

gbd20022002

你可以去看看fluent帮助，看看网格拓扑
你定义的是三维的压力入口，必然是一个面——不论是平面还是曲面——下面的循环便是给此面中所有的面元赋值

花花家园

还有一个问题，我把这段程序放进c程序中调试，总是出现[Error] C:\Documents and Settings\zhou\My Documents\C-Free\Temp\未命名1.cpp:1: udf.h: No such file or directory，这个事怎么回事？会不会影响我的下一步进程啊！

gbd20022002

请以英文命名，后缀取.c
直接在fluent编译即可……

[ 本帖最后由 gbd20022002 于 2009-9-4 11:26 编辑 ]

花花家园

哦，这样子啊！那按照你的意思，如果要用到cell就是说要定义一个体的压力了，而非单纯的入口或出口。明白了！

花花家园

哦，果然，我起初没有保存这段程序，直接就把它写在上面，它就有了个”未命名“三个中文字，怪不得，现在没有错误了！谢谢了！

花花家园

我又遇到了一个新的问题，在我迭代的时候它出来一个TIME的窗口，这里的Time Step Size 和Number of Time Steps是怎么设定的，是随意设定吗？