Corrigé du TP1.
Données
On commence par calculer un choc stationnaire réaliste dans l’air (gaz idéal avec
γ = Cp /Cv = 1.4). Par exemple
ρL = 1 kg/m3 uL = 1000 m/s pL = 100000 Pa
60 850 2450000 (1)
ρR = kg/m3 uR = m/s pR = Pa
17 3 3
Ce choc est stationnaire et vérifie la condition d’entropie (voir feuille Maple
jointe). Si on échange la gauche et la droite, c’est encore un choc stationnaire,
mais il n’est plus entropique.
On peut alors calculer d’autres quantités (utilisées dans FLUENT), les tem-
pératures TL et TR ainsi que les températures totales T0L et T0R . Pour l’air,
Cp = 1006.43 et Cv = Cp /γ
p = (γ − 1)ρε = (γ − 1)ρCv T
p (2)
T =
(γ − 1)ρCv
Pour un gaz idéal, l’enthalpie vaut
p (γ − 1)ρCv T
h(T ) = ε + = Cv T + = γCv T = Cp T (3)
ρ ρ
et la température totale est définie par
u2
h(T0 ) = h(T ) + (4)
2
On trouve
TL = 347.7638783 TR = 804.6869739
(5)
T0L = T0R = 844.5694187
Maillage GAMBIT
Il faut ensuite faire le maillage d’un tube de 20 m de long sur 1 m de haut. Dé-
clarer les frontières supérieures et inférieures comme des ”wall”. La frontière de
gauche est un ”mass-flow-inlet” et la frontière de droite un ”pressure-outlet”. Dé-
clarer les zones gauche et droite comme des ”fluid”. Sauver le maillage (”export
mesh...” en 2D) au format ”fluent 5/6”.
Calcul FLUENT
Quelques indications pour le calcul FLUENT (voir la documentation de FLU-
ENT pour plus de détails):
1
• lire le fichier .msh (file –> read –> case...)
• solver coupled explicit, unsteady
• équation d’énergie, inviscid
• choisir le matériau ”air” en gaz idéal
• dans operating conditions mettre la pression de référence à zéro !
• condition aux limites à gauche: cliquer sur ”set...” choisir ”mass flux”
comme méthode et mettre ρL uL dans ”mass flux”, T0L dans ”total tem-
perature” et pL dans ”supersonic/initial pressure”
• condition aux limites à droite: cliquer sur ”set...” mettre T0R dans ”back-
flow total temperature” et pR dans ”gauge pressure”
• Initialiser le domaine fluide avec pL , uL et TL puis ”patcher” le domaine
de droite avec pR , uR et TR
• définir une animation
• ”iterate...”
Le choc reste immobile: aucune des variables n’est modifiée.
Recommencer en échangeant la droite et la gauche: des ondes apparaissent.
Les fichiers de calcul FLUENT sont joints.
