Problem in keps turbulence code in complex geometry
Dear everyone,
I've written a code for solving turbulent flow problem using keps model in complex geometry (boundaryfitted coordinate system). The code when run as a laminar producing acceptable result. But if I run it (in liddriven cavity problem) for turbulent cases (Re=480000). The result is not correct. I took a rectangular cavity (angle of inclination =90deg) and for the mentioned Re, the program is generating a result what one would get in the middle of a of result that an unsteady code would generates (for the same problem) before it (unsteady code) converges. To be precise,the vortex after generating just below the lidcentre comes to the right wall and remained attached there just below the lid instead of going to the centre of the cavity. I have used highRe model using wall functions and presently trying to solve the flow field only. The matrix is being solved using linebyline TDMA technique and scheme is SIMPLE with staggered grids. Thanks in advance for helping me sorting out the problem. 
Have you checked the literature to ensure that the model you use accurately predicts liddriven cavity problem? Also ensure that the first grid point is in the loglaw region.

keps
http://dl.dropbox.com/u/14695550/Books/export.epsYes, the formulation working well for simple geometry both in laminar and turbulent cases. Using same scheme I've written the code for laminar case in complex geometry for liddriveninclinedcavity problem and that is working fine for square (angle=90deg) as well as inclined cavity problem. But when I changed it to solve turbulent flow problem it's not producing acceptable result. Although converging but result is wrong.
Y+ is kept in the log layer. The vortex is remained attached to the vertical wallhttp://dl.dropbox.com/u/14695550/Books/export.eps as in the attached fig. May be there is something wrong in the physics. http://dl.dropbox.com/u/14695550/Books/export.eps Thank you. 
All times are GMT 4. The time now is 15:45. 