24小时热门版块排行榜

返回列表

zq91000

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 82.2
散金: 5
帖子: 41
在线: 22.9小时
虫号: 4612611
注册: 2016-04-23
性别: GG
专业: 宏观经济学

[求助] FLUENT中UDF的编码问题

虫友们好，我这有个问题想请大家指正，，FLUENT中的UDF的代码如下面所示，
而导入到FLUENT中时出现的问题截图如图所示。
我个人感觉第四行没问题，因为这是一个宏，括号里面是两个参数，前一个是名字，后一个是FLUENT传递过来的东东，
请大家指点，不胜感激！

#include "udf.h"
#define d 0.1
#define Tm 298.15
DEFINE_ADJUST(laminar_z,d)
{
Thread *t;
real sum_diss=0.;

cell_t c;

thread_loop_c(t,d)
{
begin_c_loop(c,t)
{
real dux=C_DUDX(c,t);
real duy=C_DUDY(c,t);
real duz=C_DUDZ(c,t);
real dvx=C_DVDX(c,t);
real dvy=C_DVDY(c,t);
real dvz=C_DVDZ(c,t);
real dwx=C_DWDX(c,t);
real dwy=C_DWDY(c,t);
real dwz=C_DWDZ(c,t);
real dwz=C_DWDZ(c,t);
real V=C_VOIUME(c,t);
real mu=C_MU_L(c,t);

C_UDMI(c,t,0)=mu*(2*(dux*dux+dvy*dvy+dwz*dwz)+(duy+dvx)*(duy+dvx)+(duz+dwx)*(duz+dwx)+(dvz+dwy)*(dvz+dwy))/Tm;
sum_diss+=C_VOLUME(c,t)*C_UDMI(c,t,0);
}
end_c_loop(c,t);
}

printf("Entropy production for laminar flow: %g\n",sum_diss);
}

错误第四行.JPG

回复此楼