Des bulles haute-pression en expansion dans un liquide se rencontrent dans plusieurs problemes de surete etudies dans l'industrie nucleaire: accidents de type Borax dans des reacteurs experimentaux de type piscine, Accident Hypothetique de Perturbation du Coeur ou Hypothetical Core Disruptive Accident (HCDA) dans un reacteur a neutrons rapides refroidi au sodium. Un accident de type Borax implique une bulle haute-pression creee par l'interaction de materiaux metalliques fondus avec l'eau. Pour un reacteur a neutrons rapides refroidi au sodium, le scenario d'un HCDA suppose que le coeur du reacteur a partiellement fondu et que le combustible fondu interagit chimiquement avec le sodium liquide pour produire une large quantite de composants gazeux. L'expansion explosive de la bulle de gaz haute pression ainsi formee genere des charges significatives sur l'enceinte du reacteur. Une large gamme de nombre de Mach est observee dans un tel ecoulement diphasique compressible, avec un regime bas-Mach dans le liquide mais un regime a Mach eleve dans l'autre phase. Le calcul d'ecoulements a Mach eleve avec une methode de type Godunov en formulation volumes finis centree sur la cellule est desormais une tache de routine. Cependant, lorsque l'on traite des ecoulements bas-Mach, ces solveurs de type Godunov souffrent d'un probleme de precision et d'efficacite pour lequel des remedes ont ete proposes dans la litterature. Dans le but de de comprendre le comportement bas-Mach pour des cas avec ou sans presence d'ondes acoustiques, l'analyse numerique des schemas est tout d'abord realisee pour des ecoulements monophasiques de liquide, decrits a l'aide d'une equation d'etat de type gaz raides. La correction bas-Mach de Dellacherie et al. est etendue aux ecoulements liquides et il est montre que le probleme de precision a faible Mach peut effectivement etre resolu en centrant la pression. Cependant, l'efficacite de la simulation reste a ameliorer. Une methode originale a 2 etapes est alors proposee: dans une premiere etape, une methode de type compressibilite artificielle est mise en oeuvre pour accelerer la resolution du systeme hyperbolique decrivant la propagation d'ondes acoustiques artificielles de faible amplitude; dans une seconde etape, une transformation isentropique est appliquee pour prendre en compte les mecanismes de compression / detente en domaines fermes. La methode est etendue aux maillages mobiles a l'aide d'une approche de type ALE et egalement etendue a un modele diphasique a 7 equations avec l'approche Discrete Equation Method (DEM) afin de calculer des ecoulements gaz-liquide bas-Mach. Elle est appliquee avec succes a un probleme d'expansion bas-Mach d'une bulle dans une enceinte fermee. (auteur)