本人最近做仓储粮堆得温度场模拟,想给仓壁设为温度随时间变化而变化的边界条件!
编写了UDF如下/***********************************************************************
myudf.c
1 UDF for specifying unsteady temperature profile boundary condition
************************************************************************/
#include "udf.h"
DEFINE_PROFILE(unsteady_temperature, thread, position)
{
real t,temperature;
face_t f;
begin_f_loop(f, thread)
{
t=RP_Get_Real(“flow-time”);
temperature=-0.4821*(t/2592000-11.382)*(t/2592000-11.382)+316.36;
F_PROFILE(f,thread,position)=temperature;
}
end_f_loop(f, thread)
}
但带入Fluent12.0时候总是显示错误,说的错误是在第一行,不知道第一行有什么错误!实在找不出这个UDF哪儿出了错误。
而且当我用下面这个UDF带入的时候又显示编译成功,
/***********************************************************************
myudf.c
1 UDF for specifying unsteady temperature profile boundary condition
************************************************************************/
#include "udf.h"
DEFINE_PROFILE(unsteady_temperature, 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)=-0.4821*(y/2592000-11.382)*(y/2592000-11.382)+316.36;
}
end_f_loop(f, thread)
}
坐等高手来帮我看看下!
我实在是整不明白为什么上面一个编译不了,下面一个就可以编译了!
刚接触CFD,菜鸟啊!请高手帮我指点一二啊!非常感谢,金币大大的送上! |