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

y+=1 vs Wall Function

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

Like Tree12Likes
  • 1 Post By LuckyTran
  • 2 Post By Far
  • 2 Post By Far
  • 1 Post By LuckyTran
  • 2 Post By Far
  • 2 Post By LuckyTran
  • 2 Post By Far

Reply
 
LinkBack Thread Tools Display Modes
Old   November 16, 2012, 03:34
Default y+=1 vs Wall Function
  #1
New Member
 
Join Date: Nov 2012
Posts: 24
Rep Power: 4
paduchev is on a distinguished road
Does anybody knows how fluent chooses when to use wall function and when to model boundary layer directly when I'm using turbulenxe model K-W-SST?
paduchev is offline   Reply With Quote

Old   November 16, 2012, 06:38
Default
  #2
Senior Member
 
Shamoon Jamshed
Join Date: Apr 2009
Location: Karachi
Posts: 114
Rep Power: 8
Shamoon Jamshed is on a distinguished road
Send a message via Skype™ to Shamoon Jamshed
There is no separate option to use wall f/n in kw SST case. Just keep near wall cells fine enough to get y+=1 in FLUENT.
Shamoon Jamshed is offline   Reply With Quote

Old   November 16, 2012, 13:57
Default
  #3
Senior Member
 
Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 594
Rep Power: 11
LuckyTran is on a distinguished road
Quote:
Originally Posted by paduchev View Post
Does anybody knows how fluent chooses when to use wall function and when to model boundary layer directly when I'm using turbulenxe model K-W-SST?
Are you asking about any equation in particular? Wall functions for velocity, the k-equation or w-equation? The velocity boundary condition uses the identical enhance wall treatment found on the realizable-k-epislon model if you are familiar with that. The k-equation uses the identical wall boundary condition from the k-epsilon model. There is a blended-two-layer wall function for the w equation as well (similar blending to get the enhanced wall treatment).

Also, a bit of a technicality, but there are always wall functions used to model the wall boundary condition. Hence, it's not good practice to say wall function and modelling the boundary layer as you are saying.

Quote:
Originally Posted by Shamoon Jamshed View Post
There is no separate option to use wall f/n in kw SST case. Just keep near wall cells fine enough to get y+=1 in FLUENT.
There is no option but it is there in the background. In other words, you don't have a choice but the wall functions are always there. You can use the kw SST model on a coarse grid and it will automatically "switch" to the log-law. There actually is no switch, the blended wall function is just a smoothed version of linear and log laws.
joy2000 likes this.
LuckyTran is offline   Reply With Quote

Old   November 18, 2012, 01:19
Default
  #4
New Member
 
Join Date: Nov 2012
Posts: 24
Rep Power: 4
paduchev is on a distinguished road
Quote:
Originally Posted by LuckyTran View Post
Also, a bit of a technicality, but there are always wall functions used to model the wall boundary condition. Hence, it's not good practice to say wall function and modelling the boundary layer as you are saying.

There is no option but it is there in the background. In other words, you don't have a choice but the wall functions are always there. You can use the kw SST model on a coarse grid and it will automatically "switch" to the log-law. There actually is no switch, the blended wall function is just a smoothed version of linear and log laws.
Thanks very much. Have I understood correctly - fall function is always using, when I'm use k-w SST, but if i have fine near-wall mesh, boundary layer will be modeled by logarithmic law and if i will have coarse mesh the smoothed version of logarithmic law (wall function) will be using? But why in fluents tutorial posted: (13.2.1.8.*Grid Resolution RANS Models) "However, for simulations with high accuracy demands on the wall boundary layer (especially for heat transfer predictions) near wall meshes with y+*~1 are recommended. When wall functions are used, it is essential to avoid meshes with * values lower than y+ ~30 as the wall shear stress and the wall heat transfer can and will seriously deteriorate under such conditions."
Does it means that when i'm using k-w-SST model i must avoid meshes with y+ from 2 to 30?
And the last. I have modeled experimental data - flow around conical nozzle and the "wall function mesh" y+ near 60 showed better match with test data, than mesh with y+ = 1
paduchev is offline   Reply With Quote

Old   November 18, 2012, 11:57
Default
  #5
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,906
Blog Entries: 6
Rep Power: 38
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
In older versions there was switch which changes the formulation from integration to wall (ITW) approach to wall function (WF) at Y+ (Rather Y*) = 11.06. Below 11.06 ITW was used and above 11.06 WF was used irrespective the fact that WF were not valid in the buffer zone (from 5-30) and similarly ITW was used irrespective the fact the viscous sub layer requires Y+ = 1 or at-least Y +< 5.

Now Fluent is using the hybrid wall functions (similar to the apporach in reference give below and technology comes from CFX due to F. Menter) which is combination of ITW and WF. Upto Y+ = 6, it uses the ITW apporach and above Y+ > 30 it uses the Wall function (yes standard wall functions) and in the region between Y + 6-30 it uses the combination of both.

