Incompressible flow calculations with a consistent physical interpolation finite volume approach

The computation of incompressible three-dimensional viscous flow is discussed. A new physically consistent method is presented for the reconstruction for velocity fluxes which arise from the mass and momentum balance discrete equations. This closure method for fluxes allows the use of a cell-centered grid in which velocity and pressure unknowns share the same location, while circumventing the occurrence of spurious pressure modes. The method is validated on several benchmark problems which include steady laminar flow predictions on a two-dimensional cartesian (lid driven 2D cavity) or curvilinear grid (circular cylinder problem at Re = 40), unsteady three-dimensional laminar flow predictions on a cartesian grid (parallelopipedic lid driven cavity) and unsteady two-dimensional turbulent flow predictions on a curvilinear grid (vortex shedding past a square cylinder at Re = 22,000).