[sarah l]

Hi,
I am using non-equilibrium model in porous zone in Ansys Fluent 14.5. for specifying heat transfer coefficient for the fluid / solid interface, I have written the following udf

#include "udf.h"
DEFINE_PROFILE(Heat_trans_coeff,t,i)
{
cell_t c;
real Nu,Re,Pr;
real dens, visc; /*Fluid*/
real cond, cp; /* Fluid*/
real d_p;
real por;/*Porosity of the bed*/
d_p=0.06;
por=0.36

begin_c_loop(c,t)
{
dens = C_R(c,t); /*Density of fluid*/
visc = C_MU_L(c,t); /*Viscosity fluid*/
cond = C_K_L(c,t); /*Conductivity fluid*/
cp = C_CP(c,t); /*Specific heat fluid*/
Re = ND_MAG(C_U(c,t),C_V(c,t),C_W(c,t))*d_p*dens*por/visc ; /*Reynolds number has been calculated based on pore velocity*/
Pr = cp*visc/cond ;
Nu = 2.+ 1.1 * pow(Re,0.6) * pow(Pr,1./3.);
F_PROFILE(c,t,i) = Nu*cond/d_p ;
}
end_c_loop(c,t)
}
However this code cannot be interpreted in FLUENT. The error message is "c_loop_last: undeclared variable"
I would be grateful if somebody help me!

[anon_m]

Hello:

Were you able to solve your problem?

[sarah l]

well, I changed some of the statements to constatnt parameter and the it run!

[anon_m]

Sarah:
Thanks for your reply. Do you by any chance know how one can get the solid temperature at a specific point in the media by udf?

Thank you.

[sarah l]

you're welcome
while you use non equilibrium model , the software creats solid zone itself so you can observe the variation of Temp on that ! also you can create new plane with specific points to see more details.

regards

[othmann]

Dear Pr Sarah

I work on the heat transfer in porous media, and I need a UDF for the heat transfer coefficient fluid/solid that is a function of the solid temperarue and the fluid temperature

Thank's to answer me

Best regards

[junianatalia]

Good night,

I'm working with heat transfer and I've had serious problems. Can you help me?

[Bruno Machado]

Originally Posted by junianatalia

Good night,

I'm working with heat transfer and I've had serious problems. Can you help me?

explain your problem, post the code you wrote and we can try to help you.

[junianatalia]

heat_udf_1.h

It's my code. I'm working with tho phases in Eulerian Model. I'm having trouble describing each stage in my UFD.

The error is pointed out in this line:

heat_kmiec_1980 (c cell_t, Thread * ti, Thread * tj)

[junianatalia]

Originally Posted by Bruno Machado

explain your problem, post the code you wrote and we can try to help you.

Good morning.
Fallow my udf.

Thanks

[guilbemloco]

Dear members of the forum

I need a UDF for the heat transfer coefficient, considering the properties change with the temperature.
My function is very similar to the one Sarah posted here, if anyone knows how to correct that function that would be great!

Thank you,

G

[guilbemloco]

Hello Sarah,

Could you please show me the UDF you created after you changed some parameters to constant?

Thank you

[guilbemloco]

Well, I found that in this code there is a ";" missing after the declaration of "por=0.36".

I added that and I was able to interpret the function normally.

[junianatalia]

Originally Posted by guilbemloco

Well, I found that in this code there is a ";" missing after the declaration of "por=0.36".

I added that and I was able to interpret the function normally.

Hello,

Can you send me your UDF?

Thanks

[guilbemloco]

Originally Posted by junianatalia

Hello,

Can you send me your UDF?

Thanks

Sure, it is this:

#include "udf.h"
DEFINE_PROFILE(Heat_trans_coeff,t,i)
{
cell_t c;
real Nu,Re,Pr;
real dens, visc; /*Fluid*/
real cond, cp; /* Fluid*/
real d_p;
real por;/*Porosity of the bed*/
d_p=0.018;
por=0.3;

begin_c_loop(c,t)
{
dens = C_R(c,t); /*Density of fluid*/
visc = C_MU_L(c,t); /*Viscosity fluid*/
cond = C_K_L(c,t); /*Conductivity fluid*/
cp = C_CP(c,t); /*Specific heat fluid*/
Re = ND_MAG(C_U(c,t),C_V(c,t),C_W(c,t))*d_p*dens*por/visc ; /*Reynolds number has been calculated based on pore velocity*/
Pr = cp*visc/cond ;
Nu = 2.+ 1.1 * pow(Re,0.6) * pow(Pr,1./3.);
F_PROFILE(c,t,i) = Nu*cond/d_p ;
}
end_c_loop(c,t)
}

I was able to interpret in windows and to compile in linux systems, the reason is not clear thought.

[junianatalia]

Originally Posted by guilbemloco

Sure, it is this:

#include "udf.h"
DEFINE_PROFILE(Heat_trans_coeff,t,i)
{
cell_t c;
real Nu,Re,Pr;
real dens, visc; /*Fluid*/
real cond, cp; /* Fluid*/
real d_p;
real por;/*Porosity of the bed*/
d_p=0.018;
por=0.3;

begin_c_loop(c,t)
{
dens = C_R(c,t); /*Density of fluid*/
visc = C_MU_L(c,t); /*Viscosity fluid*/
cond = C_K_L(c,t); /*Conductivity fluid*/
cp = C_CP(c,t); /*Specific heat fluid*/
Re = ND_MAG(C_U(c,t),C_V(c,t),C_W(c,t))*d_p*dens*por/visc ; /*Reynolds number has been calculated based on pore velocity*/
Pr = cp*visc/cond ;
Nu = 2.+ 1.1 * pow(Re,0.6) * pow(Pr,1./3.);
F_PROFILE(c,t,i) = Nu*cond/d_p ;
}
end_c_loop(c,t)
}

I was able to interpret in windows and to compile in linux systems, the reason is not clear thought.

Thanks! guilbemloco