| 查看: 1867 | 回复: 1 | ||
| 【悬赏金币】回答本帖问题，作者951223将赠送您 20 个金币 | ||
951223
|
[求助]
大家有没有用fluent作吸附模拟的
|
|
|
#include"udf.h" enum {xifu,d_xifu}; DEFINE_ADJUST(xifu,domain) { Thread *t; cell_t c; real n_m; real E_m; real x_p; real A; real x_b; real n; thread_loop_c (t,domain) { if (NULL!= THREAD_STORAGE(t,SV_UDS_I(xifu))&& NULL!=T_STORAGE_R_NV(t,SV_UDSI_M1(xifu))) { begin_c_loop (c,t) { real tem=C_T(c, t); n_m=25.4; E_m=11850+7.05*tem; x_p=1.85E9/ C_UDSI(c,t,p); A=8.31429*tem*log(x_p); x_b=-A*A/E_m/E_m; C_UDSI(c,t,xifu)=n_m*exp(x_b); } end_c_loop (c,t) } } } DEFINE_ADJUST(d_xifu, domain) { Thread *t; cell_t c; thread_loop_c (t,domain) { if(NULL!=THREAD_STORAGE(t,SV_UDS_I(d_xifu))&& NULL!=T_STORAGE_R_NV(t,SV_UDSI_M1(d_xifu))) { begin_c_loop (c,t) { real physical_dt; physical_dt=RP_Get_Real("physical-time-step"  ; C_UDSI(c,t,d_xifu)=-10.56*(C_UDSI(c,t,xifu)-C_UDSI_M1(c,t,xifu))/physical_dt;} end_c_loop (c,t) } } } DEFINE_SOURCE(mass_source,c,t,dS,eqn) { real m_source; m_source=C_UDSI_M1(c,t,d_xifu); dS[eqn] = 0; return m_source; } DEFINE_SOURCE(energy_source,c,t,d S,eqn) { real e_source; e_source=-C_UDSI_M1(c,t,d_xifu)*84143E3; dS[eqn] = 0; return e_source; } 这是程序 下面是问题 大家能不能帮我看一下，谢谢各位  捕获.PNG  捕获.PNG |
回复此楼» 猜你喜欢
24年 申博 化学/材料 一作6篇sci
已经有7人回复
-
博士招生
已经有6人回复
-
普通院校药学硕士，做合成的，感觉找不到工作
已经有4人回复
-
大家访学都是怎么找的啊？
已经有3人回复
-
两类问题算是白选了~
已经有9人回复
-
前同事拿我以前未中的国自然标书申请今年的国自然，现在基金委发来调查函
已经有47人回复
-
考研，求职还是考编?
已经有19人回复
-
期刊推荐
已经有7人回复
-
真的好想读博！
已经有10人回复
-
博士论文被抄袭
已经有31人回复
 |
2楼2018-07-19 09:35:16
10