Singularité de Prandtl-Glauert

La singularité de Prandtl-Glauert est une singularité mathématique intervenant dans le calcul du coefficient de pression aérodynamique sur un corps élancé fait par Hermann Glauert<ref>Modèle:Article</ref>. L'expression obtenue montre l'existence d'une indétermination lorsque le nombre de Mach vaut l'unité. Le problème est lié à l'utilisation de simplifications du calcul hors de leur champ de validité, l'expression étant valide dans le domaine des écoulements faiblement transsoniques, typiquement Mach < 0,8. Au-delà apparaissent des fortes détentes qui peuvent être matérialisées sur un avion, lorsque les conditions s'y prêtent, par un nuage de condensation de vapeur d'eau.

Hypothèses, équations de l'écoulement

L'écoulement est supposé descriptible par les équations de Navier-Stokes pour un écoulement faiblement compressible et irrotationnel d'un gaz parfait. Par suite la vitesse dérive d'un potentiel <math>\textstyle \phi</math> :

<math>\vec{V}=\nabla \phi</math>

Le corps est supposé élancé : <math>\textstyle \phi_x \approx V_\infty>>\phi_y,\phi_z</math> partout.

L'équation de conservation de cette quantité s'écrit<ref>Modèle:Ouvrage</ref><ref name=Shapiro1>Modèle:Ouvrage</ref> :

<math>\beta^2\phi_{xx} + \phi_{yy} + \phi_{zz} = 0\,,\quad \beta=\sqrt{1-M_\infty^2}</math>

où <math>\textstyle \phi_{\alpha\alpha}</math> est la dérivée seconde par rapport à α et β le coefficient de Prandl-Glauert. Cette équation est obtenue par linéarisation des équations dans l'hypothèse :

<math>\beta^2>>(\gamma-1)M_\infty^2\frac{\phi_x}{V_\infty}</math>

où γ est l'indice adiabatique.

Dans l'équation le terme en <math>\textstyle \beta^2</math> représente l'effet de la compressibilité.

La vitesse à l'infini amont étant portée par l'axe x ainsi que la composante u de la vitesse, la condition d'entrée est :

<math>u_\infty=\phi_{x_\infty}=V_\infty\,,\quad v_\infty=\phi_{y_\infty}=0 \quad \Rightarrow \quad \phi_\infty=V_\infty x </math>

Les conditions aux limites sur l'objet sont, en remplaçant <math>\textstyle \phi_x</math> par <math>\textstyle V_\infty</math> :

<math>V_\infty \vec{n}_x+\phi_y \vec{n}_y+\phi_z \vec{n}_z=0</math>

où <math>\textstyle \vec{n}</math> est la normale au corps.

Résolution, transformation de Prandl-Glauert

Pour la résolution on utilise la transformation suivante, attribuée à Ernst Mach et Hermann Glauert :

<math>\begin{align}

\bar{x} &= x \\ \bar{y} &= \beta y \\ \bar{z} &= \beta z \\ \bar{\phi} &= \beta^2 \phi \end{align}</math>

Et par suite :

<math>\begin{align}

\bar{n}_{\bar{x}} &= \beta n_x \\ \bar{n}_{\bar{y}} &= n_y \\ \bar{n}_{\bar{z}} &= n_z \end{align}</math>

L'équation de conservation devient une équation de Laplace :

<math>\bar{\phi}_{\bar{x}\bar{x}} + \bar{\phi}_{\bar{y}\bar{y}} + \bar{\phi}_{\bar{z}\bar{z}} = 0</math>

et la condition aux limites est inchangée :

<math>V_\infty \bar{n}_{\bar{x}} + \bar{\phi}_{\bar{y}} \bar{n}_{\bar{y}} + \bar{\phi}_{\bar{z}} \bar{n}_{\bar{z}} = 0</math>

L'équation est résolue par une méthode quelconque et l'on calcule le coefficient de pression transformé <math>\textstyle \bar{C}_p</math>. On remonte ensuite à l'expression dans le domaine physique en utilisant la règle de Göthert<ref name=Shapiro1/><ref>Modèle:Article</ref> :

<math>C_p= -2\frac{\phi_x}{V_\infty} = -\frac{2}{\beta^2} \frac{\bar{\phi}_{\bar{x}}}{V_\infty} = \frac{1}{\beta^2} \bar{C}_p</math>

Cette expression illustre la singularité de Prandl-Glauert, en effet :

<math>\lim\limits_{M_\infty\to 1}C_p=\infty</math>