| 查看: 352 | 回复: 0 | |||
[求助]
fluentUDF传质设置时出错
|
|
udf代入FLUENT中解释时没有错误，但设置传质时出现错误？ D:\PROGRA~2\ANSYSI~1\v145\fluent\fluent14.5.0\win64\2d\fl1450s.exe received fatal signal (ACCESS_VIOLATION) 1. Note exact events leading to error. 2. Save case/data under new name. 3. Exit program and restart to continue. 4. Report error to your distributor. #include"udf.h" #define T_SAT 373 /*定义饱和温度*/ #define lam 100 /*定义松弛因子*/ #define LAT_H 2448 /*定义气液相变潜热*/ DEFINE_PROPERTY(density_tem,c,t) { real density,k1,cp,v; density=0.34710*pow(0.27400,-(pow(1-C_T(c,t)/647.13,0.28571))); return density; /*物性参数随温度变化*/ } DEFINE_PROPERTY(thermal_conductivity_tem,c,t) { real thermal_conductivity; thermal_conductivity=-0.2758+0.0046120*C_T(c,t)-0.000005391*pow(C_T(c,t),2); return thermal_conductivity; } DEFINE_PROPERTY(cp_tem,c,t) { real cp; cp=(92.053-0.039953*C_T(c,t)-0.00021103*pow(C_T(c,t),2)-0.00000053469*pow(C_T(c,t),3))*(1000/18); return cp; } DEFINE_PROPERTY(viscosity_tem,c,t) {real viscosity; viscosity=pow(10,-10.2158+1792.5/C_T(c,t)+0.017730*C_T(c,t)-0.000012631*pow(C_T(c,t),2)); return viscosity; } DEFINE_PROPERTY(VCp_tem,c,t) {real VCp; VCp=33.933-0.0084186*(C_T(c,t))+2.9906e-05*pow(C_T(c,t),2)-1.7825e-08*pow(C_T(c,t),4)+3.6934e-12*pow(C_T(c,t),4)*(1000/18); return VCp; } DEFINE_PROPERTY(VViscosity_tem,c,t) { real VViscosity; VViscosity=-36.826+4.2900e-01*(C_T(c,t))-1.6200e-05*pow(C_T(c,t),2); return VViscosity; } DEFINE_PROPERTY(VThermal_conductivity_tem,c,t) {real VThermal_conductivity; VThermal_conductivity=-0.2758+4.6120e-03*(C_T(c,t))-5.5391e-06*pow(C_T(c,t),2); return VThermal_conductivity; } DEFINE_SOURCE(liq_src,cell,first,ds,eqn) { Thread *mixer,*second; real m_dot_first; mixer=THREAD_SUPER_THREAD(first); second=THREAD_SUB_THREAD(mixer,1); if(C_T(cell,mixer)>=T_SAT) { m_dot_first=-lam*C_VOF(cell,first)*C_R(cell,first)*fabs(C_T(cell,first)-T_SAT)/T_SAT; ds[eqn]=-lam*C_VOF(cell,first)*fabs(C_T(cell,first)-T_SAT)/T_SAT; } else { m_dot_first=lam*C_R(cell,first)*fabs(T_SAT-C_T(cell,first))/T_SAT; ds[eqn]=0; } return m_dot_first; } DEFINE_SOURCE(vap_src,cell,second,ds,eqn) { Thread *mixer,*first; real m_dot_second; mixer=THREAD_SUPER_THREAD(second); first=THREAD_SUB_THREAD(mixer,0); if(C_T(cell,mixer)<=T_SAT) { m_dot_second=lam*C_VOF(cell,first)*C_R(cell,first)*fabs(T_SAT-C_T(cell,mixer))/T_SAT; ds[eqn]=0; } else { m_dot_second=-lam*C_VOF(cell,first)*C_R(cell,first)*fabs(T_SAT-C_T(cell,mixer))/T_SAT; ds[eqn]=-lam*C_R(cell,first)*fabs(T_SAT-C_T(cell,mixer))/T_SAT; } return m_dot_second; } DEFINE_SOURCE(enrg_src,cell,mixer,ds,eqn) { Thread *first,*second; real m_dot; first=THREAD_SUB_THREAD(mixer,0); second=THREAD_SUB_THREAD(mixer,1); if(C_T(cell,mixer)>=T_SAT) { m_dot=-lam*C_VOF(cell,first)*C_R(cell,first)*fabs(C_T(cell,first)-T_SAT)/T_SAT; ds[eqn]=-lam*C_VOF(cell,first)*(C_T(cell,first)-T_SAT)/T_SAT; } else if(C_T(cell,mixer)<=T_SAT) { m_dot=lam*C_VOF(cell,first)*C_R(cell,first)*fabs(T_SAT-C_T(cell,mixer))/T_SAT; ds[eqn]=lam*C_VOF(cell,first)*(T_SAT-C_T(cell,mixer))/T_SAT;} return LAT_H*m_dot; } |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有11人回复
-
救命帖
已经有5人回复
-
限项规定
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
招博士
已经有3人回复
-
存款400万可以在学校里躺平吗
已经有35人回复
-
最失望的一年
已经有18人回复
-
求推荐博导
已经有4人回复
-
求推荐英文EI期刊
已经有5人回复
-
疑惑？
已经有5人回复
