[soleimanian]

Hello
I want to define a velocity profile with UDF in ANSYS Fluent.
velocity profile is according to law Logarithm.
and I should use of the if expression in UDF code.
but when I interpreted this UDF code I received parse error.
I do not know why and I need to help me.
Thanks a lot.

my Code :

#include "udf.h"

/*define constant parameter*/
#define zmin 5.0;
#define ustar 0.086;
#define z0 0.3;
#define uuniform 25.0;

DEFINE_PROFILE(inlet_x_velocity, thread, position)
{
real x[ND_ND];
real Uz, Z;
face_t f;

begin_f_loop(f, thread)
{
F_CENTROID(x, f, thread);

Z= x[2];

if(Z< zmin) Uz= zmin;
if(Z>zmin) Uz= (ustar / 0.4)*(double log (double (Z / z0))*uuniform;
F_PROFILE(f, thread, position)=Uz;
}
end_f_loop(f, thread)
}

[AlexanderZ]

you forgot } at the end

next time put your error message here, so people may understand what is the problem

best regards

[soleimanian]

Originally Posted by soleimanian

Hello
I want to define a velocity profile with UDF in ANSYS Fluent.
velocity profile is according to law Logarithm.
and I should use of the if expression in UDF code.
but when I interpreted this UDF code I received parse error.
I do not know why and I need to help me.
Thanks a lot.

my Code :

#include "udf.h"

/*define constant parameter*/
#define zmin 5.0;
#define ustar 0.086;
#define z0 0.3;
#define uuniform 25.0;

DEFINE_PROFILE(inlet_x_velocity, thread, position)
{
real x[ND_ND];
real Uz, Z;
face_t f;

begin_f_loop(f, thread)
{
F_CENTROID(x, f, thread);

Z= x[2];

if(Z< zmin) Uz= zmin;
if(Z>zmin) Uz= (ustar / 0.4)*(double log (double (Z / z0))*uuniform;
F_PROFILE(f, thread, position)=Uz;
}
end_f_loop(f, thread)
}

I'm facing to this error.
Best regards

[AlexanderZ]

#include "udf.h"

/*define constant parameter*/
#define zmin 5.0
#define ustar 0.086
#define z0 0.3 
#define uuniform 25.0

DEFINE_PROFILE(inlet_x_velocity, thread, position)
{
	real x[ND_ND]; /* this will hold the position vector */
	real Uz, Z;
	face_t f;

	begin_f_loop(f,thread)
	{
		F_CENTROID(x, f, thread);
		Z= x[2];
		if (Z <= zmin) Uz = zmin;
		if (Z > zmin) Uz= (ustar / 0.4)*((double)log((double) (Z / z0))*uuniform);
		F_PROFILE(f, thread, position) = Uz;
	}
	end_f_loop(f, thread)
}

best regards