Recueil Examen ECS 1415 PDF

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 71

ECOULEMENTS COMPRESSIBLES ET SUPERSONIQUES

Contrôle de connaissance

Marc Massot

Durée 3 heures (Les notes de cours et les documents distribués en cours


sont autorisés ainsi que l’usage de calculatrices).

Le contrôle de connaissance se compose de deux problèmes dont certaines parties sont des
exercices de cours, qui peuvent être traitées de manière indépendante. Dans les deux problèmes, on
considère des écoulements de gaz parfaits idéaux polytropiques (pas de phénomènes de dissipation,
loi d’état de gaz idéal, p = ρ r T , r = R/M où M est la masse molaire du gaz ; cp ainsi que γ
sont des constantes données). Le contrôle de connaissance est volontairement long et il n’est pas
forcément nécessaire de résoudre l’intégralité du contrôle pour obtenir une très bonne évalution.
Cette évaluation reposera aussi sur la clarté de la solution et des arguments proposés.

1 Problème I : Evaluation du Mach amont par Tube de Pitot


On considère un écoulement plan permanent autour d’un dièdre symétrique d’angle de déflection
ϕ =15 degrés, représenté sur la Figure 1. Au bord d’attaque, on trouve un choc oblique d’angle σ
par rapport à l’écoulement amont de nombre de Mach M1 , de vitesse V1 = U1 ex . La zone en amont
du choc oblique est notée 1 et les grandeurs qui s’y réfèrent possèdent un indice 1.
La pression mesurée par le tube de Pitot est de 2.596 atm. Un capteur positionné au point A,
dans la zone 3, derrière l’onde détente, détecte une pression de 0.1 atm. Le but du problème est la
détermination du nombre de Mach amont M1 .

1.1 Onde de choc oblique

(a) Rappeler la structure d’une onde de choc oblique ainsi que la relation qui relie l’angle de
déflection ϕ, l’angle du choc et la donnée de l’écoulement amont. Combien de solutions de cette
relation sont envisageables, l’angle de déflection ϕ étant donné ? Qu’observe-t-on physiquement ?
(b) Donner, vu l’angle de déflection proposé, la gamme de nombre de Mach amont permettant
l’existence d’un choc oblique. Préciser comment obtenir dans ce cas la valeur de σ, le nombre de
Mach amont étant donné.
(c) Dans le cas du dièdre étudié dans ce problème, du fait de la présence d’un second changement
de direction en O ′ , préciser le domaine spatial dans lequel la solution de type onde de choc oblique
est valable.
(d) Pour un nombre de Mach amont fixé, donner l’expression de M2 en fonction de M1 , et des
angles de l’onde de choc oblique σ et de déflection ϕ.

1
C+

C+
α3
2
σ Tube de Pitot
1 O′ 3
A
V2 4
M1 > 1
V1
O ϕ = 15◦ −ϕ = −15◦

Figure 1 – Schéma de principe de la mesure de vitesse amont par utilisation d’un Tube de Pitot.

1.2 Onde de détente de Prandtl-Meyer


(a) 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.
(b) Présenter précisément la structure de la détente centrée de Prandtl-Meyer et donner l’éten-
due spatiale de validité de l’onde centrée dans le cas du dièdre considéré.
(c) Si l’on suppose connu le nombre de Mach M3 , comment trouve-t-on α3 , angle de Mach de
la dernière C+ caractéristique qui délimite l’éventail de détente centrée ?
(d) Par conservation de l’invariant de Riemann, quel est l’angle de Mach α2 , angle de Mach de
la première C+ caractéristique qui délimite l’éventail de détente centrée en amont.

1.3 Onde de choc devant le Tube de Pitot


L’écoulement dans la zone 3 doit être supersonique, et par conséquent il se forme une onde
de choc normale devant le tube de Pitot. L’écoulement juste derrière l’onde de choc est noté avec
l’indice 4 comme indiqué sur la Figure 1. Le tube de Pitot permet la mesure de la pression d’arrêt
ou pression génératrice p04 = 2.596 atm.
(a) Justifier pourquoi l’écoulement dans la zone 3 doit être supersonique dans la configuration
choisie (indépendemment de la présence du Tube de Pitot).
(b) En utilisant les tables de choc normal, et connaissant la pression p3 au point A ainsi que la
pression d’arrêt ci-dessus, donner la valeur du nombre de Mach M3 (interpoler au besoin pour être
précis).

2
(c) Expliciter les relations de saut impliquées pour tabuler p04 /p3 en fonction de M3 dans les
tables de choc normal que vous venez d’utiliser. On donnera une formule finale.

1.4 Résolution du problème de départ


(a) Calculer les valeurs successives de α3 , α2 et finalement M2 en utilisant les questions pré-
cédentes. On dispose d’une tabulation en degrés de la fonction de Prandtl-Meyer à la fin de ce
document.
(b) En déduite l’équation qui lie M1 , M2 et l’angle de l’onde de choc oblique σ en utilisant les
valeurs numériques évaluées.
(c) Démontrer que cela définit une solution unique pour le nombre de Mach amont.
(d) Résoudre l’équation implicite en M1 en faisant des itérations successives à partir de M1 = 4,
et expliquer comment vous faites varier M1 et pourquoi.
(e) Que conclure sur l’existence d’un choc oblique attaché en O ?

2 Problème II : Piston en mouvement dans un tube


Ce problème est très proche des questions traitées en cours. On s’intéresse à un tube infini à
symétrie de révolution autour de son axe (x′ Ox) de section constante, dans l’approximation d’un
écoulement instationnaire monodimensionnel. On considère un état initial, à l’instant t=0, constitué
de deux gaz au repos de chaque côté d’un piston auquel on va imposer un mouvement à partir de
l’instant initial t = 0 (Voir Figure 2). Ces gaz sont des fluides parfaits idéaux polytropiques notés
“état 4” en ce qui concerne l’état à gauche du piston (de température T4 , de vitesse u4 = 0 de masse
molaire M4 et de coefficient γ4 ) et “état 1” à droite du piston (de température T1 = T4 , de vitesse
u1 = 0, de masse molaire M1 et de coefficient γ1 ). Le mouvement imposé au piston est décrit par
sa position X(t) et le piston se trouve en x = 0 à l’instant t = 0 (on fait l’hypotèse simplificatrice
que le piston est un piston d’épaisseur infiniment fine). L’expérimentateur applique à l’instant t > 0
une force sur le piston de manière à lui imposer la trajectoire X(t) = a tn , avec a < 0 et n un entier
supérieur ou égal à un (n ≥ 1).

2.1 Cas des deux gaz dans un tube infini avec accéleration progressive
(1) Exprimer dans chaque gaz au repos les célérités du son c40 et c10 en fonction des grandeurs
de l’énoncé.
On considère dans un premier temps le cas où n = 2.
(2) Effectuer, comme dans le cours, la résolution des propriétés de l’écoulement pour les deux
gaz, pour tout temps t < T ∗ , pour lequel l’écoulement reste homoentropique dans chaque gaz.
Proposer une synthèse des principes de base permettant de trouver l’état des deux gaz. Sur quelles
grandeurs effectue-t-on la résolution ? Jusqu’à quel temps ce type de résolution basé sur le caractère
homoentropique des deux écoulements est-il valable ? Que se passe-t-il ensuite ?

3
Gaz 4 piston en x=0 a t=0 Gaz 1
x’ x

Gaz 4 Gaz 1
x’ x
O
piston en X(t) a l’instant t

Figure 2 – Tube rempli de deux gaz séparés par un piston à deux instants successifs.

(3) Tracer le diagramme des caractéristiques dans le plan (x, t). On fera une représentation
précise jusqu’au temps T ∗ puis une représentation qualitative pour les temps supérieurs. Identifier les
diverses zones d’écoulement de manière précise dans le diagramme. Une figure claire sera appréciée.
(4) Donner les expressions de la vitesse dans la direction x, de la vitesse du son c, de la pression
et de la température de chaque côté du piston. Quelle est la force totale appliquée sur le piston
par le gaz ? Si le piston est de masse M , trouver quelle est la force instantanée F (t) appliquée par
l’expérimentateur en appliquant la loi fondamentale de la dynamique.
(5) Préciser l’expression du temps et de la position de formation du choc et donner l’expression
de la vitesse en tout point (x, t) pour tout temps inférieur au temps de formation du choc.

2.2 Cas de la détente dans une partie semi-infinie avec accélération immédiate
On se place alors dans le cas n = 1 et l’on se focalise uniquement sur le gaz 1 dans cette sous-
partie. Comme dans le cas de l’onde simple de détente traité en cours, on considère une vitesse
négative uf > −2c0 /(γ − 1) et l’on considère le cas limite où l’accélération du piston est infinie en
t = 0 de telle manière que le piston se met immédiatement à la vitesse uf pour tout t > 0.
(1) Montrer que l’écoulement dans le gaz 1, sous les hypothèses données, est homoentropique
pour tout temps t > 0.
Dans la suite, on admet que le diagramme des caractéristiques a la forme donnée sur la Figure
3 de manière cohérente avec les propriétés d’invariance des équations d’Euler. Les domaines 1 et 3
sont des zones d’écoulement uniforme et les droites OA et OI les droites frontières de la région 2.
(2) Donner les propriétés de la région 2. En déduire les équations des deux droites OA et OI.
(3) Donner (u, c/c0 , p/p0 ) dans les domaines 1 et 3.
(4) Soit M (t, x) un point de la région 2. Déterminer (u, c) en ce point en fonction de (t, x).
(5) Pour un point P (τ, X) dans la région 2, vérifier que la C− de cette région qui passe par P
admet pour équation :
! " ! " 3−γ
2 c0 t 2 c0 τ t γ+1
x=− + X+ . (1)
γ−1 γ−1 τ

4
Trajectoire du piston
x = uf t

t 3 2
A

1
x

Figure 3 – Onde de détente centré au point O.

(Indice : pour cela on cherchera une équation sur la nouvelle variable f = x/t comme une fonction
du temps t et on écrira une équation différentielle sur cette variable. On montrera alors que la forme
proposée est bien solution de cette équation différentielle du premier ordre.)

2.3 Bilan sur la présence des deux gaz (Cette sous-partie aura un poids significatif
dans le barême)

On revient sur l’étude de la dynamique du piston en présence des deux gaz 1 et 4 de chaque côté
du piston. On se place dans le cas où X(t) = a tn , n = 1, d’une accéleration immédiate du piston
et mise en vitesse à l’instant initial t = 0.
(1) Expliquer ce qui se passe dans le gaz 4 et effectuer une représentation précise et justifiée
du diagramme des caractéristiques associé pour les deux gaz à partir des résultats de la sous-partie
précédente.
(2) Ce type de diagramme est proche de celui du tube à choc. Expliquer quelle est la différence
fondamentale entre le cas étudié dans cet exercice et le cas du tube à choc. Dans quel cas les deux
diagrammes sont-ils equivalents ? Représenter dans un diagramme de phase pression-vitesse les deux
cas, afin d’expliquer la différence de manière claire quand elle est présente. Quand les deux cas sont
équivalents, que dire la force appliquée par l’expérimentateur ?

