CFD Online Logo CFD Online URL
Home > Forums > General Forums > Main CFD Forum

Need urgent help Regarding Drag coefficient

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

LinkBack Thread Tools Display Modes
Old   September 16, 2006, 12:44
Default Need urgent help Regarding Drag coefficient
Posts: n/a
hi sir, can u help me to write UDF for drag coefficient. i have written UDF but its not working... i am attaching my UDF here. can u find where is the mistake in UDF. Its very urgent.

My drag coefficient function is:

Cd = max{ 24(1+0.15*pow(Re,0.687))/24, 8Em/(3(Em+4))}

Re = density(l)*Dia*|Vg-Vl|/viscosity;

Dia = Diameter of bubble

|Vg-Vl|---> Relative slip velocity

Em = g*density diff*pow(Dia,2)/surface tension;

#include #define diam2 0.005

DEFINE_EXCHANGE_PROPERTY(custom_drag,cell,mix_thre ad,s_col,f_col) { Thread *thread_g, *thread_l;

real x_vel_g, x_vel_l, y_vel_g, y_vel_l, z_vel_g, z_vel_l, U_vel, abs_v,

rho_g, rho_l, mu_l, reyp, k_g_s,sre,psv,sv,Em,f,vasu;

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

thread_g = THREAD_SUB_THREAD(mix_thread, f_col);/* gas phase */ thread_l = THREAD_SUB_THREAD(mix_thread, s_col);/* liquid phase*/

/* find phase velocities and properties*/

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

x_vel_l = C_U(cell, thread_l); y_vel_l = C_V(cell, thread_l); z_vel_l = C_W(cell, thread_l); U_vel = (x_vel_l + y_vel_l + z_vel_l)/3;

rho_g = 1.225; /* density of gas */ rho_l = 998; /* density of liquid */ mu_l = 0.0001; /* viscosity of liquid */

/*compute slip*/ abs_v = sqrt(x_vel_g*x_vel_g + y_vel_g*y_vel_g + z_vel_g*z_vel_g);

/*compute Reynold's number*/ f=(U_vel-abs_v); if (f>0)

f=f; else

f=-f; reyp = rho_l*f*diam2/mu_l;/*Reynold number*/

sre = pow(reyp,-0.313); psv = pow(reyp,-1); sv = 24*psv + 3.6*sre; Em = 133946.8*diam2*diam2; vasu = 8*Em/(3*Em + 12); if(sv < vasu)

k_g_s = vasu; else

k_g_s = sv;

return k_g_s; }

  Reply With Quote

Old   September 16, 2006, 14:15
Default Re: Need urgent help Regarding Drag coefficient
Posts: n/a
1- is the first line part of the function (missing semi colon)

2- what the third and fourth lines are? are they statements definitions or what

3- Sixth line, you cannot have and include and define predirectives on the same line.

The best you can do is to get a C language primer
  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
Drag Coefficient Convergence Problem John FLUENT 16 September 4, 2009 02:44
Drag coefficient for parcels in dieselFoam sebastian_vogl OpenFOAM Running, Solving & CFD 5 December 31, 2008 13:19
Automotive test case vinz OpenFOAM Running, Solving & CFD 98 October 27, 2008 09:43
Need urgent help Regarding Drag coefficient srinivas FLUENT 3 September 20, 2006 10:10
Need urgent help Regarding Drag coefficient srinivas FLUENT 0 September 16, 2006 12:25

All times are GMT -4. The time now is 16:21.