[Yadu krishnan]

I have used following UDF for neumann boundary condition as heat flux profile. How can I validate if it is working fine?


#include "udf.h"
DEFINE_PROFILE(Heat_flux,thread,position)
{
Thread *t;
Thread *tf;
cell_t c;
face_t f;
real x[ND_ND];
real cell[ND_ND];
int n;
int k = 1;
real wf;
real nc;
real temp;
real tempgrad ;
real Twall=100;
real delta;
begin_c_loop_all(c,t)
{
c_face_loop(c,t,n)
{
tf = C_FACE_THREAD(c,t,n);
if(THREAD_TYPE(tf) == THREAD_F_WALL)
{
f = C_FACE(c,t,n);
F_CENTROID(x,f,t);
wf = x[1];
C_CENTROID(cell,c,t);
nc = cell[1];
delta = wf-nc;
temp = C_T(c,t);
tempgrad= (temp - Twall)/delta ;
F_PROFILE(f, thread, position) = k * tempgrad;
}
}


}
end_c_loop_all(c,t)


}

[pakk]

Apply it in a simple geometry, and see if the results are the same as what you calculate by hand.

[AlexanderZ]

you have a typo
to be

F_CENTROID(x,f,tf);