2.4 Cas d’un tube de longueur fini L dans le gaz 1 (Optionnel - difficile)
On considère une zone de gaz 1 dans un tube de longueur L représenté sur la Figure 4. On
veut alors résoudre la question de la propagation et de la réflexion sur le bout du tube de l’onde de
détente initiée par le mouvement du piston.
Ayant résolu le problème associé à l’onde simple de détente avant toute interaction, on se propose
maintenant de considérer le problème de la réflexion de l’onde au bout du tube. Pour des temps
relativement faibles, on suppose que le diagramme des caractéristiques prend la forme donnée sur
la Figure 5.

5
t=0

piston Gaz
x

Figure 4 – Piston et gaz immobiles à l’instant t = 0 dans un tube de longueur L.

(a) Donner les propriétés de la région 5 de l’écoulement, quelle est la nature des courbes repré-
sentées dans cette zone de l’écoulement ? Donner leurs propriétés. Quelle est la nature des courbes
AI et IB.
(b) En admettant éventuellement la relation (1), donner les coordonnées (xI , tI ) du point I.
(c) Donner la valeur de la vitesse sur les demi-droites BB ′ et BB ′′ . En remarquant que deux
courbes caractéristiques d’une même famille ne peuvent pas se couper, démontrer que le domaine
6 est une zone d’écoulement uniforme. Donner dans cette zone les grandeurs (u, c, p, T ) en fonction
des grandeurs de l’énoncé.
(d) On souhaite à présent donner l’évolution de la pression le long de la paroi au cours du temps.
Si h désigne l’enthalpie spécifique du fluide et φ(t, x) un potentiel des vitesses de l’écoulement, on
peut toujours introduire une fonction auxiliaire g, définie tout comme φ à une contante additive
près, telle que
u2
! "
g(u, h) = φ − u x + t h + , (2)
2
ce qui permet d’obtenir les relations importantes :
! " ! " ! "
∂g ∂g ∂g
t= , x=u − . (3)
∂h u ∂h u ∂u h
Montrer que dans le domaine 2, la fonction g est constante ; elle sera alors choisie nulle. Donner les
conditions vérifiées dans le domaine 4 par g sur le tronçon AB et la courbe frontière AI. Le domaine
4 est-il un domaine d’onde simple ?
(e) Nous cherchons dans ce domaine 4 la fonction g. Nous admettons, juste pour cette dernière
question que γ = 5/3 (gaz monoatomique) et que la fonction g prend, sous cette hypothèse, la forme
générale suivante : # ! " ! "$
1 2c 2c
g(u, h) = F u+ +G u− . (4)
c γ−1 γ−1
En utilisant la question précédente, montrer que :
%! "2 &
L γ−1 2
g(u, h) = g(u, c) = u − c − c0 , (5)
(γ − 1)c 2

dans le domaine 4. En déduire sur le tronçon AB l’équation reliant t/tA à c/c0 . Donner alors les
coordonnées du point B et expliquer comment trouver p(t)/p0 en x = L au cours du temps t.

6
B ′′ B′
6

I 4
3 2
t
A

O L x

Figure 5 – Structure des caractéristiques dans le cas de la réflexion de l’onde de détente.

