[haghshenasfard]

Dear All,
This is a udf I've used to compute the particle deposition rate case. But it doesn't work , I have a problem in build files of C++ in Fluent which is:

e:\fluent.inc\fluent6.3.26\src\machine.h(114) : warning C4005: 'stdout' : macro redefinition
C:\Program Files (x86)\Microsoft Visual Studio 9.0\VC\INCLUDE\stdio.h(161) : see previous definition of 'stdout'
e:\fluent.inc\fluent6.3.26\src\machine.h(115) : warning C4005: 'stderr' : macro redefinition
C:\Program Files (x86)\Microsoft Visual Studio 9.0\VC\INCLUDE\stdio.h(162) : see previous definition of 'stderr'
e:\fluent.inc\fluent6.3.26\src\machine.h(116) : warning C4005: 'stdin' : macro redefinition
C:\Program Files (x86)\Microsoft Visual Studio 9.0\VC\INCLUDE\stdio.h(160) : see previous definition of 'stdin'
..\..\src\aa.c(35) : error C2275: 'Thread' : illegal use of this type as an expression
e:\fluent.inc\fluent6.3.26\src\mem.h(653) : see declaration of 'Thread'
..\..\src\aa.c(35) : error C2065: 'tcelll' : undeclared identifier
..\..\src\aa.c(36) : error C2275: 'cell_t' : illegal use of this type as an expression
e:\fluent.inc\fluent6.3.26\src\mem.h(173) : see declaration of 'cell_t'
..\..\src\aa.c(36) : error C2146: syntax error : missing ';' before identifier 'c'
..\..\src\aa.c(36) : error C2065: 'c' : undeclared identifier
..\..\src\aa.c(39) : error C2275: 'cell_t' : illegal use of this type as an expression
e:\fluent.inc\fluent6.3.26\src\mem.h(173) : see declaration of 'cell_t'
..\..\src\aa.c(39) : error C2146: syntax error : missing ';' before identifier 'c'
..\..\src\aa.c(39) : error C2065: 'c' : undeclared identifier
..\..\src\aa.c(40) : error C2275: 'Domain' : illegal use of this type as an expression
e:\fluent.inc\fluent6.3.26\src\mem.h(555) : see declaration of 'Domain'
..\..\src\aa.c(40) : error C2065: 'd' : undeclared identifier
..\..\src\aa.c(41) : error C2275: 'real' : illegal use of this type as an expression
e:\fluent.inc\fluent6.3.26\src\global.h(158) : see declaration of 'real'
..\..\src\aa.c(41) : error C2146: syntax error : missing ';' before identifier 'normal'
..\..\src\aa.c(41) : error C2065: 'normal' : undeclared identifier
..\..\src\aa.c(41) : error C2109: subscript requires array or pointer type
..\..\src\aa.c(42) : error C2143: syntax error : missing ';' before 'type'
..\..\src\aa.c(43) : error C2275: 'real' : illegal use of this type as an expression
e:\fluent.inc\fluent6.3.26\src\global.h(158) : see declaration of 'real'
..\..\src\aa.c(43) : error C2146: syntax error : missing ';' before identifier 'Wa'
..\..\src\aa.c(43) : error C2065: 'Wa' : undeclared identifier
..\..\src\aa.c(44) : error C2065: 'd' : undeclared identifier
..\..\src\aa.c(44) : warning C4047: '=' : 'int' differs in levels of indirection from 'Domain *'
..\..\src\aa.c(45) : error C2275: 'real' : illegal use of this type as an expression
e:\fluent.inc\fluent6.3.26\src\global.h(158) : see declaration of 'real'
..\..\src\aa.c(45) : error C2146: syntax error : missing ';' before identifier 'x'
..\..\src\aa.c(45) : error C2065: 'x' : undeclared identifier
..\..\src\aa.c(45) : error C2109: subscript requires array or pointer type
..\..\src\aa.c(46) : error C2065: 'i' : undeclared identifier
..\..\src\aa.c(46) : error C2065: 'i' : undeclared identifier
..\..\src\aa.c(46) : error C2065: 'idim' : undeclared identifier
..\..\src\aa.c(46) : error C2065: 'i' : undeclared identifier
..\..\src\aa.c(47) : error C2065: 'normal' : undeclared identifier
..\..\src\aa.c(47) : error C2065: 'i' : undeclared identifier
..\..\src\aa.c(47) : error C2109: subscript requires array or pointer type
..\..\src\aa.c(47) : error C2065: 'i' : undeclared identifier
..\..\src\aa.c(48) : error C2065: 'tcell' : undeclared identifier
..\..\src\aa.c(48) : error C2223: left of '->storage' must point to struct/union
..\..\src\aa.c(48) : error C2065: 'c' : undeclared identifier
..\..\src\aa.c(52) : error C2065: 'i' : undeclared identifier
..\..\src\aa.c(52) : error C2065: 'i' : undeclared identifier
..\..\src\aa.c(52) : error C2065: 'idim' : undeclared identifier
..\..\src\aa.c(52) : error C2065: 'i' : undeclared identifier
..\..\src\aa.c(54) : error C2065: 'i' : undeclared identifier
..\..\src\aa.c(54) : error C2065: 'normal' : undeclared identifier
..\..\src\aa.c(54) : error C2065: 'i' : undeclared identifier
..\..\src\aa.c(54) : error C2109: subscript requires array or pointer type
..\..\src\aa.c(55) : error C2065: 'i' : undeclared identifier
..\..\src\aa.c(59) : error C2065: 'tcell' : undeclared identifier
..\..\src\aa.c(59) : error C2223: left of '->storage' must point to struct/union
..\..\src\aa.c(59) : error C2065: 'c' : undeclared identifier
..\..\src\aa.c(77) : error C2065: 'tcell' : undeclared identifier
..\..\src\aa.c(77) : error C2223: left of '->storage' must point to struct/union
..\..\src\aa.c(77) : error C2065: 'c' : undeclared identifier
..\..\src\aa.c(78) : error C2065: 'tcell' : undeclared identifier
..\..\src\aa.c(78) : error C2223: left of '->storage' must point to struct/union
..\..\src\aa.c(78) : error C2065: 'tcell' : undeclared identifier
..\..\src\aa.c(78) : error C2223: left of '->head' must point to struct/union
..\..\src\aa.c(78) : error C2065: 'tcell' : undeclared identifier
..\..\src\aa.c(78) : error C2223: left of '->storage' must point to struct/union
..\..\src\aa.c(78) : error C2065: 'c' : undeclared identifier
..\..\src\aa.c(78) : fatal error C1903: unable to recover from previous error(s); stopping compilation



