Di Felice Drag Model for Spouted/Fluidised Beds

Submitted by ragingseahorse on Wed, 07/18/2018 - 11:35

Hi everyone!

I have a question about implementing Di Felice as a UDF drag model (in TFM) in FLUENT.

Would I be giving FLUENT the expression for gas-solid coefficient, K_gs or do I give it the expression for calculating the drag force?

/* compute reynolds no */
re = rho_g*abs_v*diam*vof_g/mu_g;

/* cd and voidage function */
cd = pow(0.63+4.8*sqrt(1/res),2); /* dallavalle */
beta = 4.8-0.65*exp(-0.5*pow(1.5-log10(re),2)); /* which one? */
ge = pow(vof_g,-beta);

/* CHECK: fluent wants fd or kgs? */
/* fd = 0.125*cd*rho_g*abs_v*abs_v*pi*d*d*ge; */

Also for beta, how could I check whether to use 4.8 or 3.7-0.65*exp(......?

Sorry if they're stupid questions, I'm really new to this.

Thanks so much
L.E.O.