UDF for 3D circular inlet

Wong Wai Loong

I had tried with the below UDF progamme but it cannot produce a parabolic profile at the 3D inlet. So i need some advise on my programme. thanx

#include "udf.h"

DEFINE_PROFILE(inlet_y_velocity, thread, position) { real y[ND_ND]; /* this will hold the position vector */ real x, z, r; face_t f; r=x=z;

begin_f_loop(f, thread)

{

F_CENTROID(y,f,thread);

x = y[1];

z = y[2];

r = y[3];

F_PROFILE(f, thread, position) = .07667 - r*r/(.0075*.0075)*.07667;

} end_f_loop(f, thread) }