Une méthode pseudo-spectrale 3D a été développée pour la résolution des équations de Navier-Stokes instationnaires en approximation faible nombre de Mach. Ce modèle régissant les écoulements de fluides dilatables permet de considérer de fortes variations de densité en négligeant les échelles acoustiques. Le problème de Stokes issu d'une discrétisation Fourier-Galerkin/Chebyshev-Collocation, et d'un shéma semi-implicite du second ordre en temps, est résolu avec un algorithme itératif de type Uzawa préconditionné. Le code est validé numériquement sur une solution analytique et physiquement sur deux benchmarks en régime stationnaire et instationnaire. Appliqué à un problème de convection naturelle dans une cavité différentiellement chauffée de rapport de forme 8 et pour un grand écart de température, cet algorithme a permis de mettre en évidence des modes d'instabilité 3D stationnaires et instationnaire mettant à défaut l'hypothèse 2D. La seconde application concerne un problème de combustion isobare en domaine ouvert où une première validation 2D de l'instabilité hydrodynamique de Darrieus-Landau a été entreprise suivie d'une simulation exploratoire d'un régime 3D.