Ecoulements Compressibles et Supersoniques

Feuille d’exercices IV – Compléments sur les ondes de choc obliques

1 Réflection d’une onde de choc oblique

On considère la réflection d’un choc oblique permanent pour un gaz idéal polytropique de
constante γ = 1.4 initialement généré par la déflection d’un angle ϕA = 16◦ d’un écoulement à
Mach M1 = 2.8. Voir Figure 1.

trajectoire Ch
oc 3
1 2 Ré
t

fle
en

ch
id

i
c
In
oc
Ch

Figure 1 – Réflection d’un choc oblique.

1. À l’aide de l’abaque de choc fournie, calculer l’angle de choc oblique : σA .

2. Calculer les propriétés de l’écoulement dans la zone 2, zone aval du choc incident. On calculera
en particulier le nombre de Mach M2 , les rapports de pressions et de températures.
3. On suppose que le choc initié en A atteint la paroi supérieure au point B et se réflechit au sens
où l’écoulement dans la partie aval, la zone 3, est redevenu parallèle à la paroi. Déterminer
σB et en déduire Φ, l’angle que forme le choc réflechi avec la paroi.
4. De la même manière que pour 2), calculer les propriétés de l’écoulement dans la zone 3.
5. Représenter graphiquement la réflection à l’aide de la polaire de choc dans le plan (ϕ, log(p/p1 )).

1
2 Étude qualitative de l’interaction de deux ondes de choc obliques

h oc
C
Trajectoire
3 4
tinuité
Discon
1
4'
2
Trajectoire Cho
c'

Figure 2 – Intersection de deux chocs obliques.

Sur la figure 2, on considère maintenant l’interaction de deux chocs obliques issus de la déflection
d’un écoulement à Mach M1 > 1 d’un angle ϕ2 pour la partie inférieure et d’un angle ϕ3 pour la
partie supérieure. On suppose que ces deux chocs collisionnent en E et que leur interaction donne
naissance à deux chocs réflechi et une discontinuité de contact de type ligne de glissement faisant
un angle Φ avec l’horizontale. On a également représenté deux trajectoires issues de deux points de
l’écoulement amount qui passent de chaque côté du point E. On note ϕ4 la déflection de l’écoulement
de la zone 3 vers la zone 4 et ϕ04 la déflection de la zone 2 à la zone 4’. On rappelle que par définition,
le long de la ligne de glissement on a continuité de l’angle de l’écoulement et continuité des pressions.
Question : En vous aidant du dernier résultat de l’exercice précédent, faire l’analyse qualitative
de la situation dans un diagramme (ϕ, log(p/p1 )).

3 Problème : Plaque infiniment fine dans un écoulement superso-

nique
On considère un écoulement plan permanent autour d’une plaque plane de longueur finie L,
dans la limite d’une épaisseur infiniment faible. Cet écoulement est représenté sur la Figure 3. Au
bord d’attaque, la plaque engendre deux écoulements classiques, une détente de Prandtl-Meyer et

2
une onde de choc oblique, que nous avons étudiées en cours. La structure de l’écoulement engendre
alors une force sur la plaque, induisant portance et trainée ; c’est l’objet de la première sous-partie.
Le but de la seconde sous-partie est le raccord des écoulements au bord de fuite. Tout au long de
cette étude, on se considerera suffisamment près de la plaque pour négliger les possibles interactions
entre les différentes ondes provenant des bords d’attaque et de fuite.

O
V2 2 4
V1

V3

O’
5

Figure 3 – Ondes de choc obliques et détentes de Prandtl-Meyer sur une plaque de longueur finie et
d’épaisseur nulle faisant un angle de 20 degrés avec l’écoulement amont.

3.1 Forces s’exerçant sur la plaque

A l’amont du dispositif, l’écoulement est supersonique, de nombre de Mach M1 , de vitesse V1 =
U1 ex . Dans cette question, l’écoulement est supposé partout supersonique, et pour les applications
numériques, on prendra γ = 1.4 et M1 = 3. On va analyser les principales lignes caractéristiques de
l’écoulement autour de la plaque.

3.1.1 Détente de Prandtl-Meyer

1. Rappeler brièvement comment s’effectue la résolution de la détente centrée de Prandtl-Meyer
à l’aide de la fonction du même nom.
2. L’angle ϕ vaut 20◦ . En utilisant la tabulation de la fonction de Prandtl-Meyer donnée en
annexe et issue du livre de Anderson, calculer le nombre de Mach aval M2 .

3
 3. En déduire le rapport de pression p2 /p1 en utilisant le caractère homoentropique de l’écoule-
