关于圆柱简谐运动

詹皇

最近在做关于管道振动的仿真，管道里流动的是液固两相流，udf程序如下：
#include "udf.h"

/**计算*/
DEFINE_GRID_MOTION(pipe,domain,dt,time,dtime)
{
    Thread *tf = DT_THREAD(dt);
    face_t f;
    Node *v;
   real scale;
    real sign;
        real diam=0;
        real x=0;
    real z=0;
    int n;
    /* set deforming flag on adjacent cell zone */
       

    SET_DEFORMING_THREAD_FLAG(THREAD_T0(tf));
     sign = 0.004 * sin (10 * time);
   
    begin_f_loop(f,tf)
    {
       f_node_loop(f,tf,n)
       {
                  
        v = F_NODE(f,tf,n);

                 if(NODE_Y(v)<=0.15)diam=0.03;
                 else  diam=-(NODE_Y(v)-0.2)*0.6;
                  
                scale=(diam+sign)/diam;/**直线振动变化*/

        
        if (NODE_POS_NEED_UPDATE (v))
        {x=NODE_X(v);
                 z=NODE_Z(v);
         
         NODE_POS_UPDATED(v);
                 NODE_X(v)=diam*scale*(x/pow((x*x+z*z),0.5));
                 NODE_Z(v)=diam*scale*(z/pow((x*x+z*z),0.5));
                 /*Y轴固定*/
           }
     }
  }
  end_f_loop(f,tf);
}
请问大神程序有没有问题，dynamic mesh里该如何设置？为什么我的仿真结果里固体相在出口端呈阶梯式减少，结果与不震动情况下相差太大

望山跑

你的圆柱是柔性圆柱么？

打俩鸡蛋

这个模型及代码我见过，是收缩运动，但我没算