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

in k-epsilon wall function approach high Re turbulence models: question of velocity

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

Like Tree9Likes
  • 6 Post By romant
  • 3 Post By romant

Reply
 
LinkBack Thread Tools Display Modes
Old   June 27, 2013, 04:04
Default in k-epsilon wall function approach high Re turbulence models: question of velocity
  #1
Senior Member
 
romant's Avatar
 
Roman Thiele
Join Date: Aug 2009
Location: Stockholm, Sweden
Posts: 359
Rep Power: 11
romant is on a distinguished road
Hej,

I looked through the different wall function for k-epsilon turbulence models, however I feel I am missing something. In Wilcox (Turbulence modeling for CFD 2006) as well as many other publications, in wall functions the velocity seems to be set in the first grid point. In OpenFOAM, the following properties are changed with wall functions in the example of the k-epsilon turbulence model,

  1. epsilon: based on the logarithmic law of the wall, epsilon is set in the first cell node (epsilonWallFunction)
  2. G: the turbulence production term is determined from the law of the wall and also set in the first node from the wall (epsilonWallFunction)
  3. k: is only set to zeroGradient (kqRWallFunction)
  4. velocity is not touched in any of these functions?
  5. the turbulent viscosity (nut, mut) are set based on the law of the wall (mutUWallFunction)
The question is, why?
__________________
~roman
romant is offline   Reply With Quote

Old   July 2, 2013, 05:18
Default took a while but here is a solution to the problem
  #2
Senior Member
 
romant's Avatar
 
Roman Thiele
Join Date: Aug 2009
Location: Stockholm, Sweden
Posts: 359
Rep Power: 11
romant is on a distinguished road
Hej,

I think I found the solution to this. I found it in Pope, Turbulent Flow, 2000 under the way that wall functions are created. The original question was why the velocity is not set by any of the wall function. The answer is, that wall functions are apparently more robust by setting the shear stress and not the velocity, which makes sense when looking at separation points and other nasty flow phenomena, in which the velocity can become essentially 0. The way this is for the k-epsilon model done is

nutkWallFunction
  1. determine the nominal value for the shear velocity, u_shear
  2. find y+ using the shear velocity
  3. find the nominal values for the velocity in the first grid point (based on the log law in this case), U_p,nominal
  4. set the shear stresses -<u_i u_j> (not the normal stresses in the Reynolds stresses) to u_shear * U_p,real / U_p,nominal , In OF this happens through the turbulent viscosity (nu_t)
__________________
~roman
romant is offline   Reply With Quote

Old   July 2, 2013, 13:22
Default
  #3
Senior Member
 
cfdonline2mohsen's Avatar
 
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 117
Rep Power: 7
cfdonline2mohsen is on a distinguished road
Thanks roman
it was very helpful.
__________________
“If you have an apple and I have an apple and we exchange these apples then you and I will still each have one apple. But if you have an idea and I have an idea and we exchange these ideas, then each of us will have two ideas.”
cfdonline2mohsen is offline   Reply With Quote

Old   July 2, 2013, 19:02
Default
  #4
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
hi
Whould you explain more and Why stage 4 is done?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   July 3, 2013, 04:26
Default
  #5
Senior Member
 
romant's Avatar
 
Roman Thiele
Join Date: Aug 2009
Location: Stockholm, Sweden
Posts: 359
Rep Power: 11
romant is on a distinguished road
Quote:
Originally Posted by immortality View Post
hi
Whould you explain more and Why stage 4 is done?
Step 4 is introduced in order to have a more stable solver. Instead of setting the velocity in the first point away from the wall, a force (the wall shear force) is introduced in the equation. This way, the solver won't have any 0 velocities introduced by accident (think of stagnation points), this could destabilize the solver as you would be introducing 0's into the matrix.

The easiest way to introduce this source term into the equation is by modifying the turbulent viscosity, see Pope, Turbulent Flows, 2000, pp 442-443 and Versteeg and Malalasekera, An introduction to computational fluid dynamics: The finite volume method, pp. 276-277
__________________
~roman
romant is offline   Reply With Quote

Old   March 15, 2015, 22:12
Default
  #6
Senior Member
 
Huang Xianbei
Join Date: Sep 2013
Location: CAU,China
Posts: 266
Rep Power: 4
huangxianbei is on a distinguished road
Quote:
Originally Posted by romant View Post
Hej,

  1. set the shear stresses -<u_i u_j> (not the normal stresses in the Reynolds stresses) to u_shear * U_p,real / U_p,nominal , In OF this happens through the turbulent viscosity (nu_t)
Hi,romant:
Thank you for your explain. I have a question about the wall function for non-linear models. As the non-linear model can't be depicted by a simple turbulent viscosity, how can we set the shear stress?
say we obtained the shear stress as in your form: u_shear * U_p,real / U_p,nominal, so we should set the Reynolds shear stress by u_shear * U_p,real / U_p,nominal - nu*dU/dy, here, y represents the wall-normal direction. However, the non-linear Reynolds stress is not symmetric, it has 6 shear stress components, which one should be modified? I'm confused by this problem for some days. I'm looking forward for your reply.

Thank you

Xianbei
huangxianbei is offline   Reply With Quote

Reply

Tags
turbulence model, velocity, 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
mesh airfoil NACA0012 anand_30 OpenFOAM Meshing & Mesh Conversion 12 December 12, 2011 05:16
Discussion: Reason of Turbulence!! Wen Long Main CFD Forum 3 May 15, 2009 09:52
Variables Definition in CFX Solver 5.6 R P CFX 2 October 26, 2004 02:13
Quick Question - Wall Function D.Tandra Main CFD Forum 2 March 16, 2004 05:29


All times are GMT -4. The time now is 17:51.