24小时热门版块排行榜

返回列表

sltwk

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 185.8
红花: 1
帖子: 53
在线: 36.6小时
虫号: 2670058
注册: 2013-09-21
性别: GG
专业: 工程热物理相关交叉领域

[求助] fluentUDF传质设置时出错

udf代入FLUENT中解释时没有错误，但设置传质时出现错误？

D:\PROGRA~2\ANSYSI~1\v145\fluent\fluent14.5.0\win64\2d\fl1450s.exe received fatal signal (ACCESS_VIOLATION)
1. Note exact events leading to error.
2. Save case/data under new name.
3. Exit program and restart to continue.
4. Report error to your distributor.

#include"udf.h"
#define T_SAT  373  /*定义饱和温度*/
#define lam 100 /*定义松弛因子*/
#define LAT_H  2448 /*定义气液相变潜热*/
DEFINE_PROPERTY(density_tem,c,t)
{
      real density,k1,cp,v;
   density=0.34710*pow(0.27400,-(pow(1-C_T(c,t)/647.13,0.28571)));
      return density;                                                                /*物性参数随温度变化*/
}
DEFINE_PROPERTY(thermal_conductivity_tem,c,t)
{
  real thermal_conductivity;

  thermal_conductivity=-0.2758+0.0046120*C_T(c,t)-0.000005391*pow(C_T(c,t),2);
  return thermal_conductivity;
}
DEFINE_PROPERTY(cp_tem,c,t)
{
real cp;
  cp=(92.053-0.039953*C_T(c,t)-0.00021103*pow(C_T(c,t),2)-0.00000053469*pow(C_T(c,t),3))*(1000/18);
   return cp;

}
DEFINE_PROPERTY(viscosity_tem,c,t)
{real viscosity;
viscosity=pow(10,-10.2158+1792.5/C_T(c,t)+0.017730*C_T(c,t)-0.000012631*pow(C_T(c,t),2));
   return viscosity;

}
DEFINE_PROPERTY(VCp_tem,c,t)
{real VCp;
VCp=33.933-0.0084186*(C_T(c,t))+2.9906e-05*pow(C_T(c,t),2)-1.7825e-08*pow(C_T(c,t),4)+3.6934e-12*pow(C_T(c,t),4)*(1000/18);
return VCp;
}
DEFINE_PROPERTY(VViscosity_tem,c,t)
{ real VViscosity;
VViscosity=-36.826+4.2900e-01*(C_T(c,t))-1.6200e-05*pow(C_T(c,t),2);
return VViscosity;
}
DEFINE_PROPERTY(VThermal_conductivity_tem,c,t)
{real VThermal_conductivity;
  VThermal_conductivity=-0.2758+4.6120e-03*(C_T(c,t))-5.5391e-06*pow(C_T(c,t),2);
return VThermal_conductivity;
}

DEFINE_SOURCE(liq_src,cell,first,ds,eqn)
{
   Thread *mixer,*second;
   real m_dot_first;
   mixer=THREAD_SUPER_THREAD(first);
   second=THREAD_SUB_THREAD(mixer,1);
   if(C_T(cell,mixer)>=T_SAT)
   {
      m_dot_first=-lam*C_VOF(cell,first)*C_R(cell,first)*fabs(C_T(cell,first)-T_SAT)/T_SAT;
      ds[eqn]=-lam*C_VOF(cell,first)*fabs(C_T(cell,first)-T_SAT)/T_SAT;
   }
   else
   {
      m_dot_first=lam*C_R(cell,first)*fabs(T_SAT-C_T(cell,first))/T_SAT;
      ds[eqn]=0;
   }
      return m_dot_first;
}
DEFINE_SOURCE(vap_src,cell,second,ds,eqn)
{
   Thread *mixer,*first;
   real m_dot_second;
   mixer=THREAD_SUPER_THREAD(second);
   first=THREAD_SUB_THREAD(mixer,0);
   if(C_T(cell,mixer)<=T_SAT)
   {
      m_dot_second=lam*C_VOF(cell,first)*C_R(cell,first)*fabs(T_SAT-C_T(cell,mixer))/T_SAT;
      ds[eqn]=0;
   }
   else
   {
   m_dot_second=-lam*C_VOF(cell,first)*C_R(cell,first)*fabs(T_SAT-C_T(cell,mixer))/T_SAT;
   ds[eqn]=-lam*C_R(cell,first)*fabs(T_SAT-C_T(cell,mixer))/T_SAT;
   }
   return m_dot_second;

}
DEFINE_SOURCE(enrg_src,cell,mixer,ds,eqn)
{
Thread *first,*second;
real m_dot;
first=THREAD_SUB_THREAD(mixer,0);
second=THREAD_SUB_THREAD(mixer,1);
if(C_T(cell,mixer)>=T_SAT)
  {

      m_dot=-lam*C_VOF(cell,first)*C_R(cell,first)*fabs(C_T(cell,first)-T_SAT)/T_SAT;
      ds[eqn]=-lam*C_VOF(cell,first)*(C_T(cell,first)-T_SAT)/T_SAT;
  }
   else if(C_T(cell,mixer)<=T_SAT)
  {
   m_dot=lam*C_VOF(cell,first)*C_R(cell,first)*fabs(T_SAT-C_T(cell,mixer))/T_SAT;
   ds[eqn]=lam*C_VOF(cell,first)*(T_SAT-C_T(cell,mixer))/T_SAT;}
   return LAT_H*m_dot;

}

回复此楼

» 猜你喜欢

gogogo

1楼 2014-07-05 10:16:39

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 sltwk 的主题更新

返回列表