CFD Online Logo CFD Online URL
Home > Forums > Fluent UDF and Scheme Programming

Di Felice drag model FLUENT UDF

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   September 11, 2013, 01:51
Default Di Felice drag model FLUENT UDF
Musango Lungu
Join Date: Dec 2011
Location: China
Posts: 73
Rep Power: 7
Musa is on a distinguished road
Sponsored Links
I am modeling gas solid fluidized beds using different drag models. The Di Felice drag model I used compiles properly but it cannot predict fluidization in the bed. The bed does not fluidize. The code is given below. Any suggestions or advice?

#include "udf.h"

#define pi 4.*atan(1.)

DEFINE_EXCHANGE_PROPERTY(EMMS_Bubbling,cell,mix_th read,s_col,f_col)
Thread *thread_g, *thread_s;
real x_vel_g, x_vel_s, y_vel_g, y_vel_s, abs_v, slip_x, slip_y,
rho_g, rho_s, mu_g, reyp, Res,Cd,X,vof_g, vof_s,taup, k_g_s, diam2;

/* find the threads for the gas (primary) */
/* and solids (secondary phases) */

thread_g = THREAD_SUB_THREAD(mix_thread, s_col);/* gas phase */
thread_s = THREAD_SUB_THREAD(mix_thread, f_col);/* solid phase*/

/* find phase velocities and properties*/

x_vel_g = C_U(cell, thread_g);
y_vel_g = C_V(cell, thread_g);

x_vel_s = C_U(cell, thread_s);
y_vel_s = C_V(cell, thread_s);

slip_x = x_vel_g - x_vel_s;
slip_y = y_vel_g - y_vel_s;

rho_g = C_R(cell, thread_g);
rho_s = C_R(cell, thread_s);

mu_g = C_MU_L(cell, thread_g);

diam2 = C_PHASE_DIAMETER(cell, thread_s);

/*compute slip*/
abs_v = sqrt(slip_x*slip_x + slip_y*slip_y);

/*compute Reynold's number*/

reyp = rho_g*abs_v*diam2/mu_g;

/* compute particle relaxation time */

vof_g = C_VOF(cell, thread_g);/* gas vol frac*/
vof_s = C_VOF(cell, thread_s);/* particle vol frac*/

Res = vof_g * reyp;


X = 3.7-0.65*exp(-pow(1.5-log10(Res),2)/2);

/* Di Felice Model */
k_g_s = 3/4*Cd*vof_g*vof_s*rho_g*abs_v/diam2*pow(vof_g,-X);

return k_g_s;
Musa is offline   Reply With Quote
Sponsored Links

Old   November 25, 2013, 04:20
New Member
Prashant Gupta
Join Date: Mar 2011
Location: Edinburgh
Posts: 29
Rep Power: 8
Prash is on a distinguished road
Hi Musa,

Did you find the minimum fluidization velocity using De Felice ?
I suggest take up a case from literature, with a known minimum fluidization velocity. and compare with the solution you are getting.

Alternatively, you could simple test set ups like single particle sedimentation, and block of constant particles falling to check against analytical solutions.

I am also testing these, with EDEM-Fluent.

Kind regards
Prash is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
problem with UDF for (v2-f) turbulence model in fluent artemiss1984 Fluent UDF and Scheme Programming 6 January 17, 2014 06:50
new turbulence model in Fluent C.C FLUENT 0 August 29, 2013 10:44
fluent udf problem: write specific data for every iteration in a file. nnvoro Fluent UDF and Scheme Programming 1 May 27, 2013 15:26
drag coefficient in ansys fluent 12.0 krishna FLUENT 18 May 2, 2013 13:29
Advanced Turbulence Modeling in Fluent, Realizable k-epsilon Model Jonas Larsson FLUENT 5 March 13, 2000 04:27

Sponsored Links

All times are GMT -4. The time now is 05:01.