In this paper, we revisit the simultaneous numerical solution method for the steady-state Navier-Stokes equations through Picard’s iteration. In a simultaneous method, all the primary variables at all the grid points are solved simultaneously, as opposed to the sequential method, where the primary variables are solved sequentially. This is a consequence of retaining the coupling between the primary variables. After a brief review of the previous studies where the simultaneous solution is attempted in the finite volume context, the discretization of the governing equations and the idea of Picard’s iteration are discussed. We employ second-order accurate central difference schemes to approximate the first and second-order spatial derivatives. We present numerical solutions of 2D lid-driven cavity flow and flow past a backward-facing step for low and high Reynolds numbers and compare the solutions and the performance of the method to that of an OpenFOAM solver. We conclude by discussing the advantages of the simultaneous solution method and the major hurdles in its widespread deployment.