Dear all, I have a parabolic inlet velocity profile in a 3d laminar flow. The boundary face is on x_y plane ad the duct extend to z. The velocity profile function is z = - 0.0000001446271106·(x^2 + y^2 ) + 0.0000005640230104·(x + y) - 0.0000005498116698 I need an UDF specifing this velocity profile. Many Thank's

I try to use this udf and works fine. #include "udf.h" DEFINE_PROFILE(vel_profile, thread, position) { real x[ND_ND]; /* this will hold the position vector */ real y; real k;

face_t f; begin_f_loop(f, thread) { F_CENTROID(x,f,thread); y = x[1]; k= x[0]; F_PROFILE(f, thread, position) =- 0.0000001446271106*(pow(k*1000.,2.) +pow(y*1000.,2.) ) + 0.0000005640230104*(k + y)*1000. - 0.0000005498116698; } end_f_loop(f, thread) }