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

Why y+ < 1

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

Like Tree2Likes
  • 1 Post By LuckyTran
  • 1 Post By LuckyTran

Reply
 
LinkBack Thread Tools Display Modes
Old   January 17, 2013, 11:15
Default Why y+ < 1
  #1
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
Hi all,

I allways wondered about why - for resolving the viscous sublayer - the threshold for y+ is supposed to be "1". Why not anywhere in the linear region, thus "<5" ?
As far as I understood it is needed to get correct boundary conditions for turbulent quantities (such as k) and for the velocity gradient at the wall.
Now, since the viscous layer is linear, the velocity gradient would be the same for y+ = 2. Which of the other values would be incorrect for higher y+ ?

Thanks in advance!
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   January 17, 2013, 12:03
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,662
Rep Power: 23
FMDenaro will become famous soon enough
Quote:
Originally Posted by RodriguezFatz View Post
Hi all,

I allways wondered about why - for resolving the viscous sublayer - the threshold for y+ is supposed to be "1". Why not anywhere in the linear region, thus "<5" ?
As far as I understood it is needed to get correct boundary conditions for turbulent quantities (such as k) and for the velocity gradient at the wall.
Now, since the viscous layer is linear, the velocity gradient would be the same for y+ = 2. Which of the other values would be incorrect for higher y+ ?

Thanks in advance!

The key, is to see y+ as a local Reynold number based on the computational grid length. y+ = O(1) means that you are able to solve in a region where the diffusive (and dissipative) effetcs are of the same order of the convective ones. That means also that a fluid particle the is advected in the computational cell has the same characteristic diffusive time.
FMDenaro is offline   Reply With Quote

Old   January 17, 2013, 12:37
Default
  #3
New Member
 
Guan
Join Date: Dec 2011
Location: US
Posts: 24
Rep Power: 5
foolboy007 is on a distinguished road
the fluctuation velocity close to the wall is not linear. it is proportional to y+^2
foolboy007 is offline   Reply With Quote

Old   January 17, 2013, 14:47
Default
  #4
Senior Member
 
Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 696
Rep Power: 13
LuckyTran will become famous soon enough
Quote:
Originally Posted by RodriguezFatz View Post
Hi all,

I allways wondered about why - for resolving the viscous sublayer - the threshold for y+ is supposed to be "1". Why not anywhere in the linear region, thus "<5" ?
As far as I understood it is needed to get correct boundary conditions for turbulent quantities (such as k) and for the velocity gradient at the wall.
Now, since the viscous layer is linear, the velocity gradient would be the same for y+ = 2. Which of the other values would be incorrect for higher y+ ?

Thanks in advance!
Your logic is pretty much spot on. y+ ~ 1 is a general rule (a very conservative one). You really only need y+<5 to apply the laws in the laminar sublayer. The additional restrictions come from calculation of gradients near the wall (where you now need a few points in the linear region). So if you wanted say 2-3 points in the linear region to resolve the gradients there then you would need correspondingly y+ ~ 1.5 - 2.5. If you are not worried about near wall gradients then you can use y+ ~ 5 and be okay.

Actually, it is much more important to have sufficient cells across the entire inner layer to resolve the variables changes within the inner layer. You can actually achieve solutions superior to y+<0.5 but not enough cells in inner layer, using y+ ~ 5 and with enough cells to completely resolve the inner layer. This is actually the mesh criteria for LES and DNS style grids where y+ ~ 3 is considered sufficient as long as the most important regions are resolved.
LuckyTran is offline   Reply With Quote

Old   January 18, 2013, 06:40
Default
  #5
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
Thank you all for your replies!


Quote:
Originally Posted by FMDenaro View Post
The key, is to see y+ as a local Reynold number based on the computational grid length. y+ = O(1) means that you are able to solve in a region where the diffusive (and dissipative) effetcs are of the same order of the convective ones. That means also that a fluid particle the is advected in the computational cell has the same characteristic diffusive time.
But all this is also the case for y+ = 2.


Quote:
Originally Posted by foolboy007 View Post
the fluctuation velocity close to the wall is not linear. it is proportional to y+^2
I should have said that I mean resolved sheaths of RANS turbulence models.


Quote:
Originally Posted by LuckyTran View Post
Your logic is pretty much spot on. y+ ~ 1 is a general rule (a very conservative one). You really only need y+<5 to apply the laws in the laminar sublayer. The additional restrictions come from calculation of gradients near the wall (where you now need a few points in the linear region). So if you wanted say 2-3 points in the linear region to resolve the gradients there then you would need correspondingly y+ ~ 1.5 - 2.5. If you are not worried about near wall gradients then you can use y+ ~ 5 and be okay.
What are these gradients needed for?
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   January 18, 2013, 06:45
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,662
Rep Power: 23
FMDenaro will become famous soon enough
Actually, I wrote O(1) ... resolving with some grid point below y+=1 is a sort of guarantee that you are able to solve the boundary layer profile taking into account for an accurate evaluation of the wall stress.