7
M f (M ) M f (M ) M f (M ) M f (M ) M f (M ) M f (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 f (M ). ATTENTION : le résultat


est donné en degré et non en radian.

8
Figure 6 – Abaque de choc oblique.
9
Tables du Livre de ANDERSON

Table A.1 - Table d’écoulement isentropique (p. 691–695)

p ρ
Les colonnes sont : M, p0 , ρ0 , TT0 , A
A

Table A.2 - Table de choc normal (p. 696–699)

p ρ T2 p02 p02
Les colonnes sont : M1 , p2 , ρ2 , T , , p , M2
1 1 1 p01 1

10
2
Cours Option MAE / Branche AE et Master MAS

ECOULEMENTS COMPRESSIBLES ET
SUPERSONIQUES

M. Massot
Professeur, Ecole Centrale Paris
Laboratoire EM2C - UPR CNRS 288
[email protected]
Documents autorisés - On proposera des réponses aux exercices en essayant d’être concis, en
s’efforçant de donner des réponses aussi complètes que possible et de montrer une compréhension
des questions traitées. En particulier, pour les questions de cours, on ne demande pas de recopier
les notes de cours mais de proposer une synthèse permettant de répondre à la question. Les
quatre parties du contrôle sont indépendantes; par contre dans le problème, il est important de
bien lire l’ensemble du problème et de progresser par étapes.

1
2
1 Modélisation du trafic routier (d’après Leveque 1992
noté sur 19)
On considère l’équation de conservation sur la densité de voiture ρ(t, x) :
∂t ρ + ∂x (f (ρ)) = 0,
où f (ρ) désigne le flux de voiture. En général, une modélisation standard montre que le
flux s’écrit comme le produit de la densité locale par une vitesse u, qui va ici, pour ne pas
complexifier le modèle, être donnée comme une fonction de ρ. Une possibilité permettant de
modéliser l’arrêt du trafic quand la densité de voiture atteint un maximum ρmax ainsi que la
limitation de vitesse umax quand la densité de voiture est faible, est donnée par
! "
ρ
u = umax 1 − .
ρmax
La vitesse umax représente donc la vitesse maximum autorisée et ρmax une densité maxi-
mum correspondant à un bouchon où les voitures sont pare-choc contre pare-choc, à l’arrêt.
L’expression du flux de départ est donc :
! "
ρ
f (ρ) = ρ u = ρ umax 1 − .
ρmax
La forme du flux est telle qu’il est concave (contrairement à Burgers où il est convexe) et
les discontinuités admissibles sont telles que ρG < ρD . Pour simplifier la formulation du
problème et les calculs, on se placera dans le cas où ρmax = 1 et umax = 1, ce qui revient à
adimensionner les variables.
Exercice 1 : Résolution du problème de Riemann (6/19).
I. Rappeler l’équation des courbes caractéristiques et expliquer quelle est la vitesse de
propagation de l’information dans cette équation non-linéaire.

II. Résoudre le problème de Riemann pour ρG = 15 ρmax = 51 , ρD = ρmax = 1, qui cor-


respond au feu rouge. Expliquer ce qui se passe en traçant tant le diagramme des
caractéristiques dans ce cas, que l’allure de la solution à deux instants successifs. Com-
parer la dynamique des solutions suivant la valeur de ρG en référence à la feuille d’exer-
cice 1.
III. Pour le problème de Riemann résolu, on représentera le diagramme dans le plan (x, t)
des trajectoires des véhicules. Ces trajectoires seront notées Y (t) ; elles sont telles que

⎨ dY (t) = u(ρ(t, Y (t))) = 1 − ρ(t, Y (t))



dt (1)
0
Y (0) = x .

On montrera, en particulier que les trajectoires des véhicules, dans le cas de la présence
d’une onde de choc, traversent cette onde.

3
Considérons le problème Cauchy pour la donnée initiale particulière que l’on appellera
le petit ”accordéon” du fait de la structure de sa donnée initiale par rapport au problème
traité dans le cours :

⎨ 1/2, x ∈] − ∞, −1] ∪ [1, +∞[,
0
ρ (x) = 1 + x/2, x ∈] − 1, 0], (2)
1 − x/2, x ∈ [0, 1[,

et l’on veut construire la solution en tout temps de l’équation du trafic pour cette donnée
initiale. Ce cas correspond physiquement à la présence d’une zone très localisée de petit
“bouchon” dans un trafic régulier et on se propose de déterminer la dynamique qui en
résulte.

Exercice 2 : Résolution par la méthode des caractéristiques sur [0, T ∗ [ (4/19).


I. Effectuer la résolution de l’évolution de la solution par la méthode des caractéristiques
et déterminer le temps d’existence d’une solution régulière T ∗ . On présentera la solution
du problème sur cet intervalle [0, T ∗ [ en faisant un graphique de la solution ρ(t, x) pour
deux instants T ∗ /2 et T ∗ et on reliera ces solution au diagramme des caractéristiques.
II. Que se passe-t-il à l’instant T ∗ ? Expliquer.

Exercice 3 : Résolution de la dynamique de la solution faible pour t ≥ T ∗ (9/19).


Après le temps T ∗ , un choc s’est formé qui se propage. Notons σ(t) sa vitesse pour t ≥ T ∗
et ϕ(t) sa position. Nous allons établir une équation différentielle ordinaire permettant de
caractériser les deux fonctions du temps σ(t) et ϕ(t), pour t > T ∗ et la dynamique qui en
résulte.
I. La vitesse de ce choc est donnée par

dϕ(t)
= σ(t, ϕ(t)) = 1 − ρ− (t, ϕ(t)) + ρ+ (t, ϕ(t) ,
$ %
dt
où ρ− et ρ+ sont les valeur de la solution de chaque coté de la discontinuité. Justifier
cette relation qui provient de la relation de Rankine-Hugoniot appliquée à l’onde de
choc qui permet de relier la vitesse propagation du choc à l’amplitude du saut de la
variable.
II. Quelles sont les valeurs de ϕ et de σ à l’instant T ∗ qui permettent de donner la
condition initiale de l’équation différentielle ordinaire sur ϕ(t) ?
III. En utilisant l’équation des caractéristiques dans chacun des domaines à gauche et à
droite de la discontinuité (localisée en ϕ(t) à l’instant t), domaines où la solution est
régulière, ainsi que la propriété que vérifie la solution le long des caractéristiques qui
sont des droites, trouver les valeurs de ρ− (t, ϕ(t)) et ρ+ (t, ϕ(t)) comme des fonctions
simple de t (indication : on montrera que la caractéristique est forcément issue pour

4
l’état droit de ξ appartenant au segment [0, 1] et, si elle passe par le point (t, ϕ(t)),
alors on doit avoir (ϕ(t) − ξ) = (−1 + ξ)t que l’on justifiera et qui permet de trouver
le pied de la caractéristique ξ en fonction du point (t, ϕ(t)).
IV. Nous avons donc obtenu une équation différentielle ordinaire sur ϕ(t) pour t ∈
[T ∗ , +∞] avec une donnée initiale en T ∗ ; résoudre cette équation différentielle en
séparant les variables.
V. Présenter sur deux figures l’allure de la solution aux instants suivant le temps de
formation du choc ainsi que le diagramme des caractéristiques. Que font les droites
caractéristiques par rapport à la dynamique de l’ondre de choc ? Le choc est-il entro-
pique ? Expliquer le comportement qualitatif de la solution, en particulier on pourra
évaluer l’évolution temporelle de l’amplitude du saut de la solution à travers la discon-
tinuité et de sa vitesse.

2 Exercice : Discontinuité dans une tuyère de Laval


(noté sur 9)
On considère l’écoulement dans une tuyère de Laval amorcée représenté sur la Figure 1
dans les cas où la pression de sortie est adaptée et vaut donc p2 pour une pression génératrice
de p0 . Le but de cet exercice est de montrer ce qui se passe quand la pression de sortie décrit
un intervalle [p′1 , p1 ] et génére l’apparition d’un choc dans le divergent tout en restant dans
le contexte d’un écoulement quasi-monodimensionnel.
I. Justifiez le fait que le choc stationnaire ne peut apparaı̂tre que dans le divergent.
II. Donner une description précise de ce qui se passe physiquement quand on observe un
choc dans le divergent. En particulier on indiquera le sens de variation des variables
aérothermiques dans le divergent, d’une part avant le choc et d’autre après le choc en
justifiant la réponse précisément.
III. La présence d’un choc dans le divergent change-t-il le débit de la tuyère amorcée ?
L’enthalpie totale est-elle toujours conservée ? Pourquoi ? La témpérature génératrice
saute-t-elle à travers l’onde de choc ? Expliquer.
Introduisons deux positions de choc x1 et x2 correspondant aux deux valeurs de la pression
de sortie p1(s) et p2(s) respectivement. On veut montrer dans la suite de l’exercice que x2 > x1
implique p2(s) < p1(s) .
i
IV. Si l’on note M(a) le nombre de Mach amont du choc dans la position xi , quelle est la
1 2
relation d’ordre entre M(a) et M(a) ?
V. Quelle relation d’ordre y a-t-il entre la pression génératrice après le choc p10 et p20 ?
En déduire, en utilisant la forme du débit dans la tuyère la relation d’ordre entre les
nombres de Mach en sortie dans les deux cas. Enfin en se rappelant que la la pression

5
Choc droit

A
A

p
p p0
p0

(a)
p1
(b) p0
(c)
0.528
0.528
?
(d)
p2
p0

x
x
M
M

(d)
M=1 ?
M=1

(c)
(b)

(a)

Figure 1 – Tuyere amorcée isentropique à gauche et choc dans le divergent à droite.

6
de sortie sur la pression génératrice est une fonction croissante du nombre de Mach en
sortie, en déduire la relation d’ordre entre les deux pressions de sortie.
VI. En déduire les diverses formes du profil de pression et de Mach pour plusieurs lieu
de l’onde de choc dans le divergent entre le col et la section de sortie. Quel inter-
valle de pression est compatible en sortie pour que l’hypothèse d’écoulement quasi-
monodimensionel reste adaptée ?

3 Problème : Conception d’une tuyère supersonique


par la méthode des caractéristiques (noté sur 22)
On se propose dans ce Problème d’utiliser la théorie des caractéristiques pour les écoule-
ments bi-dimensionnels permanents homentropiques homoénergétiques afin de concevoir
deux types de tuyères supersoniques dont l’écoulement en sortie de tuyère est quasi-1D à
Mach 2.4. La première est utilisée dans le cadre des souffleries supersoniques, pour lesquelles
on doit avoir un écoulement quasi-monodimensionnel d’excellente qualité et que l’on obtient
en faisant varier très lentement la section sur une relative grande distance (plus longue que
ne le montre la figure 2-haut). La seconde est une tuyère pour la propulsion spatiale pour la-
quelle on doit concevoir une tuyère de petite taille et pour laquelle la qualité de l’écoulement
aval est moins importante ; on se permet donc d’avoir une variation de section discontinue,
même si la section varie de manière continue comme indiqué sur la Figure 2-bas.
Pour des raisons de symétrie, on ne va travailler que dans la demie-tuyère supérieure
dans les deux cas. On note Θmax dans le premier cas l’angle que la fait la tuyère avec l’axe
horizontal au point d’inflexion C, c’est-à-dire l’angle maximum de l’écoulement au bord par
rapport à l’écoulement aligné avec l’axe de symétrie au centre. Pour la second tuyère, l’angle
maximum est atteint au point A pour lequel il vaut Θmax−C (C pour Court). En ce point, on
observe une détente de Prandtl-Meyer au point A. Dans le second cas, dans le triangle ABC,
les caractéristiques représentées ne sont pas rectilignes (même si elles peuvent le sembler sur
la représentation), mais le sont ailleurs. Le second cas peut être vu comme la limite du
premier quand la distance entre les points A et C du premier cas devient très petite.
Dans les deux premières sections, on suppose que la solution est connue, c’est-à-dire que
l’on a la forme de la tuyère avec un écoulement sonique au col et supersonique uniforme en
sortie à Mach 2.4. On va analyser ces écoulements puis concevoir graphiquement la tuyère
courte.

3.1 Tuyère longue (8/22)


I. Identifier les zones d’ondes simples et celles qui ne le sont pas (on se contentera de
travailler sur la zone au dessus de la ligne de symétrie). Justifier votre réponse.
II. En utilisant la conservation des invariants de Riemann et en suivant la ligne pleine
dans le cas de la première tuyère (qui se trouve être alternativement une C− ou une

7
D

C Theta_max

A
Ecoulement uniforme
E

Mach 2.4
B

F
D

Theta_max_C
A
Ecoulement de sortie

E
uniforme à Mach 2.4
B

Figure 2 – Conception de tuyère longue à écoulement très contrôlé pour une soufflerie superso-
nique à gauche et de tuyère courte pour la propulsion spatiale à droite. Les lignes sont des C+ ou
des C− ; dans la figure du bas, les C+ sont en trait plein et les C− en trait pointillé.

8
C+ ), évaluer l’évolution du nombre de Mach à chaque point de rebroussement, qu’il
soit au centre ou sur le bord de la tuyère (indication : on commencera pas le point D,
pour lequel l’écoulement est horizontal, puis on arrivera à E, puis à C et ainsi de suite
en utilisant les invariants).
III. On suppose que la forme de la partie expansion AC pour la première tuyère est
fixée. Expliquer comment on peut alors résoudre l’ensemble de l’écoulement dans la
zone ACE.
IV. Représenter les C+ dans le triangle curviligne DEC et proposer une façon de construire
le bord de la tuyère, l’écoulement étant connu sur la courbe C− : CE, par approximation
successive en discrétisant la courbe CE.

3.2 Tuyère courte (6/22)


On travaille dans un voisinage du point A et on considère ici la détente de Prandtl-Meyer
centrée au point A qui fait passer l’écoulement sonique de vitesse horizontale à un écoulement
d’angle Θmax−C positif.
I. En utilisant les invariants de Riemann, donner la valeur de la fonction de Prandtl-
Meyer juste en aval de la détente, c’est-à-dire pour un écoulement faisant un angle de
Θmax−C avec l’horizontale.
II. En supposant que la C− issue de A qui passe par E aboutit à un écoulement uniforme
horizontal à Mach 2.4 en E et dans le dièdre DEF , déterminer l’angle Θmax−C (la
valeur de la fonction de Prandtl-Meyer à Mach 2.4 est 36.75 degrés).
III. Identifier les zones d’onde simple dans cette nouvelle tuyère et indiquer précisément
où sont les caractéristiques qui sont des droites et les caractéristiques qui sont des
courbes.

3.3 Construction graphique de la tuyère courte (8/22)


On se propose de construire graphiquement la tuyère à partir d’une discrétisation de la
détente de Prandtl-Meyer en 7 C− , une première correspondant à un angle d’écoulement de
θ = 0.375 degré, qui aboutit au point numéroté 1 sur la ligne de symétrie t pour laquelle on
vérifiera que l’angle avec la verticale est de −73, 725 degré, puis 6 autre faisant un sixième
de l’angle restant pour arriver à l’angle Θmax−C calculé ci-dessus.
I. Evaluer le nombre de Mach en ce point 1.
II. On prendra ensuite une discrétisation du reste de la détente de Prandtl-Meyer et on
déterminera successivement les intersections de ces C− avec la C+ issue du point 1,
donnant ainsi les point 2 à 7 sur la figure 3.
III. La C− partant de A et aboutissant au point 7 marque la fin de la détente de Prandtl-
Meyer. Expliquer pourquoi l’écoulement au point 7 et 8 sont identiques.

9
? 35 = D

15
A 8

7
56
4
3
2

1 9 34 = E
Figure 3 – Début de construction graphique de la tuyère courte.

Le point 8 est obtenu comme intersection de la C+ rectiligne issue de 7 et de la droite de


pente 1/2(Θmax−C + θ8 ) = 1/2(Θmax−C + θ7 ) issue de A.
IV. On construit ainsi le bord de la tuyère. Expliquer comment on va itérer le processus
et aboutir au point 34, qui correspond à une approximation du point E, puis au point
35, qui correspond à une approximation du point D.
V. En fait au point 1, il y a une légère inconsistance car on devrait avoir un angle de
l’écoulement de zéro, alors que l’on a choisit 0.375 degrés dans le but de l’approximation
de la solution. Proposer une explication pour justifier cette très légère inconsistance.

10
ÉCOULEMENTS COMPRESSIBLES ET SUPERSONIQUES
Contrôle de connaissance - 16 mars 2012

Adam Larat

Durée 3 heures (Tous les documents papier sont autorisés, ainsi que l’usage de calculatrices. Tous
autres appareils électroniques sont interdits).

Le contrôle de connaissance se compose de deux questions de cours, de deux exercices portant sur
les ondes de choc normales et sur les chocs obliques et d’un problème long. Pour toutes les questions,
le fluide considéré est un gaz idéal polytropique et = 1.4. Tous les angles sont donnés en degré,
n’oubliez pas de changer le mode de calcul des fonctions trigonométriques de vos calculatrices.

1 Questions de Cours
1.1 Question 1 : le tube à choc
On s’intéresse à un tube à choc, dans l’approximation d’un écoulement instationnaire monodi-
mensionnel. On considère un état initial, à l’instant t=0, constitué de deux gaz au repos de chaque
coté d’une membrane que l’on retire à l’instant initial. Ces gaz sont des fluides parfaits idéaux po-
lytropiques non miscibles notés “état WG ”, en ce qui concerne l’état à gauche de la membrane (de
pression pG , de température TG , de vitesse uG = 0, de masse molaire MG et de coefficient G ) et
“état WD ”, à droite de la membrane (de pression pD pG , de température TD = TG , de vitesse
uD = 0, de masse molaire MD MG et de coefficient D ). Cette donnée initiale correspond au
problème de Riemann traité en cours.
1) Rappeler la structure du diagramme des caractéristiques dans le plan (x, t) et justifier sa
structure centrée. Rappeler par quels types d’ondes sont connectés les divers états constants
obtenus dans la résolution de ce problème de Riemann. Dans quelle direction se déplace l’in-
terface entre les deux fluides non-miscibles ?

1.2 Question 2 : le Bell XS1


Le premier vol supersonique - officiellement reconnu comme tel - remonte au 14 octobre 1947.
Ce jour-là, le Bell XS1 piloté par Chuck Yeager atteignait M = 1.06 et ouvrait la voie à un nouveau
régime de vol, dans lequel les ondes de choc dominent l’écoulement. À M = 1.06, on a alors observé
un arc de choc détaché du fuselage en amont du nez de l’avion profilé comme une balle de fusil. Un
peu plus tard dans le courant de l’année suivante, le même pilote dans le même avion atteignait
M = 1.45. Dans cette configuration de vol, un choc oblique attaché au nez de l’avion se formait
(voir Figure 1).
2) Expliquer la différence entre la position des ondes de chocs dans les deux régimes de vol.

1
Figure 1 – Illustration des deux configurations de vol du même avion Bell XS1 entre son premier
vol à Mach 1.06 (gauche) et un vol l’année suivante à Mach 1.45 (droite).

2 Exercice : Mesure de vitesse supersonique par un tube de Pitot


On considère une navette spatiale en phase de rentrée atmosphérique. L’écoulement autour de
la navette est supersonique. Sur le nez de cette navette se trouve un tube de Pitot qui mesure la
pression d’arrêt pA
0 au point A (voir Figure 2). On suppose en outre qu’un thermocouple permet de
mesurer la température T0A en ce même point (température d’arrêt). En régime supersonique, un
choc détaché se forme en amont du tube ; on approchera cet discontinuité en arc par un choc droit.
Un capteur de pression situé en B permet de mesurer la pression statique pB avant ce choc. Nous
allons montrer que ce dispositif permet de connaître la vitesse de la navette en suivant la démarche
suivante :

Figure 2 – Mesure de vitesse en régime supersonique par un tube de Pitot. La courbe de choc
devant le tube de Pitot est assimilée à un choc droit.

1) En introduisant l’état B et l’état A respectivement en amont et en aval du choc, poser le problème


en termes de grandeurs connues et inconnues.
2) On souhaite montrer que la vitesse de la navette peut s’exprimer en fonction du nombre de Mach
amont MB comme unique inconnue.
a) Justifier le fait que la température d’arrêt se conserve à travers le choc. Rappeler alors la
relation entre la température d’arrêt et la température statique.

2
b) En déduire la célérité du son cB en amont du choc et exprimer la vitesse de la navette en
fonction de MB comme unique inconnue.
3) Montrer que le rapport entre la pression totale en aval du choc et de la pression statique en
amont du choc peut être exprimé comme fonction du Mach amont uniquement. On mettra cette
relation sous la forme suivante :

 1
pA
0 +1 2 ( + 1)2 MB2 1
= MB
pB 2 2 MB2 ( 1)

Note : dans le cas où MB > 1, on montre que la relation ci-dessus est inversible. La table 1
ci-dessous fournit des valeurs tabulées des nombres de Mach MA et MB en fonction du rapport
(pression totale en aval du choc) / (pression statique en amont du choc).
4) Calculer la vitesse de la navette pour les données mesurées suivantes : pression statique en B
égale à 2.83 Pa, pression d’arrêt en A égale à 92.42 Pa, température d’arrêt en A égale à 1200 K.
À cet instant du vol, la constante r de l’air, définie comme le rapport de la constante universelle
des gaz parfaits, R, et de la masse molaire de l’air dans les conditions considérées, est égale à
280 J.kg 1 K 1 .

