||February 23, 2010 13:17
2D Lid Driven Cavity Flow simulation using MATLAB
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!