关于根据位置改变流体粘度的udf

gloomyjcjd

自己写了段udf  想根据流体在通道中不同位置改变其粘度   但是编译通过后 ，在初始化时总会报错
希望哪位大神来指点一二

yxyzr

发过来瞧瞧

gloomyjcjd

#include "udf.h"
DEFINE_PROPERTY(cell_viscosity, cell, thread)
{
real mu_lam;
Thread *t;
face_t f;
real x[ND_ND];
real y;

F_CENTROID(x,f,t);
y=x[1];
if(y>=0.0018)
mu_lam=0.0009;
else if(y<0.0016)
mu_lam=0.0001;
else
mu_lam=4*y-0.0063;

return mu_lam;
}

cfd-soft

F_CENTROID(x,f,t);改成C_CENTROID(x,f,t);

gloomyjcjd

多谢 现在能编译通过计算了   但是   条件语句似乎没有生效    粘度只是选择了一个条件然后取值。。。


#include "udf.h"
DEFINE_PROPERTY(cell_viscosity, cell, thread)
{
real mu_lam;
cell_t c;
real x[ND_ND];
real y;

C_CENTROID(x,c,thread);
y=x[1];
if(y<0.0016)
mu_lam=1;
else
mu_lam=11;

return mu_lam;
}