# Problem in k-eps turbulence code in complex geometry

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

 December 23, 2010, 23:46 Problem in k-eps turbulence code in complex geometry #1 New Member   Shantanu Pramanik Join Date: Dec 2010 Location: Durgapur Posts: 2 Rep Power: 0 Dear everyone, I've written a code for solving turbulent flow problem using k-eps model in complex geometry (boundary-fitted coordinate system). The code when run as a laminar producing acceptable result. But if I run it (in lid-driven 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 lid-centre 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 high-Re model using wall functions and presently trying to solve the flow field only. The matrix is being solved using line-by-line TDMA technique and scheme is SIMPLE with staggered grids. Thanks in advance for helping me sorting out the problem.

 December 27, 2010, 16:10 #2 Senior Member   N/A Join Date: Mar 2009 Posts: 189 Rep Power: 16 Have you checked the literature to ensure that the model you use accurately predicts lid-driven cavity problem? Also ensure that the first grid point is in the log-law region.

 December 28, 2010, 00:44 k-eps #3 New Member   Shantanu Pramanik Join Date: Dec 2010 Location: Durgapur Posts: 2 Rep Power: 0 Yes, 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 lid-driven-inclined-cavity 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 wall 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.