Great code, but check your boundary conditions. The zero-velocity BC along the top wall is more of a finite-difference BC. Since you're doing a finite-volume formulation, shouldn't those boundary cell velocities be the negatives of those in the first active cell? That way you have zero velocities AT THE WALL. You do have the symmetry plane BC coded correctly for a finite-volume formulation.
Great code, but check your boundary conditions. The zero-velocity BC along the top wall is more of a finite-difference BC. Since you're doing a finite-volume formulation, shouldn't those boundary cell velocities be the negatives of those in the first active cell? That way you have zero velocities AT THE WALL. You do have the symmetry plane BC coded correctly for a finite-volume formulation.
Comment only