[sariug]

Hello everyone. For my wind analysis I am trying to implant an udf for inlet boundry condition. However its not working as i expected. I expected a log profile but its constant. Can anyone see the bug ? Thanks in advance.





#include "udf.h"
/* BOUNDARY PROFILES FOR WIND SPEED Uref = 10 m/s */
/* ROUGHESS LENGTH yo= 0.5 m */
real v = 10;
real uf = 1.38;
real yo = 0.5;
/* CALCULATION OF THE PROFILE FOR HORIZONTAL WIND SPEED */
DEFINE_PROFILE(vel, thread, nv) /* function name, thread and variable number */
{ face_t f;
real x[ND_ND];
real y;
begin_f_loop (f,thread)
{
F_CENTROID(x,f,thread);
y=x[1];
F_PROFILE(f,thread,nv) = (uf/(0.42))*log((y+11+yo)/yo);
}
end_f_loop(f,thread)
}
/* CALCULATION OF THE PROFILE FOR TURBULENT KINETIC ENERGY */
DEFINE_PROFILE(tke, thread, nv) /* function name, thread and variable number */
{
face_t f;
real x[ND_ND];
real y;
begin_f_loop (f,thread)
{
F_CENTROID(x,f,thread);
y=x[1];
F_PROFILE(f,thread,nv) = 3.33*(uf*uf);
}
end_f_loop(f,thread)
}
/* CALCULATION OF THE PROFILE FOR TURBULENCE DISSIPATION RATE */
DEFINE_PROFILE(eps, thread, nv) /* function name, thread and variable number */
{
face_t f;
real x[ND_ND];
real y;
begin_f_loop (f,thread)
{
F_CENTROID(x,f,thread);
y=x[1];
F_PROFILE(f,thread,nv) = (uf*uf*uf)/(0.42*(yo+x[2]+11));
}
end_f_loop(f,thread)
}

[vasava]

Can you replace

real v = 10;
real uf = 1.38;
real yo = 0.5;

With the following and try again

#define v  10
#define uf 1.38
#define yo 0.5

[sariug]

It did not work. Do you have any other suggestion ?