24小时热门版块排行榜

返回列表

gky227227

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 10
散金: 81
红花: 3
帖子: 304
在线: 90小时
虫号: 807273
注册: 2009-07-11
专业: 化学反应工程

[求助] 相变udf求助

最近在做相变的模拟，udf使用的是论坛上大家通用的，我只是改了一下指针，结果就老是初始化报错，把udf放在下面，希望大家给看一下啊，谢谢各位，我个人感觉是温度的指针出了问题。
DEFINE_SOURCE(liq_src, cell, thread, dS, eqn)
{

  /*Thread *mix_th, *sec_th;*/

Thread *mix_thread, *thread_gas, *thread_liq;
  real m_dot_l;

  real T_SAT;             /*定义临界温度*/
  real tem_pressure;             /*定义单元压力*/

  real vof_gas,vof_liq,roul_gas,roul_liq,t_gas,t_liq,t_mix;

   mix_thread= THREAD_SUPER_THREAD(thread);
thread_gas= THREAD_SUB_THREAD(mix_thread, 1);
thread_liq= THREAD_SUB_THREAD(mix_thread, 0);
vof_gas=C_VOF(cell, thread_gas);
vof_liq=C_VOF(cell, thread_liq);
roul_gas=C_R(cell, thread_gas);
roul_liq=C_R(cell, thread_liq);
t_gas=C_T(cell, thread_gas);
t_liq=C_T(cell, thread_liq);
t_mix=C_T(cell, mix_thread);

  /*mix_th = THREAD_SUPER_THREAD(pri_th);
  sec_th = THREAD_SUB_THREAD(mix_th, 1);*/

  tem_pressure=(C_P(cell,mix_thread)+101325)/1000000;                         /*读取单元压力*/
  T_SAT=3826.36/(9.3876-log(tem_pressure))+45.47;    /*定义饱和温度与压力的关系*/

  if(t_liq>=T_SAT)
  {
m_dot_l = -0.1*vof_liq*roul_liq*fabs(t_liq - T_SAT)/T_SAT;
dS[eqn] = -0.1*roul_liq*fabs(t_liq - T_SAT)/T_SAT;
  }


  {
m_dot_l = 0.1*vof_gas*roul_gas*fabs(t_liq - T_SAT)/T_SAT;
dS[eqn] =0;
}

return m_dot_l;
  }
DEFINE_SOURCE(vap_src, cell, thread, dS, eqn)

{
  Thread *mix_thread, *thread_gas, *thread_liq;
  real m_dot_v;
  real T_SAT;             /*定义临界温度*/
  real tem_pressure;             /*定义单元压力*/
real vof_gas,vof_liq,roul_gas,roul_liq,t_gas,t_liq,t_mix;

mix_thread= THREAD_SUPER_THREAD(thread);
thread_gas= THREAD_SUB_THREAD(mix_thread, 1);
thread_liq= THREAD_SUB_THREAD(mix_thread, 0);
vof_gas=C_VOF(cell, thread_gas);
vof_liq=C_VOF(cell, thread_liq);
roul_gas=C_R(cell, thread_gas);
roul_liq=C_R(cell, thread_liq);
t_gas=C_T(cell, thread_gas);
t_liq=C_T(cell, thread_liq);
t_mix=C_T(cell, thread);

tem_pressure=(C_P(cell,mix_thread)+101325)/1000000;                         /*读取单元压力*/
  T_SAT=3826.36/(9.3876-log(tem_pressure))+45.47;    /*定义饱和温度与压力的关系*/

  if(t_gas>=T_SAT)
  {
m_dot_v = 0.1*vof_liq*roul_liq* fabs(t_gas - T_SAT)/T_SAT;
dS[eqn] = 0;  /*1*/
                           }
  else {
m_dot_v = -0.1*vof_liq*roul_liq*fabs(T_SAT-t_gas)/T_SAT;

dS[eqn] = -0.1*roul_liq*fabs(t_gas - T_SAT)/T_SAT;

   }
return m_dot_v;

}

DEFINE_SOURCE(enrg_src, cell, thread, dS, eqn)
{
Thread *mix_thread, *thread_gas, *thread_liq;
  real m_dot_l;
real vof_gas,vof_liq,roul_gas,roul_liq,t_gas,t_liq,t_mix;

  real T_SAT;             /*定义临界温度*/
  real tem_pressure;             /*定义单元压力*/

  mix_thread= THREAD_SUPER_THREAD(thread);
thread_gas= THREAD_SUB_THREAD(mix_thread, 1);
thread_liq= THREAD_SUB_THREAD(mix_thread, 0);
vof_gas=C_VOF(cell, thread_gas);
vof_liq=C_VOF(cell, thread_liq);
roul_gas=C_R(cell, thread_gas);
roul_liq=C_R(cell, thread_liq);
t_gas=C_R(cell, thread_gas);
t_liq=C_R(cell, thread_liq);
t_mix=C_T(cell, thread);

tem_pressure=(C_P(cell,mix_thread)+101325)/1000000;                         /*读取单元压力*/
  T_SAT=3826.36/(9.3876-log(tem_pressure))+45.47;    /*定义饱和温度与压力的关系*/

  if(t_mix>=T_SAT)
  {
m_dot_l = -0.1*vof_liq*roul_liq*
            fabs(t_mix - T_SAT)/T_SAT;
dS[eqn] =-0.1*vof_liq*roul_liq/T_SAT;
                           }
  else {
m_dot_l = 0.1*vof_gas*roul_gas*
            fabs(t_mix - T_SAT)/T_SAT;

dS[eqn] =0.1*vof_gas*roul_gas/T_SAT; /*1*/

   }

  return LAT_HT*m_dot_l;
}

回复此楼