24小时热门版块排行榜    

查看: 518  |  回复: 0

annie19880

新虫 (初入文坛)

[求助] 动网格不会压缩

这动网格是用C编写的  但是网格不会压缩  请问哪边出了问题呢?

#include "udf.h"

DEFINE_GRID_MOTION(wind001,domain,dt,time,dtime)
{
  Thread *tf = DT_THREAD(dt);
  face_t f;
  Node *v;
  real NV_VEC(omega),NV_VEC(axis),NV_VEC(dx);
  real NV_VEC(origin),NV_VEC(rvec);
  real sign;
  int n;

  /*set deforming flag on adjacent cell zone */
  SET_DEFORMING_THREAD_FLAG(THREAD_TO(tf));

  sign = -5.0 * sin (26.178 * time);

  Message ("time = %f, omega = %f\n", time, sign);
  NV_S(omega, =, 0.0);
  NV_D(axis, =, 0.0,1.0,0.0);
  NV_D(origin, =, 0.0, 0.0, 0.512);

  begin_f_loop(f,tf)

  {
   f_node_loop(f,tf,n)
   {
    v = F_NODE(f,tf,n);

        /* update node if x position is greater than 0.02
           and that the current node has not been previously
           visited when looping through previous faces */
   if (NODE_X(v) > 0.02 && NODE_POS_NEED_UPDATE (v))
   {
           /* indicate that node position has been update
              so that it's not updated more than once */
           NODE_POS_UPDATED(v);

           omega[1] =  sign * pow (NODE_X(v)/0.230, 0.5);
           NV_VV(rvec, =,NODE_COORD(v), -,origin);
           NV_CROSS(dx, omega, rvec);
           NV_S(dx, *=, dtime);
           NV_V(NODE_COORD(v), +=, dx);
   }
   }
  }
  end_f_loop(f,tf);
}
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 annie19880 的主题更新
信息提示
请填处理意见