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

problem with p_rgh BC and recirculation inside my domain

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 25, 2013, 11:48
Default problem with p_rgh BC and recirculation inside my domain
  #1
New Member
 
Kevin Bob
Join Date: Jan 2013
Posts: 1
Rep Power: 0
kevin is on a distinguished road
Hi everybody,

I'm trying to replicate a wind tunnel experiment by simulating the propagation of Helium in a building model using OF 2.1.1 and the firefoam solver.

My problem consist of an inlet source from which an air/helium mixture enter the fluid domain with a specified velocity profile and an outlet which should represent the ambient air.

I use the homogeneousDynOneEqEddymodel. Everything is desactivated in firefoam except the specie transport equations and I use the thermoType hsPsiMixtureThermo<multiComponentMixture<constGasT hermoPhysics>>.

I have read many posts here about p_rgh saying that totalPressure should be used for atmospheric BC and buyoantPressure everywhere else and everything works fine except that most of my mixture stay inside the domain. The flow only exit the domain at the bottom of the sides patches and nothing cross the top BC.

I'm quite new to OF and i'm a litte bit loss between all the BC availables. I have already dig into the source code and found many interesting things but can't find a way to get it working.

fixedValue for the outlet give me the same result whereas specifying a gradient BC (zeroGradient, buyoantPressure or fixedFluxPressure) let the fluid get out but gives a wrong field as expected since the pressure is specified nowhere. OutletInlet seems a good choice for me but it crash the calculation.

Any advice ?

Thanks in advance for your answers.

0/p_rgh
Code:
dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 101325;

boundaryField
{
    walls
    {
        type            buoyantPressure;
        value           $internalField;
    }
    inlet
    {
        type            buoyantPressure;
        value           $internalField;
    }
    pressure_outlets
      {
        type            totalPressure;
        p0              $internalField;
        U               U;
        phi             phi;
        rho             rho;
        psi             none;
        gamma           0;
        value           $internalField;
      }
}
0/U
Code:
dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    walls
    {   
        type            fixedValue;
        value           uniform (0 0 0);
    }
    inlet
    {
        type            timeVaryingMappedFixedValue;
        setAverage      off;
    }
    pressure_outlets
    {
        type            pressureInletOutletVelocity;
        value           uniform (0 0 0);
        phi             phi;
    }
}
kevin is offline   Reply With Quote

Old   August 23, 2016, 11:29
Default
  #2
Senior Member
 
Lieven
Join Date: Dec 2011
Location: Leuven, Belgium
Posts: 299
Rep Power: 23
Lieven will become famous soon enough
Hi Kevin and others,

did you ever figure out how to solve this problem?

Kind regards,

Lieven
Lieven is offline   Reply With Quote

Reply


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



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