Cette thèse traite de la modélisation et de l'approximation numérique des écoulements liquide-gaz compressibles. La difficulté essentielle réside dans la modélisation et l'approximation de l'interface liquide-gaz. Schématiquement, deux types de méthodes permettent l'étude de la dynamique de l'interface : l'approche eulérienne, aussi dite de capture de front ("front capturing method") et l'approche lagrangienne, de suivi de front ("front tracking method"). Nos travaux sont plutôt basés sur la méthode de capture de front. Le modèle bifluide est constitué d'un système de lois de conservation du premier ordre traduisant le bilan de masse, de quantité de mouvement et d'énergie du système physique. Ce système doit être fermé par une loi de pression du mélange gaz-liquide pour que sa résolution soit possible. Cette loi de comportement doit être choisie soigneusement, puisqu'elle conditionne les bonnes propriétés du système comme l'hyperbolicité ou l'existence d'une entropie de Lax. Les méthodes d'approximation doivent permettre de traduire au niveau discret ces propriétés. Les schémas conservatifs classiques de type Godunov peuvent être appliqués au modèle bifluide. Ils conduisent cependant à des imprécisions qui les rendent inutilisables en pratique. Enfin, l'existence de solutions discontinues rend difficile la construction de schémas d'ordre élevé. La structure complexe des solutions nécessite alors des maillages très fins pour une précision acceptable. Il est donc indispensable de proposer des algorithmes performants pour les calculateurs parallèles les plus récents. Au cours de cette thèse, nous allons aborder partiellement chacune de ces problématiques : construction d'une "bonne" loi de pression, construction de schémas numériques adaptés, programmation sur calculateur massivement multicoeur.