小木虫 --- 600万学术达人喜爱的学术科研平台

首页 >> 能源 >>查看话题

求大神们帮忙看下我的UDF

我在做水平管内含不凝气蒸汽的流动模拟,要用UDF加入质量、动量、能量源项。现在画的是二维网格,长20mm,直径18mm。网格数量1000000。但是结果得到的冷凝和潜热都是0。这是我的UDF,请大神看看有没有问题?
#include <udf.h>
#define LAT_HT 2357
#define Mv 18.01534
#define Ma 28.966
  int i=0;
  Thread *t;
  cell_t c;
  real V;
  real Vx;
  real Vy;
  real m_dot_l;
  real Tw;
  real Ps;
  real Pv;
  real Wv;
  real P;
  Thread *mix_th, *pri_th;
  V=C_VOLUME(c,t);
  Vx=C_U(c,t);
  Vy=C_V(c,t);
  
  mix_th = THREAD_SUPER_THREAD(t);
  pri_th = THREAD_SUB_THREAD(mix_th, 0);
  
  Tw=C_T(c,t);
  Wv=C_YI(c,t,i);
  P=C_P(c,t)+30000;
  Pv=Ma*Wv*P/(Ma*Wv+Mv*(1-Wv));
  Ps=2000*exp(18.5916-3991.11/(Tw-39.31))/15;
  if(Pv>Ps)
   {
     if(C_VOF(cell,pri_th)==1||C_VOF(cell,pri_th)==0)
       {
         m_dot_l=0;
        }
     else
       {
        m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
                fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;
        }
      }
  else
    {
     m_dot_l = 0;
      }
      C_UDMI(c,t,0)=m_dot_l*V;
      C_UDMI(c,t,1)=m_dot_l*LAT_HT*V;
      C_UDMI(c,t,2)=m_dot_l*Vx*V;
      C_UDMI(c,t,3)=m_dot_l*Vy*V;
   return m_dot_l;
}
DEFINE_SOURCE(vapour_species_source, c, t , dS, eqn)
{   
  real source=0;
  real V;
  V=C_VOLUME(c,t);
  source=C_UDMI(c,t,0)/V;
  dS=0.0;
  return source;
}
DEFINE_SOURCE(vapour_h_source, c, t , dS, eqn)
{   
  real source=0;
  real V;
  V=C_VOLUME(c,t);
  source=C_UDMI(c,t,1)/V;
  dS=0.0;
  return source;
}
DEFINE_SOURCE(lateral_mom_source, c, t , dS, eqn)
{   
  real source=0;
  real V;
  V=C_VOLUME(c,t);
  source=C_UDMI(c,t,2);
  dS=0.0;
  return source;
}
DEFINE_SOURCE(longi_mom_source, c, t , dS, eqn)
{   
  real source=0;
  real V;
  V=C_VOLUME(c,t);
  source=C_UDMI(c,t,3);
  dS=0.0;
  return source;
}

相关话题

用户评论

蒸发公式是什么

2楼: Originally posted by BeObjective at 2016-12-29 23:19:34
蒸发公式是什么
我的是水平管内含空气的水蒸汽冷凝,不是蒸发。冷凝的质量源项公式是    -0.1*蒸汽体积分数*密度* (T -T_SAT)/T_SAT;

3楼: Originally posted by 风去云随 at 2016-12-30 09:41:46
我的是水平管内含空气的水蒸汽冷凝,不是蒸发。冷凝的质量源项公式是    -0.1*蒸汽体积分数*密度* (T -T_SAT)/T_SAT;...
使用的哪个文献的公式

5楼: Originally posted by BeObjective at 2016-12-30 10:41:58
使用的哪个文献的公式
...
是Lee的相变系数模型求质量源项,在刘泉的   纯蒸气及含不凝气蒸气冷凝的数值研究  硕士论文中提到

5楼: Originally posted by BeObjective at 2016-12-30 10:41:58
使用的哪个文献的公式
...
是Lee的相变系数模型求质量源项,在刘泉的   纯蒸气及含不凝气蒸气冷凝的数值研究  硕士论文中提到