CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   Fluent UDF and Scheme Programming (https://www.cfd-online.com/Forums/fluent-udf/)
-   -   Probleme to interpret 3D temperature profil UDF (https://www.cfd-online.com/Forums/fluent-udf/96207-probleme-interpret-3d-temperature-profil-udf.html)

christophe06 January 16, 2012 09:23

Probleme to interpret 3D temperature profil UDF
 
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 January 16, 2012 10:58

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

Aerogirl February 19, 2012 15:08

Adapt UDF
 
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 February 20, 2012 04:49

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 February 20, 2012 08:06

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


All times are GMT -4. The time now is 15:15.