24小时热门版块排行榜

>论坛更新日志 (839)
>虫友互识 (86)
>休闲灌水 (46)
>导师招生 (21)
>论文投稿 (18)
>公派出国 (15)
>考研 (15)
>考博 (14)
>硕博家园 (11)
>基金申请 (5)
>文献求助 (4)
>教师之家 (3)
>有奖起名 (2)
>能源 (1)
>育儿交流 (1)
>出国考试 (1)

返回列表

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;
}

回复此楼