thanks in advance for any response......

the UDF is:

#include "udf.h"
#include "dpm.h"
#include "mem.h"
#include "sg.h"
#define nu_s 0.27 /*Poisson ratio of the surface*/
#define nu_p 0.27 /*Poisson ratio of the particle*/
#define VISC 1.7894e-05 /*Viscosity*/
/*#define RGAS (UNIVERSAL_GAS_CONSTANT/MW)*/
#define Tdatum 288.15
#define NUM_UDM 6 /* No. of user-defined memory locations8/
real ParticleTotalMass;
real P_Mass[6];
real P_Impact_Mass[6];
real P_Stick_Mass[6];
Domain *domain; *Get domain pointer and assign later to domain*/
/* Boundary condition macro for the deposition model*/
DEFINE_DPM_BC(udfdeposition8,p,t,f,f_normal,dim)
{
#if !RP_HOST /* Used only for 3D cases*/
real crit_vel,alpha,Tavg,Ep,vcr=0.,Es,k1,k2,calc,E,utau 1,dudz,wall_shear_stress;
real Cu,cbar,lamda,kn,val,kc,ucws,utc,h,del,ufr,ff,utau ,area,wall_shear_force,utaunew;
real vn=0.;
real utau2, utc1,vpabs=0.,MassImpact;
real nor_coeff=0.8;
real tan_coeff=0.4;
real R=287.;
real tem_Mass=0.;
real tem_Particle_Dia=0.;
real A[ND_ND],ds,es[ND_ND],A_by_es,dr0[ND_ND],ivu=0.,jvv=0.,kvw=0.,du=0.;
real tauwall1=0.,tauwall2=0.,tauwall3=0.,wallfricv1=0., wallfricv2=0.,wallfricv3=0.;
real yplus,wallfricv4=0.,wallfricv5=0.,uvel,vvel,wvel,d umag=0.;
FILE *fp;*
Thread *tcelll=P_CELL_THREAD(p);
cell_t c = P_CELL(p); /* Exact thread location of particle*/
/*cell_t c=RP_CELL(&p->cCell); */
/*Thread *tcell=RP_THREAD(&p->cCell);*/
cell_t c=P_CELL(p); /* Exact cell location of particle*/
Domain *d;
real normal[3];
int i,idim=dim;
real Wa = 0.039;
d=Get_Domain(1);
real x[ND_ND];
for (i=0.; i<idim; i++)
normal[i] = f_normal[i];
C_UDMI(c,tcell,0) += 1.0;
if(p->type==DPM_TYPE_INERT) /* Checks if inert particle is used*/
{
/*computing normal velocity*/
for(i=0.;i<idim;i++)
{
vn += p->state.V[i]*normal[i];
vpabs += pow(p->state.V[i],2.);
}
vpabs = pow(vpabs,.5);
/*computing critical velocity*/
Tavg = (F_T(f,t)+C_T(c,tcell))/2.;
/*Message("Avg temp is %g\n",Tavg);*/
/*Message("particle temp is %g\n",P_T(p));*/
/*Message("face temp from tcell is %g\n",F_T(f,tcell));*/
/*Message("normal face temp is %g\n",F_T(f,t));*/
/*NEW CORRELATION*/
Ep=(3.*(pow(10.,20.)))*exp(-0.02365*Tavg);
Es=Ep;
/*Message("Young's mod is %g\n",Ep);*/
k1 = (1.-(nu_s*nu_s))/(3.14*Es);
k2 = (1.-(nu_p*nu_p))/(3.14*Ep);
/*Message("particle density is %g\n",P_RHO(p));*/
/*Message("Ai particle density is %g\n",p->init_state.rho);*/
calc = (5.*3.14*3.14*(k1+k2))/(4.*(pow(P_RHO(p),1.5)));
E = 0.51*(pow(calc,(2./5.))); /* El-Batsh parameter*/
vcr = pow(((2.*E)/P_DIAM(p)),(10./7.));
/*Message("capture velocity is %g\n",vcr);*/
MassImpact=P_FLOW_RATE(p)*pow(10,9);
C_UDMI(c,tcell,1)+=MassImpact;
cbar = sqrt(((8.*C_RGAS(c,tcell)*C_T(c,tcell))/3.14));
lamda = (2.*C_MU_EFF(c,tcell))/(C_R(c,tcell)*cbar);
kn = (2.*lamda)/P_DIAM(p); /* Knudsen number*/
Cu = 1.+kn*(1.2+(0.41*exp(-0.88/kn))); /*Cunningham Correction Factor*/
Tavg = (F_T(f,t)+C_T(c,tcell))/2.; /*Avg of particle & surface temperature*/
/*Calculating El-Batsh parameter*/
val = ((1.-(pow(nu_s,2.)))/Es)+((1.-(pow(nu_p,2.)))/Ep);
kc = (4./3.)/val;
/*Calculating critical wall shear velocity*/
ucws = ((Cu*Wa)/(C_R(c,tcell)*P_DIAM(p)))*(pow((Wa/(P_DIAM(p)*kc)),(1./3.)));
utc = sqrt(ucws);
/* Message("wall shear vel is %g\n", utc);*/
utc1=pow(ucws,0.5);
/*Calculating wall friction velocity*/
/* Alternate Method 1*/
/*dudz = C_DUDZ(c,tcell);*/
/*utau1 = sqrt((VISC*dudz)/C_R(c,tcell));*/
/*Message("Friction vel by regular formula is %g\n",utau1);*/
yplus=F_STORAGE_R(f,t,SV_WALL_YPLUS_UTAU);
/*Alternate Method 2*/
ivu=C_U(c,tcell)*normal[0];
jvv=C_V(c,tcell)*normal[1];
kvw=C_W(c,tcell)*normal[2];
du =ivu+jvv;
Thread *t0=THREAD_T0(t);
cell_t c0=F_C0(f,t);
BOUNDARY_FACE_GEOMETRY(f,t,A,ds,es,A_by_es,dr0);
/* Message("wall distance is %g\n",ds);*/
tauwall1=C_MU_EFF(c,tcell)*(du/ds);
wallfricv1=sqrt(tauwall1/C_R(c,tcell));
/* Message("wallfrictionvelocity 1 is %g\n",wallfricv1);*/
if (yplus<11.25)
{
tauwall2=C_MU_EFF(c,tcell)*C_STRAIN_RATE_MAG(c,tce ll);
wallfricv2=sqrt(tauwall2/C_R(c,tcell));
}
else
{
wallfricv2 = (1./0.41)*log(yplus*9.);
}
/* Message("wallfrictionvelocity 2 is %g\n",wallfricv2);*/
/* Alternate Method 3*/
tauwall3=C_MU_EFF(c,tcell)*sqrt(C_DUDX(c,tcell)*(C _DUDX(c,tcell)+C_DUDX(c,tcell))+C_DUDY(c,tcell)*(C _DUDY(c,tcell)+C_DVDX(c,tcell))+C_DUDZ(c,tcell)*(C _D

UDZ(c,tcell)+C_DWDX(c,tcell))+C_DVDX(c,tcell)*(C_D VDX(c,tcell)+C_DUDY(c,tcell))+C_DVDY(c,tcell)*(C_D VDY(c,tcell)+C_DVDY(c,tcell))+C_DVDZ(c,tcell)*(C_D VDZ(c,tcell)+C_DWDY(c,tcell))+C_DWDX(c,tcell)*(C_D WDX(c,tcell)+C_DUDZ(c,tcell))+C_DWDY(c,tcell)*(C_D WDY(c,tcell)+C_DVDZ(c,tcell))+C_DWDZ(c,tcell)*(C_D WDZ(c,tcell)+C_DWDZ(c,tcell)));
wallfricv3=sqrt(tauwall3/C_R(c,tcell));
/*ACTUAL WALL FRICTION VELOCITY CALCULATION*/
/* Message("wall yplus is %g\n",yplus);*/
ds=C_WALL_DIST(c,tcell);
wallfricv4=(VISC*yplus)/(ds*C_R(c,tcell));
/* Alternate Method 5*/
uvel=C_U(c0,t0);
vvel=C_V(c0,t0);
wvel=C_W(c0,t0);
dumag=sqrt((uvel*uvel)+(vvel*vvel));
if(yplus>10)
wallfricv5=sqrt((C_MU_EFF(c,tcell)*dumag)/(C_R(c,tcell)*ds));
else
wallfricv5=sqrt((VISC*dumag)/(C_R(c,tcell)*ds));
/* Message("wallfrictionvelocity 5 is %g\n",wallfricv5);*/
fp=fopen("Impact2.txt","a");
fprintf(fp,"%6.3f %6.2f %6.2f %f %f %f %f\n", P_DIAM(p),vn,vcr,Ep,MassImpact,ucws,utaunew);
tem_Particle_Dia=P_DIAM(p)*pow(10,6);
if(tem_Particle_Dia<1)
P_Impact_Mass[0]+=P_FLOW_RATE(p)*pow(10,9);
else if(tem_Particle_Dia>1 && tem_Particle_Dia<3)
P_Impact_Mass[1]+=P_FLOW_RATE(p)*pow(10,9);
else if(tem_Particle_Dia>3 && tem_Particle_Dia<5)
P_Impact_Mass[2]+=P_FLOW_RATE(p)*pow(10,9);
else if(tem_Particle_Dia>5 && tem_Particle_Dia<7)
P_Impact_Mass[3]+=P_FLOW_RATE(p)*pow(10,9);
else if(tem_Particle_Dia>7 && tem_Particle_Dia<10)
P_Impact_Mass[4]+=P_FLOW_RATE(p)*pow(10,9);
else
P_Impact_Mass[5]+=P_FLOW_RATE(p)*pow(10,9);
/* PARTICLE DOES NOT STICK*/
if(abs(vn)>vcr)
{
C_UDMI(c,tcell,2) += 1.0;
C_UDMI(c,tcell,3) += P_FLOW_RATE(p)*pow(10,9);
for(i=0;i<idim;i++)
p->state.V[i]-=vn*normal[i];
/*Apply tangential coefficient of restitution. */
for (i=0; i<idim; i++)
p->state.V[i]*=tan_coeff;
/*Add reflected normal velocity. */
for (i=0; i<idim; i++)
p->state.V[i]-=nor_coeff*vn*normal[i];
/* Store new velocity in state0 of particle*/
for (i=0; i<idim; i++)
p->state0.V[i]=p->state.V[i];
return PATH_ACTIVE;
}
/* PARTICLE DEPOSITS*/
else
if(wallfricv4 < utc)
{
/*num of particles deposited or num of hits*/
C_UDMI(c,tcell,4) += 1.0;
/* mass of particles deposited*/
C_UDMI(c,tcell,5) += P_FLOW_RATE(p)*pow(10,9);
tem_Mass=P_FLOW_RATE(p)*pow(10,9);
if(tem_Particle_Dia<1.)
P_Stick_Mass[0]+=P_FLOW_RATE(p)*pow(10,9);
else if(tem_Particle_Dia>1 && tem_Particle_Dia<3)
P_Stick_Mass[1]+=P_FLOW_RATE(p)*pow(10,9);
else if(tem_Particle_Dia>3 && tem_Particle_Dia<5)
P_Stick_Mass[2]+=P_FLOW_RATE(p)*pow(10,9);
else if(tem_Particle_Dia>5 && tem_Particle_Dia<7)
P_Stick_Mass[3]+=P_FLOW_RATE(p)*pow(10,9);
else if(tem_Particle_Dia>7 && tem_Particle_Dia<10)
P_Stick_Mass[4]+=P_FLOW_RATE(p)*pow(10,9);
else
P_Stick_Mass[5]+=P_FLOW_RATE(p)*pow(10,9);
fp=fopen("Stick2.txt","a");
fprintf(fp,"%f %f %6.2f %f %6.2f %6.2f %f\n", tem_Particle_Dia,tem_Mass,vn,vcr,Ep,ucws,utaunew);
fclose(fp);
}
}
return PATH_ABORT;
#endif /* Only for 3D cases*/
}
/*Setting UDM locations and initializing them to 0*/
DEFINE_ON_DEMAND(reset_UDMsnew)
{
int i=0;
Thread *t;
Domain *d;
cell_t c;
face_t f;
d=Get_Domain(1);
Message("Setting UDMs \n");
thread_loop_c(t,d)
{
begin_c_loop(c,t)
{
for(i=0;i<6;i++)
C_UDMI(c,t,i)=0.0;
}
end_c_loop(c,t)
}
}

