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

k-w Turbulence model implementation

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

LinkBack Thread Tools Search this Thread Display Modes
Old   November 20, 2005, 11:29
Default k-w Turbulence model implementation
Posts: n/a
Dear Friends,

I am working on the implementation of k-w turbulence model for 2d flows and I have a few problems, which I would like to clear. The turbulence model of Wilcox(1988) has been implemented by me in a coupled fasion, with implicit time stepping.

1. The code blows off after reaching a residue of 1e-4 on fine meshes.

2. k and w are becoming negative.

There are a few references that talk about limiting k and w, but the procedure is not well explained. I believe that the blow up of the code and the simultaneous degradation in the contours is because of improper limiting. It would be a big help if anyone could clarfiy these doubts and/or point to references that contain substantial discussions in this regard.

Thanks in advance, Suneesh
  Reply With Quote

Old   November 21, 2005, 03:15
Default Re: k-w Turbulence model implementation
Posts: n/a
normally the k starts attaining -ve value first.then obviously solution blows up. adopt linearisation of k -eqn. also ensure attainable convergence of the continuity eqn.
  Reply With Quote

Old   November 23, 2005, 11:12
Default Re: k-w Turbulence model implementation
Posts: n/a
There are a lot of lower limits on K and omega to prevent them below unphysical value. For example, From, CFL3d imposes the below limit on omega and K: If either nondimensional omega or non-dimensional K drops equal to or below 1.0e-20, they will be reset to 1.0e-20. From, WIND-US software sets kmin and εmin to 10'10 and 10'12, respectively. You can also find a more complicate formula to limit omega and K from Which way is the best rests on the problem you solve
  Reply With Quote

Old   November 23, 2005, 12:19
Default Re: k-w Turbulence model implementation
Posts: n/a
Dear Peter Zhao and pkm,

Thanks for the suggestions. I shall check them out. I would also like to know if some upper limits are provided for k and w, and also like to know if there are limits on the eddy viscosity.

Thanks in advance,

Regards, Suneesh
  Reply With Quote

Old   November 23, 2005, 18:35
Default experimetal or numerical data channel flow
Posts: n/a
Dear Friends

I am working on the simulation of turbulent mixed convection in vertical channel flow , I have validated my code for forced convection and the natural convection in a tall and square cavity, I want to validate my code for mixed o natural convection in the vertical channel flow, Please anyone can send to me experimental or numerical data for the validation code.

sincerely thank a lot

  Reply With Quote


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
Superlinear speedup in OpenFOAM 13 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 06:36
Centrifugal Pump and Turbulence Model Michiel CFX 12 January 25, 2010 04:20
turbulence model equation Andy Chen FLOW-3D 4 January 1, 2010 22:45
SA Turbulence model implementation ganesh Main CFD Forum 0 March 6, 2006 13:23
A reference on implementation of Spalart-Alam. Turbulence Model? Mohammad Kermani Main CFD Forum 2 December 26, 1999 03:56

All times are GMT -4. The time now is 05:07.