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.