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

OpenFOAM wall functions

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

Like Tree1Likes
  • 1 Post By fredo490

Reply
 
LinkBack Thread Tools Display Modes
Old   May 6, 2013, 21:34
Default OpenFOAM wall functions
  #1
New Member
 
Hossein
Join Date: Feb 2012
Posts: 13
Rep Power: 4
Gildeh is on a distinguished road
Hello All,

I am working with pisoFoam solver and running several RANS turbulence models (k-e, RNG, ...). I believe that these turbulence models use the standard wall functions for near wall regions. For instance, when I look at inside the nut dictionary at 0 in k-e turbulence modeling, it uses nutkWallFunction. I have two questions:

1. What is this nutkWallFunction exactly (the concept, the formulation, and where can I find the implementation in OF)?

2. How can I change the roughness of the wall function within OpenFOAM? I believe that wall functions must have a roughness coefficient for the model calculation that we should be able to change it. More specifically, I want to see the effect of changing the roughness in my model! How is it possible to change this coefficient?

Thank you for reading this thread and thanks for your time.

Gildeh
Gildeh is offline   Reply With Quote

Old   May 7, 2013, 01:14
Default
  #2
Senior Member
 
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 233
Rep Power: 7
fredo490 is on a distinguished road
1) your answer is in the source code:
OpenFOAM-2.1.x / src / turbulenceModels / incompressible / RAS / derivedFvPatchFields / wallFunctions / nutWallFunctions / nutkWallFunction / nutkWallFunctionFvPatchScalarField.H ,gives you:
Code:
Description
    Boundary condition for turbulent (kinematic) viscosity when using wall
    functions, based on turbulence kinetic energy.
and: OpenFOAM-2.1.x / src / turbulenceModels / incompressible / RAS / derivedFvPatchFields / wallFunctions / nutWallFunctions / nutkWallFunction / nutkWallFunctionFvPatchScalarField.C gives you:
Code:
  forAll(nutw, faceI)
    {
        label faceCellI = patch().faceCells()[faceI];

        scalar yPlus = Cmu25*y[faceI]*sqrt(k[faceCellI])/nuw[faceI];

        if (yPlus > yPlusLam_)
        {
            nutw[faceI] = nuw[faceI]*(yPlus*kappa_/log(E_*yPlus) - 1.0);
        }
    }

    return tnutw;
The source can also be found here:
https://github.com/OpenFOAM/OpenFOAM...hScalarField.H

2) check the relative wall function for rough wall. You will find the description in the .H file and the formula used in the .C file.
mgg likes this.
fredo490 is offline   Reply With Quote

Old   May 7, 2013, 01:53
Default
  #3
Senior Member
 
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 233
Rep Power: 7
fredo490 is on a distinguished road
I just saw your answer while I was doing another research:
http://foam.sourceforge.net/docs/cpp/a01446.html
fredo490 is offline   Reply With Quote

Reply

Tags
pisofoam, roughness, turbulence models, wall function

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
Enhanced Wall Treatment paduchev FLUENT 18 April 11, 2014 07:37
Wall functions - questions about implementation AlmostSurelyRob OpenFOAM 33 April 8, 2014 02:42
[ICEM] Export ICEM mesh to Gambit / Fluent romekr ANSYS Meshing & Geometry 1 November 26, 2011 12:11
Automatic/hybrid/compound wall functions Far OpenFOAM Running, Solving & CFD 0 November 12, 2011 03:50
Wall Functions pierre OpenFOAM Running, Solving & CFD 0 October 1, 2005 13:13


All times are GMT -4. The time now is 14:58.