[martin986]

hey!

I have different particle diameters in my simulation and want to use the stokes cunningham draglaw but therefore you can only set one value, so for only one diameter. I made an UDF for different correction factors (because of different diameters) as the following:

DEFINE_DPM_DRAG(particle_drag_force, Re, p)
{
double cu, kn, drag_force;
real g1 = 1.231;
real g2 = 0.4695;
real g3 = 1.1783;
real l = 68*pow(10,-9);

kn = 2*l/P_DIAM(p);
cu = 1+kn*(g1+g2*exp(-g3/kn));

drag_force = 18/cu+2.7*pow(Re,0.687)/cu;

return (drag_force);
}

It would be helpful if one can inform me if I can model that in this way and I am not sure in which range of Reynoldsnumber the cunningham correction is allowed, thanks