24小时热门版块排行榜

返回列表

lvhenghua

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 349
帖子: 26
在线: 21.9小时
虫号: 3554833
注册: 2014-11-23
专业: 化工系统工程

[交流] udf中的面宏如何使用，例如F_UDMI(f,t,i),和C_UDMI(c,t,i)的使用有什么区别？

请教一下各位fluent大神，udf中的面宏如何使用，例如F_UDMI(f,t,i),和C_UDMI(c,t,i)的使用有什么区别？附上自己模拟时的一段程序，一直运行不通，请大家多多指教。
DEFINE_INIT(my_init_func2,d)
{
face_t f;
Thread *t;
real xf[ND_ND];
int id_f=5;
t=Lookup_Thread(d,id_f);

thread_loop_f(t,d)
{
  begin_f_loop(f,t)
  {
F_CENTROID(xf,f,t);
F_UDMI(f,t,0)=(arrhenius_rate(tw0))*(c0-cs(tw0))*(c0-cs(tw0))+2*mw_h2o/mw_caso4*(arrhenius_rate(tw0))*(c0-cs(tw0))*(c0-cs(tw0));/*deposite rate md*/
F_UDMI(f,t,1)=0;/*length of foul*/
F_UDMI(f,t,2)=0;/*resistence of fouling layer*/
F_UDMI(f,t,3)=tw0;/*t=0,tf=320k*/
F_UDMI(f,t,4)=(arrhenius_rate(tw0))*(c0-cs(tw0))*(c0-cs(tw0))+2*mw_h2o/mw_caso4*(arrhenius_rate(tw0))*(c0-cs(tw0))*(c0-cs(tw0));/*mf*/
F_UDMI(f,t,5)=0;/*mr*/
  }
  end_f_loop(f,t)
}
}
DEFINE_EXECUTE_AT_END(execute_at_end2)
{

Domain *d;
Thread *t;
face_t f;
d=Get_Domain(1);

int id_f=5;

t=Lookup_Thread(d,id_f);
thread_loop_f(t,d)
{
begin_f_loop(f,t)
{
F_UDMI(f,t,1)+=F_UDMI(f,t,0)*delta_t/rho_f;

F_UDMI(f,t,2)=F_UDMI(f,t,1)/lambda_f;

F_UDMI(f,t,3)=F_T(f,t)+(tw0-F_T(f,t))/(1+h*F_UDMI(f,t,2));

F_UDMI(f,t,0)=(arrhenius_rate(F_UDMI(f,t,3)))*(F_YI(f,t,caso4)*F_R(f,t)-cs(F_T(f,t)))*(F_YI(f,t,caso4)*F_R(f,t)-cs(F_T(f,t)))+2*mw_h2o/mw_caso4*(arrhenius_rate(F_UDMI(f,t,3)))*(F_YI(f,t,caso4)*F_R(f,t)-cs(F_T(f,t)))*(F_YI(f,t,caso4)*F_R(f,t)-cs(F_T(f,t)));

F_UDMI(f,t,4)=F_UDMI(f,t,0)-F_UDMI(f,t,5);/*mf*/

F_UDMI(f,t,5)=0.012*pow(F_U(f,t),1.46)*F_UDMI(f,t,4)*(1+beta*(320-F_UDMI(f,t,3)))*dp*pow(F_R(f,t)*F_R(f,t)*0.0005806*9.81,1/3);/*mr*/
}
end_f_loop(f,t)
}

}
DEFINE_ADJUST(adjust_wall_temp,d)
{
Thread *t;
face_t f;
real time;
int id_f=5;
t=Lookup_Thread(d,id_f);
time=RP_Get_Integer("time_step"

;
thread_loop_f(t,d)
  {
begin_f_loop(f,t)
{
if (time<2)
   F_UDMI(f,t,6)=tw0;
else
   F_UDMI(f,t,6)=F_T(f,t)+(tw0-F_T(f,t))/(1+h*F_UDMI(f,t,2));/*adjust wall temperature*/
}
end_f_loop(f,t)
  }
}
DEFINE_PROFILE(wall_temp,t,position)
{
face_t f;
Domain *d;
int id_f=5;
d=Get_Domain(1);
t=Lookup_Thread(d,id_f);
begin_f_loop(f,t)
{
  F_PROFILE(f,t,position)=F_UDMI(f,t,6);
}
end_f_loop(f,t)
}

回复此楼

» 猜你喜欢

现代”学阀”该如何界定已经有8人回复
课题组招2027级博士上海工程技术大学激光智能制造方向已经有3人回复
上海工程技术大学激光智能制造课题组 2027级博士研究生招生已经有3人回复
实名举报：华南理工大学物理与光电学院副院长李志远婚内两次出轨女学生（博士）已经有16人回复
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告已经有11人回复
上海工程技术大学激光智能制造课题组招收博士研究生已经有10人回复
申博发邮件已经有11人回复
各位大神，目前国内有哪些比较好用的逆合成软件？已经有12人回复
我的奶奶已经有3人回复
国社科系统bug了，是不是要放榜了？已经有9人回复

1楼 2015-12-25 09:43:59

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

相关版块跳转我要订阅楼主 lvhenghua 的主题更新

返回列表