Ouvrir le menu principal

Les relations de Rankine-Hugoniot expriment la discontinuité de diverses quantités au travers d'une onde de choc ou d'une ligne de glissement dans un gaz. Elle a été ainsi dénommée[1] en l'honneur de Pierre-Henri Hugoniot[2] et de William Rankine[3].

Sommaire

Le cas généralModifier

On s'intéresse aux équations aux dérivées partielles en dimension 1 du type :

 

Ces équations sont dites hyperboliques lorsque la matrice jacobienne de f est diagonalisable et a des valeurs réelles positives, ce que l'on suppose vérifié ici. De telles équations admettent des solutions régulières et des solutions discontinues auxquelles on s'intéresse.

On intègre l'équation de conservation ci-dessus au voisinage de l'abscisse xc de la discontinuité :

 

En utilisant la règle d'intégration de Leibniz il vient :

 

On a noté   la vitesse de propagation de la discontinuité.

En faisant   on obtient la relation de saut de w :

 

Afin de simplifier les notations on écrira

 

Équation de BurgersModifier

Un exemple simple est l'équation de Burgers qui correspond à la définition ci-dessus avec   et  .

Dans ce cas l'équation de saut s'écrira :

 
soit
 

On note qu'à la différence des équations d'Euler étudiées ci-dessous, il ne peut pas exister de choc stationnaire.

Équations d'EulerModifier

Problème instationnaireModifier

On applique la relation de saut pour chacune des équations d'Euler :

Continuité      
Quantité de mouvement      
Énergie      

On a noté :

  •   la masse volumique,
  •   la vitesse,
  •   la pression,
  •   l'énergie totale par unité de masse,
  •   l'énergie interne par unité de masse.

La discontinuité peur être de deux ordres :

  • le choc où toutes les quantités sont discontinues,
  • la discontinuité de contact où la vitesse et la pression sont continues. Cela correspond à deux lignes de courant glissant l'une à côté de l'autre sans se pénétrer   et ayant la même pression (conservation de la quantité de mouvement).

Choc droit stationnaireModifier

 
Valeurs des sauts au travers d'un choc droit.

Dans le cas d'un choc stationnaire   tel que rencontré en aérodynamique les relations de saut deviennent :

 

En simplifiant l'équation de saut d'énergie par la relation de saut sur le débit massique on obtient la conservation de l'enthalpie totale :

 

 

Dans le cas d'un gaz parfait   on extrait les relations liant les « sauts », c'est-à-dire le rapport des valeurs en aval et en amont du choc en introduisant le nombre de Mach à l'amont, supposé être une donnée du problème

 

où r est la constante spécifique du gaz.

Ces relations sont les suivantes[4]

 
 
 
 

Le rapport des masses volumiques est limité lorsque  

 

Par exemple pour l'air, formé de molécules diatomiques pour lesquelles  , la limite du rapport des masses volumiques est égal à 6.

On peut également calculer la variation d'entropie réduite   , où CV est la capacité thermique volumique

 

Cette valeur est nulle pour Mg = 1 et croît avec le nombre de Mach. Les valeurs de Mg < 1 conduisent à des valeurs négatives interdites par la thermodynamique. Il n'y a pas de choc pour  .

Choc oblique stationnaireModifier

 
Choc oblique.
 
Loi angle dièdre - angle choc- nombre de Mach pour un choc oblique.

Un choc oblique stationnaire peut être présent sur une configuration géométrique de type dièdre (voir figure). Dans le cadre des équations d'Euler on ignore les phénomènes complexes d'interaction visqueuse au voisinage immédiat de la paroi.

Écrivons la conservation du tenseur de densité de quantité de mouvement au travers du choc[5]

 

  est la normale au choc,   et    les composantes de V normale et parallèle, respectivement. D'où les relations de conservation

 
 

La conservation de la masse s'écrit

 

Des deux dernières équations on déduit la conservation de la vitesse parallèle  , que l'on notera donc  

Le système est donc identique à celui du choc droit pour les vitesses normales     et   , les vitesses parallèles     et   , les nombres de Mach de ces composantes   et  , cette dernière quantité étant donc inférieure à l'unité. Cela ne préjuge en rien du fait que l'écoulement en aval du choc soit ou subsonique ou supersonique. À partir des relations pour le choc droit on peut donc donner

  • l'expression liant l'angle inconnu θ au nombre de Mach et à l'angle du dièdre β sous forme implicite
 
qui se simplifie en
 
  • Les courbes correspondantes (voir figure) montrent que :
    • à un angle de dièdre donné peuvent correspondre deux angles de choc possibles et donc des chocs de force différente (la « force » étant mesurée par exemple par le saut de pression) ; pour un choc faible (généralement rencontré) une augmentation du nombre de Mach entraîne une diminution de l'angle de choc ;
    • pour M donné il existe un maximum de β autorisant un choc droit, au-delà une configuration différente apparaît : un choc courbe à l'avant du pied du dièdre.
  • le saut de pression
 
  • le saut de masse volumique
 
  • le nombre de Mach aval
 


RéférencesModifier

  1. (en) M. D. Salas, « The Curious Events Leading to the Theory of Shock Waves », 17th Shock Interaction Symposium, Rome,‎ (lire en ligne)
  2. Pierre-Henri Hugoniot, « Mémoire sur la propagation des mouvements dans les corps et spécialement dans les gaz parfaits (deuxième partie) », Journal de l'École Polytechnique, vol. 58,‎ , p. 1–125 (lire en ligne)
  3. (en) William Rankine, « On the thermodynamic theory of waves of finite longitudinal disturbances », Philosophical Transactions of the Royal Society of London, vol. 160,‎ , p. 277–288 (lire en ligne)
  4. Marc Buffat, « Relations à travers un choc droit »
  5. (en) John D. Anderson, Jr., Fundamentals of Aerodynamics, McGraw-Hill Education, (ISBN 0-07-001679-8)

Liens externesModifier