Wall Normal Temperature Gradient

Daniel Tanner · April 30, 2009, 13:59

I would like to access the temperature gradients normal to my wall boundary. Consider steady laminar flow in a cylindrical pipe. I solve the energy equation using a constant temperature at the wall and wish to access the wall normal temperature gradients at each face of the wall boundary (i.e., dT/dn).

How can I do this? (I am familiar with UDF's). There appears to be no F_T_G() macro and I am not sure how to calculate the face value from the cell centroid (C_F_G()).

Help would be appreciated.

contrarian · May 2, 2009, 06:14

You can assume that the cell gradient is equal to the face gradient and dot it with the face area.
You can also get the face normal temperature gradient from the flux: divide the face heat flux with the thermal conductivity.

Daniel Tanner · May 2, 2009, 07:53

Thanks for the reply. I understand how to do it conceptually but I'm having problems implementing it. I have attached a basic UDF below which I cannot get to work! I just want to save the wall temperature to UDM. It compiles and executes without error but the UDM values are zero everywhere (I have set the UDM in Fluent).

++++++++++++++++++++++++
#include "udf.h"

DEFINE_ON_DEMAND(Test)
{
Domain *domain=Get_Domain(1);
Thread *t;
face_t f;

int Zone_ID=5; /* Found in Fluent BC panel and identifies Wall */

t=Lookup_Thread(domain,Zone_ID);

begin_f_loop(f,t)
{
F_UDMI(f, t, 0) = F_T(f,t);
}
end_f_loop(f,t)
}
++++++++++++++++++++++++

Any ideas why this does not work?

contrarian · May 2, 2009, 09:51

F_T may not be populated on load. Could you try after doing one iteration?

Daniel Tanner · May 2, 2009, 11:57

Sorry and thanks for the help. I found the problem. The UDF was working fine, but Fluent cannot display the F_UDMI data as it is stored at the face, i.e., it can only display cell centered data.

I'll have to think of another way to do it.

contrarian · May 3, 2009, 00:01

You can may be store the F_UDMI at C_UDMI(c0,t0).

Maryam-A · September 20, 2015, 13:10

Hi
I would appreciate if anyone can help me.
I wanted to return derivative of UDS as wall boundary value, but I face with "fatal signal" that I found this because of boundary condition UDS.

this is part of my udf for B.C of UDS( volume fraction of nanoparticle in my case)

DEFINE_PROFILE(alpha_bc1,thread,position)
{
Thread *t;
cell_t c;
.
.
.

begin_f_loop(f,thread)
{

C_UDSI_G(c,t,0)[0]=-(D_T/D_B)*C_T_G(c,t)[0];
F_PROFILE(f,thread,position)=C_UDSI_G(c,t,0)[0];
}
end_f_loop(f,t)
}

April 30, 2009, 13:59	Wall Normal Temperature Gradient	#1
Daniel Tanner Member Daniel Tanner Join Date: Apr 2009 Posts: 54 Rep Power: 17	I would like to access the temperature gradients normal to my wall boundary. Consider steady laminar flow in a cylindrical pipe. I solve the energy equation using a constant temperature at the wall and wish to access the wall normal temperature gradients at each face of the wall boundary (i.e., dT/dn). How can I do this? (I am familiar with UDF's). There appears to be no F_T_G() macro and I am not sure how to calculate the face value from the cell centroid (C_F_G()). Help would be appreciated.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
surface tension as function of normal gradient	tom	FLUENT	2	September 7, 2016 03:48
How to input the wall temperature b.c.?	Jimmy	Siemens	9	March 14, 2008 09:31
How to input the wall temperature b.c.?	Jimmy	Siemens	16	March 7, 2008 14:11
Variable temperature wall boundary condition	Virginie	FLUENT	2	October 2, 2007 09:23
momentum equation (Normal to wall_	kk	Main CFD Forum	0	July 2, 2006 23:10

May 2, 2009, 06:14		#2
contrarian New Member Contrarian Join Date: Apr 2009 Posts: 5 Rep Power: 16	You can assume that the cell gradient is equal to the face gradient and dot it with the face area. You can also get the face normal temperature gradient from the flux: divide the face heat flux with the thermal conductivity.

May 2, 2009, 07:53		#3
Daniel Tanner Member Daniel Tanner Join Date: Apr 2009 Posts: 54 Rep Power: 17	Thanks for the reply. I understand how to do it conceptually but I'm having problems implementing it. I have attached a basic UDF below which I cannot get to work! I just want to save the wall temperature to UDM. It compiles and executes without error but the UDM values are zero everywhere (I have set the UDM in Fluent). ++++++++++++++++++++++++ #include "udf.h" DEFINE_ON_DEMAND(Test) { Domain domain=Get_Domain(1); Thread t; face_t f; int Zone_ID=5; /* Found in Fluent BC panel and identifies Wall */ t=Lookup_Thread(domain,Zone_ID); begin_f_loop(f,t) { F_UDMI(f, t, 0) = F_T(f,t); } end_f_loop(f,t) } ++++++++++++++++++++++++ Any ideas why this does not work?

May 2, 2009, 09:51		#4
contrarian New Member Contrarian Join Date: Apr 2009 Posts: 5 Rep Power: 16	F_T may not be populated on load. Could you try after doing one iteration?

May 2, 2009, 11:57		#5
Daniel Tanner Member Daniel Tanner Join Date: Apr 2009 Posts: 54 Rep Power: 17	Sorry and thanks for the help. I found the problem. The UDF was working fine, but Fluent cannot display the F_UDMI data as it is stored at the face, i.e., it can only display cell centered data. I'll have to think of another way to do it.

May 3, 2009, 00:01		#6
contrarian New Member Contrarian Join Date: Apr 2009 Posts: 5 Rep Power: 16	You can may be store the F_UDMI at C_UDMI(c0,t0).

September 20, 2015, 13:10		#7
Maryam-A Member Join Date: Nov 2014 Posts: 42 Rep Power: 11	Hi I would appreciate if anyone can help me. I wanted to return derivative of UDS as wall boundary value, but I face with "fatal signal" that I found this because of boundary condition UDS. this is part of my udf for B.C of UDS( volume fraction of nanoparticle in my case) DEFINE_PROFILE(alpha_bc1,thread,position) { Thread t; cell_t c; . . . begin_f_loop(f,thread) { C_UDSI_G(c,t,0)[0]=-(D_T/D_B)C_T_G(c,t)[0]; F_PROFILE(f,thread,position)=C_UDSI_G(c,t,0)[0]; } end_f_loop(f,t) }