CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

Turbulence Modelling

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

LinkBack Thread Tools Display Modes
Old   April 5, 2003, 04:57
Default Turbulence Modelling
Posts: n/a

I am at present working on implementation of high reynold's number K-epsilion turbulence model and k-omega turbulence model.

While refering to the book by Versteeg and Malalsekra, I have the wall functions for y+ between 30 and 500, but what should be the wall function for low values of y+ say around 1.

In case of k-omega, what should be the value of omega.

Also how should be the boundary conditions be implemented in case of both.

Thanks and Regards

  Reply With Quote

Old   April 5, 2003, 07:43
Default Re: Turbulence Modelling
Jonas Larsson
Posts: n/a
If you are going to use k-epsilon down to y+ = 1 you need to use a low-Re model. A good review of several classical low-Re k-epsilon models is:

Virendra C. Patel, Wolfgang Rodi, and Georg Scheuerer. Turbulence models for near-wall and low Reynolds number flows: A review. AIAA Journal, 23(9):1308-1319, 1985.

The k-omega model can be used down to y+ = 1 without any modification. Strictly speaking omega goes to infinity on the wall. However, I would recommend that you use Wilcox's sruface rougness model to set a wall omega (see Wilcox's book Turbulence Modeling for CFD). The surface rougness model sets:

omega_wall = friction_velocity * S_R / nu


S_R = (50/k_R+)^2 when k_R+ < 25

S_R = 100/K_R+ when k_R+ >= 25


k_R+ = friction_velocity * average_height_of_roughness / nu

If you set k_R+ to 5 (or lower) you will get a smooth wall.

This surface wall-model is much more stable than having omega go to infinity on the wall and it works very well also for smooth walls.
  Reply With Quote

Old   April 6, 2003, 03:21
Default Re: Turbulence Modelling
Posts: n/a
Hi, There are a lot of issues for the use of k-epsilon model and I have a report which I prepared after I wrote the code for k-epsilon model. If you want I can email it to you. About using the wall functions for y+ < 11.63, we do not use the equillibrium condition where G- density * epsilon = 0 is used. Instead we do this : Calculate tau_wall = mu_eff * up/yp calculate G = tau_wall * up/yp calculate epsilon = (2 * nu * Kp) /(yp^2)

About the boundary conditions implementation, I have explained in that report of mine , if you want ...

Abhijeet Vaidya IIT - Delhi
  Reply With Quote


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
Turbulence postprocessing Mohsin FLUENT 2 October 3, 2016 14:18
Discussion: Reason of Turbulence!! Wen Long Main CFD Forum 3 May 15, 2009 09:52
Turbulence modelling validation : Orifice flow Joe Main CFD Forum 2 July 26, 2006 16:41
Turbulence modelling Panos Siemens 6 July 8, 2004 02:40

All times are GMT -4. The time now is 06:29.