Important points to remember:

1. Keep Y+ below 10 for viscous dominate flows with separation.
2. Use Y+ < 1 for the transition prediction and heat transfer calculations.
3. Use Y+ < 10 for SST model (omega based models in general)
4. Use K-Epsilon model with scalable wall functions (default option in new version of Fluent). It will keep the Y+ above 11.06 so that it does not suffer the limitation of old standard wall functions. Moreover now you can carry the successive refinment of mesh, which was not possible with standard WFs.

http://num.math.uni-goettingen.de/ba...ings/knopp.pdf
Cloud and mrenergy like this.
Far is offline   Reply With Quote

Old   November 18, 2012, 12:01
Wink Matching the Experimental data is not the goal !
  #6
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,906
Blog Entries: 6
Rep Power: 38
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Quote:
And the last. I have modeled experimental data - flow around conical nozzle and the "wall function mesh" y+ near 60 showed better match with test data, than mesh with y+ = 1
First make the mesh refinement study. once you choose the fine enough mesh, then half the first cell height (in other words you are reducing Y+) and check the results. Do they change? or double the first cell height and compare the results. Again do they change?


Moreover CFD is not done to match the experimental data, it can be done in MS Excel . You must ensure the best practices of CFD and there are many source of errors e.g. turbulence modelling, grid, boundary conditions, flow scheme and it is also possible that some source of error is cancelling other factor and you are getting good results.
joy2000 and mrenergy like this.
Far is offline   Reply With Quote

Old   November 18, 2012, 22:30
Default
  #7
New Member
 
Join Date: Nov 2012
Posts: 24
Rep Power: 4
paduchev is on a distinguished road
Quote:
Originally Posted by Far View Post
3. Use Y+ < 10 for SST model (omega based models in general)
It's turbulence models requirements or Fluent's?
Why I could not use k-omega SST model with wall-function mesh with y+>30?
paduchev is offline   Reply With Quote

Old   November 18, 2012, 22:32
Default
  #8
Senior Member
 
Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 594
Rep Power: 11
LuckyTran is on a distinguished road
Quote:
Originally Posted by paduchev View Post
It's turbulence models requirements or Fluent's?
Why I could not use k-omega SST model with wall-function mesh with y+>30?
You can use a larger mesh with y+ > 30 even with k-w SST.
joy2000 likes this.
LuckyTran is offline   Reply With Quote

Old   November 19, 2012, 02:57
Default
  #9
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,906
Blog Entries: 6
Rep Power: 38
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Quote:
Originally Posted by paduchev View Post
It's turbulence models requirements or Fluent's?
Why I could not use k-omega SST model with wall-function mesh with y+>30?
strictly speaking, SST model is low Reynolds number model. Which means you need the Y+ ~ 1. But in commercial codes (as they know that it is not possible to get the Y+ 1 for every problem by every one), therefore they have implemented WF. So you can use the WF and it is automatically turned on when you have Y+ > 30 (Which is not always possible to ensure due to presence of separation, stagnation zone and your yplus can go to zero!) .

You must also understand that in the log region, you are solving the transformed K-Epsilon model and it is good idea to use the original K-epsilon (which is designed for WFs) in combination with scalable wall functions for the reasons mentioned above and in previous post.

In my limited experience K-epsilon model with scalable wall function is better choice than the SST model with hybrid functions for the Y+> 30.
joy2000 and mrenergy like this.
Far is offline   Reply With Quote

Old   November 19, 2012, 03:15
Default
  #10
New Member
 
Join Date: Nov 2012
Posts: 24
Rep Power: 4
paduchev is on a distinguished road
Quote:
Originally Posted by Far View Post
strictly speaking, SST model is low Reynolds number model. Which means you need the Y+ ~ 1. But in commercial codes (as they know that it is not possible to get the Y+ 1 for every problem by every one), therefore they have implemented WF. So you can use the WF and it is automatically turned on when you have Y+ > 30 (Which is not always possible to ensure due to presence of separation, stagnation zone and your yplus can go to zero!) .

You must also understand that in the log region, you are solving the transformed K-Epsilon model and it is good idea to use the original K-epsilon (which is designed for WFs) in combination with scalable wall functions for the reasons mentioned above and in previous post.

In my limited experience K-epsilon model with scalable wall function is better choice than the SST model with hybrid functions for the Y+> 30.

Thanks a lot, I have writed you a privat message, could you answer on it?)
paduchev is offline   Reply With Quote

Old   November 19, 2012, 03:19
Default
  #11
