| 查看: 959 | 回复: 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). |
回复此楼» 猜你喜欢
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有217人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有20人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有1人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有0人回复
