rezazarghani
Subscriber

#include "udf.h"


 


DEFINE_DPM_DRAG(particle_drag_force,Re,p)


{


  real w, drag_force;


 


  if (Re < 0.01)


    {


    drag_force=18.0;


    return (drag_force);


    }


  else if (Re < 20.0)


    {


    w = log10(Re);  


    drag_force = 18.0 + 2.367*pow(Re,0.82-0.05*w) ; 


    return (drag_force);


    }


  else


    /* Note: suggested valid range 20 < Re < 260 */


    {


    drag_force = 18.0 + 3.483*pow(Re,0.6305) ;


    return (drag_force);


    }


}