pA
0
MB MA
pB
4.50 0.42355 26.539
4.60 0.42168 27.710
4.70 0.41992 28.907
4.80 0.41826 30.130
4.90 0.41670 31.379
5.00 0.41523 32.653
5.10 0.41384 33.954
5.20 0.41252 35.280
5.30 0.41127 36.631
5.40 0.41009 38.009
5.50 0.40897 39.412

Table 1 – Table de ratio entre la pression totale en aval et la pression statique en amont du choc
pour = 1.4

3 Exercice : Optimisation d’une entrée d’air en écoulement super-


sonique
On peut montrer que la poussée d’un réacteur dépend fortement des pertes de pression généra-
trice p0 occasionnées par le système de chocs présents dans l’entrée d’air. On introduit de ce fait un
paramètre ⌘ = p00 /p0 , rapport entre la pression génératrice à la fin de l’entrée d’air et la pression
génératrice initiale. On souhaite maintenant maximiser la valeur de ⌘. Dans les questions suivantes,

3
l’écoulement amont sera caractérisé par un nombre de Mach de M = 2 et l’écoulement à l’intérieur
de la conduite sera toujours supposé subsonique.

Figure 3 – Trois configurations pour une entrée d’air dans un écoulement supersonique. Sur chacune
de ces figures, les lignes pleines bleues désignent la structure, les pointillés longs rouges les chocs
normaux et les petits pointillés rouges les chocs obliques. À gauche, cas du choc droit. Au milieu,
avec une rampe de pente de manière à obtenir un choc faible à l’avant de l’entrée. À droite, une
proue à deux pentes 0 et 1 .

Questions :
1) On se place dans le cas d’une entrée d’air de type Pitot (cf. Figure 3, en haut), où l’on observe
simplement un choc normal dans la section d’entrée du dispositif. Ce choc sépare deux états W
et W 0 de pression statique, pression d’arrêt et nombre de Mach respectifs p, p0 , M et p0 , p00 , M 0 .
a) Montrer que ⌘ s’exprime en fonction du Mach M uniquement. On mettra si possible cette
relation sous la forme :

✓ ◆ 1
( + 1) +1 M 2 1
⌘(M ) = (1)
(2 + ( 1)M 2 ) (2 M 2 ( 1))

b) Que vaut ⌘ dans le cas où M = 2 ?


2) On souhaite diminuer la violence du choc. Pour ce faire, on utilise un système à proue émergente
déflectant l’écoulement amont d’un angle , voir Figure 3 au milieu. On suppose que l’onde de
choc oblique ainsi obtenue vient "taper" l’entrée d’air pile en son sommet, comme illustré sur la
Figure 3.
a) En utilisant l’abaque de choc oblique donnée en annexe, lire l’angle du choc oblique obtenu
pour des angles de déflexion = 10, 16 et 20 (on se contentera d’une précision au degré
près). Pourquoi avoir choisi cette valeur plutôt que l’autre valeur possible ?
b) Pour chacun de ces angles de déflexion, calculer M ⇤ , le nombre de Mach de l’état intermédiaire
compris entre le choc oblique et le choc normal d’entrée du dispositif.
c) Pour chacune de ces configurations calculer alors ⌘. Conclure.
3) On souhaite améliorer encore le dispositif. On adopte une double proue émergente, avec 0 = 14
et 1 = 8 , voir Figure 3, en bas. On obtient ainsi deux états intermédiaires W ⇤ et W ⇤⇤ de
nombres de Mach respectifs M ⇤ et M ⇤⇤ séparés par un choc oblique. W ⇤ est séparé de W par
un premier choc oblique, tandis que W ⇤⇤ est séparé de W 0 par un choc normal.

4
a) À nouveau en lisant sur l’abaque donnée en annexe, donner les angles des chocs obliques 0
et 1 , ainsi que les nombres de Mach M ⇤ et M ⇤⇤ .
Note : Pour 1 on considérera une interpolation linéaire entre les valeurs obtenues pour les
deux abaques les plus proches de celle de M ⇤ .
b) Cette nouvelle configuration conduit-elle à un gain sur ⌘ ?
4) Comment pourrait-on améliorer encore le dispositif ? Quel point a-t-on négligé ? Qui peut réduire
considérablement les gains de ce système ?

4 Problème : Optimisation d’un profil d’aile dans un écoulement


supersonique
Considérons un profil en losange tel que celui défini par la figure 4 ci-dessous. Dans toute la
suite du problème on supposera A1 fixe en (0, 0), A2 fixe en (c/2, t/2). On prendra A3 mobile le
long de ( 1 ) et A4 fixe en (c/2, t/2).

Figure 4 – Profil d’aile losangique.

Plongeons maintenant ce profil dans un écoulement supersonique de fluide parfait. Il se forme :


• un choc attaché au bord d’attaque A1 du profil,
• une détente supersonique à partir du sommet A2 et une autre à partir du sommet A4 ,
• un choc attaché au bord de fuite A3 .

Question 1 : Justifier la structure de l’écoulement décrit précédemment et illustré sur la Fi-


gure 5. Pourquoi de telles ondes se forment ?

On rappelle que la résultante des forces aérodynamiques exercées sur le profil est donnée par :
Z
~
F = p.~ndl
profil

où p désigne la pression, ~n le vecteur unitaire normal au profil et pointant vers l’extérieur de celui-ci
et dl est un élément de longueur infinitésimal. F~ se décompose classiquement sur la base orthonormée
(~i, ~j) en une composante de traînée D et une composante de portance L :

F~ = D~i + L~j

5
Figure 5 – Écoulement supersonique autour du profil.

Ces grandeurs sont généralement normalisées pour obtenir les coefficients de portance et de traînée
sans dimension :
L D
CL = et CD =
q1 c q1 c
avec q1 = pM12 /2 et = 1.4 désigne le rapport des chaleurs spécifiques. Optimiser les performances
aérodynamiques du profil va alors consister à rendre le rapport CD /CL le plus proche possible de
zéro, de façon à minimiser la résistance à l’avancement exprimée par la traînée, tout en maximisant
la portance favorable au vol. Notons que dans le cadre du présent problème où les effets visqueux
sont négligés (fluide parfait), la traînée du profil est réduite à la traînée d’onde.

Question 2 : Dans le cas où A3 est en (c, 0) et A4 en (c/2, t/2), calculer l’expression exacte
de CD et CL en fonction de p2 , p3 , q1 et (t/c). Compte tenu de la valeur de CD /CL , commenter
les performances aérodynamiques de ce profil.

Question 3 : On cherche maintenant à modifier le profil de façon à améliorer le résultat précé-


dent. On déplace le point A3 le long de ( 1 ) jusqu’à la position (c, ↵t), où ↵ 2 [0; 1], comme
indiqué sur la Figure 6 ci-dessous.

Figure 6 – Profil Modifié. A3 peut descendre jusqu’à la coordonnée (c, t), c’est-à-dire qu’A3 peut
se retrouver plus bas qu’A4 .

6
a) On définit les rapports de pression : r = p2 /p3 , s = p03 /p2 . Montrer que :
CD t 2r (2↵ + 1) rs(1 2↵)
= '(r, s, ↵) ⇥ , avec '(r, s, ↵) =
CL c rs 1
b) Montrer que, dans le cas particulier où ↵ = 1, on a exactement : '(r, s, ↵) = 3.
c) Justifier le fait qu’à travers les deux détentes en A2 et A4 , on peut écrire les relations :
! !
1 1
1 1+ (M2 )2 1 1+ (M20 )2 1
= 2
1 et s= 2
1
0 2
(2)
r 1+ 2 (M3 )2 1+ 2 (M3 )

d) Décrire dans le cas général la procédure de calcul de r et s en supposant M1 , p1 , c, t et ↵


connus. On donnera toutes les formules employées, même lorsque celles-ci sont implicites
en la quantité cherchée. On détaillera successivement :
– l’obtention de p2 , M2 en fonction de l’état 1,
– l’obtention de M3 en fonction de l’état 2, puis celle de r donc de p3 ,
– l’obtention de M30 en fonction de l’état 20 puis celle de s donc de p03 .
Note : dans ce dernier cas, on fera bien attention à la nature des caractéristiques
formant l’éventail de détente en A4 .
e) On suppose maintenant : M1 = 1.7, p1 = 1, c = 1 et t = 7/100 (on considère des quantités
p
adimensionnées). Calculer '(r, s, ↵) pour ces conditions et ↵ = 1/2 puis ↵ = 1/ 2. À
chaque étape intermédiaire, on maintiendra une precision à 4 chiffres après la virgule.
Note : on utilisera l’abaque de choc oblique et la table de la fonction de Prandtl-Meyer
pour = 1.4, fournies en annexe. Pour la tabulation de la fonction de Prandtl-Meyer, on
n’hésitera pas à faire une interpolation linéaire autour des deux valeurs les plus proches
de la valeur cherchée.
Compte tenu de b), que peut-on supposer quant à l’évolution du rapport CD /CL ?

