Relations de Rankine-Hugoniot

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

Le cas général

modifier

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 Burgers

modifier

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

Un choc stationnaire implique donc nécessairement que .

Équations d'Euler

modifier

Problème instationnaire

modifier

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é peut ê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 stationnaire

modifier
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 stationnaire

modifier
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 nommées polaires de choc (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

On notera que le système dégénère pour un angle égal à l'angle de Mach tel que pour lequel la déviation est nulle et toutes les quantités sont continues.

Références

modifier
  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, New York/St. Louis/Paris etc., McGraw-Hill Education, , 772 p. (ISBN 0-07-001679-8)

Liens externes

modifier