[cfduser5226]

I am trying write a UDF to be used as an inlet BC. I have read through the tutorial in Fluent and it doesn't provide much help for a 3-d application.

My inlet velocity is composed of velocity magnitudes (x,y,z velocities). The velocity only varies in Y-direction and I have derived an equation for each velocity(u=f(y),v=f(y),w=f(y)). I am just curious how to write the UDF that will be able to define velocity the entire span of the inlet.

Thanks

[diamondx]

I am no expert in UDF, I just wrote my UDF last week; the UDF manual and the tutorials will help you understand the functions, but if you wanna learn , look for existing pieces of codes, understand them and modify them so it can suit your need... 2-D should not be that different from 3-D you just add the z component... I guess

[macfly]

Hi cfduser,

Here's your udf, I just copied 3 times an example from the udf manual version 15.0: 2.3.19.3. Example 1 - Pressure Profile

It doesn't matter 2D or 3D in this case, the DEFINE_PROFILE loops over all faces of the boundary where you apply it.

 #include "udf.h"
 DEFINE_PROFILE(vx,t,i)
 {
    real x[ND_ND];    /* this will hold the position vector */
    real y;
    face_t f;
    begin_f_loop(f,t)
      {
       F_CENTROID(x,f,t);
       y = x[1];
       F_PROFILE(f,t,i) = 1.1e5 - y*y/(.0745*.0745)*0.1e5;
      }
    end_f_loop(f,t)
 } 
 
  DEFINE_PROFILE(vy,t,i)
 {
    real x[ND_ND];    /* this will hold the position vector */
    real y;
    face_t f;
    begin_f_loop(f,t)
      {
       F_CENTROID(x,f,t);
       y = x[1];
       F_PROFILE(f,t,i) = 1.1e5 - y*y/(.0745*.0745)*0.1e5;
      }
    end_f_loop(f,t)
 } 
 
  DEFINE_PROFILE(vz,t,i)
 {
    real x[ND_ND];    /* this will hold the position vector */
    real y;
    face_t f;
    begin_f_loop(f,t)
      {
       F_CENTROID(x,f,t);
       y = x[1];
       F_PROFILE(f,t,i) = 1.1e5 - y*y/(.0745*.0745)*0.1e5;
      }
    end_f_loop(f,t)
 }