Question 4 : On se place maintenant dans le cadre de la théorie des profils minces pour préciser
le comportement de CD /CL en fonction de ↵. Cette théorie donne une approximation des
coefficients de portance et de traînée quand l’epaisseur t du profil est très petite devant la
corde c. Si yU (x) (resp. yL (x)) désigne l’équation de la partie supérieure (resp. inférieure) du
profil, c’est-à-dire de la ligne A1 A2 A3 (resp. A1 A4 A3 ), on a alors dans l’approximation des
profils minces :
8 Z ✓ ◆
> 2 1 c dyU dyL
>
> C = p + dx
> L
> M12 1 c 0 dx dx
<
(3)
>
> Z c ✓ ◆2 ✓ ◆2 !
>
> 2 1 dy U dy L
>
: CD = p 2 + dx
M1 1 c 0 dx dx

a) On considère à nouveau le profil de la figure 6. Exprimer à l’aide des relations (3) ci-
dessus, CD /CL en fonction de ↵ et de (t/c).
p
b) Montrer que, pour (t/c) fixé, CD /CL est minimal pour ↵ = 1/ 2. Comparer alors la
valeur de CD /CL à celle obtenue en 2). Commenter.

7
Figure 7 – Abaque de chocs oblique donnant l’angle de déflexion en fonction de l’angle du choc.
8
M ⌫(M ) M ⌫(M ) M ⌫(M ) M ⌫(M ) M ⌫(M ) M ⌫(M )
1.00 0.0000 1.30 6.1703 1.60 14.8604 1.90 23.5861 2.40 36.7465 3.00 49.7573
1.01 0.0447 1.31 6.4449 1.61 15.1561 1.91 23.8687 2.42 37.2289 3.02 50.1417
1.02 0.1257 1.32 6.7213 1.62 15.4518 1.92 24.1506 2.44 37.7077 3.04 50.5231
1.03 0.2294 1.33 6.9995 1.63 15.7473 1.93 24.4318 2.46 38.1831 3.06 50.9016
1.04 0.3510 1.34 7.2794 1.64 16.0427 1.94 24.7123 2.48 38.6551 3.08 51.2771
1.05 0.4874 1.35 7.5607 1.65 16.3379 1.95 24.9920 2.50 39.1236 3.10 51.6497
1.06 0.6367 1.36 7.8435 1.66 16.6328 1.96 25.2711 2.52 39.5886 3.12 52.0195
1.07 0.7973 1.37 8.1276 1.67 16.9276 1.97 25.5494 2.54 40.0503 3.14 52.3864
1.08 0.9680 1.38 8.4130 1.68 17.2220 1.98 25.8269 2.56 40.5085 3.16 52.7505
1.09 1.1479 1.39 8.6995 1.69 17.5161 1.99 26.1037 2.58 40.9633 3.18 53.1118
1.10 1.3362 1.40 8.9870 1.70 17.8099 2.00 26.3798 2.60 41.4147 3.20 53.4703
1.11 1.5321 1.41 9.2756 1.71 18.1034 2.02 26.9295 2.62 41.8628 3.22 53.8261
1.12 1.7350 1.42 9.5650 1.72 18.3964 2.04 27.4762 2.64 42.3074 3.24 54.1791
1.13 1.9445 1.43 9.8553 1.73 18.6891 2.06 28.0197 2.66 42.7488 3.26 54.5294
1.14 2.1600 1.44 10.1464 1.74 18.9814 2.08 28.5600 2.68 43.1868 3.28 54.8770
1.15 2.3810 1.45 10.4381 1.75 19.2732 2.10 29.0971 2.70 43.6215 3.30 55.2220
1.16 2.6073 1.46 10.7305 1.76 19.5646 2.12 29.6308 2.72 44.0529 3.32 55.5643
1.17 2.8385 1.47 11.0235 1.77 19.8554 2.14 30.1613 2.74 44.4810 3.34 55.9040
1.18 3.0743 1.48 11.3169 1.78 20.1458 2.16 30.6884 2.76 44.9059 3.36 56.2411
1.19 3.3142 1.49 11.6109 1.79 20.4357 2.18 31.2121 2.78 45.3275 3.38 56.5756
1.20 3.5582 1.50 11.9052 1.80 20.7251 2.20 31.7325 2.80 45.7459 3.40 56.9075
1.21 3.8060 1.51 12.1999 1.81 21.0139 2.22 32.2494 2.82 46.1611 3.42 57.2369
1.22 4.0572 1.52 12.4949 1.82 21.3021 2.24 32.7629 2.84 46.5731 3.44 57.5639
1.23 4.3117 1.53 12.7901 1.83 21.5898 2.26 33.2730 2.86 46.9820 3.46 57.8883
1.24 4.5694 1.54 13.0856 1.84 21.8768 2.28 33.7796 2.88 47.3877 3.48 58.2102
1.25 4.8299 1.55 13.3812 1.85 22.1633 2.30 34.2828 2.90 47.7903 3.50 58.5298
1.26 5.0931 1.56 13.6770 1.86 22.4492 2.32 34.7825 2.92 48.1898 3.52 58.8469
1.27 5.3590 1.57 13.9728 1.87 22.7344 2.34 35.2787 2.94 48.5863 3.54 59.1616
1.28 5.6272 1.58 14.2686 1.88 23.0190 2.36 35.7715 2.96 48.9796 3.56 59.4739
1.29 5.8977 1.59 14.5645 1.89 23.3029 2.38 36.2607 2.98 49.3700 3.58 59.7838

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

9
Formulaire
Équations d’Euler :
0 1 0 1
⇣! ⌘ ⇢ ⇢ui
@U !
+ div F (U) = 0, U = @ ⇢~u A et F = (F1 , F2 ) , avec Fi = @ ⇢ui ~u + p i A
@t
⇢E (⇢E + p)ui

Équation d’état des Gaz ✓ Parfaits


◆ r : p =p⇢rT . Gaz idéaux polytropiques : p = ( 1)⇢e
@p p
Vitesse du son : c2 = = = rT Nombre de Mach : M = u/c.
@⇢ s ⇢
2 2 1 2 2M⇤2 ( + 1)M 2
Variables soniques : c2⇤ = c + u ; M2 = 2
; M⇤2 =
+1 +1 +1 ( 1)M⇤ 2+( 1)M 2
1
Relation de Prandtl : (M⇤ )1 = ;
(M⇤ )2
✓ ◆ ✓ ◆
T0 1 2 p0 T0 1 ⇢0
Grandeurs d’arrêt : =1+ M ; = =
T 2 p T ⇢
Relation de saut dans un choc droit :
1
1+ M12 p2 2 ⇢2 u1 ( + 1)M12
M22 = 2
1; =1+ (M 2 1); = =
M12 2
p1 +1 1 ⇢1 u2 2+( 1)M12
✓ ◆
T2 2( 1)(M12 1)( M12 + 1) p2 ⇢1
=1+ ; s2 s1 = CV ln ; s 2 > s1
T1 ( + 1)2 M12 p1 ⇢2
dA du
Équation de Hugoniot : + (1 M 2 ) =0
A u
2c
Écoulement 1D instationnaires. Invariants de Riemann : J0 = s, J± = v ±
0 1 1
✓ ◆ G 1
2cG @ p1 2 G
1-Détente : v1 > vG , p1 < pG , v1 = vG 1A
G 1 pG
p
2cG p 1 pG
1-Choc : v1 < vG , p1 > pG , v1 = vG p p
pG ( + 1)p1 + ( 1)p1G
0
✓ ◆ D 1
2cD @ p1 2 D
3-Détente : v1 < vD , p1 < pD , v1 = vD + 1A
D 1 p D
p
2cD p1 pD
3-Choc : v1 > vD , p1 > pD , v1 = vD + p p
pD ( + 1)p1 + ( 1)pD
2-Onde : v1 = v2 , p1 = p2
Écoulement 2D permanents :
Chocs Obliques :

