[etudiant_IITB]

Hi all.
I am really in need of help regarding the following problem. Your concern is appreciated.
I am trying to solve a 2D solid rectangular plate conduction problem (length in x and height in y with bottom left corner of plate as the origin). The plate is insulated on the left, right and bottom edge but is subjected to a spatially (in x) varying HEAT FLUX on the top edge. For the top edge, I have made the following UDF code but it is not giving the symmetric solution. The function that you see below for the heat flux is a spatially symmetric function with respect to the centre of the plate. I have paid attention in defining the function correctly with respect to the coordinate system in use.
I believe there should be some problem with my UDF C language code below (as I have not clearly understood the use of the x[ND_ND] array).

I would be grateful if someone can please put in a small amount of their time and help me out with this problem.
//***********//
#include "udf.h"
DEFINE_PROFILE(heat_flux,thread,i)
{
real x[ND_ND];
face_t f;
real k;
begin_f_loop(f,thread)
{
k=F_CENTROID (x,f,thread);

F_PROFILE(f,thread,i) = (-3.073e+20)*k^8 + (6.145e+19)*k^7+ (-4.966e+18)*k^6 + (2.072e+17)*k^5 -(4.718e+15)*k^4 + (5.647e+13)*k^3 + (-3.084e+11)*k^2 + (6.515e+08)*k + (-2.148e+05);
// this is the function for the heat flux //
}
end_f_loop(f,thread);
}

[brunoc]

F_CENTROID(x,f,thread) will write the coordinates of the centroid on your x array:

x[0], x[1], x[2] = (x, y and z coordinates)

Use those to access each coordinate of the centroid in your equation.