SUMMARY An implicit fractional-step method for the numerical solution of the time-dependent incompressible Navier–Stokes equations in primitive variables is studied in this paper. The method, which is first-orderaccurate in the time step, is shown to converge to an exact solution of the equations. By adequately splitting the viscous term, it allows the enforcement of full Dirichlet boundary conditions on the velocity in all substeps of the scheme, unlike standard projection methods. The consideration of this method was actually motivated by the study of a well-known predictor–multicorrector algorithm, when this is applied to the incompressible Navier–Stokes equations. A new derivation of the algorithm in a general setting is provided, showing in what sense it can also be understood as a fractional-step method; this justifies, in particular, why the original boundary conditions of the problem can be enforced in this algorithm. Two different finite element interpolations are considered for the space discretization, and numerical results obtained with them for standard benchmark cases are presented. © 1998 John Wiley & Sons, Ltd.