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

what is correct boundary condition for K-w sst model at y+<1?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 4, 2013, 14:35
Default what is correct boundary condition for K-w sst model at y+<1?
  #1
Senior Member
 
Join Date: Jun 2011
Posts: 141
Rep Power: 6
mechy is on a distinguished road
Hi

I found some thread about k-w sst boundary conditions
unfortunately they are not clear in low y+

please let me know what is correct boundary conditions for k, w and nut
for y+<5

Best Regards
mechy is offline   Reply With Quote

Old   August 4, 2013, 15:42
Default
  #2
New Member
 
Andre Weiner
Join Date: Aug 2012
Posts: 27
Rep Power: 4
andre.weiner is on a distinguished road
Hello!

omega:
Menter's automatic wall funcion works for all parts of the momentum boundary layer,
in openfoam it's the omegaWallFunction http://www.cfd-online.com/Wiki/Near-...k-omega_models
or you use the solution for the viscos sublayer to set the boundary value
http://turbmodels.larc.nasa.gov/sst.html (there is an equation for omega_wall)
for y^+<5 both give the same results

k:
fixed value 0 or
zero gradient
if you look at same results for k near the wall (from other simulations for example), you will see that both are valid (again for y^+<5)

nut:
calculated or
fixed value 0 or
nutLowReWallFunction (sets the value to zero and provides a function for calculating y^+ which is needed for the yPlusRAS tool)

This is what i tried so far and it worked for mi :-)

Best regards, Andre
andre.weiner is offline   Reply With Quote

Old   August 5, 2013, 03:51
Default
  #3
Senior Member
 
Join Date: Jun 2011
Posts: 141
Rep Power: 6
mechy is on a distinguished road
Quote:
Originally Posted by andre.weiner View Post
Hello!

nut:
calculated or
fixed value 0 or
nutLowReWallFunction (sets the value to zero and provides a function for calculating y^+ which is needed for the yPlusRAS tool)

This is what i tried so far and it worked for mi :-)

Best regards, Andre
thanks so much for your answer

what function should be set for calculating y+ ? and how set it ?
can you attached an example for k,w,nut (initial file at folder 0)?

Best Regards
mechy is offline   Reply With Quote

Old   August 5, 2013, 10:29
Default
  #4
Senior Member
 
Joachim
Join Date: Mar 2012
Location: Atlanta
Posts: 139
Rep Power: 6
Joachim is on a distinguished road
Actually, I think that nutLowReWallFunction is a wall function for low-Re turbulence models. If you use LaunderSharmaKE and have y+ ~ 60, you can use it for example.

Also, try using 1e-9 (or any small value) rather than 0 when imposing the BC.

Good luck!
Joachim is offline   Reply With Quote

Old   August 5, 2013, 14:17
Default
  #5
New Member
 
Andre Weiner
Join Date: Aug 2012
Posts: 27
Rep Power: 4
andre.weiner is on a distinguished road
for k:

Code:
    myWall
    {
        type            zeroGradient;
    }
// or
    myWall
    {
        type            fixedValue;
        value           uniform 0;
    }
for omega:

Code:
    myWall
    {
        type            omegaWallFunction;
    }
for nut:

Code:
    myWall
    {
        type            fixedValue;
        value           uniform 0;
    }
// or
    myWall
    {
        type            calculated;
        value           uniform 0;    // optional value entry
    }
// what i prefer
    myPatch
    {
        type            nutLowReWallFunction;
    }
I don't know for what the nutLowReWallFunction schould or schloud not be, but if you take little look at the code:

it sets nut to zero:
Code:
tmp<scalarField> nutLowReWallFunctionFvPatchScalarField::calcNut() const
{
    return tmp<scalarField>(new scalarField(patch().size(), 0.0));
}
and provides a function for y+:

Code:
tmp<scalarField> nutLowReWallFunctionFvPatchScalarField::yPlus() const
{
....
    return y*sqrt(nuw*mag(Uw.snGrad()))/nuw;
}
Regards, Andre
andre.weiner 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
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 07:00
CFX13 Post Periodic interface EtaEta CFX 7 December 8, 2011 18:15
Epsilon boundary condition for walls in LamBremhorst LowRe model maruthamuthu_venkatraman OpenFOAM 0 March 16, 2010 04:27
SST model and boundary layers Rodrigo Basniak CFX 7 March 27, 2008 15:19


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