已知入口处风速随高度变化如函数所示:u=4.57*(0.1*x)^0.1677。希望好心人写下代码。
#include"udf.h"
DEFINE_PROFILE(velocity,t,i)
{
real x[ND_ND];
real y;
face_t f;
begin_f_loop(f,t)
{
F_CENTROID(x,f,t);
y=x[1];(x[0]x[1]x[2]分别表示xyz轴,你自己根据模型看)
F_PROFILE(f,t,i)=4.57*(0.1*y)^0.1677;(y单位是m)
}
end_f_loop(f,t)
}
加载到fluent,在速度上选择这个udf就可以了。
DEFINE_PROFILE(velocity,t,i)
{
real x[ND_ND];
real y;
face_t f;
begin_f_loop(f,t)
{
F_CENTROID(x,f,t);
y=x[1];(x[0]x[1]x[2]分别表示xyz轴,你自己根据模型看)
F_PROFILE(f,t,i)=4.57*(0.1*y)^0.1677;(y单位是m)
}
end_f_loop(f,t)
}
加载到fluent,在速度上选择这个udf就可以了。
温馨提示:答案为网友推荐,仅供参考
第1个回答 2021-05-20