| 查看: 1498 | 回复: 0 | ||
[求助]
UDF编写化学反应速率加载到fluent内出现问题
|
|
最近用udf编写化学反应速率导入到fluent内,load后一直出现Error: The UDF library you are trying to load (libudf) is not compiled for 2ddp on the current platform (win64)。UDF如下: #include "udf.h" #define CH4 0 #define O2 1 #define CO2 2 #define H2O 3 #define N2 4 DEFINE_VR_RATE(vol_reac_rate,c,t,r,wk,yk,rate,rr_t) { int i,zone_ID=4; real CH4_mole_fract,O2_mole_fract,T,P_total,P_CH4,P_O2,k1,k2,k3,k4; real total_mole; Thread *t1; Material *mix_mat,*sp; *t1=Lookup_Thread(domain,zone_ID); mix_mat=mixture_material(t1); total_mole=0; mixture_species_loop(mix_mat,sp,i) { wk=MATERIAL_PROP(sp,PROP_wk); total_mole+=C_YI(c,t,i)/wk; } CH4_mole_fract=yk[CH4]/wk[CH4])/total_mole; O2_mole_fract=yk[O2]/wk[O2]/total_mole; P_total=ABS_P(C_P(c,t),op_pres); P_CH4=P_total*CH4_mole_fract; P_O2=P_total*O2_mole_fract; T=C_T(c,t); k1=5.852e17*exp(-204000/8.314/T); k2=4.02e5*exp(103500/8.314/T); k3=5.08e5*exp(66200/8.314/T); k4=1+k2*P_CH4+k3*pow(P_O2,0.5); if (FLUID_THREAD_P(t)) *rate = k1*P_CH4*pow(P_O2,0.5)/pow(k4,2); else *rate = 0.; *rr_t = *rate; } 在帮助文档内找了一个简单的编译型UDF能成功load,初步判定是UDF程序编写有问题,还望有好心人能给看下是什么个问题,在此提前拜谢了! |
» 猜你喜欢
青椒八年已不青,大家都被折磨成啥样了?
已经有5人回复
青年基金C终止
已经有3人回复
26申博求博导推荐-遥感图像处理方向
已经有4人回复
限项规定
已经有7人回复
救命帖
已经有8人回复
招博士
已经有4人回复
西南交通大学国家级人才团队2026年博士研究生招生(考核制)—机械、材料、力学方向
已经有3人回复
英文综述是否需要润色及查重
已经有5人回复
为什么nbs上溴 没有产物点出现呢
已经有9人回复













回复此楼