ment et la pression génératrice (pression d’arrêt).
4. Evaluer les angles que font les deux C+ caractéristiques, qui délimitent la structure d’onde
simple de détente, avec la direction de l’écoulement amont.

3.1.2 Onde choc oblique

5. En utlisant les abaques de choc oblique (courbes correspondant à un nombre de Mach amont
donné dans le diagramme ϕ − α) donnés en annexe, évaluer l’angle de l’onde de choc oblique
α correspondant à l’angle de déflection ϕ.
6. Exprimer le nombre de Mach normal amont et en déduire le Mach aval M3 ainsi que le rapport
de pression p3 /p1 .
7. Représenter qualitativement les points correspondant aux zones 2 et 3 sur un diagramme
(ϕ, log(p/p1 )) en traçant la polaire de choc oblique et la courbe correspondant à la détente de
Prandtl-Meyer.

3.1.3 Trainée et portance d’onde (“wave drag and lift”)

8. Donner la force de trainée T et de portance P par unité d’envergure de la plaque (c’est-à-dire
par unité de longueur dans la direction orthogonale au plan de la Figure 3) issue du différentiel
de pression entre p2 et p3 sur la plaque de longueur L.
9. Calculer les coefficients de trainée cT et de portance cP définis par
T P
cT = , cP = , (1)
L q1 L q1

où la pression dynamique q1 vaut γ2 p1 M12 .

3.2 Raccords des écoulements en O0

On fait l’hypothèse qu’en O0 , la structure de l’écoulement retrouve sa symétrie, c’est à dire que
Φ vaut 0 et les vitesses V4 et V5 sont colinéaires à ex .
10. Expliquez pourquoi, dans la mesure où l’écoulement en centré, il est nécessaire d’avoir la
présence d’une onde de choc oblique et d’une onde de détente de Prandtl-Meyer qui relient
les états 2 et 3 aux états 4 et 5 respectivement.
11. Quelle doit être la structure de l’écoulement dans les zones 4 et 5 ainsi que les conditions de
saut à l’interface entre ces deux zones ? Quelle est la nature de cette interface ?
12. Dans la mesure où l’on se donne la direction de l’écoulement aval via l’angle Φ, donnez la
valeur du champ de pression dans les zones 4 et 5. Faites l’application numérique pour Φ = 0.
Qu’en déduisez-vous ?
13. L’angle Φ est-il un degré de liberté du problème ou bien est-il déterminé par l’écoulement
amont ? Dans le second cas, expliquez comment il doit être obtenu.
14. Calculer une approximation de Φ. On procèdera par approximation successives en partant
de Φ = 0, on expliquera si l’on doit augmenter ou diminuer Φ pour converger vers la bonne
solution. On calculera en particulier la solution pour laquelle la valeur absolue de Φ vaut 1

4
 degré et on en déduira par interpolations successives une valeur raisonnablement précise de
Φ. Commentez le résultat obtenu.
15. Calculer les rapports de pressions p4 /p1 et p5 /p1 dans les zones 4 et 5, puis donner l’angle µ
que fait la seconde onde de choc oblique issue de O0 avec l’écoulement amont, ainsi que les
nombres de Mach dans ces zones.
16. Tracer la structure des différentes ondes obtenues dans un diagramme (ϕ, log(p/p1 )) de polaire
de choc et situer en particulier les zones 4 et 5 sur ce diagramme.

5
 M ν(M ) M ν(M ) M ν(M ) M ν(M ) M ν(M ) M ν(M )
