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

Why do we need to specify k omega on wall surface by using komega sst model

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 22, 2020, 13:51
Post Why do we need to specify k omega on wall surface by using komega sst model
  #1
New Member
 
Hui Zhang
Join Date: Jun 2020
Posts: 14
Rep Power: 5
Isaiah_HZ is on a distinguished road
Hi foamers,

I am a fresh guy in openfoam. Recently, i started with flow past cylinder. And when I set the parameters of k and omega for the inlet and wall surface, some confusions came into my mind. why we need to set wall surface omega and k? They should be calculated by openfoam. Also, for the pressure boundary condition, near wall, I set zerogradient by following the tutorials, but for stagnation point, there is a nonzero gradient near wall surface. I do not understand, can you guys help me explain? Thank you a lot.
Isaiah_HZ is offline   Reply With Quote

Old   June 25, 2020, 08:00
Default
  #2
Senior Member
 
Join Date: Dec 2019
Location: Cologne, Germany
Posts: 355
Rep Power: 8
geth03 is on a distinguished road
can you share your 0-dict files?
geth03 is offline   Reply With Quote

Old   June 29, 2020, 11:55
Default
  #3
New Member
 
Hui Zhang
Join Date: Jun 2020
Posts: 14
Rep Power: 5
Isaiah_HZ is on a distinguished road
Quote:
Originally Posted by geth03 View Post
can you share your 0-dict files?
Thanks for your reply!
The following is my 0-dict files
PHP Code:
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  7
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    
version     2.0;
    
format      ascii;
    class       
volScalarField;
    
location    "0";
    
object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// k = 3/2*(Uref * I)^2 = 0.002159
// I = 0.16*Re^(-1/8)= 0.03794
dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0.002159;

boundaryField
{
    
inlet
    
{
        
type            fixedValue;
        
value           0.002159;
    }
    
outlet
    
{
        
type            zeroGradient;
    }
    
topwall
    
{
        
type            symmetryPlane;
    }
    
bottomwall
    
{
        
type            symmetryPlane;
    }
    
cylinder
    
{
        
type            fixedValue;
        
value           uniform 0.0;
    }
    
front
    
{
        
type            empty;
    }

    
back
    
{
        
type            empty;
    }
}


// ************************************************************************* // 
PHP Code:
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  7
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    
version     2.0;
    
format      ascii;
    class       
volScalarField;
    
location    "0";
    
object      nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{
    
inlet
    
{
        
type            calculated;
        
value           uniform 0;
    }
    
outlet
    
{
        
type            calculated;
        
value           uniform 0;
    }
    
topwall
    
{
        
type            symmetryPlane
        
value           uniform 0;
    }
    
bottomwall
    
{
        
type            symmetryPlane
        
value           uniform 0;
    }
    
cylinder
    
{
        
type            calculated
        
value           uniform 0;
    }
    
front
    
{
        
type            empty;
    }

    
back
    
{
        
type            empty;
    }
}


// ************************************************************************* // 
PHP Code:
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  7
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    
version     2.0;
    
format      ascii;
    class       
volScalarField;
    
object      omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


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

internalField   uniform 11765;

boundaryField
{
    
inlet
    
{
        
type            fixedValue;
        
value           uniform 11765;
    }

    
outlet
    
{
        
type            zeroGradient;

    }
    
topwall
    
{
        
type            symmetryPlane
    }
    
bottomwall
    
{
        
type            symmetryPlane
    }
    
cylinder
    
{
        
type            fixedValue;
        
value           uniform 238884;
    }
    
front
    
{
        
type            empty;
    }

    
back
    
{
        
type            empty;
    }
}


// ************************************************************************* // 
PHP Code:
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  7
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    
version     2.0;
    
format      ascii;
    class       
volScalarField;
    
object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{
    
inlet
    
{
        
type           zeroGradient;
    }

    
outlet
    
{
        
type           fixedValue;
    
value           uniform 0;
    }
    
    
topwall
    
{
    
type               symmetryPlane;
    }

    
bottomwall
    
{
    
type           symmetryPlane;
    }
    
    
cylinder
    
{
    
type           zeroGradient;
    }
    
front
    
{
        
type            empty;
    }

    
back
    
{
        
type            empty;
    }
}

// ************************************************************************* // 
PHP Code:
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  7
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    
version     2.0;
    
format      ascii;
    class       
volVectorField;
    
object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform (0 0 0);

