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

Wall functions

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 29, 2002, 13:51
Default Wall functions
  #1
John Grimm
Guest
 
Posts: n/a
As I mentioned in my previous message I am new to CFD modelling. At the moment I am still getting to grips with what I should be doing and at the moment I am looking at wall functions. Wall functions are used to model the near wall region, but which bit exactly? I use 3 layers to define the near wall region: The viscous sublayer where the flow is almost laminar, a buffer layer where viscous forces and turbulence are equally important and a fully turbulent or log-law layer where the velocity is a logarithmic function of the distance from the wall. After the log-law layer the flow is fully turbulent and the wall has no effect (I call this the core). Some books/manuals, such as Fluent, say that wall fuctions are used to model the viscous sublayer and buffer layer (from this I must assume that turbulence equations are used to model the log law and core). Others say that the wall functions model the viscous sublayer, buffer layer and log-law layer. Which is correct?

Thanks for your help, John
  Reply With Quote

Old   January 29, 2002, 15:03
Default Re: Wall functions
  #2
Anup
Guest
 
Posts: n/a
Hi, You may want to refer the book, "Turbulence Modeling for CFD" by David C. Wilcox. The book has an accompanyung software that you may find useful.

Abup.
  Reply With Quote

Old   January 29, 2002, 17:19
Default Re: Wall functions
  #3
Steve
Guest
 
Posts: n/a
Wall functions are used because some turbulence variables such as epsilon or omega are singular or unknown at a wall boundary. Consequently it is difficult to prescribe a boundary condition at the wall and integrating through the sublayer is computationally expensive.

Wall functions are really boundary conditons placed at the first grid point out from the wall (so we are not really modelling the sublayer and buffer layer, we are omitting the near-wall). They use the log-law of the wall to estimate what the velocities and turbulent quantities should be at the first grid point. Therefore to obtain sensical solutions the first grid point should be in the log-law region. So wall functions do not replace the log-law region they are merely an off-wall boundary condition. Ideally the first grid point should be placed as close to the lower bound of the log-law region so that we have as many points in the boundary layer as possible. Also the first grid space should be large to place the first grid point in the log-layer, with subsequent spacings much smaller to increase resolution of the boundary layer.

Beware, the log-law of the wall only holds for simple attached flows - as soon as the pressure gradient becomes significantly adverse, or separated they fail to provide accurate solutions. I never use them.
  Reply With Quote

Old   January 30, 2002, 04:31
Default Re: Wall functions
  #4
sylvain
Guest
 
Posts: n/a
To answer to your question, the post from Steve is quiet good. The aim of my post is just to correct a point in yours.

In the viscous sublayer the flow is NOT laminar, the fact is that turbulence in the boundary layer takes its origin from this region. But, roughtly speaking, almost all the "turbulent eddies" are align with the wall. So the effectif turbulence - perpendicular to the wall - is almost zero, then from a RANS point de vue, the flow looks laminar.

Regards,

Sylvain

  Reply With Quote

Old   January 30, 2002, 06:03
Default Re: Wall functions
  #5
John Grimm
Guest
 
Posts: n/a
Thanks everyone for responding to my messages, it is really appreciated. I just require one point of clarification (I am now assuming that wall functions only operate in the FIRST cell from the fixed boundary). In your message Steve you mentioned that wall fuctions use the log law relationship to estimate velocity and turbulence quantities. You also say that they should be close to the lower bound of the log law layer. If this is the case:

1) Why use wall functions at the lower end of the log law relationship (which they can model) when trying to model the viscous sublayer and buffer layer (which they cannot)? Is there not a better relationship and can the k-e turbulence model correctly model the log law layer without modification?

2) If wall functions can correctly model the log law layer why not save computational effort and place them towards the end of this layer?
  Reply With Quote

Old   January 30, 2002, 07:30
Default Re: Wall functions
  #6
sylvain
Guest
 
