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) } |
I find the solution, I just restarted from the beginning,it was fun
|
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:) |
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. |
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. |