Cette thèse s'intéresse au modèle diphasique de Baer-Nunziato. L'objectif de ce travail est de proposer quelques techniques de prise en compte de la disparition de phase, régime occasionnant d'importantes instabilités au niveau du modèle et de sa simulation numérique. Par des méthodes d'analyse et de simulation reposant sur les techniques d'approximation par relaxation à la Suliciu, on montre que dans ces régimes, on peut stabiliser les solutions en introduisant une dissipation de l'entropie totale de mélange. Dans une première approche dite approche Eulerienne directe, la résolution exacte du problème de Riemann pour le système relaxé permet de définir un schéma entropique extrêmement précis, et qui se révèle bien plus économique en terme de coût CPU (à précision donnée) que le schéma classique très simple de Rusanov. De plus, nous montrons que ce schéma permet de simuler avec robustesse des régimes de disparition de phase. Le schéma est développé en 1D puis étendu en 3D et intégré à un prototype de code industriel développé par EDF. La deuxième approche, dite approche par splitting acoustique, propose une séparation des ondes acoustiques rapides et des ondes de transport lentes. L'objectif est d'éviter la résonance due à l'interaction entre ces deux types d'ondes, et de permettre à long terme un traitement implicite de l'acoustique, et explicite du transport. Le schéma, très simple, permet la prise en compte simple de la disparition de phase. La nouveauté est ici l'exploitation de fermetures dissipatives nouvelles du couple vitesse et pression d'interface, qui permettent le contrôle des solutions du problème de Riemann associé à l'étape acoustique.

Dans ce travail, on s'intéresse à plusieurs aspects de l'approximation numérique des systèmes hyperboliques de lois de conservation. La première partie est dédiée à la construction de schémas d'ordre élevé sur des maillages 2D non structurés. On développe une nouvelle technique de reconstruction de gradients basée sur l'écriture de deux schémas MUSCL sur deux maillages imbriqués. Cette procédure augmente le nombre d'inconnues numériques, mais permet d'approcher la solution avec une grande précision. Dans la deuxième partie, on étudie la stabilité des schémas d'ordre élevé. On montre dans un premier temps que les inégalités d'entropie discrètes usuelles vérifiées par les schémas d'ordre élevé ne sont pas pertinentes pour assurer le bon comportement dans le régime de convergence. On propose alors une extension des techniques de limitation {\it a posteriori} pour forcer la vérification des inégalités d'entropie discrètes requises. Dans la dernière partie, on s'intéresse à la construction de schémas well-balanced pour le modèle de Saint-Venant, le modèle de Ripa et les équations d'Euler avec gravité. On propose plusieurs stratégies permettant d'obtenir des schémas numériques capables de préserver tous les régimes stationnaires au repos. On développe également des extensions d'ordre élevé.