
[Sponsors] 
July 21, 2017, 18:48 
How to set Variable Thermal Conductivity?

#1 
Member
Join Date: May 2017
Posts: 47
Rep Power: 6 
Hello,
I'm simulating a flow inside a tube, and the boundary condition at the wall is "Heat Transfer Coefficient", i.e, the thermal conductivity of the material of the tube (I will call it k from now on) over the wall thickness. The problem is that k is a function of the wall temperature (k = f(T)). Is there anyway I could get the temperature at each cell face of the wall to write an expression at the "Setup", e.g., Code:
k = const x Temperature@wall Could anybody help me? Thanks in avance. 

July 22, 2017, 06:10 

#2 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 16,840
Rep Power: 132 
Heat transfer coefficient is not usually defined from the conductivity and thickness of the wall material. The HTC is usually defined as the heat transfer due to fluid motion next to a wall. The wall conductivity and thickness are usually accounted for differently. Are you sure you are talking about the right thing?


July 22, 2017, 13:52 

#3 
Member
Join Date: May 2017
Posts: 47
Rep Power: 6 
I'll try to make myself more clear. As you can see from the attached figure, I need the Heat Flux (q") to calculate the temperature at the outside surface of the tube (T2) and the Temperature T1 to calculate my heat transfer coefficient.
The tube walls are not part of the computational domain. Capturar.PNG 

July 23, 2017, 03:15 

#4 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 16,840
Rep Power: 132 
Is this steady state or transient?


July 23, 2017, 06:58 

#5 
Member
Join Date: May 2017
Posts: 47
Rep Power: 6 

July 23, 2017, 18:33 

#6 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 16,840
Rep Power: 132 
In steady state you can combine the convective outer boundary and the conduction in the tube to a combined boundary on the fluid domain. One approach could be to add the thermal resistance of the convective section and the conduction in the pipe, and apply that as a heat flux boundary condition to the fluid. You can do this using CEL expressions and you will be able to variable k (or anything else for that matter).


July 23, 2017, 18:51 

#7  
Member
Join Date: May 2017
Posts: 47
Rep Power: 6 
Quote:


July 24, 2017, 02:39 

#8 
Senior Member
urosgrivc
Join Date: Dec 2015
Location: Slovenija
Posts: 365
Rep Power: 8 
And what is your boundary condition at your most outer wall?
I dont understand, do you just want (T) dependant material properties? if so you can do it via user function create a table of T and conductivity than set it in material properties as a form of expresion conductivity(T) 

July 24, 2017, 06:47 

#9 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 16,840
Rep Power: 132 
You don't need to calculate T1, the solver does it. If you need T1 to evaluate material properties then just use a CEL expression/interpolation function which evaluates material properties as a function of T  it will use T from the previous coefficient loop, so will converge on the correct value of T as your simulation converges.


July 25, 2017, 06:07 

#10  
Member
Join Date: May 2017
Posts: 47
Rep Power: 6 
Quote:


July 25, 2017, 06:56 

#11 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 16,840
Rep Power: 132 
I don't normally do this but I was in the mood for a bit of light entertainment so here goes:
Let the temperature at the inside of the wall be T1, outside T2 and ambient air temp T3. The wall is thickness L, thermal conductivity k and the heat transfer coefficient from T2 to T3 is h. In the wall, q"=k dT/dX which we will approximate as = k (T1T2)/L In the convection q"=h dT = h(T2T1). Equate the T2 terms and write in terms of q" and you get q" = k (T1  T3)/(L + k/h) Now you have a function for heat flux at the inside wall as a function of T1 and known variables. You can use this as the heat flux for the wall boundary condition and it will model a combined conducting wall and convection. You can make any of the constants functions of temperature, time, position or anything. I recommend you check my maths, but that's the general idea. 

July 27, 2017, 12:59 

#12 
Member
Join Date: May 2017
Posts: 47
Rep Power: 6 
ghorrocks, thank you so much for your attention, but I think I didn't make myself clear yet.
The point is that I don't know how to tell CFX that I want the value of T1. If I create an expression named "T1", I'll just have to type Code:
T Code:
areaAve(T)@wall 

July 27, 2017, 19:04 

#13 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 16,840
Rep Power: 132 
Yes, for the variable T1 just use the "T" variable in CFX. This will get the local wall temperature.


July 28, 2017, 11:12 

#14 
Member
Join Date: May 2017
Posts: 47
Rep Power: 6 
Thanks. And is there anyway I could get the q" on each cell face? I just managed to get the area average, using
Code:
areaAve(Wall Heat Flux)@wall 

July 29, 2017, 06:27 

#15 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 16,840
Rep Power: 132 
The variable Wall Heat Flux is the local q".


Tags 
coefficient, conductivity, temperature, variable 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Constant velocity of the material  Sas  CFX  15  July 13, 2010 08:56 
Install Openfoam in Suse 100 setup environment  h34gk1  OpenFOAM Installation  1  September 26, 2006 20:08 
Problem installing on 64bit with ver13  jonititan  OpenFOAM Installation  0  April 28, 2006 05:45 
Replace periodic by inletoutlet pair  lego  CFX  3  November 5, 2002 20:09 
Info: Short Course On Thermal Design of Electronic Equipment  Arnold Free  Main CFD Forum  0  August 10, 1999 10:18 