# Boundary condition for transient natural convection

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

June 2, 2013, 09:03
Boundary condition for transient natural convection
#1
New Member

Eshita Pal
Join Date: Mar 2012
Posts: 6
Rep Power: 7
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
{
}
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;
}

}
Attached Images
 tank_with_inlet_outlets.jpg (31.2 KB, 116 views)

 June 2, 2013, 10:58 #2 Senior Member   HECKMANN Frédéric Join Date: Jul 2010 Posts: 237 Rep Power: 10 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: 7 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: 10 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: 7 Thank you. I shall try if this works.

 July 4, 2013, 10:24 #6 New Member   Korichi Abdelkader Join Date: Jan 2013 Posts: 11 Rep Power: 6 Dear fredo490, I ahve similar problem, natural convection in channel (with inlet and oultlet), I have a problem to specify the boundary conditions for p and p-rgh (inlet and outlet) I need help

 August 14, 2013, 12:04 #7 Senior Member     Jose Rey Join Date: Oct 2012 Posts: 131 Rep Power: 10 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/...aeKim-OFW8.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. huadian001 likes this. 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

Join Date: May 2011
Posts: 231
Rep Power: 9
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:
``` dimensions      [0 2 -2 0 0 0 0];internalField   uniform 0;boundaryField{    inlet    {        type            zeroGradient;    }    outlet    {                type fixedValue;        value \$internalField;                  }    wall    {        type            zeroGradient;    }}  ```
and for p_rgh

PHP Code:
``` dimensions      [0 2 -2 0 0 0 0];internalField   uniform 0;boundaryField{    wall{type  fixedFluxPressure;rho rhok;value uniform 0;}outlet{type  fixedFluxPressure;rho rhok;value uniform 0;}inlet{type fixedValue;value uniform 0;} }  ```
for T

PHP Code:
``` dimensions      [0 0 0 1 0 0 0];internalField   uniform 363.15;boundaryField{    wall    {         type            fixedValue;        value           uniform 343.15;    }    inlet    {        type            fixedValue;        value           uniform 363.15;    }    outlet    {        type            zeroGradient;    }}  ```
and alpha

PHP Code:
``` dimensions      [0 2 -1 0 0 0 0];internalField   uniform 0;boundaryField{    wall    {        type            alphatJayatillekeWallFunction;        Prt             0.85;        value           uniform 0;    }    inlet    {        type            calculated;            }    outlet    {       type            calculated;    }}  ```
I have a laminar flow and i am getting results but they seems like wrong...
Can you see any mistake in the BC?

Quote:
 Originally Posted by JR22 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/...aeKim-OFW8.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.

 July 16, 2016, 04:06 #9 Senior Member   Nima Samkhaniani Join Date: Sep 2009 Location: Tehran, Iran Posts: 1,216 Blog Entries: 1 Rep Power: 17 may be its late , however look the this document to set boundary condition __________________ Telegram channel (https://telegram.me/openfoam4Iranian) My Weblog in Persian(http://openfoam.blogfa.com/) My Personal Website (http://nimasamkhaniani.ir/)

November 8, 2016, 12:40
#10
New Member

Join Date: Nov 2016
Posts: 1
Rep Power: 0
Quote:
 Originally Posted by EshitaPal Thank you. I shall try if this works.
Did you solve the problem?

 Tags boundary condition, natural convection, p_rgh

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post peob OpenFOAM Running, Solving & CFD 3 February 3, 2017 11:54 shenying0710 CFX 7 March 26, 2013 05:13 thomasyangfly FLOW-3D 3 September 11, 2012 10:14 Ank OpenFOAM 16 July 30, 2012 04:26 modrio FLUENT 2 August 11, 2005 12:29