| 查看: 962 | 回复: 0 | ||
[求助]
fluent 多相流源项UDF求助？？？
|
|
#include "udf.h" DEFINE_SOURCE(source_char_com,c,t,dS,eqn) { Thread **pt = THREAD_SUB_THREADS(t); Thread *tp = pt[0]; /* gas phase */ Thread *ts = pt[1]; /* solid phase */ real diam,Ts,phi,p,Vp,Re,Sc,Sh,kr,km,Dg,epsilon,kc,r; diam=C_PHASE_DIAMETER(c,ts); Ts=C_T(c,ts); p=2500*exp(-5.19E04/(UNIVERSAL_GAS_CONSTANT*Ts)); Vp=pow(C_V(c,ts)*C_V(c,ts)+C_U(c,ts)*C_U(c,ts),0.5); epsilon=C_VOF(c,tp); if(diam<0.00005) phi=(2*p+2.0)/(p+2.0); else if(diam>0.001) phi=1; else phi=(2*p+2.0-p*(100000*diam-0.005)/0.095)/(p+2); Re=C_R(c,tp)*Vp*diam/C_MU_L(c,tp); Dg=3.13e-04*pow((C_T(c,tp)/1500),1.75)*(101325.0/C_P(c,tp)); Sc=C_MU_L(c,tp)/C_R(c,tp)/Dg; Sh=2*epsilon+0.69*pow(Sc,1.0/3.0)*pow(Re/epsilon,0.5); km=phi*Sh*Dg*12.0/diam*UNIVERSAL_GAS_CONSTANT* C_T(c,tp) ; kr=8910*(exp(-1.4974e8/(UNIVERSAL_GAS_CONSTANT*Ts))); kc=UNIVERSAL_GAS_CONSTANT*Ts/12.0/(1.0/km+1.0/kr); r=6*(1-epsilon)*0.2505/diam/C_R(c,tp)*kc*1.43; return -r; dS[eqn]=0; } 请问这段UDF哪里出问题了啊？运行就会出现Error: received a fatal signal (Segmentation fault). |
回复此楼» 猜你喜欢
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有242人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有21人回复
-
欢迎加入课题组
已经有0人回复
