| 查看: 1243 | 回复: 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号啊？请大家帮我看看 |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
今天基金会出结果吗？20260819
已经有12人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有121人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有6人回复
-
今天放榜吗？
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教使用旋转蒸发仪的问题!
已经有11人回复
2楼2013-11-04 10:23:10
3楼2013-11-04 10:57:31
4楼2018-05-18 16:17:39
