Abstract
A scheme for the numerical simulation of incompressible flows is presented. The modeling equations are written in conservation law form. The algorithm, written in a delta form, is very robust without resorting to any degree of relaxation. The well-known cumbersome numerical implementation related to staggered grids is totally removed. The scheme is second-order-accurate in space and first-order in time. An important novelty is the separation of the pressure terms in the physical equations in a pressure flux, facilitating the tackling of different fluids. The algorithm has been applied to laminar and turbulent flows. The results are truly encouraging.