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

simpleFoam: switching flow direction with kOmegaSST RAS Model

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 2, 2014, 16:25
Default simpleFoam: switching flow direction with kOmegaSST RAS Model
  #1
Member
 
Jason G.
Join Date: Sep 2009
Location: St. Louis, IL
Posts: 89
Rep Power: 13
JasonG is on a distinguished road
I have been looking at pressure drops for internally bounded, steady-state, incompressible fluid flow cases with the simpleFoam solver. Typically, I specify a velocity at the entrance to the fluid domain and a reference pressure of 0 at the exit of the fluid domain.

A new case I am working on requires me to specify a velocity at the exit and a reference pressure at the entrance. I have been able to successfully do this with laminar assumptions, but as soon as I attempt to switch to the kOmegaSST turbulence solver I run into convergence issues. It appears the kinetic turbulence values become discontinuous near the entrance.

Any help and/or recommendations are greatly appreciated. Below are the current BC's I am implementing:

"intialConditions"
Code:
flowVelocity         (0 0 165.760612715333);
pressure             0;
turbulentKE          194.091657099187;
turbulentOmega       1731.03910939215;
kinematicvis         nu [ 0 2 -1 0 0 0 0 ] 0.0398447454916618;   // in^2/s
density             0.00007839;           //  Value of the density(sslug/in^3) pstatic = lbf/in^2 

//RAS PROPERTIES

RAS_MODEL         kOmegaSST; //kOmegaSST;    //laminar;
turb_on_off         on;    //on


#inputMode           merge
"p"

Code:
#include        "initialConditions"

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform $pressure;

boundaryField
{
 
    outlet_1          
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }



    inlet_1      
      {
        type            fixedValue;
        value           $internalField;
    }
   

    boundary       
    {
        type            zeroGradient;
    }

    symmetry      
    {
        type        symmetryPlane;
    }


}
"u"

Code:
#include        "initialConditions"

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    outlet_1   
    {
        type            outletInlet;
        outletValue     $flowVelocity;
        value           $flowVelocity;
    }




    inlet_1        
    {
        type            outletInlet;
        outletValue     $internalField;
        value           $internalField;
    }



    boundary       
    {
        type            fixedValue;
        value           $internalField;
    }

    symmetry       
    {
        type        symmetryPlane;
    }


}
"k"

Code:
#include        "initialConditions"

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform $turbulentKE;

boundaryField
{

    outlet_1
    {
       type  fixedValue;
       value $internalField;
    }



    inlet_1
    {
        type            outletInlet;
        outletValue     $internalField;
        value           $internalField;
    }



    
    boundary
    {
        type            kqRWallFunction;
        value           $internalField;
    }

    symmetry       
    {
        type        symmetryPlane;
    }

}
"nut"
Code:
dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet_1
    {
        type            calculated;
        value           uniform 0;
    }



    outlet_1
    {
        type            calculated;
        value           uniform 0;
    }



  
     boundary
    {
        type            nutLowReWallFunction;
        value           uniform 0;
    } 



    symmetry       
    {
        type        symmetryPlane;
    }
}
"omega"
Code:
 
#include        "initialConditions"

dimensions      [0 0 -1 0 0 0 0];

internalField   uniform $turbulentOmega;

boundaryField
{


    outlet_1
    {
       type  fixedValue;
       value $internalField;
    }



    inlet_1
    {
        type            outletInlet;
        outletValue     $internalField;
        value           $internalField;
    }



    boundary
    {
        type            omegaWallFunction;
        value           $internalField;
    }



    symmetry       
    {
        type        symmetryPlane;
    }

}
JasonG is offline   Reply With Quote

Old   September 25, 2014, 08:26
Default
  #2
Senior Member
 
Joachim Herb
Join Date: Sep 2010
Posts: 601
Rep Power: 18
jherb is on a distinguished road
You might try the fixedMean boundary condition for the velocity, k and omega at the outlet (assuming the the specified flow is going "outwards").
Then you do not enforce a constant value over the whole outlet but just the (area averaged) mean for the three quantities. So the turbulent block profile within your flow domain is not destroyed at the outlet.
jherb 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
Use of k-epsilon and k-omega Models Jade M Main CFD Forum 39 March 19, 2020 19:27
interFoam vs. simpleFoam channel flow comparison DanM OpenFOAM Running, Solving & CFD 12 January 31, 2020 16:26
About Some Concepts:Laminar flow, turbulent flow, steady flow and time-dependent flow Jing Main CFD Forum 8 October 5, 2018 18:02
Eul-Eul flow, k-e-kp-ep-Theta Turbulence model us FLUENT 5 April 5, 2011 03:29
Ho to model the flow in a Regenerator alexandrepereira OpenFOAM Running, Solving & CFD 2 August 25, 2010 18:35


All times are GMT -4. The time now is 22:32.