Le modèle Euler compressible bifluide ne présente pas de difficultés théoriques supplémentaires comparé au cas monofluide. Mais sa résolution numérique est notoirement plus difficile à cause du phénomène d'oscillations de pression à l'interface entre fluides. Nous présentons une approche basée sur un échantillonnage aléatoire "à la Glimm" à l'interface, qui permet de s'affranchir de ce défaut. Le schéma obtenu est applicable à des maillages non structurés, il a d'excellentes propriétés de robustesse et de convergence. Nous l'appliquons à des cas de déferlement.