| 查看: 1090 | 回复: 3 | ||
[求助]
蒸发的程序出错，顺便请大神指导下C_YI的用法
|
|
#include "udf.h" //包括常规宏 #include "sg_mphase.h" //包括体积分数 #include "mem.h" //包括温度及组分质量分数 #include "sg.h" #include "flow.h" #define T_SAT 373 //定义蒸发温度100℃ #define LAT_HT 1.e3 DEFINE_SOURCE(liq_src,cell,sec_th,dS,eqn) { Thread *mix_th,*pri_th; //定义计算区域指针 real m_dot_l; //定义液相质量转移kg/(m2*s) mix_th=THREAD_SUPER_THREAD(sec_th); pri_th=THREAD_SUB_THREAD(mix_th,2); if(C_T(cell,mix_th)>=T_SAT) { m_dot_l=-0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=-0.1*C_VOF(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; } else { m_dot_l=0.1*C_YI(cell,pri_th,2)*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=0; } return m_dot_l; } DEFINE_SOURCE(vap_src,cell,pri_th,dS,eqn) { Thread *mix_th,*sec_th; //定义计算区域指针 real m_dot_v; //定义气相质量转移kg/(m2*s) mix_th=THREAD_SUPER_THREAD(pri_th); sec_th=THREAD_SUB_THREAD(mix_th,3); if(C_T(cell,mix_th)>=T_SAT) { m_dot_v=0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=0; } else { m_dot_v=-0.1*C_YI(cell,pri_th,2)*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=-0.1*C_YI(cell,pri_th,2)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; } return m_dot_v; } DEFINE_SOURCE(enrg_src,cell,mix_th,dS,eqn) { Thread *pri_th,*sec_th; //定义计算区域指针 real m_dot; //定义液相质量转移kg/(m2*s) pri_th=THREAD_SUPER_THREAD(mix_th,2); sec_th=THREAD_SUB_THREAD(mix_th,3); if(C_T(cell,mix_th)>=T_SAT) { m_dot=-0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=-0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)/T_SAT; } else { m_dot=0.1*C_YI(cell,pri_th,2)*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=0.1*C_YI(cell,pri_th,2)*C_VOF(cell,pri_th)*C_R(cell,pri_th)/T_SAT; } return LAT_HT*m_dot; } DEFINE_SOURCE(h2o_src,cell,pri_th,dS,eqn) { Thread *mix_th,*sec_th; //定义计算区域指针 real m_dot_vw; //定义气相质量转移kg/(m2*s) mix_th=THREAD_SUPER_THREAD(pri_th); sec_th=THREAD_SUB_THREAD(mix_th,3); if(C_T(cell,mix_th)>=T_SAT) { m_dot_vw=0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=0; } else { m_dot_vw=-0.1*C_YI(cell,pri_th,2)*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=-0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; } return m_dot_vw; } 这是我的程序的内容，一初始化就出错，不知道是不是我的C_YI事业有问题，还有第一相是0还是查的ID号啊？请大家帮我看看 |
回复此楼» 猜你喜欢
一志愿211，化学学硕，310分，本科重点双非，求调剂
已经有13人回复
-
285求调剂
已经有8人回复
-
22408 266求调剂
已经有6人回复
-
264求调剂
已经有8人回复
-
材料工程专业日语生求调剂
已经有9人回复
-
297求调剂
已经有15人回复
-
277、学硕，求调剂 数一104，
已经有12人回复
-
机械专硕273请求调剂
已经有3人回复
-
080100力学316求调剂
已经有3人回复
-
288环境专硕,求调材料方向
已经有39人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教使用旋转蒸发仪的问题!
已经有11人回复
2楼2013-11-04 10:23:10
3楼2013-11-04 10:57:31
4楼2018-05-18 16:17:39
