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

pressure in buoyantBoussinesqSimpleFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 27, 2017, 05:36
Default pressure in buoyantBoussinesqSimpleFoam
  #1
New Member
 
Hugo
Join Date: Jul 2017
Posts: 14
Rep Power: 8
Gohu8 is on a distinguished road
Hello everyone,

I use buoyantBoussinesqSimpleFoam solver. My case is close to the windaroundBuildings tutorial. Each parameter seems to work fine, but when I want to read the pressure, I select p and the scale starts from approximately -712 until 707. So I can read that the average pressure around the buildings is about -360. The units are supposed to be in m²/s². So to have the real pressure in Pascal, I am supposed to multiply it by rho. But here rho = 1.2 kg/m^3 ~ 1 kg/m^3, so it should be a pressure of -360 Pa which is completely wrong !

PS: for the same time, when I read p_rgh, the scale starts from -100 until 17... and the average pressure around buildings seem to be close to 0....

I used simflow for this case. In initial conditions: p-rho*g*h=0.
In boudary conditions: for the inlet: p-rho*g*h is zero gradient
for the outel: p-rho*g*h is fixed value =0
Gohu8 is offline   Reply With Quote

Old   July 27, 2017, 15:30
Default
  #2
Senior Member
 
M. C.
Join Date: May 2013
Location: Italy
Posts: 286
Blog Entries: 6
Rep Power: 16
student666 is on a distinguished road
Have you set internal field for p and p_rgh equal to 0 or 10e5 Pa?

Sent from my ASUS_X008D using CFD Online Forum mobile app
student666 is offline   Reply With Quote

Old   August 8, 2017, 06:38
Default
  #3
New Member
 
Hugo
Join Date: Jul 2017
Posts: 14
Rep Power: 8
Gohu8 is on a distinguished road
Quote:
Originally Posted by student666 View Post
Have you set internal field for p and p_rgh equal to 0 or 10e5 Pa?

Sent from my ASUS_X008D using CFD Online Forum mobile app

EDIT: I found in the folders that p and p_rgh are both defined as:

internalField uniform 0.0


Here is some screenshots:

for p:




for p_rgh:

Gohu8 is offline   Reply With Quote

Old   August 8, 2017, 08:10
Default
  #4
New Member
 
Hugo
Join Date: Jul 2017
Posts: 14
Rep Power: 8
Gohu8 is on a distinguished road
Quote:
Originally Posted by student666 View Post
Have you set internal field for p and p_rgh equal to 0 or 10e5 Pa?

Sent from my ASUS_X008D using CFD Online Forum mobile app
In p and p_rgh I have:

internalField unfinorm 0.0


Here some screenshots:

for p:




for p_rgh:

Gohu8 is offline   Reply With Quote

Old   June 13, 2022, 11:50
Default
  #5
Member
 
Mahmoud
Join Date: Nov 2020
Location: United Kingdom
Posts: 43
Rep Power: 5
Mahmoud Abbaszadeh is on a distinguished road
Dear Foamers,

It’s been a while that I’ve had the same problem in setting the correct BC for similar problems. Here is the review of how I solved this issue:

Problem: steady state open channel flow simulation. In this scenario, normally the simpleFoam is employed. The BCs are given by:


BCs: inlet, outlet, wall, topWall

U:

inlet
{
type fixedValue;
value uniform (0 0.1589 0);
}

outlet
{
type zeroGradient;
}

wall
{
type noSlip;
}

atmosphere
{
type slip;
}

p:

atmosphere
{
type zeroGradient;
}

inlet
{
type zeroGradient;
}
outlet
{
type fixedValue;
value uniform 0;
}

wall
{
type fixedFluxPressure; //remember that we cannot use Drichlet and Neuman BC at the same time
value uniform 0;
}


As one can see, selecting fixedValue BC for the outlet pressure is physically incorrect even though one can get result. To resolve this issue, I decided to use the hydrostatic pressure as the outlet BC for the pressure. To this end, one need to either define gravity in the governing equation or else use a steady state solver that has gravity in the equations. I went for the second approach. The only available solver for incompressible fluids that has this feature is buoyantBoussinesqSimpleFoam. When this solver is selected, the BCs needs to be defined for the U, p, and p_rgh. Generally, p_rgh=p-rho*g*h. I have searched a lot and finally could find the correct way of setting BC for the foregoing parameters.


U:

inlet
{
type fixedValue;
value uniform (0 0.1589 0);
}

outlet
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
inletValue uniform (0 0 0);
}

wall
{
type noSlip;
}

atmosphere
{
type slip;
}


P:

atmosphere
{
type zeroGradient;
}

inlet
{
type zeroGradient;
}

outlet
{
type calculated;
value $internalField;
}

wall
{
type zeroGradient;

}



P_rgh

inlet
{
type fixedFluxPressure;
}

outlet
{
type fixedValue;
value uniform 0;
}

wall
{
type fixedFluxPressure;
}

atmosphere
{
type fixedFluxPressure;
}


let me know if it is not clear.

Keep Foaming
Mahmoud Abbaszadeh is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Periodic flow using Cyclic - comparison with Fluent nusivares OpenFOAM Running, Solving & CFD 30 December 12, 2017 06:35
Pressure loss Velocity coupling CFXMUFFIN CFX 1 February 6, 2016 05:43
Discharge of Pressure Vessel into Pipe with Regulator gajowni2 System Analysis 0 October 31, 2015 19:57
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 07:27
Unsteady pressure differential between inlet and outlet of the pipe for single phase joshi20h FLUENT 0 September 26, 2012 13:41


All times are GMT -4. The time now is 06:35.