modifying radial distribution function
 

hey friends,
i have to write udf for modifying radial distribution function for eulerian granular simulation for fluidised bed reactor.i dont know haw to do ....please help me out.....

this is the udf which i hv written ......but there is error while iterating.

#include "udf.h"

DEFINE_PROPERTY(cell_radial,c,t)
{
real rad_dis;
real a_a;

a_a=1-pow(((1.-C_VOF(c,t))/0.63),(1/3));
rad_dis=(3/5)*pow(a_a,-1);

return rad_dis;
}

Hello, ravindra.
Were you able to solve this problem? I have the exact same problem and I can't write a right UDF for it.
Thanks in advance!