| 查看: 359 | 回复: 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; } |
回复此楼» 猜你喜欢
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有6人回复
-
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有6人回复
-
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有10人回复
-
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）
已经有9人回复
-
有院领导为了换新车，用横向课题经费买了俩车
已经有10人回复
-
遇见不省心的家人很难过
已经有24人回复
-
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性
已经有6人回复
-
酰胺脱乙酰基
已经有13人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有10人回复
