Modélisation de l'écoulement diphasique dans les injecteurs Diesel

CHAPITRE 6 : VALIDATION DU CODE CAVIF 125fl uide sur la surface du piston est bien celle du piston, validant ainsi la condition limite à la paroi mobile.La vitesse de déplacement de l’onde est égale à la vitesse du son dans le gaz, et le profil de vitesse correspondbien à l’expression donnée par Landau et Lifchitz. La vitesse est dirigée partout dans le sens de01.2e−03−501.0e−03w (m/s)−100ρ (kg/m 3 )8.0e−04−150Solution numeriqueSolution analytique6.0e−04−2000e+00 2e−02 4e−02 6e−02 8e−02z (m)4.0e−040e+00 2e−02 4e−02 6e−02 8e−02z (m)FIG. 6.3: Compression d’un gaz isentropique. A gauche : profils de vitesse aux temps t =2 5 µs, 5 0µs, 7 5 µs, 100µs, 12 5 µs, 15 0µs, 17 5 µs, 2 00µs, 2 2 5 µs, 2 5 0µs. A droite : profils demasse volumique aux mêmes temps200.0Solution numeriqueSolution analytique150.05.0e−04w (m/s)100.0ρ (kg/m 3 )4.0e−0450.03.0e−040.00e+00 2e−02 4e−02 6e−02 8e−02 1e−01z (m)2.0e−040e+00 2e−02 4e−02 6e−02 8e−02 1e−01z (m)FIG. 6.4: Détente d’un gaz isentropique. A gauche : profils de vitesse aux temps t =2 5 µs, 5 0µs, 7 5 µs, 100µs, 12 5 µs, 15 0µs, 17 5 µs, 2 00µs, 2 2 5 µs, 2 5 0µs. A droite : profils demasse volumique aux mêmes tempsdéplacement du piston et décroît de façon monotone (en valeur absolue) entre la paroi du piston et le basdu domaine. Notons que la pression et la masse volumique décroissent de façon monotone également, etdans le même sens. La compression étant isentropique, nous avons vérifié la constance du terme P V γ , avecγ = 1.4 tout au long de la compression. Nous avons arrêté le calcul avant l’apparition de l’onde de choc.Un test de détente isentropique a également été effectué, en commençant le calcul avec les mêmesconditions initiales que celles du test précédent. Maintenant a = 8.0 × 10 7 . En comparant les résultatsnumériques avec la solution analytique, on voit un très bon accord à l’instar de la compression. Nous avons