| 查看: 953 | 回复: 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). |
回复此楼» 猜你喜欢
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有69人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
-
英国全奖博士招聘-深度学习与量子物理
已经有0人回复
-
间接带隙半导体有效质量求助
已经有0人回复