1 +1 M12 2 2 2+( 1)M12 sin2
= tan 1 ; M 2 sin ( ') =
tan ' 2 M12 sin 1 2 M12 sin2 + 1
p2 M2 p1 v2 v1 '
Ondes de Mach : ⇡ p 2 1 '; ⇡ tan ↵ ' = p 2
p1 M1 1 v 1 M1 1
Angle de Mach : ↵(M ) = arcsin(1/M ); Pente des C± caractéristiques : ✓ ± ↵(M ).

10
Tableau Récapitulatif :
9
@C0 p + ⇢U @C0 U = 0 sur C0 trajectoire = ✓ ◆
dy v
=
; dx u
@ C0 s = 0 sur C0 trajectoire C0

p
M2 1
@ C± ✓ ± @ C± p = 0
⇢U 2
ou
@ C± U T @ C± s @ C± H
± tan(↵) @C± ✓ = 0,
U U2
sur les C± caractéristiques de pente locale a± faisant un angle ±↵ avec la
direction de l’écoulement.

Écoulement homoentropique et homoénergétique :

J± = ±✓ (U ) = constante, , R± = ±✓ f (M ) = constante, sur C±

Fonction de Prandtl-Meyer :
r r hp i
+1 1
f (M ) = arctan (M 2 1) arctan M2 1
1 +1
r (4)
r
⇡ +1 1
= ↵(M ) + arctan (M 2 1)
2 1 +1

11
2
ECOULEMENTS COMPRESSIBLES ET SUPERSONIQUES
Contrôle de connaissance

Marc Massot, Adam Larat

Durée 3 heures (Les notes de cours, y compris sur support informatique, et les documents
distribués en cours sont autorisés ainsi que l’usage de calculatrices).

Le contrôle de connaissance se compose de deux exercices et d’un problème ; ces trois parties
sont indépendantes. Le contrôle est volontairement long et permet de couvrir l’ensemble des points
traités en cours ; il ne sera pas nécessaire de tout faire pour obtenir un très bon résultat. La qualité
de la rédaction sera un point important pour l’évaluation. Dans les exercices et le problème, on
considère des écoulements de gaz parfait idéaux polytropiques (pas de phénomènes de dissipation,
loi d’état de gaz idéal, p = ρ r T , r = R/M où M est la masse molaire du gaz ; cp ainsi que γ sont
des constantes données) sans force massique extérieure.

1 Exercice de cours : Piston en mouvement dans un tube rempli


par deux gaz
Cet exercice est très proche des questions traitées en cours. On s’intéresse à un tube infini à
symétrie de révolution autour de son axe (x′ Ox) de section constante, dans l’approximation d’un
écoulement instationnaire monodimensionnel. On considére un état initial, à l’instant t=0, constitué
de deux gaz au repos de chaque coté d’un piston auquel on va imposer un mouvement à partir de
l’instant initial t = 0. Ces gaz sont des fluides parfait idéaux polytropique notés “état 4” en ce qui
concerne l’état à gauche du piston (de température T4 , de pression p4 , de vitesse u4 = 0 de masse
molaire M4 et de coefficient adiabatique γ4 ) et “état 1” à droite du piston (de température T1 = T4 ,
de pression p1 , de vitesse u1 = 0, de masse molaire M1 et de coefficient adiabatique γ1 ).
Le mouvement imposé au piston est décrit par sa position X(t) et le piston se trouve en x = 0
à l’instant t = 0 (on fait l’hypotèse simplificatrice que le piston est un piston d’épaisseur infiniment
fine). L’expérimentateur applique pour tout instant t > 0 une force sur le piston de manière à lui
imposer la trajectoire X(t) = a tn , avec a > 0 et n un entier supérieur ou égal à un (n ≥ 1).
(1) Exprimer dans chaque gaz au repos les célérités du son c40 et c10 en fonction des grandeurs
de l’énoncé.
(2) Dans le cas n = 2 à quoi est homogène la constante a ? Et dans le cas n = 1 ?
On considère dans un premier temps le cas où n = 2

1
Gaz 4 piston en x=0 a t=0 Gaz 1
x’ x

Gaz 4 piston en X(t) a l’instant t Gaz 1


x’ x
O

Figure 1 – Tube rempli de deux gaz séparés par un piston à deux instants successifs.

(3) Effectuer, comme dans le cours, la résolution des propriétés de l’écoulement pour les deux
gaz, pour tout temps t < T ∗ , pour lequel l’ecoulement reste homoentropique dans chaque gaz.
Proposer une synthèse des principes de base permettant de trouver l’état des deux gaz. Sur quelles
grandeurs effectue-t-on la résolution ? Jusqu’à quel temps ce type de résolution basé sur le caractère
homoentropique des deux écoulements est-il valable ? Que se passe-t-il ensuite ?
(4) Tracer le diagramme des caractéristiques dans le plan (x, t). On fera une représentation
précise jusqu’au temps T ∗ puis une représentation qualitative pour les temps supérieurs. Identifier les
diverses zones d’écoulement de manière précise dans le diagramme. Une figure claire sera appréciée.
Proposer une représentation des solutions pour les deux gaz dans un diagramme de phase de type
pression-vitesse. Sur quelles courbes se trouvent les deux solutions ? Justifier votre réponse.
(5) Donner les expressions de la vitesse dans la direction x, de la vitesse du son c, de la pression et
de la température d’un côté et de l’autre au voisinage du piston. Quelle est la force totale appliquée
sur le piston par le gaz ? Si le piston est de masse M , trouver quelle est la force instantanée F (t)
appliquée par l’expérimentateur en appliquant la loi fondamentale de la dynamique.
(6) Dans ce cas n = 2, préciser l’expression du temps et de la position de formation du choc.
Montrer que cette position et ce temps correspondent à la première apparition d’une onde de choc
et au premier point de croisement des caractéristiques avec la caractéristique issue de (0, 0). Donner
l’expression de la vitesse en tout (x, t) pour un temps inférieur au temps de formation du choc.

(7) On se place alors dans le cas n = 1. Effectuer une représentation précise et justifiée du
diagramme des caractéristiques associé.
(8) Ce type de diagramme est proche de celui du tube à choc. Expliquer quelle est la différence
fondamentale entre le cas étudié dans cet exercice et le cas du tube à choc. Dans quel cas les deux
diagrammes sont-ils equivalents ? Représenter dans un diagramme de phase pression-vitesse les deux
cas, afin d’expliquer la différence de manière claire quand elle est présente. Quand les deux cas sont
équivalents, que dire la force appliquée par l’expérimentateur ?

2 Exercice : Optimisation d’un moteur et de sa poussée


On considère la tuyère convergente-divergente d’un moteur de type moteur de fusée. Cette tuyère
est telle que sa section au col est Ac = 1 cm2 . Les conditions en amont de la section d’entrée sont

2
les conditions d’arrêt : p0 = 100 atm, T0 = 2000 K. Le gaz de combustion est un gaz parfait tel que
r = 415.7 J.kg−1 K −1 et γ = 1.25. On supposera dans toute la suite que la tuyère est amorcée.
(1) Calculer ṁ = ρuA, le débit-masse de la tuyère.

La poussée du moteur fusée est donnée par la formule :


P = ṁue + Ae (pe − pa )
où ue désigne la vitesse dans la section d’éjection, pe est la pression dans la section d’éjection et pa
la pression ambiante. On admettra ici que la poussée est maximale quand la tuyère est adaptée à
la pression ambiante, c’est-à-dire quand la pression dans la section d’éjection est égale à la pression
ambiante, avec un écoulement isentropique dans toute la tuyère, y compris dans la section d’éjection.
Au cours d’un vol, la pression ambiante va naturellement varier avec l’altitude et il n’est donc pas
possible d’avoir une tuyère adaptée à tout moment. On se pose ici la question de savoir s’il est
préférable d’adapter la tuyère aux conditions au sol ou à grande altitude.
Supposons la pression ambiante connue et la tuyère adaptée.
(2) Donner les relations qui permettent de calculer le Mach Me , la température Te et la vitesse
ue dans la section d’éjection. Préciser également la relation permettant de calculer l’aire de la section
d’éjection.

(3) On considère tout d’abord le cas d’une tuyère adaptée à la pression ambiante au sol :
pa = 1 atm.
3-(a) Déterminer l’aire de la section d’éjection d’une telle tuyère.
3-(b) Calculer la poussée développée au sol et à une altitude de 23.5 km, à laquelle la pression
ambiante vaut pa = 0.03 atm.

(4) On considère maintenant le cas d’une tuyère adaptée à la pression ambiante à l’altitude de
23.5 km : pa = 0.03atm.
4-(a) Déterminer l’aire de la section d’éjection d’une telle tuyère.
4-(b) Calculer la poussée développée au sol et à une altitude de 23.5 km, à laquelle la pression
ambiante vaut pa = 0.03 atm.
4-(c) Montrer que pour une valeur de la pression ambiante très proche de 1atm, un choc droit
apparaît dans la section d’éjection de la tuyère.
4-(d) Calculer la poussée développée dans ces conditions (on assimilera le résultat obtenu à la
valeur de la poussée exactement au sol).

(5) Comparer les résultats obtenus en 3) et 4) - que l’on pourra regrouper dans un tableau du
type ci-dessous - et tirer une conclusion argumentée quant au choix de telle ou telle tuyère (tenant
compte des poussées prévues et de la dimension de la tuyère).

Tuyère Poussée au sol (N) Poussée à h = 23.5km (N) Aire de Ae (cm2 )


Adaptée au sol
Adaptée à h = 23.5km

3
Figure 2 – Dispositif de statoréacteur (par raison de symétrie, on ne présente que la moitié des lignes
caractéristiques de l’écoulement extérieur au dispositif).

3 Problème : Le statoréacteur
On considère un écoulement plan permanent dans un statoréacteur selon le dispositif présenté
sur la Figure 2.
A l’amont du dispositif, l’écoulement est supersonique, de nombre de Mach M∞ , de vitesse U∞ ex
parallèle à l’axe de symétrie du dispositif.
La partie interne du dispositif débute par une tuyère lentement divergente, que l’on appelle aussi
diffuseur. Sa forme a été calculée à l’aide de la méthode des caractéristiques pour qu’elle fournisse
en régime supersonique un écoulement parfaitement plan uniforme dans une section (de la section
d’entrée E à la section interne C à partir de laquelle la chambre est de section constante AC ).
L’épaisseur des parois de la chambre sera négligée dans toute l’étude.

3.1 Partie 1 : Ecoulement externe


Dans cette question, l’écoulement est supposé partout supersonique, et pour les applications
numériques, on prendra γ = 1.4 et M∞ = 2. On va analyser les principales lignes caractéristiques
de l’écoulement extérieur autout du dispositif. Compte tenu de la symétrie, on raisonne dans le
demi-plan supérieur.
(a) L’angle φ vaut 15o . Le graphe de la polaire de choc donne pour la déviation de l’écoulement
σ derrière l’onde de choc oblique attachée en O une valeur de σ = 45, 34o . Le nombre de Mach M1

4
p
et la pression p1 dans la région 1 valent M1 = 1.446 et p 1 = 2, 196. Détailler les équations qui

permettent d’obtenir ces valeurs.
(b) Expliquer comment on détermine la valeur de M3 dans la région 3, à l’aval du faisceau de
détente de Prandtl-Meyer (domaine 2) centré sur la point O′ et qui est schématisé par le secteur
(O′ IJ) sur la Figure 2 (on donnera sans application numérique les équations nécessaires).
p3
(c) Le calcul donne M3 = 1, 962. Exprimer alors le rapport p∞ où p3 désigne la pression dans la
p1
région 3, en fonction de p∞ , M1 et M3 . Calculer ce rapport et vérifier qu’il est légèrement supérieur
à 1.
(d) Calculer le rapport pp03 entre la pression génératrice p03 de la région 3 et la pression
0∞
p3
génératrice amont p0∞ , en fonction de p∞ , M3 et M∞ . On vérifiera que ce rapport est légèrement
inférieur à 1 et on expliquera pour quelle raison.
(e) On veut démontrer que l’écoulement dans la région 3 est partout uniforme dans le triangle
O′ JO′′ . Par construction on a placé l’extrémité O′′ du dispositif de manière à ce qu’elle se trouve à
l’amont de la caractéristique C− issue du point I, intersection de la caractéristique montante initiale
du faisceau 2 avec le choc oblique issue de O. Montrer alors à l’aide des invariants de Riemann que
ce principe garantit bien que l’écoulement est effectivement uniforme dans la région délimitée par
le triangle O′ JO′′ . Pour cela on montrera d’abord que l’écoulement est uniforme le long du segment
O′ O′′ . On montrera ensuite qu’il est uniforme dans le triangle O′ JO′′ .

3.2 Partie 2 : Ecoulement interne


Dans cette partie, on supposera que dans la section S indiquée sur la Figure 2 (plan de sortie
du dispositif) l’écoulement extérieur est le même qu’à l’infini amont : Mach M∞ = 2, vitesse U∞ ex ,
pression p∞ , pression et température génératrice p0∞ et T0∞ . Ces même conditions s’appliquent
aussi à la section d’entrée E du dispositif. On prendra finalement pour les applications numériques
γ = 1.4.
(a) Donner l’expression du débit massique dans une section du diffuseur en fonction de l’aire
de la section, des conditions génératrices et du nombre de Mach M∞ dans la section d’entrée E
du diffuseur de la Figure 2 et le nombre de Mach MC dans la section d’entrée C de la chambre.
Justiifier que MC > M∞ .
(b) La tuyère a été conçue de manière à obtenir un nombre de Mach MC = 2.77 à l’entrée de la
chambre. Donner l’expression du rapport de pression pp∞
C
où pC désigne la pression dans la chambre.
On vérifiera que ce rapport est inférieur à 1.
(c) Pour qu’à l’aval de la section S et à proximité de l’orifice de sortie du dispositif de la Figure
2, l’écoulement extérieur (U∞ ex , p∞ ) et l’écoulement issu de la chambre (UC ex , pC ) s’équilibrent,
il faut qu’ils vérifient deux conditions : il faut qu’ils soient parallèles et il faut que leurs pressions
soient identiques. D’après la question précédente, il n’y a pas identité des pressions dans le plan de
sortie.
c-(i) Montrer que dans la région proche de l’orifice de sortie, l’équilibre des deux écoulements
ne peut pas être assuré par la seule présence d’un choc droit attaché en O ′′ .
c-(ii) Montrer que la présence seule d’un choc oblique attaché en O ′′ n’est pas non plus suffi-
sante.

