SUMMARYA new regularization method is proposed for the Galerkin approximation of the incompressible NavierStokes equations with Q1=P0 element, by newly introducing a square-type linear form into the variational divergence-free constraint regularized with the global pressure jump (GPJ) method. The addition of the square-type linear form is intended to eliminate the hydrostatic pressure mode appearing in conÿned ows, and to make the discretized matrix positive deÿnite and then non-singular without the pressure pegging trick. E ects of the free parameters for the regularization on the solutions are numerically examined with a 2-D driven cavity ow problem. Furthermore, the convergences in the conjugate gradient iteration for the solution of the pressure Poisson equation are compared among the mixed method, the GPJ method and the present method for both leaky and non-leaky 3-D driven cavity ows. Finally, the non-leaky 3-D cavity ows at di erent Re numbers are solved to compare with the literature data and to demonstrate the accuracy of the proposed method.