[czhao86]

Hi,

i wanna exert an external force on particles, and below is my draft. But it is not working, anyone have some idea? Thank you very much.

#include "udf.h"
 
DEFINE_DPM_BODY_FORCE(particle_body_force,p,i)
{
real bforce=0;

Domain *domain;
domain = Get_Domain(1);

Thread *t;
cell_t c;
real xc[ND_ND];

thread_loop_c(t,domain)
{
begin_c_loop_all(c, t)    
  {        
       C_CENTROID(xc,c,t);
       bforce = xc[0];
               
  }                         
end_c_loop_all(c, t)
}

return (bforce);
}

[czhao86]

Thread *t=P_CELL_THREAD(p);
cell_t c=P_CELL(p);
C_CENTROID(xc,c,t);

The code above seems to work.