5
c-(iii) On admet que l’équilibre des deux écoulements dans la région aval proche de l’orifice de
sortie nécessite la présence de trois familles de lignes caractéristiques issues de O′′ : un choc
oblique attaché, une ligne de glissement et un faisceau de détente de Prandtl-Meyer. Faire un
schéma qualitatif indiquant la position de ces trois phénomènes autour de O ′′ .

3.3 Partie 3 : Propulsion


Comme dans la partie précédente, on supposera que dans la section S indiquée sur la Figure 2
(plan de sortie du dispositif) l’écoulement extérieur est le même qu’à l’infini amont : Mach M∞ = 2,
vitesse U∞ ex , pression p∞ , pression et température génératrice p0∞ et T0∞ . Ces même conditions
s’appliquent aussi à la section d’entrée E du dispositif.
On réalise maintenant dans la chambre du dispositif une combustion grâce au mélange avec l’air
et à l’ignition d’un carburant introduit par des injecteurs (non montrés sur la Figure). Le débit
d’injection du carburant est supposé négligeable devant le débit du dispositif. On va étudier l’effet
de cette combustion sur le fonctionnement du dispositif.
On va montrer que le réchauffement du fluide consécutif aux réactions de combustion exo-
thermiques est équivalent à l’effet d’un changement de section du conduit interne du dispositif. Le
processus de combustion est une source d’entropie le long des lignes de courant et doit être modélisé.
(a) Rappeler les équations de Hugoniot établies dans le cours dans le cas d’un écoulement
quasi-monodimensionel homoentropique et dans le cas de la présence d’une source de chaleur.
(b) Quel est l’effet d’un apport de chaleur sur la vitesse de l’écoulement dans la chambre du
dispositif de la Figure 2 si l’écoulement est supersonique et continu ? Comment peut-on simuler cet
effet sans apport de chaleur.
(c) On s’intéresse au bilan des efforts sur les parois du dispositif en l’absence d’apport de chaleur
par la combustion. En considérant d’une part la pression qui règne sur la rampe extérieure OO ′ du
diffuseur, dérrière le choc oblique attaché en O et d’autre part, la pression sur la paroi interne du
diffuseur où l’écoulement est supposé continu, expliquer sans équation dans quelle direction s’exerce
la force résultante qui agit sur le dispositif en l’absence d’apport de chaleur par la combustion.
(d) Expliquer comment une combustion dans la chambre du dispositif peut changer la conclusion
de la question (c).

3.4 Partie 4 : Analogie avec un système à deux cols soniques


On considére maintenant le dispositif de la Figure 3 où la chambre de la section constante de
la Figure 2 a été remplacée par un bloc convergent de section variable, dont l’aire de la section est
comprise entre AC et AS . On fera les applications numériques avec γ = 1.4, M∞ = 2 et MC = 2.77.
Remarque : Comme schématisé sur la Figure 3, le bloc convergent est non solidaire du dispositif
et ne transmet donc par les efforts qu’il subit.
(a) Quand on ferme progressivement la section aval AS , l’expérience montre que, pour un taux
AS
de contraction A E
inférieur à une certaine valeur, un choc droit permanent s’établit dans le diffuseur.
′ AS
On note p0∞ la pression d’arrêt à l’aval de ce choc. Déterminer l’expresion de A E
en fonction de
p′
M∞ , de γ et du rapport p0∞0∞
(M 2 , γ), avec M le nombre de Mach devant le choc droit, dans le cas
où la section aval S de section AS est sonique et joue le rôle d’un “second col” dans le dispositif.

6
Figure 3 – Dispositif analogue à celui de la Figure 2 simulant l’effet d’une combustion dans la chambre
(par raison de symétrie, on ne présente que la moitié des lignes caractéristiques de l’écoulement extérieur au
dispositif).

p′
(b) Le rapport p0∞0∞
(M 2 , γ) est une fonction monotone décroissante du nombre de Mach M
devant le choc. En déduire comment se déplace le choc droit dans le dispositif quand on ferme
progressivement la section AS à partir de la valeur considérée dans la question précédente.
(c) A partir des conclusions de la question (b), décrire comment évolue, qualitativement, le bilan
des efforts (établi dans la partie précédente sans bloc convergent) en présence du bloc convergent
quand on ferme progressivement la section AS à partir de la valeur considérée dans la question (a).
(d) Dans quelle condition la résultante des efforts de pression qui s’appliquent sur la paroi
interne du diffuseur de la Figure 3 est-elle maximale ? Pourquoi peut-on affirmer que dans ce cas le
bilan des efforts est propulsif, c’est-à-dire dirigé vers l’amont ?
(e) Que se passe-t-il si on diminue encore AS ?

7
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.

8
ECOULEMENTS COMPRESSIBLES ET SUPERSONIQUES
Contrôle de connaissance

Marc Massot

Durée 3 heures (Les notes de cours et les documents distribués en cours


sont autorisés ainsi que l’usage de calculatrices).

Le contrôle de connaissance se compose d’un exercice et d’un problème. L’exercice se compose


de quatre sous-parties (les parties 1.1 et 1.2 peuvent être traitées de manière indépendante). Dans
l’exercice et le problème, on considère des écoulements de gaz parfait idéaux polytropiques (pas de
phénomènes de dissipation, loi d’état de gaz idéal, p = ρ r T , r = R/M où M est la masse molaire
du gaz ; cp ainsi que γ sont des constantes données).

1 Exercice : Piston en mouvement dans un tube


Cet exercice est très proche des questions traitées en cours. On s’intéresse à un tube infini à
symétrie de révolution autour de son axe (x′ Ox) de section constante, dans l’approximation d’un
écoulement instationnaire monodimensionnel. On considére un état initial, à l’instant t=0, constitué
de deux gaz au repos de chaque coté d’un piston auquel on va imposer un mouvement à partir de
l’instant initial t = 0 (Voir Figure 1). Ces gaz sont des fluides parfait idéaux polytropique notés
“état 4” en ce qui concerne l’état à gauche du piston (de température T4 , de vitesse u4 = 0 de masse
molaire M4 et de coefficient γ4 ) et “état 1” à droite du piston (de température T1 = T4 , de vitesse
u1 = 0, de masse molaire M1 et de coefficient γ1 ). Le mouvement imposé au piston est décrit par
sa position X(t) et le piston se trouve en x = 0 a l’instant t = 0 (on fait l’hypotèse simplificatrice
que le piston est un piston d’épaisseur infiniment fine). L’expérimentateur applique à l’instant t > 0
une force sur le piston de manière à lui imposer la trajectoire X(t) = a tn , avec a < 0 et n un entier
supérieur ou égal à un (n ≥ 1).

1.1 Cas des deux gaz dans un tube infini avec accéleration progressive
(1) Exprimer dans chaque gaz au repos les célérités du son c40 et c10 en fonction des grandeurs
de l’énoncé.
On considère dans un premier temps le cas où n ≥ 2
(2) Effectuer, comme dans le cours, la résolution des propriétés de l’écoulement pour les deux
gaz, pour tout temps t < T ∗ , pour lequel l’ecoulement reste homoentropique dans chaque gaz.
Proposer une synthèse des principes de base permettant de trouver l’état des deux gaz. Sur quelles
grandeurs effectue-t-on la résolution ? Jusqu’à quel temps ce type de résolution basé sur le caractère
homoentropique des deux écoulements est-il valable ? Que se passe-t-il ensuite ?

1
Gaz 4 piston en x=0 a t=0 Gaz 1
x’ x

Gaz 4 Gaz 1
x’ x
O
piston en X(t) a l’instant t

Figure 1 – Tube rempli de deux gaz séparés par un piston à deux instants successifs.

(3) Tracer le diagramme des caractéristiques dans le plan (x, t). On fera une représentation
précise jusqu’au temps T ∗ puis une représentation qualitative pour les temps supérieurs. Identifier les
diverses zones d’écoulement de manière précise dans le diagramme. Une figure claire sera appréciée.
(4) Donner les expressions de la vitesse dans la direction x, de la vitesse du son c, de la pression
et de la température de chaque coté du piston. Quelle est la force totale appliquée sur le piston
par le gaz ? Si le piston est de masse M , trouver quelle est la force instantanée F (t) appliquée par
l’expérimentateur en appliquant la loi fondamentale de la dynamique.
(5) Dans le cas n = 2, préciser l’expression du temps et de la position de formation du choc et
donner l’expression de la vitesse en tout point (x, t) pour tout temps inférieur au temps de formation
du choc.

1.2 Cas de la détente dans une partie semi-infinie avec accélération immédiate
On se place alors dans le cas n = 1 et l’on se focalise uniquement sur le gaz 1 dans cette sous-
partie. Comme dans le cas de l’onde simple de détente traité en cours, on considère une vitesse
négative uf > −2c0 /(γ − 1) et l’on considère le cas limite où l’accélération du piston est infinie en
t = 0 de telle manière que le piston se met immédiatement à la vitesse uf pour tout t > 0.
(1) Montrer que l’écoulement dans le gaz 1, sous les hypothèses données, est homoentropique
pour tout temps t > 0.
Dans la suite, on admet que le diagramme des caractéristiques a la forme donnée sur la Figure
2 de manière cohérente avec les propriétés d’invariance des équations d’Euler. Les domaines 1 et 3
sont des zones d’écoulement uniforme et les droites OA et OI les droites frontières de la région 2.
(2) Donner les propriétés de la région 2. En déduire les équations des deux droites OA et OI.
(3) Donner (u, c/c0 , p/p0 ) dans les domaines 1 et 3.
(4) Soit M (t, x) un point de la région 2. Déterminer (u, c) en ce point en fonction de (t, x).
(5) Pour un point P (τ, X) dans la région 2, vérifier que la C− de cette région qui passe par P
admet pour équation :
! " ! " 3−γ
2 c0 t 2 c0 τ t γ+1
x=− + X+ . (1)
γ−1 γ−1 τ

2
Trajectoire du piston
x = uf t

t 3 2
A

1
x

Figure 2 – Onde de détente centré au point O.

(Indice : pour cela on cherchera une équation sur la nouvelle variable f = x/t comme une fonction
du temps t et on écrira une équation différentielle sur cette variable. On montrera alors que la forme
proposée est bien solution de cette équation différentielle du premier ordre.)

1.3 Bilan sur la présence des deux gaz (Cette sous-partie aura un poids significatif
dans le barême)

On revient sur l’étude de la présence des deux gaz 1 et 4 de chaque côté du piston dans le cas
où n = 1 d’une accélerration immédiate du piston.
(1) Expliquer ce qui se passe dans le gaz 4 et effectuer une représentation précise et justifiée
du diagramme des caractéristiques associé pour les deux gaz à partir des résultats de la sous-partie
précédente.
(2) Ce type de diagramme est proche de celui du tube à choc. Expliquer quelle est la différence
fondamentale entre le cas étudié dans cet exercice et le cas du tube à choc. Dans quel cas les deux
diagrammes sont-ils equivalents ? Représenter dans un diagramme de phase pression-vitesse les deux
cas, afin d’expliquer la différence de manière claire quand elle est présente. Quand les deux cas sont
équivalents, que dire la force appliquée par l’expérimentateur ?

1.4 Cas d’un tube de longueur fini L dans le gaz 1