Howevewr, in our group LESinItaly, we did several comparison for LES using an unresolved grid and prescribing the no-slip condition. Apart the wall coefficient, the statistics in the inner part of the flow were acceptable
FMDenaro is offline   Reply With Quote

Old   January 18, 2013, 06:56
Default
  #7
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
Quote:
Originally Posted by FMDenaro View Post
Actually, I wrote O(1) ... resolving with some grid point below y+=1 is a sort of guarantee that you are able to solve the boundary layer profile taking into account for an accurate evaluation of the wall stress.

Howevewr, in our group LESinItaly, we did several comparison for LES using an unresolved grid and prescribing the no-slip condition. Apart the wall coefficient, the statistics in the inner part of the flow were acceptable
Ok, I got you. Now I am aware of getting higher accuracy for higher number of gridpoints. My initial question was more like if I am doing complete nonsense with y+ = 3, such as using a log-wall function in the buffer layer that actually rapes the equations.
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   January 18, 2013, 07:02
Default
  #8
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,662
Rep Power: 23
FMDenaro will become famous soon enough
Quote:
Originally Posted by RodriguezFatz View Post
Ok, I got you. Now I am aware of getting higher accuracy for higher number of gridpoints. My initial question was more like if I am doing complete nonsense with y+ = 3, such as using a log-wall function in the buffer layer that actually rapes the equations.
I think that the key is to forget for a while the theoretical results for turbulence over flat plate/channel/pipe .... in practical applications is so rare to have some boundary that can be "assumed" as a theoretical wall ... From my experience, linear law, log law, etc. are useful to understand some features of the flow but I doubt that can be used to prescribe general boundary condition when your first grid point close to wall exceeds y+=1, expecially in LES. Maybe in RANS this approach can be better justified
FMDenaro is offline   Reply With Quote

Old   January 18, 2013, 08:54
Default
  #9
Senior Member
 
Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 696
Rep Power: 13
LuckyTran will become famous soon enough
Quote:
Originally Posted by RodriguezFatz View Post
I should have said that I mean resolved sheaths of RANS turbulence models.

What are these gradients needed for?
As you said, the velocity gradients need to be resolved in order to get the proper shear stress. The shear stress is important for both the laminar (molecular diffusion) component of the shear stress. If you are doing RANS (anything that uses Boussinesq hypothesis) then the same shear stress is needed to calculate the turbulent diffusion (for example, k equation in k-epsilon or k-omega).

To properly calculate the gradient at solid boundaries, you need sufficient data points there so that the discretization error is not large for gradients (basically finite difference approximations to derivatives).

If on top of fluid flow you are also trying to simulate convective heat transfer (or mass transfer). You need to have finer grids because the thermal & concentration profile can only be properly resolved only after the hydrodynamic profile has already been resolved. i.e. the thermal & concentration gradients are guaranteed to be always less accurate than the velocity gradient.

I'm glad we were able to have this discussion as I find people nowadays that blindly accept the y+~ 1 rule, and then apply a huge stretch factor to their grids and then claiming that they have "resolved the viscous sublayer"
Kent_Neo likes this.
LuckyTran is offline   Reply With Quote

Old   January 18, 2013, 08:57
Default
  #10
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
Quote:
Originally Posted by LuckyTran View Post
To properly calculate the gradient at solid boundaries, you need sufficient data points there so that the discretization error is not large for gradients (in terms of finite difference approximations to derivatives).
What I was trying to say: When the profile is linear, it does not matter if I have y+=1 or y+=4, because due to linearity it's always the same gradient!
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   January 18, 2013, 09:14
Default
  #11
Senior Member
 
Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 696
Rep Power: 13
LuckyTran will become famous soon enough
Quote:
Originally Posted by RodriguezFatz View Post
What I was trying to say: When the profile is linear, it does not matter if I have y+=1 or y+=4, because due to linearity it's always the same gradient!
Oh now I see your point. If you know the y+ then you know what the gradient is. But the problem is you don't know the y+.

It is linear yes but it is an unknown linear gradient. That is because of the unknown friction velocity (recall the formulas for u+ and y+). Hence although you know the gradients to be linear, you do not know the value of the gradient.

