lid driven cavity primitive variables
 

I'm working on modifying the code that was written by Benjamin Seibold at MIT. The matlab code can be found here:

http://www-math.mit.edu/cse/codes/mi...navierstokes.m

I'm non dimensionalized using u* = u/U and v = v*UH/L. I am having trouble figuring out how this affects the poisson pressure solver. So far my modifications have not resulted in a divergence free velocity field.