| 查看: 1011 | 回复: 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). |
回复此楼» 猜你喜欢
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有244人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有113人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有69人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有39人回复
