udf

冷却塔模拟

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

冷却塔模拟

# include "udf.h"
const int zone_ID=3;
DEFINE_DPM_BODY_FORCE(bforce,p,i)
{      
        int n;
        real bforce,acc;
        real A[ND_ND];
        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的区域加载体积力。大神帮忙看看？