24小时热门版块排行榜    

查看: 524  |  回复: 1

chen1989

铁虫 (初入文坛)

[求助] 关于fluent+UDF的颗粒粘度编写老是出错!求各位大神帮忙!

层流流动、二维模型、欧拉混合模型,首相是水,次相是我自己定义的一种物质
代码如下:
#include "udf.h"
#include "sg_mphase.h"

DEFINE_PROPERTY(hydrate_viscosity_1,c,t)
{
  Thread *pri_th,*sec_th;
  real liquidviscosity=0.00179;
  real hydrateviscosity;
  real time=CURRENT_TIME;
  real vof;
  if(time<1.0)
  {
    return 0.002;
  }
  else
  {
          pri_th=THREAD_SUB_THREAD(t,0);
      sec_th=THREAD_SUB_THREAD(t,1);
          hydrateviscosity=C_T(c,sec_th)/100000;
          return hydrateviscosity;
  }                 
}

DEFINE_DELTAT(mydeltat,domain)
{
  real time_step;
  real flow_time=CURRENT_TIME;
  if(flow_time<10)
     time_step=0.4;
  else
     time_step=1;
  return time_step;
}

一旦计算在1s之后,也就是引入了C_VOF之后,fluent就会出现致命错误....
各位大神给帮帮忙,指点指点。
粘度方程内容并不是这样的,我是为了调试方便所以这么写,粘度和体积分数相关。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stitch928

新虫 (小有名气)

你颗粒粘度的计算公式是?
2楼2015-09-24 10:32:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chen1989 的主题更新
信息提示
请填处理意见