[christophe06]

Hi,

I have defined a UDF for 3D Inlet temperature profil and it doesn't works anymore after a crash of fluent. My profil has become flat and I couldn't have good interpretation of all my UDF know. It's really a bad point because I need to give a report about these CFD results this week. I have writing the temperature profil as followed:

#include "udf.h"
"define domain_ID 2 //

enum{OC, NUM_OF_USED_UDM};

DEFINE_PROFILE(inlet_temp, t, i)
{
real x[ND_ND];
real r;
face_t f;
real Diameter = 6.e-2;
real T_max = 5.e3;
real T_max = 3.e2;
begin_f_loop(f, t)
{
F_CENTROID(x,f,t);
y = x[1];
z = x[2];
r = sqrt(pow(y,2)+pow(z,2));
F_PROFILE(f, t, i) = (T_max-T_paroi)*(1-pow(2.0*r/Diameter,2))+T_paroi;
}
end_f_loop(f, thread)
}

[christophe06]

I find the solution, I just restarted from the beginning,it was fun

[Aerogirl]

Hi there,
I want to do the exact same UDF as you but my inlet is rectangular shaped and I was wondering if you could explain how to adapt your code to fit non-cylindrical coordinates? I am trying to model a hot spot temperature on the inlet of a gas turbine duct.

Your help would be much appreciated.
Kind Regards
AeroGirl

[christophe06]

Hi,

I think that you can define your profil along one direction, for example X and add other condition on Y.

#include "udf.h"
"define domain_ID 2 //

enum{OC, NUM_OF_USED_UDM};
DEFINE_PROFILE(inlet_x_temperature, thread, position)
{

real x[ND_ND];

real r;

face_t f;

real lenght = 6.e-2;

real T_max = 5.e3;

real T_max = 3.e2;

begin_f_loop(f, t)

{

F_CENTROID(x,f,t);

x = x[0];

F_PROFILE(f, t, i) = (T_max-T_paroi)*(1-pow(2.0*x/Diameter,2))+T_paroi;

}

end_f_loop(f, thread)

}
DEFINE_PROFILE(inlet_y_temperature, thread, position)
{

real x[ND_ND];

real r;

face_t f;

real lenght2 = 2.e-2;

real T_max = 5.e3;

real T_max = 3.e2;

begin_f_loop(f, t)

{

F_CENTROID(y,f,t);

y = x[1];

F_PROFILE(f, t, i) = (T_max-T_paroi)*(1-pow(2.0*y/Diameter,2))+T_paroi;

}

end_f_loop(f, thread)

}Maybe this script will be ok.

[Aerogirl]

Thankyou I will give this a try and see how I get on!