24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

hjnature

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 90.4
帖子: 43
在线: 40.1小时
虫号: 3610981
注册: 2014-12-24
性别: GG
专业: 多相流热物理学

[求助] 在源项中调用温度梯度已有2人参与

想在两相流源项的计算中使用温度梯度，但是为什么udf能顾编译通过，但是算了一步之后计算就报错呢，源程序如下，我是fluent新手，希望各位大神解答，
#include "udf.h"
#include "sg_mphase.h"
#define T_SAT 373
#define L 230000.0
#define d 1000
#define k_L 0.6791
DEFINE_SOURCE(liq_src, cell, thread, dS, eqn)
{
real void_s, tg,m_dot_l=0;
void_s=C_VOF(cell,thread);
if (NULL != THREAD_STORAGE(thread,SV_T_G))
{
tg=C_T_G(cell,thread)[0];
  if(C_T(cell,thread)>=T_SAT&&void_s>0&&void_s<1)
{
m_dot_l =-k_L*tg*d/L;
dS[eqn] =0;
                           }
  else if(C_T(cell,thread)<T_SAT&&void_s>0&&void_s<1)
{

m_dot_l =k_L*tg*d/L;

dS[eqn] =0.;

   }
else{
m_dot_l =0;
}
}
  return m_dot_l;
}

DEFINE_SOURCE(vap_src, cell, thread, dS, eqn)

{
real void_s, tg,m_dot_v=0;
void_s=C_VOF(cell, thread);
if (NULL != THREAD_STORAGE(thread,SV_T_G))
{
tg=C_T_G(cell,thread)[0];
  if(C_T(cell, thread)>=T_SAT&&void_s>0&&void_s<1)
{
m_dot_v =k_L*tg*d/L;
dS[eqn] = 0.;
                           }
  else if(C_T(cell, thread)<T_SAT&&void_s>0&&void_s<1)
{
m_dot_v = -k_L*tg*d/L;

dS[eqn] = 0;

   }
  else
{
  m_dot_v =0;
}
  }
return m_dot_v;
}

回复此楼