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[1]. 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 ϕ {\displaystyle \textstyle \phi }  :

V = ϕ {\displaystyle {\vec {V}}=\nabla \phi }

Le corps est supposé élancé : ϕ x V >> ϕ y , ϕ z {\displaystyle \textstyle \phi _{x}\approx V_{\infty }>>\phi _{y},\phi _{z}} partout.

L'équation de conservation de cette quantité s'écrit[2],[3] :

β 2 ϕ x x + ϕ y y + ϕ z z = 0 , β = 1 M 2 {\displaystyle \beta ^{2}\phi _{xx}+\phi _{yy}+\phi _{zz}=0\,,\quad \beta ={\sqrt {1-M_{\infty }^{2}}}}

ϕ α α {\displaystyle \textstyle \phi _{\alpha \alpha }} 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 :

β 2 >> ( γ 1 ) M 2 ϕ x V {\displaystyle \beta ^{2}>>(\gamma -1)M_{\infty }^{2}{\frac {\phi _{x}}{V_{\infty }}}}

où γ est l'indice adiabatique.

Dans l'équation le terme en β 2 {\displaystyle \textstyle \beta ^{2}} 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 :

u = ϕ x = V , v = ϕ y = 0 ϕ = V x {\displaystyle u_{\infty }=\phi _{x_{\infty }}=V_{\infty }\,,\quad v_{\infty }=\phi _{y_{\infty }}=0\quad \Rightarrow \quad \phi _{\infty }=V_{\infty }x}

Les conditions aux limites sur l'objet sont, en remplaçant ϕ x {\displaystyle \textstyle \phi _{x}} par V {\displaystyle \textstyle V_{\infty }}  :

V n x + ϕ y n y + ϕ z n z = 0 {\displaystyle V_{\infty }{\vec {n}}_{x}+\phi _{y}{\vec {n}}_{y}+\phi _{z}{\vec {n}}_{z}=0}

n {\displaystyle \textstyle {\vec {n}}} 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 :

x ¯ = x y ¯ = β y z ¯ = β z ϕ ¯ = β 2 ϕ {\displaystyle {\begin{aligned}{\bar {x}}&=x\\{\bar {y}}&=\beta y\\{\bar {z}}&=\beta z\\{\bar {\phi }}&=\beta ^{2}\phi \end{aligned}}}

Et par suite :

n ¯ x ¯ = β n x n ¯ y ¯ = n y n ¯ z ¯ = n z {\displaystyle {\begin{aligned}{\bar {n}}_{\bar {x}}&=\beta n_{x}\\{\bar {n}}_{\bar {y}}&=n_{y}\\{\bar {n}}_{\bar {z}}&=n_{z}\end{aligned}}}

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

ϕ ¯ x ¯ x ¯ + ϕ ¯ y ¯ y ¯ + ϕ ¯ z ¯ z ¯ = 0 {\displaystyle {\bar {\phi }}_{{\bar {x}}{\bar {x}}}+{\bar {\phi }}_{{\bar {y}}{\bar {y}}}+{\bar {\phi }}_{{\bar {z}}{\bar {z}}}=0}

et la condition aux limites est inchangée :

V n ¯ x ¯ + ϕ ¯ y ¯ n ¯ y ¯ + ϕ ¯ z ¯ n ¯ z ¯ = 0 {\displaystyle V_{\infty }{\bar {n}}_{\bar {x}}+{\bar {\phi }}_{\bar {y}}{\bar {n}}_{\bar {y}}+{\bar {\phi }}_{\bar {z}}{\bar {n}}_{\bar {z}}=0}

L'équation est résolue par une méthode quelconque et l'on calcule le coefficient de pression transformé C ¯ p {\displaystyle \textstyle {\bar {C}}_{p}} . On remonte ensuite à l'expression dans le domaine physique en utilisant la règle de Göthert[3],[4] :

C p = 2 ϕ x V = 2 β 2 ϕ ¯ x ¯ V = 1 β 2 C ¯ p {\displaystyle 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}}

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

lim M 1 C p = {\displaystyle \lim \limits _{M_{\infty }\to 1}C_{p}=\infty }

Références

  1. (en) H. Glauert, « The Effect of Compressibility on the Lift of an Aerofoil », Proceedings of the Royal Society A, vol. 118, no 779,‎ (lire en ligne)
  2. (en) Arnold M. Kuethe et Chuen-Yen Chow, Fundation of aerodynamics. Base of Aerodynamic Design, John Wiley and Sons, (ISBN 0-471-12919-4, lire en ligne)
  3. a et b (en) Asher H. Shapiro, The Dynamics and Thermodynamics of Compressible Fluid Flow, vol. 1, The Ronald Press Company,
  4. (de) B.H. Göthert, « Ebene und räumliche Strömung bei hohen Unterschallgeschwindigkeiten: Erweiterung der Prandtl'schen Regel », Zentrale fuer Wissenschaftliches Berichtswesen, no 127,‎ (lire en ligne)


  • (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Prandtl–Glauert transformation » (voir la liste des auteurs).
  • icône décorative Portail de la physique
  • icône décorative Portail de l’aéronautique