Cette thèse a pour but de présenter les travaux liés à l'étude numérique des écoulements critiques diphasiques dans la tuyère primaire d'un éjecteur utilisé pour améliorer le rendement d'un cycle de réfrigération. En premier lieu, les bases d'un modèle d'écoulement multiphasique ainsi que les termes de source associés sont présentés. Les développements entourant les points singuliers ainsi que les conditions critiques sont également présentés. Deux schémas numériques unidimensionnels sont proposés pour l'étude d'un écoulement critique diphasique particulaire. Les comparaisons des flux massiques critiques ainsi que des profils de pression numériques et expérimentaux montrent que le modèle représente bien l'évolution de l'écoulement diphasique. L'existence d'un point critique situé à un endroit différent du col géométrique a été vérifiée. L'étude numérique présente également une validation des conditions critiques utilisées ainsi que les effets des déséquilibres sur la position du point critique et sur l'écoulement. Des travaux supplémentaires sont cependant requis pour mieux comprendre le phénomène entourant les chocs en écoulements diphasiques et pour accroître la plage de taux de vide couverte par les schémas numériques.