当前位置: 首页 > 能源 >求大神们帮忙看下我的UDF

求大神们帮忙看下我的UDF

作者 风去云随
来源: 小木虫 350 7 举报帖子
+关注

我在做水平管内含不凝气蒸汽的流动模拟,要用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[eqn]=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[eqn]=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[eqn]=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[eqn]=0.0;
  return source;
} 返回小木虫查看更多

今日热帖
  • 精华评论
  • BeObjective

    蒸发公式是什么

  • 风去云随

    引用回帖:
    2楼: Originally posted by BeObjective at 2016-12-29 23:19:34
    蒸发公式是什么

    我的是水平管内含空气的水蒸汽冷凝,不是蒸发。冷凝的质量源项公式是    -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的相变系数模型求质量源项,在刘泉的   纯蒸气及含不凝气蒸气冷凝的数值研究  硕士论文中提到

  • 果子酒

    同志你做的怎么样了哇?

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