SUMMARYA Hermitian-Fourier numerical method for solving the Navier-Stokes equations with one non-homogeneous direction had been presented by Schiestel and Viazzo (Internat. J. Comput. Fluids 1995; 24(6):739). In the present paper, an extension of the method is devised for solving problems with two non-homogeneous directions. This extension is indeed not trivial since new algorithms will be necessary, in particular for pressure calculation. The method uses Hermitian ÿnite di erences in the non-periodic directions whereas Fourier pseudo-spectral developments are used in the remaining periodic direction. Pressure-velocity coupling is solved by a simpliÿed Poisson equation for the pressure correction using direct method of solution that preserves Hermitian accuracy for pressure. The turbulent ow after a backward facing step has been used as a test case to show the capabilities of the method. The applications in view are mainly concerning the numerical simulation of turbulent and transitional ows.