24小时热门版块排行榜

返回列表

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

bianboshen

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 15
帖子: 13
在线: 3.6小时
虫号: 2482995
注册: 2013-05-26
专业: 流体力学

[求助] UDF近壁面cell源项编写已有1人参与

各位大神，目前我想向所有贴近壁面第一层cell中添加能量源项，代码如下，但是加载后运行结果与不加载是一样的，大神能不能看看问题出现在了哪里？
DEFINE_SOURCE(liq_enrg_src, c, sec_th, dS, eqn)
{

  Thread *mix_th, *pri_th, *t;
  Domain *domain=Get_Domain(1);
  int zone_ID=18;/*wall ID*/
  real m_dot_d=0;
  real kd;
  real Cnt;
  real a;
  real Area[ND_ND];
  real b;
  real d;
  real vg;
  real Reg;
  real tr;
  real Enrg_Source=0;
  face_t f;
  mix_th = THREAD_SUPER_THREAD(sec_th);
  pri_th = THREAD_SUB_THREAD(mix_th,0);
  t=Lookup_Thread(domain,zone_ID);
  begin_f_loop(f, t)
{
if(c == F_C0(f,t))
{
  Cnt=(1.0-C_VOF(c,pri_th))*C_R(c,sec_th);
  a=Cnt/C_R(c,pri_th);
  if(a<0.3)
  {
   kd=0.18*sqrt(0.026/(C_R(c,pri_th)*0.008));
  }
  else
  {
   kd=pow(a,-0.85)*sqrt(0.026/(C_R(c,pri_th)*0.008));
  }
  m_dot_d=kd*Cnt;
  vg=sqrt(pow(C_U(c,pri_th),2)+pow(C_V(c,pri_th),2)+pow(C_W(c,pri_th),2));

  Reg=C_R(c,pri_th)*vg*0.008/C_MU_L(c,pri_th);

  d=0.00796*0.026/C_R(c,pri_th)/pow(C_VOF(c,pri_th)*vg,2)*pow(Reg,2/3)*pow(C_R(c,pri_th)/C_R(c,sec_th),-1/3)*pow(C_MU_L(c,pri_th)/C_MU_L(c,sec_th),2/3);

  tr=3.14*sqrt(C_R(c,sec_th)*d*d*d/16/0.026);

b=pow(18*pow(C_K_L(c,pri_th)*tr,3)*C_R(c,pri_th)*LAT_HT*pow(m_dot_d,5)/pow(C_R(c,sec_th),4)/C_MU_EFF(c,pri_th)/(1.0-C_VOF(c,pri_th))/(F_T(f,t)-T_SAT),0.25);
F_AREA(Area,f,t);
Enrg_Source=(F_T(f,t)-T_SAT)*b/C_VOLUME(c,mix_th)/(1-C_VOF(c, pri_th))*NV_MAG(Area)*1000000;
dS[eqn]=0;
//Enrg_Source=2;
Message("The value of Enrg_Source is: %f\n", Enrg_Source);
}
else
{
Enrg_Source=0;
dS[eqn]=0;
}
return Enrg_Source;
}
end_f_loop(f, t)
  }

回复此楼