Posts: n/a
k-e model is drawing to correctly resolved the log law with wall function - i.e y+ in [30:150]. If you want to model the viscous sublayer and the buffer layer, you have to damp the turbulent viscosity to get the correct behavior. More over, for y+ greater than 10, eps behaves like 1/y, so to the mesh has to be very fine in this region.

To answer to the second question : if we know - without curvature of the wall nor adverse pressure gradient - where the log law region begins (y+ > 30), it is very difficult to predict where it stops. Moreover, computed values of u,k and eps at y+=150 are not egal to there BC values at same y+.

Regards,

Sylvain
  Reply With Quote

Old   January 30, 2002, 15:17
Default Re: Wall functions
  #7
Steve
Guest
 
Posts: n/a
One reason for having the first grid point close to the lower bound is that you do need some resolution of the boundary layer, i.e., you don't want the first point at y+=150, 2nd at Y+=300 etc. The log law is only reliable for a small range of flows so it is best to integrate through as much of the log-layer as possible.

Steve
  Reply With Quote

Old   January 31, 2002, 16:57
Default Re: Wall functions
  #8
pepe
Guest
 
Posts: n/a
I was just wondering what you do use, steve, instead of wall profiles?
  Reply With Quote

Old   January 31, 2002, 17:06
Default Re: Wall functions
  #9
Steve
Guest
 
Posts: n/a
Omega based models resolved down to the wall, i.e, with wall boundary conditions (hydraulically smooth conditions for Omega) and a grid with y+<1. The flows I deal with always involve pressure gradient induced separation so this approach is neccessary.
  Reply With Quote

Old   January 31, 2002, 17:48
Default Re: Wall functions
  #10
uwe schrottke
Guest
 
Posts: n/a
Hello Steve,

what is an "omega based model"

Uwe
  Reply With Quote

Old   January 31, 2002, 17:54
Default Re: Wall functions
  #11
Steve
Guest
 
Posts: n/a
By Omega based models I ment models that use omega (or turbulent frequency) as the second tubulence variable. Examples of these are Wilcox k-omega, SST, BSL, Stress-Omega, EASM k-omega.
  Reply With Quote

Old   February 1, 2002, 02:56
Default Re: Wall functions
  #12
kim
Guest
 
Posts: n/a
What is the diffrence in building grid for k-epsilon and k-omega ? Is it much finer in wall region for k-omega? Is k-omega superior over k-epsilon in flows with adverse pressure gradient (b.l. separation)?
  Reply With Quote

Old   February 1, 2002, 13:33
Default Re: Wall functions
  #13
Jongdae Kim
Guest
 
Posts: n/a
The reason of using wall functions is to decrease degree of freedoms of the simulation as we discussed above.

But more important fact is the physical phenomena of the flow close to the wall.

Significant amount of turbulent production comes from around wall (related to the fact that Sylvain mentioned). Without resolving this region appropriately, we can not get resonable results from the trubulent flow analysis.

Kim.
  Reply With Quote

Old   February 1, 2002, 19:34
Default Re: Wall functions
  #14
Steve
Guest
 
Posts: n/a
Yes, a finer grid is neccessary when using k-omega (unless it is a Low Reynolds number k-e model, i.e. integrates down to the wall). K-omega is superior for Adverse Pressure gradient flows as k-epsilon typically will delay separation due to an over-estimate of the length scale in the defect layer region. SST is better again for APGs and does not have the freestream dependence of k-omega.
  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
Modelling Turbulent Wall Functions victorwfreire CFX 11 January 14, 2015 16:48
LES with wall functions, wall boundary, Re 10000 subba CFX 4 July 10, 2007 06:04
Wall Functions pierre OpenFOAM Running, Solving & CFD 0 October 1, 2005 14:13
Wall functions - iterative approach agg Main CFD Forum 1 May 13, 2005 13:25
LES wall functions barrat Main CFD Forum 10 September 1, 2004 09:29


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