# 2D Lid Driven Cavity Flow simulation using MATLAB

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

 February 23, 2010, 14:17 2D Lid Driven Cavity Flow simulation using MATLAB #1 New Member   Joseph Join Date: Feb 2010 Posts: 1 Rep Power: 0 I have been trying to simulate the 2D lid driven cavity flow using the SIMPLE algorithm on the staggered grid, in which velocity is assumed to be coupled with pressure and momentum equations are iterated and u, v and pressure p are corrected by a correction term u_corr, v_corr and p_corr respectively. The looping stops when p_corr decreases to a sufficiently small value. However, the flow obtain does not seem to look like how it should look like, i.e. a main vortex at the centre with counter vortices at the bottom corners. Instead, the plot I obtained shows the boundary conditions (Lid velocity at top, with no-slip condition at other walls), but the velocity arrows in the middle seem to be of a magnitude 100 times smaller than lid velocity and in random directions. It would appear that I have a problem with the correction terms, u_corr, v_corr and p_corr being too small, but i checked the discretized equations over and over and they seem to be correct.. Does anyone know what could be the likely reason for this? Appreciate any help or suggestions to modify my code!

June 14, 2010, 13:25
#2
New Member

Jame
Join Date: Jun 2010
Posts: 3
Rep Power: 7
I have been also facing similar kind of problem. In my case after few iterations oscillations start and the solution fails. However the same code worked really well for Coutte flow. Can ne1 help??
Quote:
 Originally Posted by josephlm I have been trying to simulate the 2D lid driven cavity flow using the SIMPLE algorithm on the staggered grid, in which velocity is assumed to be coupled with pressure and momentum equations are iterated and u, v and pressure p are corrected by a correction term u_corr, v_corr and p_corr respectively. The looping stops when p_corr decreases to a sufficiently small value. However, the flow obtain does not seem to look like how it should look like, i.e. a main vortex at the centre with counter vortices at the bottom corners. Instead, the plot I obtained shows the boundary conditions (Lid velocity at top, with no-slip condition at other walls), but the velocity arrows in the middle seem to be of a magnitude 100 times smaller than lid velocity and in random directions. It would appear that I have a problem with the correction terms, u_corr, v_corr and p_corr being too small, but i checked the discretized equations over and over and they seem to be correct.. Does anyone know what could be the likely reason for this? Appreciate any help or suggestions to modify my code!

June 23, 2010, 08:54
#3
New Member

zt
Join Date: Jun 2010
Posts: 1
Rep Power: 0
Quote:
 Originally Posted by josephlm I have been trying to simulate the 2D lid driven cavity flow using the SIMPLE algorithm on the staggered grid, in which velocity is assumed to be coupled with pressure and momentum equations are iterated and u, v and pressure p are corrected by a correction term u_corr, v_corr and p_corr respectively. The looping stops when p_corr decreases to a sufficiently small value. However, the flow obtain does not seem to look like how it should look like, i.e. a main vortex at the centre with counter vortices at the bottom corners. Instead, the plot I obtained shows the boundary conditions (Lid velocity at top, with no-slip condition at other walls), but the velocity arrows in the middle seem to be of a magnitude 100 times smaller than lid velocity and in random directions. It would appear that I have a problem with the correction terms, u_corr, v_corr and p_corr being too small, but i checked the discretized equations over and over and they seem to be correct.. Does anyone know what could be the likely reason for this? Appreciate any help or suggestions to modify my code!
Hi

Do you have to solve your system with SIMPLE algorithm? I mean you can use picard iteration which is also very simple. I just wrote a 2D Navier Stokes code for Matlab which iterates the solution of matrix K by using picard iterations. It gives reasonable results for the cavity problem.

 June 24, 2010, 01:58 #4 New Member   Dhileep Join Date: Apr 2010 Posts: 10 Rep Power: 7 Hi, i have written a matlab code for a 2D natural convection problem using malalasekara versteeg as reference. can you suggest me some reference books for lid driven cavity problem.

 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 quarkz Main CFD Forum 9 January 20, 2013 06:54 EmersonKB CFX 5 July 2, 2009 08:17 prasanthnitt Main CFD Forum 0 June 18, 2009 06:39 max91 CFX 1 July 29, 2008 20:28 ani Main CFD Forum 3 April 2, 2006 13:27

All times are GMT -4. The time now is 01:38.