| 查看: 968 | 回复: 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). |
回复此楼» 猜你喜欢
070200求调剂，一志愿某211，288分
已经有20人回复
-
法国博士后职位
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有267人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有9人回复
