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

Boundary conditions for k and epsilon on a periodic channel flow

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 22, 2019, 17:02
Default Boundary conditions for k and epsilon on a periodic channel flow
  #1
New Member
 
Maveryck Andres
Join Date: Mar 2019
Posts: 4
Rep Power: 2
Maveryck is on a distinguished road
Hi everyone,

I'm simulating a flow through a channel. I impose a mean velocity which traduces in a pressure gradient calculated by openFoam.


Here are boundary and initial conditions for Pressure:

Code:
boundaryField
{
    inlet
    {
        type            cyclic;
    }

    outlet
    {
        type            cyclic;
    }

    FixedWalls
    {
        type            zeroGradient;
    }

    frontAndBack
    {
        type            zeroGradient;
    }
}
Velocity Conditions:
Code:
boundaryField
{
     inlet 
     { 
        type            cyclic;
     } 

    outlet
    {
        type            cyclic;
    }

    FixedWalls
    {
        type            noSlip;
    }

    frontAndBack
    {
        type            noSlip;
    }
}
But it is not clear to me how to determine which boundary conditions are suitable for k and epsilon in this simulation. Until now I have tried the following:

Boundary conditions for k:

Code:
boundaryField
{
    inlet 
    { 
        type            cyclic;
    } 

    outlet
    {
        type            cyclic;
    }

    FixedWalls
    {
        type            kLowReWallFunction;
        value           uniform 0.75;
    }

    frontAndBack
    {
        type            kLowReWallFunction;
        value           uniform 0.75;
    }
}
B.C. for epsilon
Code:
boundaryField
{
    inlet 
    { 
        type            cyclic;
    } 

    outlet
    {
        type            cyclic;
    }

    FixedWalls
    {
        type            epsilonWallFunction;
        value           uniform 14;
    }

    frontAndBack
    {
        type            epsilonWallFunction;
        value           uniform 14;
    }
}
I'm getting small values for nut at the walls, which means I get higher WallShearStress values than expected and thus I don't get proper uPlus-yPlus profiles (which is what I'm looking for).

Has anyone worked with similar cases? any suggestion will be wellcome

Thanks, Maveryck.

Note: y+ at the walls have values around 5-8.
Maveryck is offline   Reply With Quote

Old   August 26, 2019, 02:58
Default
  #2
C-L
Member
 
Charlie Lloyd
Join Date: Feb 2016
Posts: 55
Rep Power: 5
C-L is on a distinguished road
Your boundary conditions are fine for those y+ values.
I expect this is because of the turbulence model you are using - k-epsilon is unsuitable for these flows unless you use highRe wall functions (i.e y+ levels are in the log-region).

Try v2-f if you want to use an epsilon-based model, or use an omega-based model if you don't want to use the four-equation one.
C-L is offline   Reply With Quote

Old   August 26, 2019, 05:08
Default
  #3
Member
 
Philipp Conen
Join Date: Apr 2019
Location: GER, NRW
Posts: 34
Blog Entries: 2
Rep Power: 2
pconen is on a distinguished road
When you have problems finding suitable values you can try out the cfd-online tools.
https://www.cfd-online.com/Tools/turbulence.php
__________________
Greetings

Philipp Conen
pconen is offline   Reply With Quote

Reply

Tags
boundary condition, cyclic boundary condition, rans model

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Radiation in semi-transparent media with surface-to-surface model? mpeppels CFX 11 August 22, 2019 07:30
Question about adaptive timestepping Guille1811 CFX 25 November 12, 2017 17:38
Basic Nozzle-Expander Design karmavatar CFX 20 March 20, 2016 08:44
An error has occurred in cfx5solve: volo87 CFX 5 June 14, 2013 17:44
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 07:00


All times are GMT -4. The time now is 04:10.