CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Boundary condition for transient natural convection

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   June 2, 2013, 09:03
Default Boundary condition for transient natural convection
  #1
New Member
 
Eshita Pal
Join Date: Mar 2012
Posts: 6
Rep Power: 5
EshitaPal is on a distinguished road
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;
}



}
Attached Images
File Type: jpg tank_with_inlet_outlets.jpg (31.2 KB, 56 views)
EshitaPal is offline   Reply With Quote

Old   June 2, 2013, 10:58
Default
  #2
Senior Member
 
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 236
Rep Power: 7
fredo490 is on a distinguished road
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).
fredo490 is offline   Reply With Quote

Old   June 2, 2013, 11:27
Default
  #3
New Member
 
Eshita Pal
Join Date: Mar 2012
Posts: 6
Rep Power: 5
EshitaPal is on a distinguished road
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.
EshitaPal is offline   Reply With Quote

Old   June 2, 2013, 11:48
Default
  #4
Senior Member
 
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 236
Rep Power: 7
fredo490 is on a distinguished road
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.
fredo490 is offline   Reply With Quote

Old   June 2, 2013, 12:15
Default
  #5
New Member
 
Eshita Pal
Join Date: Mar 2012
Posts: 6
Rep Power: 5
EshitaPal is on a distinguished road
Thank you. I shall try if this works.
EshitaPal is offline   Reply With Quote

Old   July 4, 2013, 10:24
Default
  #6
New Member
 
Korichi Abdelkader
Join Date: Jan 2013
Posts: 6
Rep Power: 4
AbdelkaderDZ is on a distinguished road
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
AbdelkaderDZ is offline   Reply With Quote

Old   August 14, 2013, 12:04
Default
  #7
Senior Member
 
JR22's Avatar
 
Jose Rey
Join Date: Oct 2012
Posts: 122
Rep Power: 8
JR22 will become famous soon enough
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.

Last edited by JR22; August 19, 2013 at 20:39. Reason: added boundary condition change for p_rgh
JR22 is offline   Reply With Quote

Old   April 29, 2015, 10:51
Default
  #8
Senior Member
 
rkhr
Join Date: May 2011
Posts: 209
Rep Power: 7
Kanarya is on a distinguished road
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?
thanks in advance!


Quote:
Originally Posted by JR22 View Post
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.
Kanarya is offline   Reply With Quote

Reply

Tags
boundary condition, natural convection, p_rgh

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


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 FLOW-3D 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


All times are GMT -4. The time now is 02:38.