[neerajpdh]

Hello all,

I intent to find out heat flux at the inlet boundary and feed the mass flux that is the function of heat flux.
I have written a UDF for mass flux. I am not sure whether its correct or not as at each face cell i am getting the same heat flux.

#include "udf.h"

DEFINE_PROFILE(mass_flux_AP, t, i)
{
face_t f;
real NV_VEC(farea);
real mflux_ap;
real area_ap;
real heat_ap;
real heat_flux_ap;


begin_f_loop(f, t)
{


F_AREA(farea, f, t);
area_ap = NV_MAG(farea);
heat_ap = BOUNDARY_HEAT_FLUX(f, t);
heat_flux_ap= fabs(heat_ap /area_ap);

if ( heat_flux_ap >= 1200000)
{
mflux_ap = (heat_flux_ap * 0.0000002082) + 0.538;
}
else
{
mflux_ap = 5.9302;
}

F_PROFILE(f, t, i) = mflux_ap;
}
end_f_loop(f, t)
}


DEFINE_PROFILE(mass_flux_binder, t, i)
{
face_t f;
real NV_VEC(farea);
real mflux;
real area;
real heat;
real heat_flux;


begin_f_loop(f, t)
{

F_AREA(farea, f, t);
area = NV_MAG(farea);
heat = BOUNDARY_HEAT_FLUX(f, t);
heat_flux= fabs(heat / area);

if ( heat_flux >= 1200000)
{
mflux = (heat_flux * 0.00005374 )+ 1.82;
}
else
{
mflux =7.4;
}

F_PROFILE(f, t, i) = mflux;
}
end_f_loop(f, t)

}


Please help me in this regard.