Cette thèse est consacrée à la simulation numérique de problèmes d'écoulements de fluides raides régis par des systèmes de lois de bilan non homogènes, dans des configurations monodimensionnelles et bidimensionnelles. La méthode numérique utilisée est une extension d'un schéma à deux pas (SRNH), comportant un paramètre \alpha^n_(j+\frac(1)(2)) ajustable, proposé par le professeur F.Benkhaldoun dans un cadre monodimensionnel. Ainsi, en un premier temps on a introduit une variante SRNHR, obtenue en remplaçant la vitesse numérique (\frac(\Delta x)(\Delta t)) par la vitesse de Rusanov locale, en vue de l'extension du schéma au cas bidimensionnel. Par la suite, une analyse de stabilité du schéma, révèle que celui-ci peut être d'ordre 1 ou 2 selon la valeur du paramètre \alpha^n_(j+\frac(1)(2)). Une stratégie de variation de ce paramètre, basée sur la théorie des limiteurs a alors été adoptée. Le schéma peut ainsi être rendu d'ordre 1 dans les zones à forte variation de l'écoulement, et d'ordre 2, là où l'écoulement est régulier. Ensuite on a établi les conditions pour que ce schéma respecte la C-propriété exacte introduite par Bermùdez et Vazquez. Une étude d'implémentation des conditions aux limites, adaptée à ce schéma, a également été menée en se basant sur les invariants de Riemann. Dans la deuxième partie de la thèse, on a appliqué ce schéma à des systèmes monophasiques homogènes et non homogènes. Par exemple on a réalisé la simulation du problème de rupture de barrage sur une marche, pour des configurations 1D et 2D, en menant en particulier une étude de convergence numérique via la détermination des courbes d'erreurs. Enfin, on a utilisé le schéma pour la simulation numérique de systèmes diphasiques (Ransom 1D et 2D).