
[Sponsors] 
February 19, 2012, 09:58 
a problem of natural convection boundary conditions

#1 
New Member
Adrian
Join Date: Feb 2012
Posts: 6
Rep Power: 7 
I want to solve a problem like this，water in cylindrical tank at 293.15K，and a single tube located in the center of the tank，I choose fixed wall temperature of tube at 363.15K，for reducing the calculation，I just modeled 1/4 of tank. My working fluid is water 97 in library of CFX 12, when i choose that，CFX choose full bouyance model，so without boussinique model. i set gravity accelerate 9.8, ref density 1000kg/m3. At first，I set heating tube wall as fixed temperature wall，and adiab wall conditions both at bottom and tank wall，also at the top free surface，but i cannot get result，Imbalance of Henergy is 100%，i confused maybe energy cannot fulfill conservation rule？but when i changed top free surface boundary condition to fixed temperature wall of 293.15K，it still cannot convergence，also because of Henergy imbalance. calculations are all in transient，about 500s. And mesh elements about 30000 or 300000 have same result，mesh quality is good. Then i choose opening boundary for top free surface，in transient calculation i got the result but temperature gradient is so small bulk region is still 293.15K，and i cannot capture the velocity vector. sincerely hope someone can give me advice，I am a new bird in CFX calculation，and my English is poor. Thank you! 

February 19, 2012, 17:26 

#2 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 13,803
Rep Power: 107 
I see many problems with what you have said so far:
* For a steady state run you need to balance the heat inputs to the system. Your hot tube is putting heat into the system, so where is it coming out? * Why did you choose a 90 degree slice? For low Rayleigh numbers (if you don't know what Rayleigh number is then google it, you will need to know it for a simulation like this) the flow is 2D so can be modelled with a very thin slice. But for moderate to high Rayleigh numbers the flow is 3D and it is unlikely you can apply symmetry planes at all. * For moderate to high Rayleigh numbers the flow will be transient, so a steady state run will not work. * Rather than trying all sorts of boundary conditions and seeing what works, how about you place boundary conditions appropriate to what you are modelling? What did the experiment you are trying to reproduce do? * I assume by water 97 you mean the IAPWS mode. Do not use these unless you need the very high levels of accuracy it will give, because they significantly increase the numerical instability are are much harder to converge. Just use standard constant properties water while you are trying things out. Only go back to IAPWS once the simulation works properly and you need the additional accuracy... and have time to sort out all the convergence issues. 

February 19, 2012, 21:27 

#3  
New Member
Adrian
Join Date: Feb 2012
Posts: 6
Rep Power: 7 
Quote:
And I want to know，tank free surface with atmosphere pressure should be simulated in what boundary condition？ 

February 20, 2012, 01:54 

#4 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 13,803
Rep Power: 107 
No, the first thing to do is to work out the Rayleigh number. Only then will you know what type of flow you are going to get, then you make the decision of a 3D or 2D mesh from that. What is the Rayleigh number of your flow? If you don't know, then look it up and work it out.
If your top surface is a free surface you have more options. Does the free surface move such that it affects the flow? Does liquid evaporate at it? What is the heat conditions on the other side of the interface  convection with the gas and/or radiation? 

February 20, 2012, 02:47 

#5  
New Member
Adrian
Join Date: Feb 2012
Posts: 6
Rep Power: 7 
Quote:
i know that the rule of laminar turn to the turbulance is about Ra=10^8~10^10, but choose different characteristic length will get very different Ra number，for vertical tube when i choose the length of tube for characteristic, Ra number reach 10^15, in some natural convection correlation it is recommended tube length will be appropriate characteristic length. in CFX which one i should choose tube diameter or tube length. finally, i want to calculate a very long tube, as long as 7 meters. i want to study the variation of boundary layer thickness，and when the water in the tank has formed thermal stratification what difference in heat transfer. after the calculation of single tube，i need calculated tube bundles horizontal and vertical. so any recommendation, thank you! 

February 20, 2012, 06:14 

#6 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 13,803
Rep Power: 107 
In that the flow is almost certainly full 3D and transient. You can model the top free surface with a slip wall with a convection boundary condition if you like, this is pretty crude but a reasonable start.


February 29, 2012, 02:19 

#7 
New Member
Adrian
Join Date: Feb 2012
Posts: 6
Rep Power: 7 
I've tried modeling full part of tank with a slender tube in it. Tank with 0.1*0.1*2, and tube with 6mm diameter, bouyance model concerned, water with constant properties, in transient calculation Henergy imbanlance as high as 100% after 30s. I find that temperature gradient appear near the wall within a very thin layer. But I' ve calculated the same problem just with a small cylingderical tank and a relative short tube in it. In transient calculation it has a good result, water in tank formed thermal stratified. I assumed that maybe mesh quality is not good. But I've tried more 。。。。mesh. It still not worked. I want to know if the natural convection heat transfer calculation between a very slender tube with tank can be modeled in CFX?
THANK YOU VERY MUCH! 

February 29, 2012, 06:57 

#8 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 13,803
Rep Power: 107 
If you are getting that high imbalances are you sure your simulation is converged?
CFX can model what you are talking about, the problem is with your set up. It is just a matter of finding the problem. 

Tags 
boundary condition, bouyancy, natural conection 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Domain Imbalance  HMR  CFX  5  October 10, 2016 05:57 
Wind turbine simulation  Saturn  CFX  45  February 8, 2016 05:42 
Implementation of boundary conditions for FVM  Tom  Main CFD Forum  7  August 26, 2014 05:58 
Natural Convection Problem  Helium  marzoa  STARCCM+  0  April 18, 2011 14:12 
Problem with Boundary conditions  Mahiboobswamulu  Main CFD Forum  10  August 26, 2003 13:24 