[valerhain]

hi dear friends
I am supposed to write a udf for heat generation rate which is a function at this form:





I prepared this code for simulating :

#include "udf.h"
#define EPS 0.01
#define SIGMA 0.004
DEFINE_SOURCE(cell_x_source, cell, thread, dS, eqn)
{
real x[ND_ND];
real source;
C_CENTROID(x,cell,thread);

source=EPS*exp(((x[0]-0.4)*(x[0]-0.4)+(x[1]-0.004)*(x[1]-0.004))/(-SIGMA*SIGMA));
dS[eqn]=0;
return source;
}


but is not work. this code dont have any effect in my results. dear friend I have less than 5 days for simulating this. please help me and I will thankfull of you...