[rmn_990]

Does anyone know the solution to this thread?

[pakk]

The solution to this thread is to ignore it.

haghshenasfard did not say if this is the first UDF that he/she tried or that other UDFs worked, haghshenasfard did not say if he/she tried to simplify the UDF.
haghshenasfard used Fluent version 6, which you should not use anymore.
haghshenasfard probably had his/her compiler set up in the wrong way. If haghshenasfard would have tried to compile one of the smaller examples from the Fluent help, his/her problem would have been much more clear.

[rmn_990]

Originally Posted by pakk

The solution to this thread is to ignore it.

haghshenasfard did not say if this is the first UDF that he/she tried or that other UDFs worked, haghshenasfard did not say if he/she tried to simplify the UDF.
haghshenasfard used Fluent version 6, which you should not use anymore.
haghshenasfard probably had his/her compiler set up in the wrong way. If haghshenasfard would have tried to compile one of the smaller examples from the Fluent help, his/her problem would have been much more clear.

Actually I was working with the same UDF and I received the same errors. So, I searched for the reason in this forum and I noticed this thread.

I will try to find the reason for this error and I will share the outcome.[/url]

[pakk]

So, do you only have problems with this specific UDF?
Have you been able to compile other UDFs?

[rmn_990]

Originally Posted by pakk

