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

totalPressure Boundary condition

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

Like Tree1Likes
  • 1 Post By klio

Reply
 
LinkBack Thread Tools Display Modes
Old   February 4, 2010, 06:39
Default totalPressure Boundary condition
  #1
Senior Member
 
BastiL
Join Date: Mar 2009
Posts: 471
Rep Power: 11
bastil is on a distinguished road
Dear Foamers,

we are currently struggling with the total pressure inlet BC (totalPressure). We do not fully understand the meaning of some of the require parameters in the "p"-File:
Code:
/* PRESSURE INLET */
    patchname_pressureinlet
    {
        type            totalPressure;
        p0              uniform 1.e5; /* Total pressure p0 = p + 1/2*rho*U^2 */
        U               U;
        phi             phi;
        rho             rho;/none;
        psi             psi;/none;
        gamma           1;
    }
We are not aware of the meaning of phi, rho, psi and gamma. Looking into the sources I found the following error:
Code:
        FatalErrorIn
        (
            "totalPressureFvPatchScalarField::updateCoeffs()"
        )   << " rho or psi set inconsitently, rho = " << rhoName_
            << ", psi = " << psiName_ << '.' << nl
            << "    Set either rho or psi or neither depending on the "
               "definition of total pressure." << nl
            << "    Set the unused variables to 'none'."
            << "\n    on patch " << this->patch().name()
            << " of field " << this->dimensionedInternalField().name()
            << " in file " << this->dimensionedInternalField().objectPath()
            << exit(FatalError);
We are unsure how to set phi, rho, psi and gamma "consitantly". We want to use it with the rhoSimpleFoam-solver.
  • Does it depend on the solver type? (pressure-based/density-based?).
  • However, we see differences when using the different formulations.
  • It can also be seen from the sources that calculation of the boundary values are depending on the formulation. So we need to know how to set these values.
It would be very helpful if somebody can comment on this. Thanks.
bastil is offline   Reply With Quote

Old   February 25, 2010, 15:35
Default totalPressure Boundary condition
  #2
New Member
 
sankarv
Join Date: Feb 2010
Posts: 26
Rep Power: 7
sankarv is on a distinguished road
Dear Bastil

Did you find answers to your question ? I also have the same question and have not received any reply yet. Please let me know if you have figured out ?

Thanks
Vaidya
sankarv is offline   Reply With Quote

Old   April 16, 2012, 23:18
Default
  #3
Member
 
Saba Saeb
Join Date: Dec 2010
Location: Erlangen, Germany
Posts: 32
Rep Power: 6
saba_saeb is on a distinguished road
Exactly same problem, any idea so far?!
saba_saeb is offline   Reply With Quote

Old   April 17, 2012, 03:06
Default
  #4
Member
 
Klio
Join Date: Mar 2012
Location: Austria
Posts: 45
Rep Power: 5
klio is on a distinguished road
Quote:
Originally Posted by bastil View Post

/* PRESSURE INLET */
patchname_pressureinlet
{
type totalPressure;
p0 uniform 1.e5; /* Total pressure p0 = p + 1/2*rho*U^2 */
U U;
phi phi;
rho rho;/none;
psi psi;/none;
gamma 1;
}


"totalPressureFvPatchScalarField::updateCoeffs ()"
<< " rho or psi set inconsitently, rho = " << rhoName_
<< ", psi = " << psiName_ << '.' << nl
<< " Set either rho or psi or neither depending on the "
"definition of total pressure." << nl
<< " Set the unused variables to 'none'."
<< "\n on patch " << this->patch().name()
<< " of field " << this->dimensionedInternalField().name()
<< " in file " << this->dimensionedInternalField().objectPath()
<< exit(FatalError);
[.
I'm really new using OpenFoam but I guess you only have to add U, phi, rho or psi if you have set or calculated initial BC for those! (Maybe with potentialFoam)
And it is also depending on incompressible or compressible simulation. Either you have a field for compressible simulations or just the value in the transport properties.
If you don't add them they will be calculated automatically.
And I also guess that psi is about particle movement depending on the velocity (streamfunction) and directly depending on rho so if you fix one of them the other will be calculated!

But may I'm wrong
saba_saeb likes this.

Last edited by klio; April 17, 2012 at 03:25.
klio is offline   Reply With Quote

Old   July 4, 2013, 09:09
Default
  #5
New Member
 
Camille Bilger
Join Date: Jul 2013
Posts: 15
Rep Power: 4
kmou is on a distinguished road
Hello,

I know this is an old thread but I was looking around for the same answer and I found in which header file OpenFoam describes these variables:
http://foam.sourceforge.net/docs/cpp/a04465_source.html
kmou is offline   Reply With Quote

Reply

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
Boundary Conditions Thomas P. Abraham Main CFD Forum 20 July 7, 2013 05:05
Transient outlet boundary condition problem jwillie2000 CFX 1 December 7, 2009 18:07
Axis Boundary Condition..what is it? CFDtoy FLUENT 6 February 13, 2007 06:51
How to set boundary condition in Fluent for the fo Peiyong FLUENT 1 November 10, 2006 12:44
How to resolve boundary condition problem? sam FLUENT 2 July 20, 2003 02:19


All times are GMT -4. The time now is 16:44.