
[Sponsors] 
September 5, 2010, 22:18 
Thin Wall Heat Transfer BC for rhoSimpleFoam

#1 
Member
Stefano Wahono
Join Date: Aug 2010
Location: Melbourne, Australia
Posts: 42
Rep Power: 6 
Good afternoon everyone!
I'm just starting to look at using OpenFOAM to solve a heat transfer problem across a 'thin' (thickness approaching zero) wall in a high Re compressible flow. I am using the rhoSimpleFoam solver, and I don't currently have a separate solid mesh for the wall. I have problem with specifying a wall BC for temperature that will solve the heat transfer across a thin wall (conjugate fluidsolid heat transfer) without having a separate mesh for the solid. Specifying zeroGradient BC at the wall led to the wall being treated as an adiabatic wall. The appropriate temperature BC at the wall should be constant heat flux across the wall. I think Fluent does the heat transfer across the wall somehow by creating the 'shadow' wall, and then mapping the temperature on one side to the other side of the wall, and solve using an outer iteration to satisfy the constant heat flux condition across the wall. Can anyone please help me if there is already an existing BC with OF that does this? Thank you very much for your help. Best regards, Stefano
__________________
Stefano Wahono Defence Science and Technology Organisation Propulsion Systems 

September 5, 2010, 22:27 

#2 
Member
Stefano Wahono
Join Date: Aug 2010
Location: Melbourne, Australia
Posts: 42
Rep Power: 6 
I think I need to clarify this further.
The amount of heat flux through the wall is dependant on the flow solution.
__________________
Stefano Wahono Defence Science and Technology Organisation Propulsion Systems 

September 6, 2010, 07:38 

#3  
Senior Member
Join Date: Jan 2010
Location: Stuttgart
Posts: 129
Rep Power: 7 
Hi Stefano,
I used "groovyBC" to solve the problem you described. Here http://openfoamwiki.net/index.php/Co...s_are_defined: you can read how to use this boundary condition. But there are also many threads around groovyBC in this forum. For a thin wall with convective and radiative heat loss, my boundary is: Quote:


September 6, 2010, 21:18 

#4 
Member
Stefano Wahono
Join Date: Aug 2010
Location: Melbourne, Australia
Posts: 42
Rep Power: 6 
Hi Chrisi,
Thank you very much for your quick reply. That looks like a complex radiation heat trf model! Yes, I have been looking in the forum, and it seems that the most commonly used BC for the thin wall heat trf is the groovyBC. May I ask you several things regarding your thin wall BC? 1) where is alphaEff value being assigned? 2) I assume Cp is of the fluid? 3) What is alphaa? In my case, Tinf differs on both sides of the wall. > Hot flow, T_HOT ~ 600 C, MACH=0.5 =============THIN=STAINLESS=STEEL=WALL========= > Cold flow, T_COLD ~ 50 C, MACH = 0.1 4) So, what do you think is the appropriate value of Tinf and Trad in this case? 5) I assumed your BC does not incorporate the effect of convective heat transfer (i.e conv heat trf coeff and local Nusselt No)? So, basically this BC is similar to the 'fixedGradient' type, however with the ability to calculate and update the temperature gradient at the two sides of the 'wall' at every flow iteration? Thank you very much, again. I think OF is missing a very important feature here, a feature which has been taken for granted when using FLUENT Best Regards, Stefano
__________________
Stefano Wahono Defence Science and Technology Organisation Propulsion Systems 

September 7, 2010, 07:06 

#5 
Senior Member
Join Date: Jan 2010
Location: Stuttgart
Posts: 129
Rep Power: 7 
Hi Stefano,
alphaEff=(alpha_lam+alphat) is derived and need not be defined. cp is of the fluid (here air) alphaa is the heat transfer coefficient of the outer surface of the pipe. The boundary incorporate the effect of convection. It is a mixed boundary out of radiation and convection. For your temperatures, I do not exactly know how your case looks like, so I cannot give any advice for them. Regards Chrisi 

September 9, 2010, 00:00 

#6 
Member
Stefano Wahono
Join Date: Aug 2010
Location: Melbourne, Australia
Posts: 42
Rep Power: 6 
Thanks for your patience and help, Chrisi!
I will try this BC and let you know. Best regards, Stefano 

September 15, 2010, 01:53 

#7 
Member
Stefano Wahono
Join Date: Aug 2010
Location: Melbourne, Australia
Posts: 42
Rep Power: 6 
Hi Chrisi,
I have tested your thin wall heat transfer BC on rhoSimpleFoam using a simple case. The resulting 'wall temperature' is different to that calculated using Fluent. I found that the calculated wall temperature is sensitive to the 'given value of Tinf and Trad. From looking at your eqn, I know that Tinf and Trad is just an 'initial guess' value used to calculate the qDOT. The wall Ts is then improved using the seven qDot equations (or iterations). ==> Is my understanding correct? If yes, then I don't understand why the resulting wall Ts is sensitive to the Tinf and Trad? Maybe Tinf and Trad should be read from internalField at every flow iteration? What determines the Tinf and Trad you use? Should we use another value of 'valueFraction'? Thank you very much. Best Regards, Stefano P.S Yes, you are right, alphaeff is calculated in the solution, and it accounts for the heat convection.
__________________
Stefano Wahono Defence Science and Technology Organisation Propulsion Systems 

