冷却塔模拟 发表于 2017-6-18 18:57:22

udf

请问如何控制离散相只在某特定区域加载体积力？谢谢

冷却塔模拟 发表于 2017-6-19 08:05:23

# include "udf.h"
const int zone_ID=3;
DEFINE_DPM_BODY_FORCE(bforce,p,i)
{      
        int n;
        real bforce,acc;
        real A;
        cell_t c;
        Thread *tc,*tf;
        Domain*d;
        d=Get_Domain(1);
        c=RP_CELL(&p->cCell);
        tc=RP_THREAD(&p->cCell);
        c_face_loop(c,tc,n)
        {
   tf=C_FACE_THREAD(c,tc,n);
       if(zone_ID==THREAD_ID(tf))
       {
      if(i==1)
                {
          bforce=1.001e-4;
                }
         
         else
          {
             bforce=0;
         }
                   }
      else
          {
             bforce=0;
         }
          acc=bforce/P_MASS(p);
return acc;
        }
}
这个是写的udf，编译出错。在编号为3的区域加载体积力。大神帮忙看看？
页: [1]
查看完整版本: udf