自然对流换热壁面UDF函数求助

jinzhp

各位大侠，小弟现在做热仿真，用到壁面边界条件为自然对流换热系数，公式为a=C1*(C2*T+C3)^0.25，C1、C2、C3为我计算中用到的系数，T为壁面温度。求哪位大侠赐程序！不胜感激！

jinzhp

#include "udf.h"
#include <math.h>
#include DEFINE_PROFILE(Coefficient,t,i)
{
real te;
te=C_T(c,t);
face_t f;
begin_f_loop(f,t)
{
F_CENTROID(x,f,t);
F_PROFILE(f,t,i)=1.127408803*pow((te-40),0.25);
}
end_f_loop(f,t)
}
请各位大侠看看程序有什么问题吗

JohnQiang

仅供参考：
x没定义，real x【ND_ND】；
貌似C_T（c,t）中的c应该定义一下吧；
你的公式中不含坐标变量，F_CENTROID(x,f,t);也就没用

jinzhp

#include "udf.h"
DEFINE_PROFILE(Coefficient,t,i)
{
real te;
face_t f;
te=F_T(f,t);
begin_f_loop(f,t)
{
F_PROFILE(f,t,i)=1.127408803*pow((te-40),0.25);
}
end_f_loop(f,t)
}