boundaryField
{
     
inlet
    
{
        
type            fixedValue;
        
value           uniform (1 0 0);
    }

    
outlet
    
{
        
type            zeroGradient;
    }
    
    
topwall
    
{
    
type               symmetryPlane;
    }

    
bottomwall
    
{
    
type               symmetryPlane;
    }
    
    
cylinder
    
{
    
type               noSlip;
    }

    
front
    
{
        
type            empty;
    }

    
back
    
{
        
type            empty;
    }

}


// ************************************************************************* // 
Isaiah_HZ is offline   Reply With Quote

Old   June 30, 2020, 06:01
Default
  #4
Senior Member
 
Join Date: Dec 2019
Location: Cologne, Germany
Posts: 355
Rep Power: 8
geth03 is on a distinguished road
at walls you can use wall function from the turbulence models:

for k for example you can use:

wallName
{
type kqRWallFunction;
value $internalField
}

there are also wall functions available for epsilon or nut ...

check out this document:
http://www.tfd.chalmers.se/~hani/kur...nfoamFinal.pdf
geth03 is offline   Reply With Quote

Old   June 30, 2020, 07:23
Default
  #5
Senior Member
 
Mikko
Join Date: Jul 2014
Location: The Hague, The Netherlands
Posts: 243
Rep Power: 12
Flowkersma is on a distinguished road
Hi Hui,

Just like for the other transport equations, you will need to specify boundary conditions for the k and omega. If you read the original papers by Menter et al, they also discuss about the boundary conditions.

Best, Mikko
Flowkersma is offline   Reply With Quote

Old   June 30, 2020, 10:18
Default
  #6
New Member
 
Hui Zhang
Join Date: Jun 2020
Posts: 14
Rep Power: 5
Isaiah_HZ is on a distinguished road
Quote:
Originally Posted by geth03 View Post
at walls you can use wall function from the turbulence models:

for k for example you can use:

wallName
{
type kqRWallFunction;
value $internalField
}

there are also wall functions available for epsilon or nut ...

check out this document:
http://www.tfd.chalmers.se/~hani/kur...nfoamFinal.pdf
thank you Geth!
But for wall functions, we still need to specify a value for wall surface. Can you explain why the fixed valued is used, can we calculated it (I mean openfoam calculated it by some equations)?
Isaiah_HZ is offline   Reply With Quote

Old   June 30, 2020, 10:21
Default
  #7
New Member
 
Hui Zhang
Join Date: Jun 2020
Posts: 14
Rep Power: 5
Isaiah_HZ is on a distinguished road
Quote:
Originally Posted by Flowkersma View Post
Hi Hui,

Just like for the other transport equations, you will need to specify boundary conditions for the k and omega. If you read the original papers by Menter et al, they also discuss about the boundary conditions.

Best, Mikko
Mikko, thank you for your help!
Could you send the link or website to download this paper?

best,
Hui
Isaiah_HZ is offline   Reply With Quote

Old   July 1, 2020, 05:05
Default
  #8
Senior Member
 
Mikko
Join Date: Jul 2014
Location: The Hague, The Netherlands
Posts: 243
Rep Power: 12
Flowkersma is on a distinguished road
Hi Hui,

For instance, this paper describes the boundary conditions. If you are looking which model exactly OpenFOAM uses then have a look at the comments in the kOmegaSST header file .
Flowkersma is offline   Reply With Quote

Old   July 1, 2020, 09:21
Default
  #9
Senior Member
 
Join Date: Dec 2019
Location: Cologne, Germany
Posts: 355
Rep Power: 8
geth03 is on a distinguished road
Quote:
Originally Posted by Isaiah_HZ View Post
thank you Geth!
But for wall functions, we still need to specify a value for wall surface. Can you explain why the fixed valued is used, can we calculated it (I mean openfoam calculated it by some equations)?
i just use it to start the simulation, in the beginning the initial values are not correct, but given a really small time step for the start of the simulation, the values will be corrected.

for my case, i did conduct some test, to see how dependent my results would be with regards to initial values of turbulence properties.
the result is, that the turbulence values will adjust themselves to correct values, if the time step initially is small enough.
geth03 is offline   Reply With Quote

Reply

Tags
boundaries condition, flow past cylinder, komega sst 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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
One more time: wall functions - SST CellZone OpenFOAM Pre-Processing 22 June 11, 2021 13:27
Wrong multiphase flow at rotating interface Sanyo CFX 14 February 7, 2017 17:19
Basic Nozzle-Expander Design karmavatar CFX 20 March 20, 2016 08:44
Superlinear speedup in OpenFOAM 13 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 05:36
SST model with rough wall? David CFX 4 November 30, 2005 11:22


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