# Problem Encountered inLid Driven Cavity in Finite Volume Method with SIMPLE Algorithm

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

July 3, 2015, 19:47
Problem Encountered inLid Driven Cavity in Finite Volume Method with SIMPLE Algorithm
#1
New Member

Join Date: Jul 2015
Posts: 2
Rep Power: 0
I have written my code in MATLAB and I need help from the forum on where things could be going wrong, my streamline plot appears not to be properly predicted and I am also not getting vortices on the corners of the cavity.

Please any help with code modification will be appreciated.

Thank you.
Attached Images
 Streamline-RE100.jpg (58.7 KB, 28 views)
Attached Files
 Lid Driven Cavity Matlab Code.txt (26.2 KB, 44 views)

Last edited by Abdels83; July 4, 2015 at 10:44. Reason: I have updated the .m file and the code is now running. I also added the streamline plot for Reynolds number of 100.

 July 4, 2015, 07:55 #2 Senior Member   Filippo Maria Denaro Join Date: Jul 2010 Posts: 6,728 Rep Power: 71 your M-file did not run, you have some errors to correct (Lx is not defined when used, normres is not defined, etc)

July 4, 2015, 10:46
#3
New Member

Join Date: Jul 2015
Posts: 2
Rep Power: 0
Quote:
 Originally Posted by FMDenaro your M-file did not run, you have some errors to correct (Lx is not defined when used, normres is not defined, etc)
Thank you Filippo, I have corrected the error and the code is now running. I also attached the streamline plot for Reynolds number of 100.

Many Thanks.

 July 4, 2015, 12:17 #4 Senior Member   Filippo Maria Denaro Join Date: Jul 2010 Posts: 6,728 Rep Power: 71 if you are using primitive variables you must debug your solution (u,v,p). The streamlines are just a post-processing that can be affected by some factors. However, the figure you attached shows a quite regular pattern for Re=100. Use a fine grid (for example 100x100) to ensure Re_h =1 and check the velocity profiles (Ghia&Ghia paper).