So, do you only have problems with this specific UDF?
Have you been able to compile other UDFs?

Not exactly this specific UDF. My problem is to find a way to write a UDF for considering stick boundary condition of particles when they hit the surface.
I think what fluent consider for each interaction between wall and particle is based on the centre of the particle and not the surface of the particle.
In other words, I want particles to attach to the wall as soon as the surface of the spherical particle hits the wall instead of the centre of particles.

[pakk]

Originally Posted by rmn_990

Not exactly this specific UDF. My problem is to find a way to write a UDF for considering stick boundary condition of particles when they hit the surface.
I think what fluent consider for each interaction between wall and particle is based on the centre of the particle and not the surface of the particle.
In other words, I want particles to attach to the wall as soon as the surface of the spherical particle hits the wall instead of the centre of particles.

That is not an answer to my question. I'll try again.

On the system that you are using now, have you ever successfully compiled any UDF? It doesn't matter what the UDF was about.

[rmn_990]

Originally Posted by pakk

That is not an answer to my question. I'll try again.

On the system that you are using now, have you ever successfully compiled any UDF? It doesn't matter what the UDF was about.

Yes, sir. There was no problem for me to compile this code.

[pakk]

Then you have completely confused me.

This thread shows compilation problems, and you said that you had the same errors. But now you say you had no problems with compiling... What? Why are you here then?

[rmn_990]

Originally Posted by pakk

Then you have completely confused me.

This thread shows compilation problems, and you said that you had the same errors. But now you say you had no problems with compiling... What? Why are you here then?

Sorry for making this confusion for you but this is the UDF that I am using and I can compile it without any error:



It would be highly appreciated if you can help me with my problem of considering trap boundary condition when the outer surface of a particle hits the wall instead of the center of the particle.

[pakk]

Originally Posted by rmn_990

Actually I was working with the same UDF and I received the same errors.

And later we find out that you were not working with the same UDF, and didn't get any errors...

Anyway: DPM assumes point particles. If your particle size is not negligible compared to the length scales in your geometry, DPM is the wrong tool. 6DOF would be better.