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

Abdels83 · July 3, 2015, 18:47

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 find code attached.

Please any help with code modification will be appreciated.

Thank you.

FMDenaro · July 4, 2015, 06:55

your M-file did not run, you have some errors to correct (Lx is not defined when used, normres is not defined, etc)

Abdels83 · July 4, 2015, 09:46

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.

FMDenaro · July 4, 2015, 11:17

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).

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[blockMesh] non-orthogonal faces and incorrect orientation?	nennbs	OpenFOAM Meshing & Mesh Conversion	7	April 17, 2013 05:42
formulation of finite element&finite volume method	ravindra	FLUENT	0	June 27, 2007 12:04
finite difference method for navier-stokes problem	dallybird	Main CFD Forum	5	February 17, 2003 22:00
tidal flow simulation using finite volume method	Jason Qiu	Main CFD Forum	0	October 20, 2002 02:34
SIMPLE Algorithm & Conjugate Gradient Method	Abhijit Tilak	Main CFD Forum	7	March 20, 2001 00:01

July 4, 2015, 06:55		#2
FMDenaro Senior Member Filippo Maria Denaro Join Date: Jul 2010 Posts: 6,768 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, 11:17		#4
FMDenaro Senior Member Filippo Maria Denaro Join Date: Jul 2010 Posts: 6,768 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).