Put another way, you don't know the value of y+ until you know the friction velocity, which you do not know until you have solved the problem. Until then, you must solve for y+ (solve for friction velocity). To solve for friction velocity, you must somehow resolve the shear stresses in that region. Hence your grid requirements are motivated by the need to solve for these gradients.
Kent_Neo likes this.
LuckyTran is offline   Reply With Quote

Old   January 18, 2013, 09:16
Default
  #12
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,662
Rep Power: 23
FMDenaro will become famous soon enough
Quote:
Originally Posted by RodriguezFatz View Post
What I was trying to say: When the profile is linear, it does not matter if I have y+=1 or y+=4, because due to linearity it's always the same gradient!

again, what you addressed is a result of a statistical theory... then, if you solve for the statistical variable, as in RANS, what you say can be somehow right. But for other formulation you do not solve for the statistically averaged velocity...
FMDenaro is offline   Reply With Quote

Old   January 18, 2013, 09:30
Default
  #13
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
Quote:
Originally Posted by FMDenaro View Post
again, what you addressed is a result of a statistical theory... then, if you solve for the statistical variable, as in RANS, what you say can be somehow right. But for other formulation you do not solve for the statistically averaged velocity...
Yes, that is what I ment.
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   January 18, 2013, 09:40
Default
  #14
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
Quote:
Originally Posted by LuckyTran View Post
Oh now I see your point. If you know the y+ then you know what the gradient is. But the problem is you don't know the y+.

It is linear yes but it is an unknown linear gradient. That is because of the unknown friction velocity (recall the formulas for u+ and y+). Hence although you know the gradients to be linear, you do not know the value of the gradient.

Put another way, you don't know the value of y+ until you know the friction velocity, which you do not know until you have solved the problem. Until then, you must solve for y+ (solve for friction velocity). To solve for friction velocity, you must somehow resolve the shear stresses in that region. Hence your grid requirements are motivated by the need to solve for these gradients.
It seems that we talk at cross-purposes. Maybe it's due to my language problems.
Imagine a standard universal wall velocity profile. Now you solve the equations numerically on a grid. Somewhere in your discretized equations the friction velocity appears. For the velocity gradient you will use the first cell's velocity divided by the distance to the wall. This is a pretty good assumption, since you know from the books, that the velocity profile is linear. Now, what I was trying to say is, that this velocity gradient will allways have - exaclty - the same value, whether it is calculated at y+=1, 2,3,4,5. Since the velocity is a linear function there, the "sample rate" of your numerical grid is not important.
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   January 18, 2013, 09:43
Default
  #15
Senior Member
 
Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 696
Rep Power: 13
LuckyTran will become famous soon enough
Quote:
Originally Posted by RodriguezFatz View Post
For the velocity gradient you will use the first cell's velocity divided by the distance to the wall. This is a pretty good assumption, since you know from the books, that the velocity profile is linear.
But you do not know this velocity in the first place. And that is why you cannot compute the y+ and velocity profile. You can only arrive at this velocity by properly accounting for the fluid stress (by solving the momentum balance).
LuckyTran is offline   Reply With Quote

Old   January 18, 2013, 09:45
Default
  #16
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
I think since we use an iterative solver, we take the gradient of the last iteration and keep doing until convergence.
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   January 18, 2013, 10:22
Default
  #17
Senior Member
 
flotus1's Avatar
 
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,107
Rep Power: 19
flotus1 will become famous soon enoughflotus1 will become famous soon enough
Let me join the discussion with a new thought:
If the y+ of the first cell is 5, the extent of the next cells can only be equal or larger than 5 for a mesh with a reasonable expansion ratio.

While the velocity is a linear function in the region of y+<5, it becomes nonlinear beyond this point and may not be resolved adequately with cells larger than y+=5.
Thus a y+ of 1 is a conservative estimate to ensure that also the nonlinear region y+>5 is resolved reasonably.
flotus1 is offline   Reply With Quote

Old   January 18, 2013, 10:33
Default
  #18
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,103
Rep Power: 16
RodriguezFatz will become famous soon enough
Flotus, I totally agree. But this is "just" an accuracy problem. You don't abuse any model equation with complete nonsense by ignoring it, right?
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz 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
Y+ < 1 for Standard K-Epsilon Leepox FLUENT 1 May 9, 2012 14:40
REy < 200 fluentguy Main CFD Forum 1 May 5, 2009 23:58
Fluent : Hydraulic Diameter < 1mm ok??? jerome-gmc FLUENT 0 July 4, 2005 05:35
< 20 species problem Tomasz Ochrymiuk FLUENT 1 September 14, 2000 09:47
FVM < FEM ? Tony Main CFD Forum 15 May 18, 1999 14:03


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