24小时热门版块排行榜    

查看: 820  |  回复: 0

fan5310xm

金虫 (正式写手)

[求助] fluent动网格模拟,udf代码可否指点一下。

本人在做动网格模拟,需要用到udf,代码遇到一些困难,能否帮我解释一下。

DEFINE_UDS_FLUX(my_flux, f, t, i)  
{
  cell_t  c0,  c1 = -1;
  Thread *t0, *t1 = NULL;

  real NV_VEC(psi_vec), NV_VEC(A), flux = 0.0;

  c0 = F_C0(f,t);
  t0 = F_C0_THREAD(f,t);
  F_AREA(A, f, t);

if (BOUNDARY_FACE_THREAD_P(t))
    {
      real dens;

      
      if (NNULLP(THREAD_STORAGE(t,SV_DENSITY)))
        dens = F_R(f,t);   
      else
        dens = C_R(c0,t0);
      NV_DS(psi_vec,  =, F_U(f,t), F_V(f,t), F_W(f,t), *, dens);

      flux = NV_DOT(psi_vec, A);
    }
  else
    {
      c1 = F_C1(f,t);      
      t1 = F_C1_THREAD(f,t);

      NV_DS(psi_vec,  =, C_U(c0,t0),C_V(c0,t0),C_W(c0,t0),*,C_R(c0,t0));
      NV_DS(psi_vec, +=, C_U(c1,t1),C_V(c1,t1),C_W(c1,t1),*,C_R(c1,t1));

      flux = NV_DOT(psi_vec, A)/2.0;
    }

  return flux;  
}
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 fan5310xm 的主题更新
信息提示
请填处理意见