CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   FLUENT (https://www.cfd-online.com/Forums/fluent/)
-   -   plz need your help plzzzzzzz plzzzzzzzzzz plzzzzzzzz plzzzzzzz (https://www.cfd-online.com/Forums/fluent/88064-plz-need-your-help-plzzzzzzz-plzzzzzzzzzz-plzzzzzzzz-plzzzzzzz.html)

kingnizarios May 7, 2011 09:27

plz need your help plzzzzzzz plzzzzzzzzzz plzzzzzzzz plzzzzzzz
 
help me with this error plzzzzz
error{cpp -I"C:\Fluent.Inc\fluent6.3.26/src" -I"C:\Fluent.Inc\fluent6.3.26/cortex/src" -I"C:\Fluent.Inc\fluent6.3.26/client/src" -I"C:\Fluent.Inc\fluent6.3.26/multiport/src" -I. -DUDFCONFIG_H="<udfconfig.h>" "C:\Users\souhail\Desktop\mass.c"
Error: C:\Users\souhail\Desktop\mass.c: line 6: parse error.
}
source:
#include "udf.h"
#include "threads.h"
#include "metric.h"
#include "mem.h"
#include "sg mphase.h"
DEFINE MASS TRANSFER(watertoair,cell,t,from index,from species index,to index,to species index)
{
/* variable declarations*/
real m lg1;
real x[ND ND];
double Ro,Ri,rad;
real rN=0.001; /*radius of nozzle=1mm*/
real v0=155;/*INPUT velocity*/
real termRo,termRi;
double
dMdrRo,dMdrRi,dMdRo1,dMdRo2,dMdRo,dMdRi1,dMdRi2,dM dRi,dRodx,dRidx;
real const1=v0*rN*rN*5.62*998.2;/*alpha × rho × vel × rN × rN × 5.62 */
m lg1=0; /*mass transfer term as reconginzed by FLUENT */
C CENTROID(x,cell,t);
Ro=0.05*x[0]+ rN;/*Ro is the width in the outer region.. Eq 4*/
Ri=0.65*sqrt(rN*x[0])+rN;/*Ri is the width in the inner region.. Eq 3*/
rad=x[1]; /*radial coordinate*/
/* de?ning the mass ?ow rate derivatives in the two regions */
dMdrRo=4.5*const1*pow((1-
pow(rad/Ro,1.5)),2)*pow(rad/Ro,0.5)*(1/(Ro*Ro*Ro));
dMdrRi=4.5*const1*pow((1-
pow(rad/Ri,1.5)),2)*pow(rad/Ri,0.5)*(1/(Ri*Ri*Ri));
dMdRo1=4.5*const1*rad*pow((1-
pow(rad/Ro,1.5)),2)*pow(rad/Ro,0.5)*(1/(Ro*Ro*Ro*Ro));
dMdRo2=-2*const1*pow((1-pow(rad/Ro,1.5)),3)/(Ro*Ro*Ro);
dMdRo=dMdRo1+dMdRo2;
dMdRi1=4.5*const1*rad*pow((1-
pow(rad/Ri,1.5)),2)*pow(rad/Ri,0.5)*(1/(Ri*Ri*Ri*Ri));
dMdRi2=-2*const1*pow((1-pow(rad/Ri,1.5)),3)/(Ri*Ri*Ri);
dMdRi=dMdRi1+dMdRi2;
dRodx=0.05;/* derivative of Eq 4 wrt x*/
dRidx=0.325*rN/sqrt(rN*x[0]);/* derivative of Eq 3 wrt x*/
termRo=dMdrRo+dMdRo*dRodx;
termRi=dMdrRi+dMdRi*dRidx;
if(rad<=Ri) /*if within the inner region of the jet.. Eq 2 */
{m lg1=termRi/4000;} /*..the number 1/4000 comes from the cell dim */
if((rad<=Ro)&&(rad>Ri)) /*if within the outer region of the jet.. Eq 2*/
{m lg1=termRo/4000;}
if(rad>Ro) /*neglect mass transfer in mist region*/
{m lg1=0;}
return m lg1;
}

laurentb May 10, 2011 09:31

try sg_mphase.h instead of sg mphase.h


All times are GMT -4. The time now is 20:34.