[Sai Krishna]

i have tried one code to give temperature variation from 500 to 800 kelvin along z direction for a wall in XZ plane. interpretation and compilation were done without errors but during post processing the temperature results, i didn't observe any variation and its showing 500K all over the plane. presently iam using ansys 18.2 version.
code is as follows.

#include"udf.h"
DEFINE_PROFILE(temperature,thread,position)
{
real x[ND_ND]; //defining position vector//
real z;
face_t f;
begin_f_loop(f,thread)
{
F_CENTROID(x,f,thread)
z=x[2];
F_PROFILE(f,thread,position)=500+3*(z);
}
end_f_loop(f,thread)
}

my doubts are:
1. is there anything wrong in the code....?
2. does the way of defining the position vector is same for 2D and 3D
models?
3. during interpretation do i need to use command name as "cpp" or "cc -E"?
4. does this code works for 3D cases, since i have seen many examples
where 2D cases have used this code...

any help will be greatly appreciated....

thanks...

[obscureed]

Hi Sai Krishna,

1. You need a semicolon at the end of the line containing F_CENTROID. I am surprised that this compiled without errors.
2. Yes, the F_CENTROID command (and many others) are the same in 2D and 3D. It would not be safe to access x[2] in 2D, but this is correct in 3D.
3. I would not recommend interpretation at all.
4. This should be OK in 3D.

One issue with your code is probably the formula "500+3*(z)". The units of z are [metre], so you would need a domain size of 100m to reach 800K. This is possible, but you should check whether this is correct for your case.

Good luck!
Ed