[LiLAJ]

dear all,
I'm working on 2d heated vertical channel,
Inlet boundary :
pressure inlet is local bernoulli equation P= -0.5*v^2
inlet velocity gradient equals zero.
the problem how to define the velocity
I used the following udf code but results were dissatisfied.

include "udf.h"
DEFINE_PROFILE(inlet_boundary, thread, position)
{
real x[ND_ND];
real y;
face_t f;
begin_f_loop(f, thread)
{
F_CENTROID(x, f, thread);
y = x[1];
if (y == 0.0)
{
cell_t c0 = F_C0(f, thread);
real v = /* calculate or obtain the velocity from your simulation */
F_U(f, thread) = 0.0;
F_V(f, thread) = v;
F_W(f, thread) = 0.0;
C_V_G(c0, thread)[1] = 0.0;
F_P(f, thread) = -0.5 * v * v;
}
}
end_f_loop(f, thread);
}