La simulation d'écoulements dans des cavités cylindriques en rotation présente une difficulté particulière en raison de l'apparition de singularités sur l'axe. Le présent travail propose une méthode collocative pseudospectrale suffisamment efficace et précise pour surmonter cette difficulté et résoudre les équations 3D de Navier-Stokes écrites en coordonnées cylindriques. Cette méthode a été développée dans le cadre des différentes études menées au laboratoire M2P2, utilisant une méthode collocative de type Chebychev dans les directions radiale et axiale et Fourier-Galerkin dans la direction azimutale [thêta]. Pour éviter de prescrire des conditions sur l'axe, une nouvelle approche a été développée. Le domaine de calcul est défini par (r,[thêta],z)∈[-1,1]×[0,2π]×[-1,1] avec un nombre N pair de points de collocation dans la direction radiale. Ainsi, r=0 n'est pas un point de collocation. La distribution de points de type Gauss-Lobatto selon r et z densifie le maillage seulement près des parois ce qui rend l'algorithme bien adapté pour simuler les écoulements dans des cavités cylindriques en rotation. Dans la direction azimutale, le chevauchement des points dû à la discrétisation est évitée par l'introduction d'un décalage égal à π/2K à [thêta]>π dans la transformée de Fourier. La méthode conserve la convergence spectrale. Des comparaisons avec des résultats expérimentaux et numériques de la littérature montrent un très bon accord pour des écoulements induits par la rotation d'un disque dans des cavités cylindriques fermées. / When simulating flows in cylindrical rotating cavities, a difficulty arises from the singularities appearing on the axis. In the same time, the flow field itself does not have any singularity on the axis and this singularity is only apparent. The present work proposes an efficient and accurate collocation pseudospectral method for solving the 3D Navier-Stokes equations using cylindrical coordinates. This method has been developed in the framework of different studies of rotor-stator flows, using Chebyshev collocation in the radial and axial directions and Fourier-Galerkin approximation in the azimuthal periodic direction [thêta]. To avoid the difficulty on the axis without prescribing any pole and parity conditions usually required, a new approach has been developed. The calculation domain is defined as (r,[thêta];,z)∈[-1,1]×[0,2π]×[-1,1] using an even number N of collocation points in the radial direction. Thus, r=0 is not a collocation point. The method keeps the spectral convergence. The grid-point distribution densifies the mesh only near the boundaries that makes the algorithm well-suited to simulate rotating cavity flows where thin layers develop along the walls. In the azimuthal direction, the overlap in the discretization is avoided by introducing a shift equal to π/2K for [thêta]>π in the Fourier transform. Comparisons with reliable experimental and numerical results of the literature show good quantitative agreements for flows driven by rotating discs in cylindrical cavities. Associated to a Spectral Vanishing Viscosity, the method provides very promising LES results of turbulent cavity flows with or without heat transfer.