The purpose of this paper is to investigate the performance of a certain family of low-order quadrilateral finite elements in the solution of the stationary Navier-Stokes equations governing incompressible fluid flows. These finite elements are derived from the ‘overlapping finite elements’, first developed for the solution of problems in solid mechanics [5,9,48] and incorporate features of both meshfree and traditional finite element methods. One of their most remarkable properties is the insensitivity to mesh distortions. Also, since the Shepard functions are replaced by a suitable interpolation, the resulting basis functions are entirely polynomial, which allows the numerical integration of the weak forms to be performed with few integration points [6,49]. We also discuss the theoretical reasons for the stability of the solutions, that is, the absence of spurious pressure modes, and show that the proposed discretization scheme passes the relevant inf-sup test.