September 15, 2010, 04:11 

#8 
Senior Member
Join Date: Jan 2010
Location: Stuttgart
Posts: 129
Rep Power: 7 
Hi!
Tinf is the air temperature of the environment which cools or heats your boundary. Trad is the radiation temperature of your environment which can be different to the air temperature. Both temperatures are not initial guesses. This are fixed values for your environment. The temperatures Tw are calculated iterative . And Tw1 is an initial guess to get a first approximation for the wall temperature. The wall temperature is needed to calculate the massflux. I hope I could hepl you. Regards Chrisi 

November 16, 2010, 10:31 

#9 
Member
Antonio Liggieri
Join Date: Aug 2010
Posts: 66
Rep Power: 5 
Hello Chrisi & Swahono,
your calculation domains contain internal walls. How do you create your mesh? I tried to convert a 3Dfluent mesh with internal walls using the option "writeSets". This is possible for the fluentMeshToFoam converter and unfortunately NOT for fluent3DMeshToFoam. But I need to convert 3D Meshes. So do you have an idea, how deal with a 3DMesh having internal walls? Second Question. I'm using groovyBC as well. Also for heatFlux problems and it looks like this: SURFACE_SOURCE { type groovyBC; value uniform 283; valueExpression "283"; gradientExpression "gradT"; fractionExpression "0"; variables "heatFlux=2000;Cp0=1005;rho0=1.18;gradT=heatFl ux/(alphaEff * Cp0 * rho0);"; timelines ( ); } Up to now I used it in models without radiation. As soon as radiation is implemented the solver gives you the following error: > FOAM FATAL ERROR: Parser Error at "1.1118" :"field alphaEff not existing or of wrong type" "heatFlux/(alphaEff * Cp0 * rho0)" " ^^^^^^^^ " From function parsingValue in file PatchValueExpressionDriver.C at line 188. Does anybody have a clue, where this problem could come from? I compared the appearance of alphaEff in buoyantSimpleFoam and in buoyantSimpleRadiationFoam. The appearance is exactly the same!! Thanx in advance for any hints! Tony 

April 2, 2013, 20:15 

#10 
Senior Member
Mohammad Shakil Ahmmed
Join Date: Oct 2012
Location: AU
Posts: 122
Rep Power: 5 
The problem is in the declaration:FOAM FATAL ERROR:
Parser Error at "1.1118" :"field alphaEff not existing or of wrong type" "heatFlux/(alphaEff * Cp0 * rho0)" " ^^^^^^^^ "  you've not declared "alphaEff" in the variables  and one parsing error of > 1.1118 to be 1.11e18 

April 3, 2013, 15:03 

#11  
Senior Member
Join Date: Jan 2010
Location: Stuttgart
Posts: 129
Rep Power: 7 
Dear all.
Due to some diverging problems, I modified the bc a little bit. It now looks that way: Quote:
alphaEff is derived from the used solver. I used that bc for the simplePorousFoam solver. Kind regards Chrisi 

September 6, 2013, 12:31 

#12  
Senior Member

Quote:
What are Tw, d and l? Why you are iterating the procedure to find heatFlux7? Thank you, Mojtaba
__________________
Complex Heat & Flow Simulation Research Group If you can't explain it simply, you don't understand it well enough. "Richard Feynman" Last edited by Mojtaba.a; October 4, 2013 at 17:30. 

October 4, 2013, 11:49 

#13 
Senior Member
Join Date: Jan 2010
Location: Stuttgart
Posts: 129
Rep Power: 7 
Hi,
Tw is the wall temperature d is the wall thickness l is the thermal conductivity of the wall. The reason for the iterative solution is that you cannot calculate the heatflux through the wall from scratch. It depends on the adjusting wall temperature, which I calculate iteratively. I am sorry but I do not have any reference for that measure. I just compared the received results with FLUENT. That showed good agreement. Kind regards Chrisi 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Conjugate heat transfer: coupled wall temperature  Sarah  FLUENT  7  August 12, 2013 22:36 
Flow around pipes  heat transfer coefficient on the wall of pipe  doodek  Main CFD Forum  2  November 23, 2009 09:48 
how to model heat transfer for the coupled wall  jing113cn  FLUENT  0  June 22, 2009 03:56 
Concentric tube heat exchanger (AirWater)  Young  CFX  5  October 6, 2008 23:17 
Wall Heat Transfer Plot  Paul  CFX  0  April 22, 2005 04:43 