1.00 0.0000 1.40 8.9870 1.80 20.7251 2.40 36.7465 3.20 53.4703 4.00 65.7848
1.01 0.0447 1.41 9.2756 1.81 21.0139 2.42 37.2289 3.22 53.8261 4.05 66.4393
1.02 0.1257 1.42 9.5650 1.82 21.3021 2.44 37.7077 3.24 54.1791 4.10 67.0820
1.03 0.2294 1.43 9.8553 1.83 21.5898 2.46 38.1831 3.26 54.5294 4.15 67.7132
1.04 0.3510 1.44 10.1464 1.84 21.8768 2.48 38.6551 3.28 54.8770 4.20 68.3332
1.05 0.4874 1.45 10.4381 1.85 22.1633 2.50 39.1236 3.30 55.2220 4.25 68.9423
1.06 0.6367 1.46 10.7305 1.86 22.4492 2.52 39.5886 3.32 55.5643 4.30 69.5406
1.07 0.7973 1.47 11.0235 1.87 22.7344 2.54 40.0503 3.34 55.9040 4.35 70.1285
1.08 0.9680 1.48 11.3169 1.88 23.0190 2.56 40.5085 3.36 56.2411 4.40 70.7062
1.09 1.1479 1.49 11.6109 1.89 23.3029 2.58 40.9633 3.38 56.5756 4.45 71.2738
1.10 1.3362 1.50 11.9052 1.90 23.5861 2.60 41.4147 3.40 56.9075 4.50 71.8317
1.11 1.5321 1.51 12.1999 1.91 23.8687 2.62 41.8628 3.42 57.2369 4.55 72.3801
1.12 1.7350 1.52 12.4949 1.92 24.1506 2.64 42.3074 3.44 57.5639 4.60 72.9192
1.13 1.9445 1.53 12.7901 1.93 24.4318 2.66 42.7488 3.46 57.8883 4.65 73.4491
1.14 2.1600 1.54 13.0856 1.94 24.7123 2.68 43.1868 3.48 58.2102 4.70 73.9701
1.15 2.3810 1.55 13.3812 1.95 24.9920 2.70 43.6215 3.50 58.5298 4.75 74.4824
1.16 2.6073 1.56 13.6770 1.96 25.2711 2.72 44.0529 3.52 58.8469 4.80 74.9863
1.17 2.8385 1.57 13.9728 1.97 25.5494 2.74 44.4810 3.54 59.1616 4.85 75.4818
1.18 3.0743 1.58 14.2686 1.98 25.8269 2.76 44.9059 3.56 59.4739 4.90 75.9691
1.19 3.3142 1.59 14.5645 1.99 26.1037 2.78 45.3275 3.58 59.7838 4.95 76.4486
1.20 3.5582 1.60 14.8604 2.00 26.3798 2.80 45.7459 3.60 60.0915 5.00 76.9202
1.21 3.8060 1.61 15.1561 2.02 26.9295 2.82 46.1611 3.62 60.3968 5.05 77.3843
1.22 4.0572 1.62 15.4518 2.04 27.4762 2.84 46.5731 3.64 60.6998 5.10 77.8409
1.23 4.3117 1.63 15.7473 2.06 28.0197 2.86 46.9820 3.66 61.0005 5.15 78.2902
1.24 4.5694 1.64 16.0427 2.08 28.5600 2.88 47.3877 3.68 61.2990 5.20 78.7324
1.25 4.8299 1.65 16.3379 2.10 29.0971 2.90 47.7903 3.70 61.5953 5.25 79.1677
1.26 5.0931 1.66 16.6328 2.12 29.6308 2.92 48.1898 3.72 61.8893 5.30 79.5962
1.27 5.3590 1.67 16.9276 2.14 30.1613 2.94 48.5863 3.74 62.1812 5.35 80.0180
1.28 5.6272 1.68 17.2220 2.16 30.6884 2.96 48.9796 3.76 62.4709 5.40 80.4332
1.29 5.8977 1.69 17.5161 2.18 31.2121 2.98 49.3700 3.78 62.7584 5.45 80.8421
1.30 6.1703 1.70 17.8099 2.20 31.7325 3.00 49.7573 3.80 63.0438 5.50 81.2448
1.31 6.4449 1.71 18.1034 2.22 32.2494 3.02 50.1417 3.82 63.3271 5.55 81.6413
1.32 6.7213 1.72 18.3964 2.24 32.7629 3.04 50.5231 3.84 63.6083 5.60 82.0319
1.33 6.9995 1.73 18.6891 2.26 33.2730 3.06 50.9016 3.86 63.8874 5.65 82.4166
1.34 7.2794 1.74 18.9814 2.28 33.7796 3.08 51.2771 3.88 64.1645 5.70 82.7956
1.35 7.5607 1.75 19.2732 2.30 34.2828 3.10 51.6497 3.90 64.4395 5.75 83.1689
1.36 7.8435 1.76 19.5646 2.32 34.7825 3.12 52.0195 3.92 64.7125 5.80 83.5368
1.37 8.1276 1.77 19.8554 2.34 35.2787 3.14 52.3864 3.94 64.9836 5.85 83.8993
1.38 8.4130 1.78 20.1458 2.36 35.7715 3.16 52.7505 3.96 65.2526 5.90 84.2565
1.39 8.6995 1.79 20.4357 2.38 36.2607 3.18 53.1118 3.98 65.5197 5.95 84.6085

Table 1 – Tabulation de la fonction de Prandtl-Meyer, notée ici ν(M ). Le résultat est donné en
degré et non en radian.