Le but de cette these est de simuler, par plusieurs methodes numeriques, un ecoulement compressible bifluide, turbulent. Dans une premiere partie du travail correspondant a une analyse bibliographique, nous etudions diverses modelisations d'ecoulements turbulents, avec une approche privilegiee pour la theorie statistique du splitting et les modeles multi-fluides. Nous batissons l'equation d'evolution de l'energie cinetique turbulente du melange, generalisant celle d'un ecoulement mono-fluide. Ceci permet l'elaboration d'un modele bifluide turbulent. La seconde partie du travail consiste a donner une resolution exacte de la partie convective du systeme dont on prouve l'hyperbolicite. Nous nous interessons particulierement aux criteres entropiques permettant de caracteriser l'unique solution. Cette analyse debouche sur la resolution du systeme pour des conditions initiales correspondant a un probleme de Riemann. Dans la troisieme partie, nous presentons les resultats numeriques de l'implementation du schema de Godunov. Plusieurs cas 2D de turbulence bifluide sont donnes, ainsi qu'une validation numerique par comparaison de la solution 1D a la solution exacte du solveur de Riemann. Pour reduire le temps de calcul, nous proposons deux solveurs approches: le schema Vfroe-nc et le schema de Rusanov. Ceux-ci montrent cependant leurs faiblesses quant a la fiabilite du resultat. Ainsi, nous avons privilegie une resolution par solveur exact du modele bifiuide turbulent puisque la robustesse du schema le necessitait. Ce travail se poursuit par la prise en compte de la partie diffusive du modele K-Epsilon complet. (auteur)