# Problem in k-eps turbulence code in complex geometry

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

 December 24, 2010, 00: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, 17:10 #2 Senior Member   N/A Join Date: Mar 2009 Posts: 189 Rep Power: 9 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, 01: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.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post bwg FLUENT 1 January 13, 2010 14:09 massimo CD-adapco 0 January 14, 2003 05:58 Jiaying Xu CD-adapco 0 August 18, 2001 10:57 Frank Muldoon Main CFD Forum 1 January 5, 1999 11:09 Lans Main CFD Forum 13 October 27, 1998 11:20

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