
[Sponsors] 
June 2, 2013, 09:03 
Boundary condition for transient natural convection

#1 
New Member
Eshita Pal
Join Date: Mar 2012
Posts: 6
Rep Power: 6 
Hello,
I am using Openfoam 2.2 and trying to simulate natural circulation in a tank with multiple inlets at the bottom and multiple outlets at the top. As I run the transient buoyantBoussinesqPimpleFoam solver, the flow reverses very quickly, and the water starts flowing out of the inlet. Are my boundary conditions appropriate? Thanks. U file dimensions [0 1 1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { pipe { type fixedValue; value uniform (0 0 0); } inlet1 { type pressureInletOutletVelocity; value uniform (0.1414 0 0); } inlet2 { type pressureInletOutletVelocity; value uniform (0.130636565 0.054111437 0); } inlet3 { type pressureInletOutletVelocity; value uniform (0.099984898 0.099984898 0); } outlet { type zeroGradient; } symmetry { type symmetryPlane; } wall { type fixedValue; value uniform (0 0 0); } } p_rgh file dimensions [0 2 2 0 0 0 0]; internalField uniform 0; boundaryField { pipe { type buoyantPressure; rho rhok; value uniform 0; } inlet1 { type totalPressure; p0 uniform 51.993; // this is equal to g*(height of the tank) U U; phi phi; rho rhok; psi none; gamma 1; value uniform 0; } inlet2 { type totalPressure; p0 uniform 51.993; U U; phi phi; rho rhok; psi none; gamma 1; value uniform 0; } inlet3 { type totalPressure; p0 uniform 51.993; U U; phi phi; rho rhok; psi none; gamma 1; value uniform 0; } outlet { type buoyantPressure; rho rhok; value uniform 0; } symmetry { type symmetryPlane; } wall { type buoyantPressure; rho rhok; value uniform 0; } } 

June 2, 2013, 10:58 

#2 
Senior Member
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 237
Rep Power: 9 
Question, why using buoyantBoussinesqPimpleFoam ? Do you have any heat transfer ? Is it to consider the Gravity ?
Is your flow really driven by the total pressure ? The pressure you have set looks quite small (especially for water). 

June 2, 2013, 11:27 

#3 
New Member
Eshita Pal
Join Date: Mar 2012
Posts: 6
Rep Power: 6 
The tubes inside the geometry has surface heat flux b.c. There is a total heat generation of 2MW in a volume of 198 m3. I do consider gravity for the natural convection flow generated. The idea is to generate natural convection flow that will exit through the outlets at the top.
The working fluid is water. I am not very sure about the pressure condition. The tank is part of a natural circulation loop. The flow should thus be generated due to the heat transfer from the hot tubes. The pressure is not the driving force. 

June 2, 2013, 11:48 

#4 
Senior Member
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 237
Rep Power: 9 
My mistake, I forgot your topic tittle while I was reading it
Well, I don't think that a pressureInletOutletVelocity is suitable in your case. Personally, I would have used a simple "zeroGradient" for the velocity boundaries. But I don't know if it is enough, you can give it a try. To me, your pressure boundaries are ok but not the velocity. If your problem persists, you can try to increase slightly the pressure at the inlet to "force" the fluid to go up and once the global flow behavior is catch, you can return to the original pressure value. 

June 2, 2013, 12:15 

#5 
New Member
Eshita Pal
Join Date: Mar 2012
Posts: 6
Rep Power: 6 
Thank you. I shall try if this works.


August 14, 2013, 12:04 

#7 
Senior Member
Jose Rey
Join Date: Oct 2012
Posts: 128
Rep Power: 9 
If anybody is still interested, I just saw that there was a discussion on natural convection in the 2013 OpenFoam Workshop:
http://www.openfoamworkshop2013.org/...aeKimOFW8.tar I also had some problems with getting buoyantBoussinesqPimpleFoam to converge. I changed buoyantPressure BCs for fixedFluxPressure and the problem started converging. Give it a try. Last edited by JR22; August 19, 2013 at 20:39. Reason: added boundary condition change for p_rgh 

April 29, 2015, 10:51 

#8  
Senior Member
rkhr
Join Date: May 2011
Posts: 231
Rep Power: 8 
Hi I am using buoyantBoussinesqSimpleFoam for heat transfer in single pipe..so I have one inlet and outlet and wall..
I have BC for pressure as following: PHP Code:
PHP Code:
PHP Code:
PHP Code:
Can you see any mistake in the BC? thanks in advance! Quote:


July 16, 2016, 04:06 

#9 
Senior Member
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,193
Blog Entries: 1
Rep Power: 16 
may be its late , however look the this document to set boundary condition
__________________
Telegram channel (https://telegram.me/openfoam4Iranian) My Weblog (http://openfoam.blogfa.com/) Training Course on OpenFOAM at (http://www.isme.ir/) 

Tags 
boundary condition, natural convection, p_rgh 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
mixed inflow/outflow downstream boundary condition question  peob  OpenFOAM Running, Solving & CFD  2  August 14, 2014 09:07 
CFX fails to calculate a diffuser pipe flow  shenying0710  CFX  7  March 26, 2013 05:13 
natural convection boundary condition  thomasyangfly  FLOW3D  3  September 11, 2012 10:14 
Boundary condition for natural convection  Ank  OpenFOAM  16  July 30, 2012 04:26 
boundary condition in natural ventilation  modrio  FLUENT  2  August 11, 2005 12:29 