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

why boudary slip velocity is so larte using udf

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 4, 2005, 00:06
Default why boudary slip velocity is so larte using udf
  #1
cxzhao
Guest
 
Posts: n/a
hello:

i am compute a nozzle now, and using udf to considerate slip velocity on the wall, but the result show that the slip velocity on the wall is larger than other position near the wall, the udf code is as following.

i added this udf into fluent on the plant(boudary condition-wall-shear condition -specified shear ) , the shear stress computed by slip veloctiy and viscosity in udf. is it right?

i think maybe somthing wrong with shear stress ,can anybody tell me how to reslove this problem??

#include "udf.h" DEFINE_PROFILE(leftwallslipx,thread,index) {

real x[ND_ND];

face_t f;

cell_t c;

double kn;

double k;

double u;

double coef;

double rho;

double D;

double T;

double P;

double l;

double miu;

double dtdxw,dudx,dudy,dudw;

Thread *t0;

FILE *fp;

Thread *cell_thread;

real time=CURRENT_TIME;

fp=fopen("test.txt","w");

begin_f_loop(f,thread)

{

if(time<0.0001)

{ F_PROFILE(f,thread,index)=0;

}

else

{

c=F_C0(f,thread);

cell_thread=THREAD_T0(thread);

T=C_T(c,cell_thread);

P=C_P(c,cell_thread)+99876;

l=0.034*0.001*T/P; //m,atm

miu=C_MU_EFF(c,cell_thread);

rho=P/(287*T); //P=rho*R*T,pa,kg/m^3,K,R=287J/kg*K

kn=l/12*100000; //knudson number // fprintf(fp,"%f,%f,%15.12f,%f,%f,%f,%f\n",T,P,l,miu ,rho,kn);

if(kn>0.001)

{

dtdxw=fabs(C_T_G(c,cell_thread)[0]);

dudx=fabs(C_DUDY(c,cell_thread)+C_DVDY(c,cell_thre ad));

dudy=fabs(C_DUDX(c,cell_thread)+C_DVDX(c,cell_thre ad));

dudw=(dudx+dudy)*1.414;

k=miu*1000000/5;

u=l*dudw+0.75*miu/rho/T*dtdxw;//velocity // fprintf(fp,"%f,%f,%15.12f,%f,%f,%f,%f,%f,%f\n",T,P ,l,miu,rho,kn,dtdxw,dudx,u);

F_PROFILE(f,thread,index)=k*u*0.707; // shear stress

}

else

F_PROFILE(f,thread,index)=0;

}

}

end_f_loop(f,thread)

fclose(fp);

}

  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
How to define a parabolic velocity inlet of particles with UDF zumaqiong FLUENT 2 February 22, 2010 10:46
UDF for velocity farne FLUENT 1 February 25, 2009 01:29
UDF problem : inlet velocity in cyl. coord. system Jongdae Kim FLUENT 0 June 15, 2004 11:21
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 09:11
slip velocity on boundary (V.P.M) Lee, Juhee Main CFD Forum 3 October 22, 1998 15:30


All times are GMT -4. The time now is 22:04.