[Neethu_N]

I need to give a parabolic heat flux at the top wall of a 2D square metal plate. I wrote the following code for it.

#include "udf.h"

DEFINE_PROFILE(heat_flux_parabola, thread, heat)
{
real P[ND_ND];
real x = P[0];
face_t f;

begin_f_loop(f, thread)
{
F_CENTROID(P, f, thread);
F_PROFILE(f, thread, heat)= (-4*pow(10,8)*pow((x-0.05),2)) + pow(10,6);
}
end_f_loop(f, thread)
}

The graph of the heat flux plotted in Matlab (1 Heat flux matlab.png) and that of wall heat flux plotted in CFD-Post (1 Heat flux cfd.png) is attached.

The matlab plot is as expected, with a peak value of 10^6, but the fluent variation is cut off at around 10^5, with a linear variation in the middle. Can somebody please explain why this is happening? Is there a a limit for heat flux in FLUENT?

[Neethu_N]

Answering my own question, there was a mistake in the code.
x should be assigned the value inside the loop only. That solved it.
Thank you!