| 查看: 1091 | 回复: 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号啊？请大家帮我看看 |
回复此楼» 猜你喜欢
求调剂
已经有6人回复
-
22408 266求调剂
已经有7人回复
-
求调剂
已经有14人回复
-
273求调剂
已经有43人回复
-
307求调剂
已经有13人回复
-
327求调剂
已经有7人回复
-
338求调剂
已经有6人回复
-
326分，一志愿沪9，求生物学调剂
已经有3人回复
-
材料科学与工程320求调剂，080500
已经有9人回复
-
0703化学调剂 348分
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教使用旋转蒸发仪的问题!
已经有11人回复
3楼2013-11-04 10:57:31
2楼2013-11-04 10:23:10
4楼2018-05-18 16:17:39