On considère une zone de gaz 1 dans un tube de longueur L représenté sur la Figure 3. On
veut alors résoudre la question de la propagation et de la réflexion sur le bout du tube de l’onde de
détente initiée par le mouvement du piston.
Ayant résolu le problème associé à l’onde simple de détente avant toute interaction, on se propose
maintenant de considérer le problème de la réflexion de l’onde au bout du tube. Pour des temps
relativement faibles, on suppose que le diagramme des caractéristiques prend la forme donnée sur
la Figure 4.
(a) Donner les propriétés de la région 5 de l’écoulement, quelle est la nature des courbes repré-
sentées dans cette zone de l’écoulement ? Donner leurs propriétés. Quelle est la nature des courbes

3
t=0

piston Gaz
x

Figure 3 – Piston et gaz immobiles à l’instant t = 0 dans un tube de longueur L.

AI et IB.
(b) En admettant éventuellement la relation (1), donner les coordonnées (xI , tI ) du point I.
(c) Donner la valeur de la vitesse sur les demi-droites BB ′ et BB ′′ . En remarquant que deux
courbes caractéristiques d’une même famille ne peuvent pas se couper, démontrer que le domaine
6 est une zone d’écoulement uniforme. Donner dans cette zone les grandeurs (u, c, p, T ) en fonction
des grandeurs de l’énoncé.
(d) On souhaite à présent donner l’évolution de la pression le long de la paroi au cours du temps.
Si h désigne l’enthalpie spécifique du fluide et φ(t, x) un potentiel des vitesses de l’écoulement, on
peut toujours introduire une fonction auxiliaire g, définie tout comme φ à une contante additive
près, telle que
u2
! "
g(u, h) = φ − u x + t h + , (2)
2
ce qui permet d’obtenir les relations importantes :
! " ! " ! "
∂g ∂g ∂g
t= , x=u − . (3)
∂h u ∂h u ∂u h

Montrer que dans le domaine 2, la fonction g est constante ; elle sera alors choisie nulle. Donner les
conditions vérifiées dans le domaine 4 par g sur le tronçon AB et la courbe frontière AI. Le domaine
4 est-il un domaine d’onde simple ?
(e) Nous cherchons dans ce domaine 4 la fonction g. Nous admettons, juste pour cette dernière
question que γ = 5/3 (gaz monoatomique) et que la fonction g prend, sous cette hypothèse, la forme
générale suivante : # ! " ! "$
1 2c 2c
g(u, h) = F u+ +G u− . (4)
c γ−1 γ−1
En utilisant la question précédente, montrer que :
%! "2 &
L γ−1
g(u, h) = g(u, c) = u − c − c20 , (5)
(γ − 1)c 2

dans le domaine 4. En déduire sur le tronçon AB l’équation reliant t/tA à c/c0 . Donner alors les
coordonnées du point B et expliquer comment trouver p(t)/p0 en x = L au cours du temps t.

4
B ′′ B′
6

I 4
3 2
t
A

O L x

Figure 4 – Structure des caractéristiques dans le cas de la réflexion de l’onde de détente.

2 Problème : Evaluation du Mach amont par Tube de Pitot


On considère un écoulement plan permanent autour d’une dièdre symétrique d’angle de déflection
ϕ =15 degrés, représenté sur la Figure 5. Au bord d’attaque, on trouve un choc oblique d’angle σ
par rapport à l’écoulement amont de nombre de Mach M1 , de vitesse V1 = U1 ex . La zone en amont
du choc oblique est notée 1 et les grandeurs qui s’y réfèrent possèdent un indice 1.
La pression mesurée par le Tube de Pitot est de 2.596 atm. Un capteur positionné au point A,
dans la zone 3, derrière l’onde détente, détecte une pression de 0.1 atm. Le but du problème est la
détermination du nombre de Mach amont M1 .

2.1 Onde de choc oblique

(a) Rappeler la structure d’une onde de choc oblique ainsi que la relation qui relie l’angle de
déflection ϕ, l’angle du choc et la donnée de l’écoulement amont. Combien de solutions de cette
relation sont envisageables, l’angle de déflection ϕ étant donné ? Qu’observe-t-on physiquement ?
(b) Donner, vu l’angle de déflection proposé, la gamme de nombre de Mach amont permettant
l’existence d’un choc oblique. Préciser comment obtenir dans ce cas la valeur de σ, le nombre de
Mach amont étant donné.
(c) Dans le cas du dièdre étudié, préciser le domaine spatial dans lequel cette solution est valable.
(d) Pour un nombre de Mach amont fixé, donner l’expression de M2 en fonction de M1 , et des
angles de l’onde de choc oblique σ et de déflection ϕ.

2.2 Onde de détente de Prandtl-Meyer


(a) Rappeler brièvement comment s’effectue la résolution de la détente centrée de Prandtl-Meyer

5
C+

C+

2
3
σ Tube de Pitot

1 O
A
V2 4
M1 > 1
V1
O ϕ = 15◦ −ϕ = −15◦

Figure 5 – Schéma de principe de la mesure de vitesse amont par utilisation d’un Tube de Pitot.

à l’aide de la fonction du même nom.


(b) Présenter précisément la structure de la détente centrée de Prandtl Meyer et donner l’étendue
spatiale de validité de l’onde centrée dans le cas du dièdre considéré.
(c) Si l’on suppose connu le nombre de Mach M3 , comment trouve-t-on ν3 , angle de Mach de
la dernière C+ caractéristique qui délimite l’éventail de détente centrée ?
(d) Par conservation de l’invariant de Riemann, quel est l’angle de Mach ν2 , angle de Mach de
la première C+ caractéristique qui délimite l’éventail de détente centrée en amont.

2.3 Onde de choc devant le Tube de Pitot


L’écoulement dans la zone 3 doit être supersonique, et par conséquent il se forme une onde
de choc normale devant le tube de Pitot ; l’écoulement juste derrière l’onde de choc est noté avec
l’indice 4 comme indiqué sur la Figure 5. Le tube de Pitot permet la mesure de la pression génératrice
p04 = 2.596 atm.
(a) Justifier pourquoi l’écoulement dans la zone 3 doit être supersonique dans la configuration
choisie (indépendemment de la présence du Tube de Pitot).
(b) En utilisant les tables de choc normal, et connaissant la pression p3 au point A ainsi que la
pression d’arrêt ci-dessus, donner la valeur du nombre de Mach M3 (interpoler au besoin pour être
précis).
(c) Expliciter les relations de saut impliquées pour tabuler p04 /p3 en fonction de M3 dans les
tables de choc normal que vous venez d’utiliser. On donnera une formule finale.

6
2.4 Résolution du problème de départ
(a) Calculer les valeurs successives de ν3 , ν2 et finalement M2 en utilisant les questions précé-
dentes.
(b) En déduite l’équation qui lie M1 , M2 et l’angle de l’onde de choc oblique σ en utilisant les
valeurs numériques évaluées.
(c) Démontrer que cela définit une solution unique pour le nombre de Mach amont.
(d) Résoudre l’équation implicite en M1 en faisant des itérations successives à partir de M1 = 4,
et expliquer comment vous faites varier M1 et pourquoi.
(e) Que conclure sur l’existence d’un choc oblique attaché en O ?

7
p = ρ r T r = R/M M cp
γ

t = 0 ] − ∞, 0]
x=0 t=0 X(τ ) τ ∈ [0, +∞[

X (τ ) > 0 X ′′ (τ ) > 0 τ ∈]0, τacc [ τacc
uf t = τacc P
(X(τ ), τ )

Gaz piston en x=0 a t=0


x’ x

Gaz detendu piston en X(t) a l’instant t


x’ x
O
vitesse finale uf

t
2
τacc
P
τ
1
x
O X(τ )

(x, t)

(x, t)

u(t, x) c(t, x) M

J+
J+ C+

(p, u)

t=0 x1 < 0

X1 (t) X1 (0) = x1 W (t) =


(p(t, X1 (t)), u(t, X1 (t), s(t, X1 (t))t W
r1
Trajectoire du piston
x = uf t

2 3 t
A

1
x

O
O

uf
t>0 (p, u)

u c x/t
P (X, τ ) 2 C−
P
! " ! " 3−γ
2 c0 t 2 c0 τ t γ+1
x= + X− .
γ−1 γ−1 τ

C+ u+c
u c x/t t d(x/t)
dt + x/t =
dx
dt
f = x/t t
B′ B ′′
6

4 I
2 3
A t

1
x
L O

AI IB
(xI , tI ) I
BB ′ BB ′′

(u, c, p, T )
C+

(M 1, gamma)
C+
V1

O
V2 2

V3 phi=20°

3
alpha

M1 V1 =
U1 ex
γ = 1.4 M1 = 3
ϕ o

M2
p2 /p1

p/pi
C+

ϕ−α
α ϕ
α

p3 /p2

(ϕ, log(p/p1 ))

O′

T P

p2 p3 L
cT cP
T P
cT = , cP = ,
L q1 L q1
γ 2
q1 2 p1 M1
C+

1
C+

O
V2 2 4
V1
mu

V3 Phi

alpha
3 O’
5

C+

C+
O′
O′

Φ
Φ=0

Φ
Φ=0 Φ
Φ
Φ

p4 /p1 p5 /p1 µ
O′

(ϕ, log(p/p1 ))
C+

1
K
C+

O
V2 2
V1
mu
4
V3
J

alpha
3 O’ 5 C−

C+ I’

C−

C+

C−

I
O′ I ′ J
I I′ J C− I
C+ O′ O
O I

O′ K
O′
2
p = ρrT cp γ

p4 T4
u4 = 0 M4 γ4
p1 << p4 T1 = T4 u1 = 0 M1 >> M4
γ1
(x, t)

(p, u)

u2 = u3
u2 = u3

∑ µ ∂∏ γ4 −1
c1 (γ4 − 1)(M12 − 1) p1 2γ1 2γ4
1− = 1+ (M1 − 1)
2
.
c4 (γ1 + 1)M1 p4 γ1 + 1
p p
c1 = γ1 RT1 /M1 c4 = γ4 RT4 /M4
γ1 = γ4 = 1.4 p4 /p1 >> 1 (M12 − 1)/M1
M1 /M4 9
M1 T2

M∞ U∞ ex

AC

γ = 1.4 M∞ = 2
φ o

σ O σ = 45, 34o
p1
M1 p1 M1 = 1.446 p∞ = 2, 196

M3

p3
M3 = 1, 962 p∞ p3
p1
p∞ M1 M3
1
p03
p0∞ p03
p3
p0∞ p∞ M 3 M∞

� ��
��

C−

� ��
� �� � ��

S
M∞ = 2 U∞ ex
p∞ p0∞ T0∞
E
γ = 1.4

M∞ E
MC C
MC > M∞
MC = 2.77
pC
p∞ pC
1
S
(U∞ ex , p∞ ) (UC ex , pC )
O��
O��

O��

O��

S
M∞ = 2
U∞ ex p∞ p0∞ T0∞
E

OO�
O

AC AS γ = 1.4 M∞ = 2 MC = 2.77
AS
AS
AE
AS
p�0∞ AE
p�0∞
M∞ γ p0∞ (M 2 , γ) M
S AS
p�0∞
p0∞ (M , γ)
2 M

AS

AS

AS

Vous aimerez peut-être aussi