Simulation in primitive variables of incompressible flow with pressure Neumann condition

We develop a velocity- pressure algorithm, in primitive variables and finitc differences, for incompressible viscous flow with a Neumann pressure boundary condition. The pressure field is initialized by least- squares and up-dated from the Poisson equation in a direct weighted manner. Simulations with the cavity problem were made for severa) Reynolds numbers. H was obtained the expected displacement of the central vortex as well as the development of secondary and terciary eddies