| 查看: 1734 | 回复: 2 | ||
| 【悬赏金币】回答本帖问题，作者无敌小张！将赠送您 5 个金币 | ||
[求助]
fluent的udf
|
||
|
大神们，有段fluent的udf，但是看不懂，有没有大神能给解答一下是什么意思 研一萌新没啥金币，直接分享，大家有需要就拿去用 #include "udf.h" #define U 18 #define I 200 #define PI 3.1415928 #define R 0.003 #define n 0.8 #define v 0.003 #define L 0.005 #define u0 1.26e-6 DEFINE_SOURCE(S_X,cell,thread,dS,eqn) { real xc[ND_ND]; real t; real source; real x,y,z; real r; C_CENTROID(xc,cell,thread); x=xc[0]; y=xc[1]; z=xc[2]; t=RP_Get_Real("flow-time"  ;r=sqrt((x-v*t)*(x-v*t)+y*y); if(r<R) { source=-1*u0*I*I/(4*PI*PI*R*R*r)*exp(-1*r*r/(2*R*R)) *(1-exp(-1*r*r/(2*R*R)))*(1-z/L)*(1-z/L)*x/L; dS[eqn]=0; } else source=dS[eqn]=0; return source; } DEFINE_SOURCE(S_Y,cell,thread,dS,eqn) { real xc[ND_ND]; real t; real source; real x,y,z; real r; C_CENTROID(xc,cell,thread); x=xc[0]; y=xc[1]; z=xc[2]; t=RP_Get_Real("flow-time" ;r=sqrt((x-v*t)*(x-v*t)+y*y); if(r<R) { source=-1*u0*I*I/(4*PI*PI*R*R*r)*exp(-1*r*r/(2*R*R)) *(1-exp(-1*r*r/(2*R*R)))*(1-z/L)*(1-z/L)*y/L; dS[eqn]=0; } else source=dS[eqn]=0; return source; } DEFINE_SOURCE(S_Z,cell,thread,dS,eqn) { real xc[ND_ND]; real t; real source; real x,y,z; real r; C_CENTROID(xc,cell,thread); x=xc[0]; y=xc[1]; z=xc[2]; t=RP_Get_Real("flow-time" ;r=sqrt((x-v*t)*(x-v*t)+y*y); if(r<R) { source=-1*u0*I*I/(4*PI*PI*L*r*r) *(1-exp(-1*r*r/(2*R*R)))*(1-z/L)*(1-z/L); dS[eqn]=0; } else source=dS[eqn]=0; return source; } DEFINE_PROFILE(heat_flux,thread,index) { real xc[ND_ND]; face_t f; real t; real heat,T0=300; real x,y,z; real r; real hc = 80, sigmab = 5.67e-8, epsilon = 0.06; t=RP_Get_Real("flow-time" ;begin_f_loop(f,thread) { F_CENTROID(xc,f,thread); x=xc[0]; y=xc[1]; z=xc[2]; r=sqrt((x-v*t)*(x-v*t)+y*y); if(r<sqrt(6)*R) heat=n*I*U/(2*PI*R*R)*exp(-1*r*r/(2*R*R)); else heat=0; F_PROFILE(f, thread, index) = heat - hc*(F_T(f,thread), - T0) - sigmab*epsilon*(pow(F_T(f,thread), 4.) - pow(T0, 4.)); } end_f_loop(f,thread) } |
回复此楼» 猜你喜欢
0857 资源与环境 285分
已经有3人回复
-
标题:捡漏预警|08工科/09农学调剂!英语要求低，过线即有机会!
已经有12人回复
-
一志愿安徽大学材料工程专硕313分，求调剂的学校
已经有8人回复
-
313求调剂
已经有4人回复
-
083000环境科学与工程调剂
已经有8人回复
-
材料调剂，307分
已经有12人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
070300化学 280 一志愿太原理工 求调剂
已经有9人回复
-
一志愿北京化工大学材料与化工296分求调剂
已经有12人回复
-
哈工大材料324求调剂
已经有8人回复
2楼2021-09-12 00:03:46
3楼2021-10-25 09:00:32