Senior Member
 
Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 594
Rep Power: 11
LuckyTran is on a distinguished road
Quote:
Originally Posted by Far View Post
strictly speaking, SST model is low Reynolds number model. Which means you need the Y+ ~ 1.
I think this is mis-leading. Strictly speaking, none of the commonly available turbulence models are restricted to any Re. The SST formulation is already a hybrid blending of a the k-epsilon and k-omega models and is not restricted to low or high Reynolds number grids. What is restricted is the wall modelling approach for the wall boundary conditions which can be high or low Re and still be consistent with the turbulence model being used.

If one desired (and it seems paduchev definitely wants to) use a grid with y+ of wall adjacent cell to be >30, one can certainly do so. Will a low Re number grid do better? Probably. But I don't think any of the previous explanations were of any help.
joy2000 and mrenergy like this.
LuckyTran is offline   Reply With Quote

Old   December 1, 2012, 10:08
Default
  #12
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,906
Blog Entries: 6
Rep Power: 38
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Please look at the document by Lars for details of SST model http://www.os-cfd.ru/cfd_docs/chalme...5ebe3f5385c13b

And don't confuse the hybrid wall treatment with the hybrid turbulence model like SST model which is combination of K-omega in inner and buffer layer. And K-epsilon model in outer and fully turbulent region.

The hybrid wall treatment is used to relax the Y+ requirements on the K-omega model in SST model. Which requires Y+<2

Moreover SST model was made by mentor in1994 and work on enhanced wall treatment (automatic wall treatment) was published in 2002-2003. Which means SST model was there and without hybrid wall functions. Now question arises what type of wall treatment was being used before the intorudction of hybrid wall functions. The answer is it is same as required by K-omega model. This is also clearly mentioned in CFX theory guide where this model was implemented first by mentor.

http://www.kxcad.net/ansys/ansys_cfx.../i1346068.html

Quote:
Automatic Near-Wall Treatment for Omega-Based Models

Automatic near-wall treatment automatically switches from wall-functions to a low-Re near wall formulation as the mesh is refined.

One of the well known deficiencies of the model is its inability to handle low turbulent Reynolds number computations. Complex damping functions can be added to the model, as well as the requirement of highly refined near-wall grid resolution (y+ < 0.2) in an attempt to model low turbulent Reynolds number flows. This approach often leads to numerical instability. Some of these difficulties may be avoided by using the model, making it more appropriate than the model for flows requiring high near-wall resolution (e.g., high wall heat transfer, transition). However, a strict low-Reynolds number implementation of the model would also require a near wall grid resolution of at least y+ < 2. This condition cannot be guaranteed in most applications at all walls. For this reason, a new near wall treatment was developed by ANSYS CFX for the based models that allows for a smooth shift from a low-Reynolds number form to a wall function formulation. This near wall boundary condition, named automatic near wall treatment in ANSYS CFX, is used as the default in all models based on the -equation (standard , Baseline , SST, -Reynolds Stress).
I am also attaching the some notes from Wilcox book on turbulence modelling.
Attached Images
File Type: jpg K-omegavs K-epsilon.jpg (97.5 KB, 61 views)
joy2000 and mrenergy like this.

Last edited by Far; December 1, 2012 at 10:24.
Far is offline   Reply With Quote

Old   December 2, 2012, 05:34
Default
  #13
Senior Member
 
Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 594
Rep Power: 11
LuckyTran is on a distinguished road
Quote:
Originally Posted by Far View Post
Please look at the document by Lars for details of SST model http://www.os-cfd.ru/cfd_docs/chalme...5ebe3f5385c13b

And don't confuse the hybrid wall treatment with the hybrid turbulence model like SST model which is combination of K-omega in inner and buffer layer. And K-epsilon model in outer and fully turbulent region.
Exactly. The SST approach for high-Re grid is equivalent to a standard k-epsilon model. Which means that the requirements should be the same as the standard k-epsilon model and not the k-omega model (for a high-Re grid). In the low-Re grid case I agree, the requirements for k-w SST should be treated as same as k-omega.
LuckyTran is offline   Reply With Quote

Old   December 2, 2012, 12:05
Default
  #14
Senior Member
 
François Grégoire
Join Date: Jan 2010
Location: Laval University, Canada
Posts: 366
Rep Power: 8
macfly is on a distinguished road
Here's the complete and possibly fresher version of Lars Davidson's document that Far is referring to in post #12, very nice : http://www.tfd.chalmers.se/~lada/pos...-modelling.pdf
macfly 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
channelFoam for a 3D pipe AlmostSurelyRob OpenFOAM 3 June 24, 2011 13:06
Wall function for velocity? johnblund OpenFOAM 0 March 10, 2011 09:50
Need some wall function approaches! yka8150 Main CFD Forum 0 September 21, 2009 23:08
Problem with compile the setParabolicInlet ivanyao OpenFOAM Running, Solving & CFD 6 September 5, 2008 20:50
Droplet Evaporation Christian Main CFD Forum 2 February 27, 2007 07:27


All times are GMT -4. The time now is 08:26.