24小时热门版块排行榜    

查看: 172  |  回复: 0

正林时代

新虫 (小有名气)

[交流] udf解释错误,麻烦大神给看看

提示第27行错误  Error: F:\example\xin\profile.c: line 27: parse error.

/*profile.c*/
/*UDF for specifying steady-state velocity profile boundary condi
#include"udf.h"
#define zo 10.0
#define uo 23.7
#define zb 5.0
#define zg 450.0
#define tur 0.31
#define cmu 0.09
#define a 0.16
DEFINE_PROFILE(inlet_z_velocity,thread,index)
{
real x[ND_ND];
real z;
face_t f;
begin_f_loop(f,thread)
{
F_CENTROID(x,f,thread);
z=x[2];
F_PROFILE(f,thread,index)=uo*pow(z/zo,a);
}
end_f_loop(f,thread)
}

/*PROFILE FOR KINETIC ENERGY*/
DEFINE_PROFILE(k_profile,thread,index)
{
real x[ND_ND];
real z,lr,um;
face_t f;
begin_f_loop(f,thread)
{
F_CENTROID(x,f,thread);
z=x[2];
if(z<=zb)
lr=tur*pow(1.0,1.);
else
lr=0.1*pow(z/zg,-0.25);
um=uo*pow(z/zo,a);
F_PROFILE(f,thread,index)=1.5*pow(um*lr,2.);
}
end_f_loop(f,thread)
}

/*PROFILE FOR DISSIPATION RATE*/
DEFINE_PROFILE(E_profile,thread,index)
{
real x[ND_ND];
real z,lr,um,l,k;
face_t f;
begin_f_loop(f,thread)
{
F_CENTROID(x,f,thread);
z=x[2];
if(z<=zb)
lr=tur*pow(1.0,1.);
else
lr=0.1*pow(z/zg,-0.25);
um=uo*pow(z/zo,a);
k=1.5*pow(um*lr,2.);
l=100*pow(z/30,0.5);
F_PROFILE(f,thread,index)=pow(cmu,3/4)*pow(k,3/2)/l;
}
end_f_loop(f,thread)
}
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 正林时代 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见