24小时热门版块排行榜

返回列表

tjj_wei

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 293
帖子: 33
在线: 30.2小时
虫号: 2358097
注册: 2013-03-18
专业: 重力学

[求助] 蒸发的程序出错，顺便请大神指导下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号啊？请大家帮我看看

回复此楼