CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > FLUENT

can you help? Its great for me.

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 25, 2006, 22:45
Default can you help? Its great for me.
  #1
Tajul
Guest
 
Posts: n/a
HI I am Tajul. I am a little bit despair about Fluent User Forum. I posted messages two times about my same problem but no one answered. This is the 3rd time I am sending the same problem. I do not think that no one does not know the answer. I am working with microchannel flow with slip wall boundary condition. The slip wall boundary condition is uw=Kn.(delu/deln)where uw is the fluid velocity at the wall, n is the normal coordinate pointing inward from the channel wall. Kn is the Knudsen number where kn=1/H.(viscosity*sqrt(pi)/sqrt(2*density*pressure)). H is the height of the channel. I wrote a UDF like this but giving floating point error and invalid number. I do not not how to set boundary conditions and hook this UDF. If any one help me I will be grateful to him and apreciate his kindness.

#include <udf.h> #include <mem.h>

#define PI 3.141592

DEFINE_PROFILE(wall_vel, t, i) {

face_t f;

cell_t c0;

Thread *t0;

double lamda, viscosity, density, Kn, R, Temperature, pressure, H, val, sigma_m, g_u, g_v, g_mag;

H = 0.0000012;

sigma_m = 1.0;

begin_f_loop(f,t)

{

c0=F_C0(f,t);

t0=t->t0;

viscosity = C_MU_L(c0,t0);

density = C_R(c0,t0);

R=C_RGAS(c0,t0);

pressure=C_P(c0,t0);

Temperature=C_T(c0,t0);

pressure = density*R*Temperature;

lamda = viscosity * sqrt(PI) / (sqrt(2.*density*pressure));

Kn = lamda / H;

g_u = C_U_G(c0,t0)[1];

g_v = C_U_G(c0,t0)[2];

g_mag=.5*(g_u+g_v);

val = ((2. - sigma_m)/sigma_m)*Kn*g_u;

F_PROFILE(f,t,i)=val;

}

end_f_loop(f,t) }

Thanking Tajul
  Reply With Quote

Old   January 26, 2006, 04:33
Default Re: can you help? Its great for me.
  #2
HVN
Guest
 
Posts: n/a
In your include you have forgetton the "udf.h"
  Reply With Quote

Old   January 26, 2006, 04:50
Default Re: can you help? Its great for me.
  #3
Tajul
Guest
 
Posts: n/a
Hi HVN, I have written that in my original UDF. It was left at the time of copying. If you have any other suggestions I'll apreciate. Thanks Tajul

#include <udf.h> #include <mem.h>

  Reply With Quote

Reply

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
Why is there a great difference between numeca and cfx???? sukiska NUMECA 0 October 12, 2011 10:21
Linux distros that work great with Fluent 6.3 on Intel PCs s195404 FLUENT 1 October 25, 2009 18:00
great guy--------------- mech FLUENT 0 May 18, 2007 09:12
fit a great number of point with Bezier spline JLUCHENG Main CFD Forum 3 April 20, 2003 20:22
great idea ! essemiani FLUENT 0 November 22, 1999 10:04


All times are GMT -4. The time now is 03:42.