[thermal_joe]

Hello,
I’m facing a problem with the implementation of a 2d gaussian heat distribution (static).
The setup is simple, a steel plate with a top surface (where the heat flux is applied and the other surfaces are treated as a wall with a constant heat transfer coefficient.


The used udf is subsequently announced:
#include "udf.h"
#define rb 0.00235
#define f1 3
#define Q 1900

DEFINE_PROFILE(gauss_heat_flux,t,i)
{
real x[ND_ND];
real y,y1;
face_t f;
begin_f_loop(f,t)
{
F_CENTROID(x,f,t);
y = x[1];
y1=x[0];
F_PROFILE(f,t,i) = (Q/rb*rb)*exp(-(f1*(y*y+y1*y1))/rb*rb);
}
end_f_loop(f,t)
}

The result should be a small heated spot (at X=0,Y=0) with a gaussian heat flux distribution but what I got is a plain heated area.

I don't know what is the problem with my udf, something should be wrong.
Maybe someone has a clue about this problem.


With best regards
Joe